形态学处理 - 图文(8)

391

392 %% 颗粒分析

393 clc

394 clear

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

396

397 sumpixels=zeros(1,36);

398 for k=0:35

399 se=strel('disk',k);

400 fo=imopen(f,se);

401 sumpixels(k+1)=sum(fo(:));

402 end

403

404 %可以看到，连续开运算之间的表面积会减少

405 plot(0:35,sumpixels),xlabel('k'),ylabel('surface area');

406 title('表面积和结构元素半径之间的关系');

407%其运算结果如下：

408

409 figure,plot(-diff(sumpixels));%diff()函数为差分或者近似倒数，即相邻2个之间的差值

410 xlabel('k'),ylabel('surface area reduction');

411 title('减少的表面积和结构元素半径之间的关系');

412%其运算结果如下：

413

414 %% 使用重构删除复杂图像的背景

415 clc

416 clear

417 f=imread('.\\images\\dipum_images_ch09\\Fig0930(a)(calculator).tif');

418 subplot(221),imshow(f);

419 title('灰度级重构原图像');

420

421 f_obr=imreconstruct(imerode(f,ones(1,71)),f);

422 subplot(222),imshow(f_obr);

423 title('经开运算重构图');

424

425 f_o=imopen(f,ones(1,71));

426 subplot(223),imshow(f_o);

427 title('经开运算后图');

428

429 f_thr=imsubtract(f,f_obr);

430 subplot(224),imshow(f_thr);

431 title('顶帽运算重构图')

432%使用重构删除复杂图像的背景1:

433

434 f_th=imsubtract(f,f_o)

435 figure,subplot(221),imshow(f_th);

436 title('经顶帽运算图');

437

438 g_obr=imreconstruct(imerode(f_thr,ones(1,11)),f_thr);

439 subplot(222),imshow(g_obr);