1、并行圖像細化算法和c代碼實現(xiàn)圖像細化一般作為一種圖像預(yù)處理技術(shù)出現(xiàn),目的是提取源圖像的骨架,即是將原圖像中線條寬度大于1個象素的線條細化成只有一個象素寬,形成骨架,形成骨架后能比較容易的分析圖像,如提取圖像的特征.細化分成串行細化和并行細化2中,串行細化即是一遍檢測滿足細化條件的點一邊刪除細化點,并行細化即是檢測細化點的時候不進行點的刪除只進行標(biāo)記,而在檢測完整幅圖像后一次性去除要細化的點.細化基本思想是層層剝奪,即從線條邊緣開始一層一層向里剝奪,直到線條剩下一個象素的為止.進行細化算法前要先對圖像進行2值化,即圖像中直包含黑和白2中顏色.細化算法:在微觀上取檢測點的8個臨域(由于是并行細化,

2、有些模板要擴展為12臨域),如下xxxxoxxxx其中o為檢測點x為其相鄰點以下用1代表黑色點,用0代表白色點,用x代表任意顏色的點,要剝奪(刪除)的點應(yīng)滿足一下8個模板中的一個.模板a(向右擴大)0x1x01110x1x模板b(向右擴大)00xx0111x11x模板c(向右擴大)x11x011100xx模板d111x1x000模板e1x01101x0模板fx00110x1x模板gx1x110x00模板h(向下擴大)000x1x111x1x符合以上8個模板的點為要剝奪的點,因為符合這8個模板的點可以確認(rèn)為線條邊沿上的點.而試事實上經(jīng)過這8個模板并行細化后還有下面2種特殊的邊沿點保留了下來,特殊

3、邊沿點1000010111特殊邊沿點2001011001造成這種2種特殊點的原因擴大后的模板a和擴大后的模板h,擴大的的本意是防止偶數(shù)列(行)的線條被完全消去(并行細化并然的).解決方法是在并行細化后再進行一次串行細化,選取縮小后的模板a和模板h模板a(縮小后)0x10110x1模板h(縮小后)000x1x111其中縮小后的模板a解決了特殊情況1,縮小后的模板h解決了特殊情況2,注意這次是串行細化了.一下根據(jù)這個原理用c+builder6.0實現(xiàn),關(guān)鍵代碼如下./-bcb6 代碼#include #pragma hdrstop#include#include unit1.h#includefi

4、le1.h#include#include#include#pragma pack(1)using namespace std;/* 程序：圖像細化 作者:sboom(lingch) 日期:05年1月18日*/bmp文件頭struct bitmapfileheader_ short type; int bfsize; short re1,re2; int offbits;/bmp信息頭struct bitmapinfo_ long size; long width,height; short planes,bitcount; long comp,sizeimg; long xpels,ypel

5、s; long used,important;/-將bmp彩色表的數(shù)據(jù)校正到bcb 的tcolor的數(shù)據(jù)。tcolor* switchcolor(unsigned char r,unsigned char g,unsigned char b) tcolor *re=new tcolor; *re=(r | g8 | boffbits,0); /-顯示一些信息 form1-edit1-text=inttostr(bmph-bfsize); form1-edit2-text=inttostr(bmpi-width); form1-edit3-text=inttostr(bmpi-height);

6、form1-edit4-text=inttostr(bmpi-comp); form1-edit5-text=inttostr(bmpi-used); int i,j,k,l,wc,pos; long n=bmph-bfsize- bmph-offbits;/象素總數(shù) unsigned char *image=new unsigned charn; /位圖矩陣 fread(image,n,1,f);/讀入位圖矩陣 int skip=0; /bmp 文件4字節(jié)對齊 if(bmpi-width%4=0) skip=0; else skip=4-bmpi-width%4; unsigned char

7、 color=0; tcolor *tc; /! /2值化 for(i=0;in;i+) if(unsigned char)imagei0xa0) imagei=(unsigned char)0; else imagei=(unsigned char)0xff; / int flag=1; long x,b; unsigned char *om=new unsigned charn; /標(biāo)記矩陣 for(i=0;in;i+) /初始化 omi=0; while(flag=1) /flag=0時迭代結(jié)束 flag=0; for( i=2;iheight-2;i+) for(j=2;jwidth-

8、2;j+) /模板a if(image(i-1)*(bmpi-width+skip)+j-1=0xff & image(i-1)*(bmpi-width+skip)+j+1=0 & image(i)*(bmpi-width+skip)+j-1=0xff & image(i)*(bmpi-width+skip)+j-0=0 & image(i)*(bmpi-width+skip)+j+1=0 & image(i)*(bmpi-width+skip)+j+2=0 & image(i+1)*(bmpi-width+skip)+j-1=0xff & image(i+1)*(bmpi-width+ski

9、p)+j+1=0) om(i)*(bmpi-width+skip)+j=0xff; flag=1; continue; /模板b if(image(i-1)*(bmpi-width+skip)+j-0=0 & image(i-1)*(bmpi-width+skip)+j+1=0 & image(i)*(bmpi-width+skip)+j-1=0xff & image(i)*(bmpi-width+skip)+j-0=0 & image(i)*(bmpi-width+skip)+j+1=0 & image(i)*(bmpi-width+skip)+j+2=0 & image(i+1)*(bmp

10、i-width+skip)+j-1=0xff & image(i+1)*(bmpi-width+skip)+j-0=0xff) om(i)*(bmpi-width+skip)+j=0xff; flag=1; continue; /模板c if(image(i-1)*(bmpi-width+skip)+j-1=0xff & image(i-1)*(bmpi-width+skip)+j-0=0xff & image(i)*(bmpi-width+skip)+j-1=0xff & image(i)*(bmpi-width+skip)+j-0=0 & image(i)*(bmpi-width+skip

11、)+j+1=0 & image(i)*(bmpi-width+skip)+j+2=0 & image(i+1)*(bmpi-width+skip)+j-0=0 & image(i+1)*(bmpi-width+skip)+j+1=0) om(i)*(bmpi-width+skip)+j=0xff; flag=1; continue; /模板d if(image(i-1)*(bmpi-width+skip)+j-1=0xff & image(i-1)*(bmpi-width+skip)+j-0=0xff & image(i-1)*(bmpi-width+skip)+j+1=0xff & imag

12、e(i)*(bmpi-width+skip)+j-0=0 & image(i+1)*(bmpi-width+skip)+j-1=0 & image(i+1)*(bmpi-width+skip)+j-0=0 & image(i+1)*(bmpi-width+skip)+j+1=0) om(i)*(bmpi-width+skip)+j=0xff; flag=1; continue; /模板e if(image(i-1)*(bmpi-width+skip)+j-1=0 & image(i-1)*(bmpi-width+skip)+j+1=0xff & image(i)*(bmpi-width+ski

13、p)+j-1=0 & image(i)*(bmpi-width+skip)+j-0=0 & image(i)*(bmpi-width+skip)+j+1=0xff & image(i+1)*(bmpi-width+skip)+j-1=0 & image(i+1)*(bmpi-width+skip)+j+1=0xff) om(i)*(bmpi-width+skip)+j=0xff; flag=1; continue; /模板f if(image(i-1)*(bmpi-width+skip)+j-0=0 & image(i)*(bmpi-width+skip)+j-1=0 & image(i)*(

14、bmpi-width+skip)+j-0=0 & image(i)*(bmpi-width+skip)+j+1=0xff & image(i+1)*(bmpi-width+skip)+j-0=0xff & image(i+1)*(bmpi-width+skip)+j+1=0xff) om(i)*(bmpi-width+skip)+j=0xff; flag=1; continue; /模板g if(image(i-1)*(bmpi-width+skip)+j-0=0xff & image(i-1)*(bmpi-width+skip)+j+1=0xff & image(i)*(bmpi-width

15、+skip)+j-1=0 & image(i)*(bmpi-width+skip)+j-0=0 & image(i)*(bmpi-width+skip)+j+1=0xff & image(i+1)*(bmpi-width+skip)+j-0=0) om(i)*(bmpi-width+skip)+j=0xff; flag=1; continue; /模板h if(image(i-2)*(bmpi-width+skip)+j-0=0 & image(i-1)*(bmpi-width+skip)+j-1=0 & image(i-1)*(bmpi-width+skip)+j-0=0 & image(i

16、-1)*(bmpi-width+skip)+j+1=0 & image(i)*(bmpi-width+skip)+j-0=0 & image(i+1)*(bmpi-width+skip)+j-1=0xff & image(i+1)*(bmpi-width+skip)+j-0=0xff & image(i+1)*(bmpi-width+skip)+j+1=0xff) om(i)*(bmpi-width+skip)+j=0xff; flag=1; continue; for(i=0;in;i+) if(omi=0xff) imagei=0xff; /-第二次串行細化 for( i=2;iheigh

17、t-2;i+) for(j=2;jwidth-2;j+) /縮小后的模板a if(image(i-1)*(bmpi-width+skip)+j-1=0xff & image(i-1)*(bmpi-width+skip)+j+1=0 & image(i)*(bmpi-width+skip)+j-1=0xff & image(i)*(bmpi-width+skip)+j-0=0 & image(i)*(bmpi-width+skip)+j+1=0 & image(i+1)*(bmpi-width+skip)+j-1=0xff & image(i+1)*(bmpi-width+skip)+j+1=0

18、) image(i)*(bmpi-width+skip)+j=0xff; flag=1; continue; /縮小后的模板h if(image(i-1)*(bmpi-width+skip)+j-1=0 & image(i-1)*(bmpi-width+skip)+j-0=0 & image(i-1)*(bmpi-width+skip)+j+1=0 & image(i)*(bmpi-width+skip)+j-0=0 & image(i+1)*(bmpi-width+skip)+j-1=0xff & image(i+1)*(bmpi-width+skip)+j-0=0xff & image(i+1)*(bmpi-width+skip)+j+1=0xff) image(i)*(bmpi-width+sk