形态学处理 - 图文(7)

317 fe=imerode(f,ones(51,1));%竖线腐蚀

318 subplot(322),imshow(fe);

319 title('使用竖线腐蚀后的结果');

320

321 fo=imopen(f,ones(51,1));%竖线做开运算

322 subplot(323),imshow(fo);

323 title('使用竖线做开运算结果');

324

325 fobr=imreconstruct(fe,f);t做标记

326 subplot(324),imshow(fobr);

327 title('使用竖线做重构开运算');

328

329 ff=imfill(f,'holes');%对f进行孔洞填充

330 subplot(325),imshow(ff);

331 title('对f填充孔洞后的图像');

332

333 fc=imclearborder(f,8);%清除边界，2维8邻接

334 subplot(326),imshow(fc);

335 title('对f清除边界后的图像');

336%图像重构过程显示如下：

337

338 %% 使用顶帽变换和底帽变换

339 clc

340 clear

341 f=imread('.\\images\\dipum_images_ch09\\Fig0926(a)(rice).tif');

342 subplot(221),imshow(f);

343 title('顶帽底帽变换原始图像');

344

345 se=strel('disk',10);%产生结构元素

346 %顶帽变换是指原始图像减去其开运算的图像

347 %而开运算可用于补偿不均匀的背景亮度，所以用一个大的结构元素做开运算后

348 %然后用原图像减去这个开运算，就得到了背景均衡的图像，这也叫做是图像的顶帽运算

349 f1=imtophat(f,se);%使用顶帽变换

350 subplot(222),imshow(f1);

351 title('使用顶帽变换后的图像');

352

353 %底帽变换是原始图像减去其闭运算后的图像

354 f2=imbothat(imcomplement(f),se);%使用底帽变换，为什么原图像要求补呢？

355 ò=imbothat(f,se);%使用底帽变换

356 subplot(223),imshow(f2);

357 title('使用底帽变换后的图像');

358

359 %顶帽变换和底帽变换联合起来用，用于增加对比度

360 f3=imsubtract(imadd(f,imtophat(f,se)),imbothat(f,se));%里面参数好像不合理？

361 subplot(224),imshow(f3);

362 title('使用顶帽底帽联合变换后图像');

363%顶帽底帽变换过程图像如下：

364

365 %%使用开运算和闭运算做形态学平滑

366 %由于开运算可以除去比结构元素更小的明亮细节，闭运算可以除去比结构元素更小的暗色细节

367 %所以它们经常组合起来一起进行平滑图像并去除噪声

368 clc

369 clear

370 f=imread('.\\images\\dipum_images_ch09\\Fig0925(a)(dowels).tif');

371 subplot(221),imshow(f);

372 title('木钉图像原图');

373

374 se=strel('disk',5);%disk其实就是一个八边形

375 fo=imopen(f,se);%经过开运算

376 subplot(222),imshow(f);

377 title('使用半径5的disk开运算后的图像');

378

379 foc=imclose(fo,se);

380 subplot(223),imshow(foc);

381 title('先开后闭的图像');

382

383 fasf=f;

384 for i=2:5

385 se=strel('disk',i);

386 fasf=imclose(imopen(fasf,se),se);

387 end

388 subplot(224),imshow(fasf);

389 title('使用开闭交替滤波后图像');

390%使用开运算和闭运算做形态学平滑结果如下：