[電腦基礎(chǔ)知識(shí)]matlab圖像處理基礎(chǔ)實(shí)例(編輯修改稿)
2024-09-17 15:01 本頁面
【文章內(nèi)容簡(jiǎn)介】 ble(rgb(:,:,2)),sob,39。replicate39。)。bx=imfilter(double(rgb(:,:,3)),sob,39。replicate39。)。by=imfilter(double(rgb(:,:,3)),sob,39。replicate39。)。r_gradiant=mat2gray(max(rx,ry))。g_gradiant=mat2gray(max(gx,gy))。b_gradiant=mat2gray(max(bx,by))。rgb_gradiant=rgb2gray(cat(3,r_gradiant,g_gradiant,b_gradiant))。gxx=rx.^2+gx.^2+bx.^2。gyy=ry.^2+gy.^2+by.^2。gxy=rx.*ry+gx.*gy+bx.*by。theta=.5*((atan(2*gxy./(gxxgyy+eps))))。g1=.5*((gxx+gyy)+(gxxgyy).*cos(2*theta)+2*gxy.*sin(2*theta))。theta=theta+pi/2。g2=.5*((gxx+gyy)+(gxxgyy).*cos(2*theta)+2*gxy.*sin(2*theta))。g1=g1.^5。g2=g2.^5。rgb_vectorgradiant=mat2gray(max(g1,g2))。diff=abs(rgb_vectorgradiantrgb_gradiant)。subplot(331)。imshow(rgb)。title(39。彩色原圖39。)subplot(332)。imshow(r_gradiant)。title(39。R分量邊緣39。)subplot(333)。imshow(g_gradiant)。title(39。G分量邊緣39。)subplot(334)。imshow(b_gradiant)。title(39。B分量邊緣39。)subplot(335)。imshow(rgb_gradiant)。title(39。分量合成邊緣39。)subplot(336)。imshow(rgb_vectorgradiant)。title(39。向量梯度邊緣39。)subplot(337)。imshow(diff)。title(39。差別39。)平滑再合成銳化rgb1=imread(39。39。)。rgb=im2double(rgb1)。fr=rgb(:,:,1)。fg=rgb(:,:,2)。 fb=rgb(:,:,3)。lapmatrix=[1 1 1。1 8 1。1 1 1]。fr_filtered=imfilter(fr,lapmatrix,39。replicate39。)。fg_filtered=imfilter(fg,lapmatrix,39。replicate39。)。fb_filtered=imfilter(fb,lapmatrix,39。replicate39。)。rgb_tmp=cat(3,fr_filtered,fg_filtered,fb_filtered)。rgb_filtered=imsubtract(rgb,rgb_tmp)。i1=(fr+fg+fb)/3。tmp1=min(min(fr,fg),fb)。tmp2=fr+fg+fb。tmp2(tmp2==0)=eps。%???????s=13.*tmp1./tmp2。tmp1=.5*((frfg)+(frfb))。tmp2=sqrt((frfg).^2+(frfb).*(fgfb))。theta=acos(tmp1./(tmp2+eps))。h1=theta。h1(fbfg)=2*pih1(fbfg)。h1=h1/(2*pi)。h1(s==0)=0。lapmatrix=[1 1 1。1 8 1。1 1 1]。i=imfilter(i1,lapmatrix,39。replicate39。)i=imfilter(i1,lapmatrix,39。replicate39。)。subplot(331)。imshow(rgb1)。title(39。pri39。)subplot(332)。imshow(fr_flitered)。title(39。fr_flitered39。)subplot(333)。imshow(fg_filtered)。title(39。fg_filtered39。)subplot(334)。imshow(fb_filtered)。title(39。fb_filtered39。)subplot(335)。imshow(h1)。title(39。h139。)subplot(336)。imshow(s)。title(39。s39。)subplot(337)。imshow(i1)。title(39。i139。)subplot(338)。imshow(rgb_tmp)。title(39。mixture39。)subplot(339)。imshow(cat(3,h1,s,i1))。title(39。hsi mixture39。)彩色圖像向量梯度的計(jì)算 rgb=imread(39。39。)。 sob=fspecial(39。sobel39。)。 rx=imfilter(double(rgb(:,:,1)),sob,39。replicate39。)。 ry=imfilter(double(rgb(:,:,1)),sob39。,39。replicate39。)。 gx=imfilter(double(rgb(:,:,2)),sob,39。replicate39。)。gy=imfilter(double(rgb(:,:,2)),sob39。,39。replicate39。)。 bx=imfilter(double(rgb(:,:,3)),sob,39。replicate39。)。by=imfilter(double(rgb(:,:,3)),sob39。,39。replicate39。)。 gxx=rx.^2+gx.^2+bx.^2。 gyy=ry.^2+gy.^2+by.^2。 gxy=rx.*ry+gx.*gy+bx.*by。 theta=.5*(atan(2*gxy./(gxxgyy+eps)))。 g1=.5*((gxx+gyy)+(gxxgyy).*cos(2*theta)+2*gxy.*sin(2*theta))。 theta=theta+pi/2。 g2=.5*((gxx+gyy)+(gxxgyy).*cos(2*theta)+2*gxy.*sin(2*theta))。 g1=g1.^.5。 g2=g2.^.5。 rgb_gradiant=mat2gray(max(g1,g2))。 subplot(121)。imshow(rgb)。title(39。pri39。) subplot(122)。imshow(rgb_gradiant)。title(39。rgb\_gradiant39。)得到圖像邊界，并計(jì)算周長(zhǎng)、面積和重心坐標(biāo) a=imread(39。39。)。 a1=bwperim(a)。 l=0。 [m,n]=size(a1)。 for i=1:m*nif(a1(i)==1)l=l+1。endend [m,n]=size(a)。s=0。 for i=1:m*nif(a(i)==1)s=s+1。endend x=0。y=0。 for i=1:mfor j=1:nif(a(i,j)==1)x=i+x。y=j+y。endendend x=x/s。y=y/s。 l,s,x,y a=imread(39。39。)%么有加；號(hào) 33 33 33 33 33 33 33 19 33 28 33 33 33 33 33 33 33 33 33 19 33 33 19 33 19 33 33 33 19 21 33 33 33 33 33 33 33 33 33 28 33 33 33 19 33 33 33 33 33 21 33 33 33 33 45 57
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號(hào)-1