Matlab实现描绘子的计算.pdfVIP

8-2 表示与描述二——描述 一、实验目的 ： 1. 理解和掌握基本的边界描绘子，并能利用 Matlab 实现描绘子的计算； 2. 掌握基本的区域描绘子； 3. 熟练使用 IPT 函数 regionprops 计算区域描绘子。 二、实验内容： 1. 对一幅含噪灰度图像，经平滑、阈值分割后，提取目标边界，并计算傅里 叶描绘子。 2. 使用 IPT 函数 regionprops 计算区域描绘子。 三、实验步骤： 1.对一幅含噪灰度图像，经平滑、阈值分割后，提取目标边界，并计算傅里叶描 绘子。 参考代码： %% 傅里叶描述子 clc clear f=imread(Fig1113(a)(chromo_original).tif); subplot(331),imshow(f); whos f;%其实是个二值图像 title(人体染色体原始图); f=im2double(f);%不再是二值图像了 h=fspecial(gaussian,15,9); g=imfilter(f,h,replicate); g=im2bw(g,0.7);%其中的 graythresh()函数为自动阈值，处理后 g 为二值 图像了 subplot(332),imshow(g); title(经阈值处理后); b=boundaries(g); b=b{1}; bim=bound2im(b,344,270);%344,270 为图片的尺寸 subplot(333),imshow(bim); title(提取的边界图); z=frdescp(b);%对边界坐标 b 进行傅里叶变换系数 (有多少个点就有多少 个系数)，将 b 的坐标点看成是复平面中的某个复数 z14=ifrdescp(z,546);%用 50%的描述子进行逆变换 z546im=bound2im(z14,344,270); subplot(334),imshow(z546im); title(546 个描述子恢复后); z110=ifrdescp(z,110);%用 50%的描述子进行逆变换 z110im=bound2im(z110,344,270); subplot(335),imshow(z110im); title(110 个描述子恢复后); z56=ifrdescp(z,56);%用 50%的描述子进行逆变换 z56im=bound2im(z56,344,270); subplot(336),imshow(z56im); title(56 个描述子恢复后); z28=ifrdescp(z,28);%用 50%的描述子进行逆变换 z28im=bound2im(z28,344,270); subplot(337),imshow(z28im); title(28 个描述子恢复后); z14=ifrdescp(z,14);%用 50%的描述子进行逆变换 z14im=bound2im(z14,344,270); subplot(338),imshow(z14im); title(14 个描述子恢复后); z8=ifrdescp(z,8);%用 50%的描述子进行逆变换 z8im=bound2im(z8,344,270); subplot(339),imshow(z8im); title(8 个描述子恢复后); 2. 使用 IPT 函数 regionprops 计算区域描绘子。 参考代码： %% 函数 regionprops 的使用 clc clear B = [1 1 1 0 0 0 0 0 1 1 1 0 1 1 0 0 1 1 1 0 1 1 0 0 1 1 1 0 0 0 0 0 1 1 1 0