• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            Make Oneself More Strong

            mcs51a( 高興就寫點:-),不高興就不寫 )

            圖像形態學膨脹的結構分解opencv實驗

            寬度W,高度H 的圖像,做膨脹操作,如果膨脹的結構元素structure element,大小為kw,kh,那么就需要做W*H*kw*kh次運算,運算量比較大。

            根據圖像形態學的理論,膨脹滿足結合律,即B和C為結構元素。假設一個結構元素S可以表示為兩個結構元素B和C的膨脹,即S=B⊕C,則A⊕S=A⊕(B⊕C)=(A⊕B)⊕C,換言之,用S膨脹A等同于用B先膨脹A,再用C膨脹前面的結果。我們稱S能夠分解成B和C兩個結構元素。結合律很重要,因為計算膨脹所需要的時間正比于結構元素中的非零像素的個數。通過上述推導,分解結構元素,然后再分別用子結構元素進行膨脹操作往往會實現很客觀的速度的增長。

            同樣,腐蝕也可以做結構分解,腐蝕滿足公式 ,B和C為結構元素,同樣如果一個結構元素S可以表示為兩個結構元素B和C的膨脹,即S=B⊕C,那么用S腐蝕A等同于用B先腐蝕A,再用C腐蝕前面的結果。公式推導省略。

            圖像形態學膨脹和腐蝕介紹,可見
                http://en.wikipedia.org/wiki/Erosion_(morphology)
                http://en.wikipedia.org/wiki/Dilation_(morphology)
              

            做結構分解后的運算量為W*H*(kw+kh)


            代碼如下

             1#include <cv.h>
             2#include <highgui.h>
             3#include <stdio.h>
             4
             5
             6int main(int argc, char** argv)
             7{
             8    if(argc<2)  {
             9        printf("has no param\n");
            10        return 0;
            11    }

            12
            13    IplConvKernel *element1 = cvCreateStructuringElementEx( 12500, CV_SHAPE_RECT, 0);
            14    IplConvKernel *element2 = cvCreateStructuringElementEx( 16100, CV_SHAPE_RECT, 0);
            15    IplConvKernel *element3 = cvCreateStructuringElementEx( 162500, CV_SHAPE_RECT, 0);
            16    
            17    IplImage* src=cvLoadImage(argv[1],1); 
            18    
            19
            20    if( src!= NULL)
            21    {
            22        
            23        IplImage* img = cvCreateImage(cvGetSize(src),IPL_DEPTH_8U,1);
            24        cvCvtColor(src,img,CV_BGR2GRAY);
            25        cvReleaseImage(&src);
            26        
            27        
            28        IplImage* tmp = cvCreateImage(cvGetSize(img),IPL_DEPTH_8U,1);
            29        cvDilate( img, tmp, element1, 1);
            30        
            31        IplImage* dst1 = cvCreateImage(cvGetSize(img),IPL_DEPTH_8U,1);            
            32        cvDilate( tmp, dst1, element2, 1);
            33        
            34        
            35        IplImage* dst2 = cvCreateImage(cvGetSize(img),IPL_DEPTH_8U,1);    
            36        cvDilate( img, dst2, element3, 1);
            37        
            38        IplImage* diff = cvCreateImage(cvGetSize(img),IPL_DEPTH_8U,1);    
            39        cvSub(dst2,dst1,tmp,NULL);
            40        cvEqualizeHist( tmp, diff);
            41        cvReleaseImage(&tmp);
            42        
            43        int noZeroCount=cvCountNonZero(diff);
            44        printf("no zerocount %d\n",noZeroCount);
            45        
            46        
            47        cvNamedWindow("img",CV_WINDOW_AUTOSIZE);
            48        cvShowImage("img",img);
            49        cvNamedWindow("dst1",CV_WINDOW_AUTOSIZE);
            50        cvShowImage("dst1",dst1);
            51        cvNamedWindow("dst2",CV_WINDOW_AUTOSIZE);
            52        cvShowImage("dst2",dst2);
            53        cvNamedWindow("diff",CV_WINDOW_AUTOSIZE);
            54        cvShowImage("diff",diff);
            55        
            56        
            57        cvWaitKey();
            58        cvDestroyAllWindows();
            59        
            60        
            61        cvReleaseImage(&img);
            62        cvReleaseImage(&dst1);
            63        cvReleaseImage(&dst2);
            64        cvReleaseImage(&diff);
            65        
            66        cvReleaseStructuringElement(&element1);
            67        cvReleaseStructuringElement(&element2);
            68        cvReleaseStructuringElement(&element3);
            69        
            70        
            71        
            72        
            73        
            74    }

            75    else
            76    {
            77        printf("error,not load\n");
            78    }

            79    return 0;
            80}
            ;
            81

             

             

            posted on 2013-02-21 10:39 mcs51a 閱讀(2304) 評論(0)  編輯 收藏 引用

            久久久久亚洲精品无码网址| 久久国产色av免费看| 精品久久久久久综合日本| 久久精品免费观看| 国内精品久久久久影院老司| 精品久久久久久国产潘金莲| 久久99精品久久久久久不卡| 中文字幕久久亚洲一区| 九九精品99久久久香蕉| 亚洲Av无码国产情品久久| 久久综合狠狠综合久久| 国产精品日韩深夜福利久久| 亚洲成色WWW久久网站| 精品久久人人做人人爽综合| 人妻无码αv中文字幕久久| 久久无码AV中文出轨人妻| 久久国产精品久久| 青青草原精品99久久精品66| 久久久黄色大片| 亚洲精品NV久久久久久久久久 | 久久毛片一区二区| 久久免费高清视频| 久久丫精品国产亚洲av| 欧美久久久久久| 性做久久久久久久久老女人| 久久久久久毛片免费看| 亚洲一区中文字幕久久| 久久99国产精品久久| 婷婷久久久亚洲欧洲日产国码AV | 久久中文字幕人妻熟av女| 久久精品国产99国产精品| 婷婷综合久久中文字幕| 欧美综合天天夜夜久久| 国产精品久久久久9999高清| 91精品国产91久久久久福利| 狠狠干狠狠久久| 国产精品一久久香蕉国产线看| 久久亚洲精品成人av无码网站| 久久精品中文騷妇女内射| 久久天天躁狠狠躁夜夜avapp| 久久精品国产亚洲AV香蕉|