【正文】 ,imshow(label2rgb(L, jet, [.])) hold on for k = 1:length(B)boundary = B{k}。noholes39。)。 figure,imshow(bw) % fill any holes, so that regionprops can be used to estimate % the area enclosed by each of the boundaries bw = imfill(bw,39。,2)。s cap se = strel(39。imshow(bw) % remove all object containing fewer than 30 pixels bw = bwareaopen(bw,30)。threshold = graythresh(I)。)。%Step 1: Read image %Step 2: Threshold the image %Step 3: Remove the noise %Step 4: Find the boundaries %Step 5: Determine which objects are round RGB = imread(39。clear。b=transpose(a)。A2:C439。c:/matlab/work/39。%將matlab數(shù)據(jù)寫(xiě)到excel中 a=ones(3)。%在調(diào)用regionprops之前必須將二值圖像轉(zhuǎn)變?yōu)闃?biāo)注矩陣 L = bwlabel(BW)。%只顯示某個(gè)下標(biāo)所對(duì)應(yīng)的物體圖像 bw2=ismember(L,N)。%創(chuàng)建一個(gè)只包含面積大于80的二值圖像 idx = find([] 80)。)。imview(I), imview(BW)%可以使用下面的方法創(chuàng)建相應(yīng)的向量： regionprops(L,39。r = [21 21 75 121 121 75]。)。 %只顯示Solidity numdown=find([] mexdata%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %roipoly函數(shù)的用法 I = imread(39。all39。 BW2 = ismember(labeled,idx)。position39。,[ ])%畫(huà)出每個(gè)物體的外接矩形 bb=[]。,[ ])% rectangle(39。,2)end%畫(huà)出1和2號(hào)物體的外接矩形% rectangle(39。,39。%直接使用打印語(yǔ)句打印序號(hào) text(centr(2*k1),centr(2*k),signal)%按序標(biāo)記物體text(centr(2*k1)30,centr(2*k)30,soli_string)%標(biāo)注每個(gè)Solidity值 endfor k=1:numobjects plot(mexdata(k).ConvexHull(:,1),mexdata(k).ConvexHull(:,2),...39。%d39。%等價(jià)于轉(zhuǎn)字符串 % signal=num2str(nums(k))。%39。for k = 1:numobjects soli=mexdata(k).Solidity。%以下內(nèi)容畫(huà)在同一figure中 centr=[]。)。 mexdata=regionprops(labeled,39。)。,39。 numobjects rgb_label=label2rgb(labeled,spring,39。figure,imshow(nobord)。%figure,imshow(nosmall)。%figure,imshow(bwero)。%figure,imshow(ifill)。holes39。)。%figure,imshow(bwsdil),title(39。,2,0)。se0=strel(39。line39。% figure,imshow(imedcanny)。canny39。% figure,imshow(imed)。figure,imshow(ig)。% figure,imshow(i)。39。close all。 max([]) std([])/mean([])%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %處理混合圖像 clear。all39。 [labeled,numobjects]=bwlabel(nobord,4)。 nobord=imclearborder(nosmall,4)。 nosmall=bwareaopen(bwero,150,4)。 bwero=imerode(ifill,[se90 se0])。)。ifill=imfill(bwsdil,39。dilated39。bwsdil=imdilate(imedcanny,[se90 se0])。line39。,2,90)。 se90=strel(39。)。imedcanny=edge(imed,39。figure,imshow(ig)。imshow(i)。39。clear。, [ ]) figure,imshow(nutdata(1).Image)%只顯示1號(hào)物體的圖像 figure,imshow(nutdata(1).ConvexImage)%畫(huà)出1號(hào)物體的凸多邊形 std([])/ mean([])std([])/ mean([])std([])/ mean([]) std([])/ mean([])std([])/ mean([])std([])/ mean([])std([])/ mean([])std([])/ mean([])std([])/ mean([])std([])/ mean([])%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %處理花椒皮 close all。, [ ])%畫(huà)矩形 rectangle(39。 min([]) rectangle(39。all39。% figure,imshow(rgb_label)。shuffle39。c39。 [labeled,numobjects]=bwlabel(nosmall,4)。nosmall=bwareaopen(imedfill,150)。)。imedfill=imfill(imedcanny,39。canny39。 imed=medfilt2(ig)。ig=rgb2gray(i)。)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %處理花椒子 i=imread(39。 nosmall=bwareaopen(ifill,150)。 ifill=igcfill|vfill。)。 vfill=imfill(vcanny,39。canny39。 v=imadjust(ig,[0 t],[],1)。 figure,imshow(nosmall)。 figure,imshow(ifill)。 figure,imshow(vfill)。hole39。)。 vcanny=edge(v,39。 t=。 figure,imshow(nosmall)。 figure,imshow(ifill)。 figure,imshow(nosmall)。 figure,imshow(imedfill)。hole39。)。%中值濾波后對(duì)圖像邊界有一定的損傷??！ imedcanny=edge(imed,39。 figure,imshow(igcfill)。hole39。,(graythresh(ig)*.1))。igcanny_thresh=edge(ig,39。hole39。)。%igcanny=edge(ig,39。 ig=rgb2gray(i)。)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%直接對(duì)灰度圖進(jìn)行canny運(yùn)算 i=imread(39。%（3）逆操作，回復(fù)矩形原來(lái)大小 bw3=imdilate(bw2,se)。bw2=imerode(bw1,se)。39。,[40 30])。%腐蝕與膨脹聯(lián)合操作 %（1）創(chuàng)建結(jié)構(gòu)元素 se=strel(39。iopen=imopen(bw,se)。disk39。%閾值處理再取反bw=~im2bw(i,graythresh(i))。)。%又叫邊界提取 %去除孤立象素點(diǎn)nosinglepixel=bwmorph(bw,39。canny39。b=bwperim(bw,8)。remove39。)。figure,imshow(BW1)BW2 = bwmorph(BW1,39。39。,Inf)。BW3 = bwmorph(BW1,39。remove39。)。 Subtract the opening from the input imageClass SupportThe input image BW1 can be numeric or must be 2D, real and output image BW2 is = imread(39。 With N = Inf, remove pixels so that an object without holes shrinks to a minimally connected stroke, and an object with holes shrinks to a ring halfway between the hold and outer boundary 39。 With N = Inf, thicken objects by adding pixels to the exterior of objects without connected previously unconnected objects 39。 Remove end points of lines without removing small objects pletely.39。 With N = Inf, remove pixels on the boundaries of objects without allowing objects to break apart 39。shrink objects with holes to connected rings 39。shrink39。 Set a pixel to 0 if its 4connected neighbors are all 139。 Perform binary opening(erosion followed by dilation)39。s 39。majority39。hbreak39。s surrounded by 139。fill39。erode39。dilate39。diag39。close39。s surrounded by 039。clean39。bridge39。bothat39。%標(biāo)記堊白米粒圖像，便于計(jì)數(shù) rate_chalky_grains=num_chalk/num*100。%顯示只含有堊白米粒的圖像 figure,imshow(result)。 [r,c]=find(pare)。% pare=(bw2)amp。%設(shè)置堊白面積的下限 pare=(l)amp。 [l,num]=bwlabel(bw2)。 bw2=imerode(bw1,se)。 bw1=imerode(bw,se)。 figure,imshow(bw)。 figure,imshow(bw_chalk)。 figure,imshow(bw_v)。 bw_chalk=im2bw(chalk,)。 bw_v=im2bw(v,)。 imhist(v,256)。 t=。 figure,imshow(j)。 j=medfilt2(i,[5 5])。 imshow(rgb)。)。 rgb=imread(39。close all。 graindata(numup,:)=[]。graindata numdown=find([] graindata(numdown,:)=[]。all39。figure,imshow(rgb_label)。shuffle39。c39。[label,numobjects]=bwlabel(imnobord,4)。figure,imshow(bw)。level=graythresh(i3)。i3=imadjust(i2,stretchlim(i2),[0 1])。i2=imsubtract(i,background)。,15))。background=imopen(i,strel(39。)。clc i=imread(39。mean(jarea)max(jarea)min(jarea)hist(jarea,255) %std([]）/(Mean([])jstd=std([])jmean=Mean([])jcv=jstd/jmean std([])/ mean([])%面積的變異系數(shù) std([])/ mean([])%橢圓的變異系數(shù) std([])/ mean([]) std([])/ mean([]) std([])/ mean([])%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % close all。)。numobjectsjdata=regionprops(labeled,39。%figure,imshow(mask)。%mask=bwsdilamp。figure,imshow(bwc)。bwc=imclose(bwnobord,se)。disk39。figure,imshow(bwnobord)。%figure,imshow(bwnobord)。%imshow(i2fill)。)。%i2fill=imfill(bwero,39。%bwero=imerode(bwsdil,[se90 se0])。)。ifill=imfill(bwsdil,39。dilated39。bwsdil=imdilate(icanny,[se90 se0])。line39。,2,90)。se90=strel(39。)。icanny=edge(i