1、-作者xxxx-日期xxxx細胞圖像處理大作業(yè)【精品文檔】 實驗報告實驗課程： 數(shù)字圖像處理 實驗名稱： 骨髓細胞計數(shù) 實驗地點： 計算機中心405 學生姓名： 學 號： 指導(dǎo)教師： 實驗時間： 年 月 一、實驗要求：將被染成紫色的白細胞從較暗的紅細胞和背景中分割出來，去除細胞粘連，用數(shù)字圖像處理技術(shù)實現(xiàn)圖像中紫色細胞的自動計數(shù)。二、實驗過程1、對原圖像進行預(yù)處理（1）由于源圖像是彩色圖像，為了便于進行仿真實驗，在進行處理之前需要將彩色圖像轉(zhuǎn)換為灰度圖像。用到的函數(shù)為rgb2gray（2）獲取圖像的全局閾值并將圖像二值化，以方便后面實驗中對圖像進行形態(tài)學修飾的進行（3）在圖像獲取過程中，由于熱

2、傳輸?shù)扔绊懀瑫a(chǎn)生高斯噪聲，較適合采用中值濾波的方法進行去噪處理。因為，中值濾波的方法能夠較好的保護細胞的邊緣信息，滿足本實驗的需求；（4）刪除小面積圖像（5）對圖像進行開閉運算。開和閉這兩種運算可以除去比結(jié)構(gòu)元素小的特定圖像細節(jié)，同時保證不產(chǎn)生全局幾何失真。開運算可以把比結(jié)構(gòu)元素小的突刺濾掉，切斷細長搭接而起到分離作用；閉運算可以把比結(jié)構(gòu)元素小的缺口或孔填充上，搭接短的間隔而起到連接作用。（6）對圖像進行孔洞填充，（7）開閉運算會留下一些偏離的孤立像素，使用bwareaopen來移除少于特定像素個數(shù)的斑點。BW2= bwareaopen(BW,P)從二值圖像中移除所以少于P像素值的連通塊，得

3、到另外的二值圖像BW3。實驗圖片顯示如下：2、分割計數(shù) 本次實驗采用分水嶺法對進行過預(yù)處理的圖像進行分割和計數(shù)，最后計數(shù)得紫色細胞個數(shù)為35三、實驗心得體會 通過這學期的學習，我獲益匪淺，學到了很多有趣的知識，知道了什么是灰度圖，二值化，閾值，分水嶺等很多有趣的圖像處理方法，知道了怎么使用matlab幫助查找自己不會的函數(shù)，也成功的自己獨立完成了細胞的計數(shù)。這門課也激起了我對圖像處理的興趣，課程結(jié)束后也會繼續(xù)學習更多的函數(shù)知識，不會停止這方面的學習。四、源代碼Image = imread('b1.bmp'); subplot(321);imshow(Image); title(

4、'原圖'); Image=rgb2gray(Image);subplot(322); imshow(Image); title('灰度圖'); Theshold = graythresh(Image);BW = im2bw(Image,1.3*Theshold); subplot(323),imshow(BW); title('二值化圖像'); BW1= medfilt2(BW,13 13);BW1=bwareaopen(BW,50,4)subplot(3,2,4); imshow(BW1); title('刪除小面積后圖像');

5、 se=strel('disk',5);fo=imopen(BW1,se);foc=imclose(fo,se)subplot(3,2,5); imshow(foc); title('開閉運算后圖像');BW2=imfill(foc,'holes');subplot(3,2,6); imshow(BW2);title('孔洞填充后圖像');BW3=bwareaopen(BW2,20);D=bwdist(BW3);D=-D;D(BW3)=-Inf;L=watershed(D);rgb=label2rgb(L,'jet'

6、;);figure;imshow(rgb,);title('分水嶺分割計算后圖像');number=max(max(L)-1;number方法一源程序：A=imread('y1.bmp');B=rgb2gray(A); h=0 1 0,1 -4 1,0 1 0; J1=conv2(B,h,'same'); K1=medfilt2(J1,3 3);¨ BW=edge(K1,'canny'); SE= strel('line',11,90);BW1=imdilate(BW,SE); BW2=bwmorph(B

7、W1,'thin',Inf); H, theta, rho = hough(BW2, 'ThetaResolution', 0.2); subplot(221),imshow(H, , 'XData', theta, 'YData', rho, 'InitialMagnification', 'fit') axis on, axis normal xlabel('theta'), ylabel('rho') ,peaks = houghpeaks(H, 5); hold on plot(theta(peaks(:, 2), rho(peaks(:, 1), .'linestyle', 'none', 'marker', 's', 'color', 'w'); lines = houghlines(BW2, theta, rho, peaks); subplot(222),imshow(BW2), hold on for k = 1:length(lines) xy