纹理 结构特征分析.doc

实验2_纹理_结构特征分析一、实验目的加深对SAR影像纹理、结构特征及相关概念的理解,能编程实现SAR影像纹理、结构特征的统计分析。二、实验软件Matlab三、、。四、,编写统计分析影像纹理特征和结构特征的程序;)纹理特征主程序CTFmain:依次实现的功能为:a打开并读取影像文件;b获取灰度影像;c扩展原始数据;d统计纹理特征频次;e对纹理特征频次排序并得到原始编码索引;f要求用户输入纹理特征分类个数目;g对纹理特征分类并得到类别编码对原始编码的索引;h生成纹理特征图。functionim=CTFmain%选择影像文件[fn,pn]=uigetfile({'*.jpg;*.tif;*.png;*.gif',...'AllImageFiles(*.jpg,*.tif,*.png,*.gif)';...'*.fig','Figures(*.fig)';'*.*','AllFiles'},...'请选择影像文件');iffn~=0%添加当前路径addpath(pn);%读取影像文件dat=imread(fn);%取出起始时间t1=clock;%获取灰度影像data0=GetGMap(dat);fprintf('%s\n','灰度影像生成完毕!');%扩展原始数据矩阵data0为datadata=EnlargeMat(data0,1,1);fprintf('%s\n','影像数据扩展完毕!');%统计每种纹理特征向量出现的频数,得到纹理特征查询矩阵stvm和原始编码图codm[stvm,codm]=StatTVm(data);fprintf('%s\n','纹理特征统计完毕!');%对各种纹理特征向量频数由大到小排序,得到排序后的纹理特征原始编码索引id0=SortGetId0(stvm);%取出终了时间t2=clock;%显示程序运行时间fprintf('%s%.1f%s\n','该阶段运行时间:',etime(t2,t1),'s');%提示用户输入纹理特征分类数目n=input('请根据统计数据输入纹理特征类别的数目n[10]:');ifisempty(n)n=10;end%取出起始时间t1=clock;%合并相似的纹理特征,得到频次和类别编码索引fid以及原始编码和类别编码索引iid[fid,iid]=IcprtTF(stvm,id0,n);fprintf('%s\n','特征类别编码完毕!');%生成纹理特征图im=CreatIM(codm,fid,iid);fprintf('%s\n','纹理特征图生成完毕!');%取出终了时间t2=clock;%显示程序运行时间fprintf('%s%.1f%s\n','该阶段运行时间:',etime(t2,t1),'s');%显示纹理特征图imshow(im,[min(min(im)),max(max(im))]);elseerror('您没有选择影像文件!');endend2)获取灰度图像子程序GetGMap:实现了根据输入数据结构采用相应算法得到灰度图像的功能,即:a如果是单层影像数据,直接认为是灰度图像,有data0=dat;b如果是3层影像数据,采用gray=*r+*g+*b公式计算出灰度图像:c如果是多余3层的影像数据,采用所有数据层平均的方法得到灰度图像。functiondata0=GetGMap(dat)ds=size(dat);data0=zeros(ds(1),ds(2));dss=size(ds);ifdss(2)==3data0=*dat(:,:,1)+*dat(:,:,2)+*dat(:,:,3);elseifdss(2)>3fori=1:dss(2)data0=data0+dat(:,:,i);endelsedata0=dat;endendend3)扩展原始数据子程序EnlargeMat:实现将数据data0扩展为data的功能,具体的方法是:a将data0的4个角点值分别赋给在其所处对角线另一端的data的角点;b将data0的初始端制定数目的行赋给data最末端制定数目的行,而data0最末端制定数目的行赋给data初始端制定数目的行;c与b相同的方法给data的初始端制定数目的列和最末端制定数目的列赋值。functiondata=EnlargeMat(data0,dr,dc)%取出原始数据的行列数[row,col]=size(data0);%扩展data0为datadata=zer