数字图像处理图像描述..docVIP

DIP-6 图像描述 一．题目要求 对下图象提取多边形对象的边缘，根据傅氏描述算子重建该图多边形的边界 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 1 1 1 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 算法设计 傅立叶描述子（Fourier Descriptor，简称FD）常用来表示单封闭曲线的形状特征，其基本思想是将目标轮廓曲线建模成一维序列，对该序列进行一维的傅立叶变换，从而获得一系列的傅立叶系数，用这些系数对该目标轮廓进行描述。傅立叶描述子方法有一系列优点，如：计算原理简单，描述清晰，具有由粗及精的特性等。计算原理简单可以使得特征提取更加稳定，因为在计算的过程中，无须设置大量控制参数就可以获得结果，计算的一致性好。傅立叶描述子具有明确的物理或几何意义，它比某些特征描述子（如Hu 不变矩）更具直观性。此外，由于任何一个序列经傅立叶变换后，其能量主要集中于少数几个低频傅立叶系数 上，因此采用极少的傅立叶系数就可以描绘该序列特征。同时，随着傅立叶系数 的增多，该序列的细节特征得以更好地描述。因此，傅立叶描述子对目标轮廓有 非常好的由粗及精的描述能力。一个傅立叶描述子的构建包括两步：首先，定义 一种好的表示（representation）方法对轮廓曲线进行描述；然后，采用傅立叶理论对该曲线进行变换。不同的曲线表示方法有不同的特性，一个好的表示方法 应该使最终获得的傅立叶描述子具有尺度、旋转、平移不变性及起始点的无关性。 傅立叶描述子，是物体形状边界曲线的傅立叶变换系数，是物体边界曲线信号的频域分析结果。它是一种描述不受起始点移动尺寸变化及旋转影响的曲线的方法。傅立叶描述子的基本思想，是把坐标的序列点看作复数： 即x 轴作为实轴，y 轴作为虚轴，边界的性质不变。这种表示方法的优点，是将 一个二维问题简化成一个一维问题。对s(k)的傅立叶变换为： 傅立叶描述子序列反映了原曲线的形状特征，同时，由于傅立叶变换具有能 量集中性，因此，少量的傅立叶描述子就可以重构出原曲线。 三．实现代码 主函数 %图像傅里叶算子边界描述 %生成图像边界矩阵 clc; r=zeros(15,15); %构造15*15的0矩阵 r(3,8:13)=1; %按要求构造“R”型图像 r(4,7)=1; r(4,13)=1; r(5,6)=1; r(5,13)=1; r(6,3:5)=1; r(6,13)=1; r(7,3)=1; r(7,13)=1; r(8,3)=1; r(8,13)=1; r(9,3)=1; r(9,13)=1; r(10,3)=1; r(10,13)=1; r(11,3)=1; r(11,12)=1; r(12,3)=1; r(12,11)=1; r(13,3:10)=1; subplot(121); imshow(r); title(构建的图像边界); %主函数 %图像傅里叶算子边界描述 i=1; for m=1:15 for n=1:15 if r(m,n)==1 %将边界坐标存入S矩阵 s(i,1)=m; s(i,2)=n; i=i+1; end end end z=frdescp(s); %图像傅里叶算子边界逆描述 z34=ifrdescp(z,34); z34=uint8(z34); I=zeros(15,15); for k=1:34 I(z34(k,1),z34(k,2))=1;