模式识别第三次作业.docVIP

模式识别作业三 一、编程要求： 编程实现KL变换，并对TM六波段图像进行演算。 KL变换的思想是：从n维特征选取m维特征，删去的n-m维特征不肯定就是无用的信息，如何在信息损失最小的状况下选取特征，在理论上就显得更严密些。通常采纳正交变换，得到新的经变换的模式， 以保证信息损失最小状况下获得有利于分类的特征。 二、编程思想： 将6副图象依次输入获得灰度值存在一个6*size(size为一副图象的像素数)的二维数组中，计算每个波段的灰度均值，然后计算协方差矩阵，得出特征值和特征向量，再重新计算要输出的处理的图象的灰度。 三、核心代码： void CKLchangeView::Oninput() //TM图象数据读入 { AfxMessageBox(留意:请一次将6副图象都读入!);//读图象的方法 int flag=0; while(flag6) { CFileDialog dlg(TRUE,bmp,.bmp); if(dlg.DoModal()==IDOK) { m_NewFileName[flag]=dlg.GetFileName(); for(int j=0;j5;j++) { if(flag==j||m_NewFileName[flag]!=m_NewFileName[j]) continue; else { AfxMessageBox(已打开过的图象!请按方法重新打开!); flag=0; for(int i=0;i6;i++) m_NewFileName[i]=; return; } } NewFileName=m_NewFileName[flag]; flag++; }//推断是否是已经打开的图象，若已经打开，则要重新按方法打开 else //不读图象，或是只是要看图象 { for(int i=0;i6;i++) m_NewFileName[i]=; return; } CDC *pDC = GetWindowDC(); CFile f(NewFileName,0); m_pDib.Read(f);//读图象 if(m_pDib.m_nColorTableEntries!=256)//推断是否是灰度图象 { AfxMessageBox(请打开灰度图象!); return; } CKLchangeDoc* pDoc = GetDocument(); pDoc-pDib=m_pDib; OnDraw(pDC);//调用OnDraw函数绘图 sizeImage = m_pDib.GetDimensions(); sizeSaveImage = m_pDib.GetDibSaveDim(); a=sizeSaveImage.cx * sizeSaveImage.cy;//象素数 if(flag==1) tm=CMatrix(6,a);//开第一副图象时为矩阵安排内存 //将每幅图象象素灰度值放入矩阵中 if(flag==1) for(int i=0;ia;i++) { tm.m_pData[i]=*(m_pDib.m_lpImage+sizeSaveImage.cx*(i/sizeImage.cx)+i%sizeImage.cx); } if(flag==2) for(int i=0;ia;i++) { tm.m_pData[a+i]=*(m_pDib.m_lpImage+sizeSaveImage.cx*(i/sizeImage.cx)+i%sizeImage.cx); } if(flag==3) for(int i=0;ia;i++) { tm.m_pData[2*a+i]=*(m_pDib.m_lpImage+sizeSaveImage.cx*(i/sizeImage.cx)+i%sizeImage.cx); } if(flag==4) for(int i=0;ia;i++) { tm.m_pData[3*a+i]=*(m_pDib.m_lpImage+sizeSaveImage.cx*(i/sizeImage.cx)+i%sizeImage.cx); } if(flag==5) for(int i=0;ia;i++) { tm.m_pData[4*a+i]=*(m_pDib.m_lpImage+sizeSaveImage.cx*(i/sizeImage.cx)+i%sizeImage.cx); } if(flag==6