第九周 项目三--稀疏矩阵的三元组表示的实现及应用

tup.h

[cpp]  view plain  copy
  在CODE上查看代码片 派生到我的代码片
  1. #ifndef TUP_H_INCLUDED  
  2. #define TUP_H_INCLUDED  
  3.   
  4. #define M 6  
  5. #define N 7  
  6. #define MaxSize  100         //矩阵中非零元素最多个数  
  7. typedef int ElemType;  
  8.   
  9. typedef struct  
  10. {  
  11.     int r;                  //行号  
  12.     int c;                  //列号  
  13.     ElemType d;             //元素值  
  14. } TupNode;                  //三元组定义  
  15.   
  16. typedef struct  
  17. {  
  18.     int rows;               //行数  
  19.     int cols;               //列数  
  20.     int nums;               //非零元素个数  
  21.     TupNode data[MaxSize];  
  22. } TSMatrix;                 //三元组顺序表定义  
  23.   
  24. void CreatMat(TSMatrix &t,ElemType A[M][N]);  //从一个二维稀疏矩阵创建其三元组表示  
  25. bool Value(TSMatrix &t,ElemType x,int i,int j);  //三元组元素赋值  
  26. bool Assign(TSMatrix t,ElemType &x,int i,int j); //将指定位置的元素值赋给变量  
  27. void DispMat(TSMatrix t);//输出三元组  
  28. void TranTat(TSMatrix t,TSMatrix &tb);//矩阵转置  
  29.   
  30. #endif // TUP_H_INCLUDED  

tup.cpp

[cpp]  view plain  copy
  在CODE上查看代码片 派生到我的代码片
  1. #include "stdio.h"  
  2. #include "tup.h"  
  3.   
  4. void CreatMat(TSMatrix &t,ElemType A[M][N])  //从一个二维稀疏矩阵创建其三元组表示  
  5. {  
  6.     int i,j;  
  7.     t.rows=M;  
  8.     t.cols=N;  
  9.     t.nums=0;  
  10.     for (i=0; i<M; i++)  
  11.     {  
  12.         for (j=0; j<N; j++)  
  13.             if (A[i][j]!=0)     //只存储非零元素  
  14.             {  
  15.                 t.data[t.nums].r=i;  
  16.                 t.data[t.nums].c=j;  
  17.                 t.data[t.nums].d=A[i][j];  
  18.                 t.nums++;  
  19.             }  
  20.     }  
  21. }  
  22.   
  23. bool Value(TSMatrix &t,ElemType x,int i,int j)  //三元组元素赋值  
  24. {  
  25.     int k=0,k1;  
  26.     if (i>=t.rows || j>=t.cols)  
  27.         return false;               //失败时返回false  
  28.     while (k<t.nums && i>t.data[k].r) k++;                  //查找行  
  29.     while (k<t.nums && i==t.data[k].r && j>t.data[k].c) k++;//查找列  
  30.     if (t.data[k].r==i && t.data[k].c==j)   //存在这样的元素  
  31.         t.data[k].d=x;  
  32.     else                                    //不存在这样的元素时插入一个元素  
  33.     {  
  34.         for (k1=t.nums-1; k1>=k; k1--)  
  35.         {  
  36.             t.data[k1+1].r=t.data[k1].r;  
  37.             t.data[k1+1].c=t.data[k1].c;  
  38.             t.data[k1+1].d=t.data[k1].d;  
  39.         }  
  40.         t.data[k].r=i;  
  41.         t.data[k].c=j;  
  42.         t.data[k].d=x;  
  43.         t.nums++;  
  44.     }  
  45.     return true;                        //成功时返回true  
  46. }  
  47.   
  48. bool Assign(TSMatrix t,ElemType &x,int i,int j)  //将指定位置的元素值赋给变量  
  49. {  
  50.     int k=0;  
  51.     if (i>=t.rows || j>=t.cols)  
  52.         return false;           //失败时返回false  
  53.     while (k<t.nums && i>t.data[k].r) k++;                  //查找行  
  54.     while (k<t.nums && i==t.data[k].r && j>t.data[k].c) k++;//查找列  
  55.     if (t.data[k].r==i && t.data[k].c==j)  
  56.         x=t.data[k].d;  
  57.     else  
  58.         x=0;                //在三元组中没有找到表示是零元素  
  59.     return true;            //成功时返回true  
  60. }  
  61.   
  62. void DispMat(TSMatrix t)        //输出三元组  
  63. {  
  64.     int i;  
  65.     if (t.nums<=0)          //没有非零元素时返回  
  66.         return;  
  67.     printf("\t%d\t%d\t%d\n",t.rows,t.cols,t.nums);  
  68.     printf("\t------------------\n");  
  69.     for (i=0; i<t.nums; i++)  
  70.         printf("\t%d\t%d\t%d\n",t.data[i].r,t.data[i].c,t.data[i].d);  
  71. }  
  72.   
  73. void TranTat(TSMatrix t,TSMatrix &tb)       //矩阵转置  
  74. {  
  75.     int p,q=0,v;                    //q为tb.data的下标  
  76.     tb.rows=t.cols;  
  77.     tb.cols=t.rows;  
  78.     tb.nums=t.nums;  
  79.     if (t.nums!=0)                  //当存在非零元素时执行转置  
  80.     {  
  81.         for (v=0; v<t.cols; v++)        //tb.data[q]中的记录以c域的次序排列  
  82.             for (p=0; p<t.nums; p++)    //p为t.data的下标  
  83.                 if (t.data[p].c==v)  
  84.                 {  
  85.                     tb.data[q].r=t.data[p].c;  
  86.                     tb.data[q].c=t.data[p].r;  
  87.                     tb.data[q].d=t.data[p].d;  
  88.                     q++;  
  89.                 }  
  90.     }  
  91. }  

main.cpp

[cpp]  view plain  copy
  在CODE上查看代码片 派生到我的代码片
  1. /*                  
  2. Copyright (c)2016,烟台大学计算机与控制工程学院                  
  3. All rights reserved.                  
  4. 文件名称:3.cpp                  
  5. 作    者:  陈朋           
  6. 完成日期:2016年11月3日                  
  7. 版 本 号:v1.0                     
  8. 问题描述:             
  9. 输入描述:无              
  10. 程序输出:若干。               
  11. */          
  12. #include <stdio.h>  
  13. #include "tup.h"  
  14. int main()  
  15. {  
  16.     TSMatrix t,tb;  
  17.     int x,y=10;  
  18.     int A[6][7]=  
  19.     {  
  20.         {0,0,1,0,0,0,0},  
  21.         {0,2,0,0,0,0,0},  
  22.         {3,0,0,0,0,0,0},  
  23.         {0,0,0,5,0,0,0},  
  24.         {0,0,0,0,6,0,0},  
  25.         {0,0,0,0,0,7,4}  
  26.     };  
  27.     CreatMat(t,A);  
  28.     printf("b:\n");  
  29.     DispMat(t);  
  30.     if (Assign(t,x,2,5)==true)  //调用时返回true  
  31.         printf("Assign(t,x,2,5)=>x=%d\n",x);  
  32.     else  //调用时返回false  
  33.         printf("Assign(t,x,2,5)=>参数错误\n");  
  34.     Value(t,y,2,5);  
  35.     printf("执行Value(t,10,2,5)\n");  
  36.     if (Assign(t,x,2,5)==true)  //调用时返回true  
  37.         printf("Assign(t,x,2,5)=>x=%d\n",x);  
  38.     else  //调用时返回false  
  39.         printf("Assign(t,x,2,5)=>参数错误\n");  
  40.     printf("b:\n");  
  41.     DispMat(t);  
  42.     TranTat(t,tb);  
  43.     printf("矩阵转置tb:\n");  
  44.     DispMat(tb);  
  45.     return 0;  
  46. }  
  47. 运行结果:

### 光流法C++源代码解析与应用 #### 光流法原理 光流法是一种在计算机视觉领域中用于追踪视频序列中运动物体的方法。它基于亮度不变性假设,即场景中的点在时间上保持相同的灰度值,从而通过分析连续帧之间的像素变化来估计运动方向和速度。在数学上,光流场可以表示为像素位置和时间的一阶导数,即Ex、Ey(空间梯度)和Et(时间梯度),它们共同构成光流方程的基础。 #### C++实现细节 在给定的C++源代码片段中,`calculate`函数负责计算光流场。该函数接收一个图像缓冲区`buf`作为输入,并初始化了几个关键变量:`Ex`、`Ey`和`Et`分别代表沿x轴、y轴和时间轴的像素强度变化;`gray1`和`gray2`用于存储当前帧和前一帧的平均灰度值;`u`则表示计算出的光流矢量大小。 #### 图像处理流程 1. **初始化和预处理**:`memset`函数被用来清零`opticalflow`数组,它将保存计算出的光流数据。同时,`output`数组被填充为白色,这通常用于可视化结果。 2. **灰度计算**:对每一像素点进行处理,计算其灰度值。这里采用的是RGB通道平均值的计算方法,将每个像素的R、G、B值相加后除以3,得到一个近似灰度值。此步骤确保了计算过程的鲁棒性和效率。 3. **光流向量计算**:通过比较当前帧和前一帧的灰度值,计算出每个像素点的Ex、Ey和Et值。这里值得注意的是,光流向量的大小`u`是通过`Et`除以`sqrt(Ex^2 + Ey^2)`得到的,再乘以10进行量化处理,以减少计算复杂度。 4. **结果存储与阈值处理**:计算出的光流值被存储在`opticalflow`数组中。如果`u`的绝对值超过10,则认为该点存在显著运动,因此在`output`数组中将对应位置标记为黑色,形成运动区域的可视化效果。 5. **状态更新**:通过`memcpy`函数将当前帧复制到`prevframe`中,为下一次迭代做准备。 #### 扩展应用:Lukas-Kanade算法 除了上述基础的光流计算外,代码还提到了Lukas-Kanade算法的应用。这是一种更高级的光流计算方法,能够提供更精确的运动估计。在`ImgOpticalFlow`函数中,通过调用`cvCalcOpticalFlowLK`函数实现了这一算法,该函数接受前一帧和当前帧的灰度图,以及窗口大小等参数,返回像素级别的光流场信息。 在实际应用中,光流法常用于目标跟踪、运动检测、视频压缩等领域。通过深入理解和优化光流算法,可以进一步提升视频分析的准确性和实时性能。 光流法及其C++实现是计算机视觉领域的一个重要组成部分,通过对连续帧间像素变化的精细分析,能够有效捕捉和理解动态场景中的运动信息
微信小程序作为腾讯推出的一种轻型应用形式,因其便捷性与高效性,已广泛应用于日常生活中。以下为该平台的主要特性及配套资源说明: 特性方面: 操作便捷,即开即用:用户通过微信内搜索或扫描二维码即可直接使用,无需额外下载安装,减少了对手机存储空间的占用,也简化了使用流程。 多端兼容,统一开发:该平台支持在多种操作系统与设备上运行,开发者无需针对不同平台进行重复适配,可在一个统一的环境中完成开发工作。 功能丰富,接口完善:平台提供了多样化的API接口,便于开发者实现如支付功能、用户身份验证及消息通知等多样化需求。 社交整合,传播高效:小程序深度嵌入微信生态,能有效利用社交关系链,促进用户之间的互动与传播。 开发成本低,期短:相比传统应用程序,小程序的开发投入更少,开发期更短,有助于企业快速实现产品上线。 资源内容: “微信小程序-项目源码-原生开发框架-含效果截图示例”这一资料包,提供了完整的项目源码,并基于原生开发方式构建,确保了代码的稳定性与可维护性。内容涵盖项目结构、页面设计、功能模块等关键部分,配有详细说明与注释,便于使用者迅速理解并掌握开发方法。此外,还附有多个实际运行效果的截图,帮助用户直观了解功能实现情况,评估其在实际应用中的表现与价值。该资源适用于前端开发人员、技术爱好者及希望拓展业务的机构,具有较高的参考与使用价值。欢迎查阅,助力小程序开发实践。资源来源于网络分享,仅用于学习交流使用,请勿用于商业,如有侵权请联系我删除!
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值