《数字摄影测量学》课堂实验报告(遥感).doc

河南理工大学测绘学院 《数字摄影测量学》教学实验报告 （专业必修课） ┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄┄ 实验成绩： 评 语: 指导老师签名： 2012年 月 日 实习报告一：点特征提取 一、、double WINAPI Det(double N[], int n) { if(n==4) return (N[0]*N[3]-N[1]*N[2]); else return 0; } //计算矩阵N之迹（即对角线元素之和）的值 //参数：大小为n的数组、数组宽度 //返回值：矩阵N之迹 double WINAPI Tr(double N[], int w) { double sum=0.0; for(int i=0;iw;i++) sum += *(N+i*w+i); return sum; } //判断该点是否为初选点 //参数：大小为n的数组、数组维数、阈值 //返回值：是，返回真；否，返回假 BOOL WINAPI Mid(int dg[], int n,int T) { BOOL flag=FALSE; int i=0,num=0; for(;in;i++) if(*(dg+i)T) num++; if(num1) flag=TRUE; return flag; } /////////////////////////////////////////////////////////////////////////// //特征提取的函数（仅对灰度图象） //参数：块图象数据指针或数组、块图象宽、块图象高、局部非最大窗口宽、局部非最大窗口高 //返回值：标记特征点位的二值化块图象指针 void WINAPI ForstnerCha(BYTE image[],BYTE* newimage, int width, int height, int win_w, int win_h) { static int T=10;//提取初选点的阈值 static double Tq=0.40;//阈值 int i,j,l,m,n;//计数器 int dg[4];//上下左右四个方向的差分绝对值 // BYTE* newimage=new BYTE[width*height];//为保存特征点开辟的数组 BYTE t,t1,t2,t3,t4;//替换值变量 //提取初选点 for(i=0;iheight;i++) for(j=0;jwidth;j++) { if(i==0 || j==0 || i==height-1 || j==width-1) *(newimage+i*width+j)=0;//窗口边缘点置零 else { t=*(image+i*width+j); t1=*(image+(i+1)*width+j); t2=*(image+i*width+(j+1)); t3=*(image+(i-1)*width+j); t4=*(image+i*width+(j-1)); dg[0]=abs(t-t1);//向下的差分绝对值 dg[1]=abs(t-t2);//向右的差分绝对值 dg[2]=abs(t-t3);//向上的差分绝对值 dg[3]=abs(t-t4);//向左的差分绝对值 //若四个方向的差分绝对值有任意两个大于阈值，为真，否则为假 if( Mid(dg,4,T) ) *(newimage+i*width+j)=255;//若为真 else *(newimage+i*width+j)=0;//若为假 } } float* w=new float[width*height];//存放各点权值的数组 int k=3/2; double N[4];//窗口的灰度协方差矩阵 double gu,gv;//像素的Robert`s梯度 float q;//兴趣值 //提取极值点 for(i=0;iheight;i++) for(j=0;jwidth;j++) { *(w+i*width+j)=0.0;//权值数组各元素置零 if(*(newimage+i*width+j)!=0) { for(l=0;l4;l++) N[l]=0;//灰度协方差矩阵置零 for(m=i-k;m=i+k-1;