形态学处理 - 图文(3)

127 title('使用square(3)开操作后的图像');

128

129 fc=imclose(f,se);

130 subplot(222),imshow(fc);

131 title('使用square闭操作后的图像');

132

133 foc=imclose(fo,se);

134 subplot(223),imshow(foc);

135 title('使用square(3)先开后闭操作后的图像')

136

137 fco=imopen(fc,se);

138 subplot(224),imshow(fco);

139 title('使用square(3)先闭后开操作后的图像');

140%指纹图像开闭操作过程结果如下：

141

142 %% bwhitmiss击中或击不中变换

143 clc

144 clear

145 f=imread('.\\images\\dipum_images_ch09\\Fig0913(a)(small-squares).tif');

146 imshow(f);

147 title('击中或不击中原始图像');

148%击中或不击中原始图像显示结果如下：

149

150 B1=strel([0 0 0;0 1 1;0 1 0]);%击中：要求击中所有1的位置

151 B2=strel([1 1 1;1 0 0;1 0 0]);%击不中，要求击不中所有1的位置

152 B3=strel([0 1 0;1 1 1;0 1 0]);%击中

153 B4=strel([1 0 1;0 0 0;0 0 0]);%击不中

154 B5=strel([0 0 0;0 1 0;0 0 0]);%击中

155 B6=strel([1 1 1;1 0 0;1 0 0]);%击不中

156

157 g=imerode(f,B1)&imerode(~f,B2)%利用定义来实现击中或击不中

158 figure,subplot(221),imshow(g);

159 title('定义实现组1击中击不中图像');

160

161 g1=bwhitmiss(f,B1,B2);

162 subplot(222),imshow(g1);

163 title('结构数组1击中击不中后的图像');

164

165 g2=bwhitmiss(f,B3,B4);

166 subplot(223),imshow(g2);

167 title('结构数组2击中击不中的图像');

168

169 g3=bwhitmiss(f,B5,B6);

170 subplot(224),imshow(g3);

171 title('结构数组3击中击不中的图像');

172%击中击不中变换后图像如下：

173

174 %%makelut

175 clc

176 clear

177

178 f=inline('sum(x(:))>=3');%inline是用来定义局部函数的

179 lut2=makelut(f,2)%为函数f构造一个接收2*2矩阵的查找表

180 lut3=makelut(f,3)

181

182 %% Conway生命游戏

183 clc

184 clear

185 lut=makelut(@conwaylaws,3);

186 bw1= [0 0 0 0 0 0 0 0 0 0

187 0 0 0 0 0 0 0 0 0 0

188 0 0 0 1 0 0 1 0 0 0

189 0 0 0 1 1 1 1 0 0 0

190 0 0 1 0 0 0 0 1 0 0

191 0 0 1 0 1 1 0 1 0 0

192 0 0 1 0 0 0 0 1 0 0