?? omlib.cpp
字號:
unsigned char **aveImage; /*圖像的固定尺寸平均值*/
unsigned char **smoothImage1; /*前一次迭代平滑的結果*/
unsigned char **smoothImage2; /*本次迭代平滑的結果*/
imageHeight=Row;
imageWidth=Col;
aveImage=(unsigned char**)fspace_2d(imageHeight,imageWidth);
smoothImage1=(unsigned char**)fspace_2d(imageHeight,imageWidth);
smoothImage2=(unsigned char**)fspace_2d(imageHeight,imageWidth);
/*start to iterately smoothing*/
for(i=0; i<imageHeight; i++)
for(j=0; j<imageWidth; j++) smoothImage1[i][j] = OrigImg[i][j];
a=0; /*標記循環是否結束？*/
while(a==0)
{
/*求取圖像的固定尺寸平均值*/
for(i=n;i<imageHeight-n;i++)
for(j=n;j<imageWidth-n;j++)
{
s=0;
for(k=i-n;k<=i+n;k++)
for(m=j-n;m<=j+n;m++)
s=s+smoothImage1[k][m];
aveImage[i][j]=s/((2*n+1)*(2*n+1));
}
for(i=0;i<n;i++)
for(j=0;j<imageWidth;j++)
aveImage[i][j]=aveImage[i+n][j];
for(i=imageHeight-n;i<imageHeight;i++)
for(j=0;j<imageWidth;j++)
aveImage[i][j]=aveImage[i-n][j];
for(j=0;j<n;j++)
for(i=0;i<imageHeight;i++)
aveImage[i][j]=aveImage[i][j+n];
for(j=imageWidth-n;j<imageWidth;j++)
for(i=0;i<imageHeight;i++)
aveImage[i][j]=aveImage[i][j-n];
/*對每一點，尋找五個區域中均勻性最好的區域，并以此區域的平均值作為平滑值*/
for(i=2*n;i<imageHeight-2*n;i++)
for(j=2*n;j<imageWidth-2*n;j++)
{
/*求五個區域的均勻性*/
h0=abs(aveImage[i+n][j-n]+aveImage[i+n][j+n]-aveImage[i-n][j-n]-aveImage[i-n][j+n])+
abs(aveImage[i-n][j+n]+aveImage[i+n][j+n]-aveImage[i-n][j-n]-aveImage[i+n][j-n]);
h1=abs(aveImage[i][j-2*n]+aveImage[i][j]-aveImage[i-2*n][j-2*n]-aveImage[i-2*n][j])+
abs(aveImage[i-2*n][j]+aveImage[i][j]-aveImage[i-2*n][j-2*n]-aveImage[i][j-2*n]);
h2=abs(aveImage[i][j]+aveImage[i][j+2*n]-aveImage[i-2*n][j]-aveImage[i-2*n][j+2*n])+
abs(aveImage[i-2*n][j+2*n]+aveImage[i][j+2*n]-aveImage[i-2*n][j]-aveImage[i][j]);
h3=abs(aveImage[i+2*n][j-2*n]+aveImage[i+2*n][j]-aveImage[i][j-2*n]-aveImage[i][j])+
abs(aveImage[i][j]+aveImage[i+2*n][j]-aveImage[i][j-2*n]-aveImage[i+2*n][j-2*n]);
h4=abs(aveImage[i+2*n][j]+aveImage[i+2*n][j+2*n]-aveImage[i][j]-aveImage[i][j+2*n])+
abs(aveImage[i][j+2*n]+aveImage[i+2*n][j+2*n]-aveImage[i][j]-aveImage[i+2*n][j]);
/*尋找五個區域中均勻性最好的區域，并以此區域的平均值作為平滑值*/
if((h0<=h1)&&(h0<=h2)&&(h0<=h3)&&(h0<=h4))
smoothImage2[i][j]=aveImage[i][j];
else if((h1<=h0)&&(h1<=h2)&&(h1<=h3)&&(h1<=h4))
smoothImage2[i][j]=aveImage[i-n][j-n];
else if((h2<=h0)&&(h2<=h1)&&(h2<=h3)&&(h2<=h4))
smoothImage2[i][j]=aveImage[i-n][j+n];
else if((h3<=h0)&&(h3<=h1)&&(h3<=h2)&&(h3<=h4))
smoothImage2[i][j]=aveImage[i+n][j-n];
else if((h4<=h0)&&(h4<=h1)&&(h4<=h2)&&(h4<=h3))
smoothImage2[i][j]=aveImage[i+n][j+n];
}
for(i=0;i<2*n;i++)
for(j=0;j<imageWidth;j++)
smoothImage2[i][j]=smoothImage2[2*n][j];
for(i=imageHeight-2*n;i<imageHeight;i++)
for(j=0;j<imageWidth;j++)
smoothImage2[i][j]=smoothImage2[imageHeight-2*n-1][j];
for(i=0;i<imageHeight;i++)
for(j=0;j<2*n;j++)
smoothImage2[i][j]=smoothImage2[i][2*n];
for(i=0;i<imageHeight;i++)
for(j=imageWidth-2*n;j<imageWidth;j++)
smoothImage2[i][j]=smoothImage2[i][imageWidth-2*n-1];
x++;
if(x==5) a=1; /*迭代5次*/
else
{
for(i=0;i<imageHeight;i++)
for(j=0;j<imageWidth;j++)
smoothImage1[i][j]=smoothImage2[i][j];
}
} /*End of while*/
/*saving the smoothed data */
for(i=0;i<imageHeight;i++)
for(j=0;j<imageWidth;j++)
SmoothImg[i][j]=smoothImage2[i][j];
/*free the unused memmory*/
dspace_2d(aveImage,imageHeight,imageWidth);
dspace_2d(smoothImage1,imageHeight,imageWidth);
dspace_2d(smoothImage2,imageHeight,imageWidth);
}
/*************************************************************************
*
* 函數名稱：
* ErosiontionDIB()
*
* 參數:
* LPSTR lpDIBBits - 指向源DIB圖像指針
* LONG lWidth - 源圖像寬度（象素數，必須是4的倍數）
* LONG lHeight - 源圖像高度（象素數）
* int nMode - 腐蝕方式，0表示水平方向，1表示垂直方向，2表示自定義結構元素。
* int structure[3][3]
- 自定義的3×3結構元素。
*
* 返回值:
* BOOL - 腐蝕成功返回TRUE，否則返回FALSE。
*
* 說明:
* 該函數用于對圖像進行腐蝕運算。結構元素為水平方向或垂直方向的三個點，中間點位于原點；
* 或者由用戶自己定義3×3的結構元素。
*
* 要求目標圖像為只有0和255兩個灰度值的灰度圖像。
************************************************************************/
BOOL WINAPI ErosionDIB(LPSTR lpDIBBits, LONG lWidth, LONG lHeight, int nMode , int structure[3][3])
{
// 指向源圖像的指針
LPSTR lpSrc;
// 指向緩存圖像的指針
LPSTR lpDst;
// 指向緩存DIB圖像的指針
LPSTR lpNewDIBBits;
HLOCAL hNewDIBBits;
//循環變量
long i;
long j;
int n;
int m;
//像素值
unsigned char pixel;
// 暫時分配內存，以保存新圖像
hNewDIBBits = LocalAlloc(LHND, lWidth * lHeight);
if (hNewDIBBits == NULL)
{
// 分配內存失敗
return FALSE;
}
// 鎖定內存
lpNewDIBBits = (char * )LocalLock(hNewDIBBits);
// 初始化新分配的內存，設定初始值為255
lpDst = (char *)lpNewDIBBits;
memset(lpDst, (BYTE)255, lWidth * lHeight);
if(nMode == 0)
{
//使用水平方向的結構元素進行腐蝕
for(j = 0; j <lHeight; j++)
{
for(i = 1;i <lWidth-1; i++)
{
//由于使用1×3的結構元素，為防止越界，所以不處理最左邊和最右邊的兩列像素
// 指向源圖像倒數第j行，第i個象素的指針
lpSrc = (char *)lpDIBBits + lWidth * j + i;
// 指向目標圖像倒數第j行，第i個象素的指針
lpDst = (char *)lpNewDIBBits + lWidth * j + i;
//取得當前指針處的像素值，注意要轉換為unsigned char型
pixel = (unsigned char)*lpSrc;
//目標圖像中含有0和255外的其它灰度值
if(pixel != 255 && *lpSrc != 0)
return FALSE;
//目標圖像中的當前點先賦成黑色
*lpDst = (unsigned char)0;
//如果源圖像中當前點自身或者左右有一個點不是黑色，
//則將目標圖像中的當前點賦成白色
for (n = 0;n < 3;n++ )
{
pixel = *(lpSrc+n-1);
if (pixel == 255 )
{
*lpDst = (unsigned char)255;
break;
}
}
}
}
}
else if(nMode == 1)
{
//使用垂直方向的結構元素進行腐蝕
for(j = 1; j <lHeight-1; j++)
{
for(i = 0;i <lWidth; i++)
{
//由于使用1×3的結構元素，為防止越界，所以不處理最上邊和最下邊的兩列像素
// 指向源圖像倒數第j行，第i個象素的指針
lpSrc = (char *)lpDIBBits + lWidth * j + i;
// 指向目標圖像倒數第j行，第i個象素的指針
lpDst = (char *)lpNewDIBBits + lWidth * j + i;
//取得當前指針處的像素值，注意要轉換為unsigned char型
pixel = (unsigned char)*lpSrc;
//目標圖像中含有0和255外的其它灰度值
if(pixel != 255 && *lpSrc != 0)
return FALSE;
//目標圖像中的當前點先賦成黑色
*lpDst = (unsigned char)0;
//如果源圖像中當前點自身或者上下有一個點不是黑色，
//則將目標圖像中的當前點賦成白色
for (n = 0;n < 3;n++ )
{
pixel = *(lpSrc+(n-1)*lWidth);
if (pixel == 255 )
{
*lpDst = (unsigned char)255;
break;
}
}
}
}
}
else
{
//使用自定義的結構元素進行腐蝕
for(j = 1; j <lHeight-1; j++)
{
for(i = 0;i <lWidth; i++)
{
//由于使用3×3的結構元素，為防止越界，所以不處理最左邊和最右邊的兩列像素
//和最上邊和最下邊的兩列像素
// 指向源圖像倒數第j行，第i個象素的指針
lpSrc = (char *)lpDIBBits + lWidth * j + i;
// 指向目標圖像倒數第j行，第i個象素的指針
lpDst = (char *)lpNewDIBBits + lWidth * j + i;
//取得當前指針處的像素值，注意要轉換為unsigned char型
pixel = (unsigned char)*lpSrc;
//目標圖像中含有0和255外的其它灰度值
if(pixel != 255 && *lpSrc != 0)
return FALSE;
//目標圖像中的當前點先賦成黑色
*lpDst = (unsigned char)0;
//如果原圖像中對應結構元素中為黑色的那些點中有一個不是黑色，
//則將目標圖像中的當前點賦成白色
//注意在DIB圖像中內容是上下倒置的
for (m = 0;m < 3;m++ )
{
for (n = 0;n < 3;n++)
{
if( structure[m][n] == -1)
continue;
pixel = *(lpSrc + ((2-m)-1)*lWidth + (n-1));
if (pixel == 255 )
{
*lpDst = (unsigned char)255;
break;
}
}
}
}
}
}
// 復制腐蝕后的圖像
memcpy(lpDIBBits, lpNewDIBBits, lWidth * lHeight);
// 釋放內存
LocalUnlock(hNewDIBBits);
LocalFree(hNewDIBBits);
// 返回
return TRUE;
}
/*************************************************************************
*
* 函數名稱：
* DilationDIB()
*
* 參數:
* LPSTR lpDIBBits - 指向源DIB圖像指針
* LONG lWidth - 源圖像寬度（象素數，必須是4的倍數）
* LONG lHeight - 源圖像高度（象素數）
* int nMode - 膨脹方式，0表示水平方向，1表示垂直方向，2表示自定義結構元素。
* int structure[3][3]
- 自定義的3×3結構元素。
*
* 返回值:
* BOOL - 膨脹成功返回TRUE，否則返回FALSE。
*
* 說明:
* 該函數用于對圖像進行膨脹運算。結構元素為水平方向或垂直方向的三個點，中間點位于原點；
* 或者由用戶自己定義3×3的結構元素。
*
* 要求目標圖像為只有0和255兩個灰度值的灰度圖像。
************************************************************************/
BOOL WINAPI DilationDIB(LPSTR lpDIBBits, LONG lWidth, LONG lHeight, int nMode , int structure[3][3])
{
// 指向源圖像的指針
LPSTR lpSrc;
// 指向緩存圖像的指針
LPSTR lpDst;
// 指向緩存DIB圖像的指針
LPSTR lpNewDIBBits;
HLOCAL hNewDIBBits;
//循環變量
long i;
long j;
int n;
int m;
//像素值
unsigned char pixel;
// 暫時分配內存，以保存新圖像
hNewDIBBits = LocalAlloc(LHND, lWidth * lHeight);
if (hNewDIBBits == NULL)
{
// 分配內存失敗
return FALSE;
}
// 鎖定內存
lpNewDIBBits = (char * )LocalLock(hNewDIBBits);
// 初始化新分配的內存，設定初始值為255
lpDst = (char *)lpNewDIBBits;
memset(lpDst, (BYTE)255, lWidth * lHeight);
if(nMode == 0)
{
//使用水平方向的結構元素進行膨脹
for(j = 0; j <lHeight; j++)
{
for(i = 1;i <lWidth-1; i++)
{
//由于使用1×3的結構元素，為防止越界，所以不處理最左邊和最右邊的兩列像素
// 指向源圖像倒數第j行，第i個象素的指針
lpSrc = (char *)lpDIBBits + lWidth * j + i;
// 指向目標圖像倒數第j行，第i個象素的指針
lpDst = (char *)lpNewDIBBits + lWidth * j + i;
//取得當前指針處的像素值，注意要轉換為unsigned char型
pixel = (unsigned char)*lpSrc;
//目標圖像中含有0和255外的其它灰度值
if(pixel != 255 && pixel != 0)
return FALSE;
//目標圖像中的當前點先賦成白色
*lpDst = (unsigned char)255;
//源圖像中當前點自身或者左右只要有一個點是黑色，
//則將目標圖像中的當前點賦成黑色
for (n = 0;n < 3;n++ )
{
pixel = *(lpSrc+n-1);
if (pixel == 0 )
{
*lpDst = (unsigned char)0;
break;
}
}
}
}
}
else if(nMode == 1)
{
//使用垂直方向的結構元素進行膨脹
for(j = 1; j <lHeight-1; j++)
{
for(i = 0;i <lWidth; i++)
{
//由于使用1×3的結構元素，為防止越界，所以不處理最上邊和最下邊的兩列像素
// 指向源圖像倒數第j行，第i個象素的指針
lpSrc = (char *)lpDIBBits + lWidth * j + i;
// 指向目標圖像倒數第j行，第i個象素的指針
lpDst = (char *)lpNewDIBBits + lWidth * j + i;
//取得當前指針處的像素值，注意要轉換為unsigned char型
pixel = (unsigned char)*lpSrc;
//目標圖像中含有0和255外的其它灰度值
if(pixel != 255 && *lpSrc != 0)
return FALSE;
//目標圖像中的當前點先賦成白色
*lpDst = (unsigned char)255;
//源圖像中當前點自身或者上下只要有一個點是黑色，
//則將目標圖像中的當前點賦成黑色
for (n = 0;n < 3;n++ )
{
pixel = *(lpSrc+(n-1)*lWidth);
if (pixel == 0 )
{
*lpDst = (unsigned char)0;
break;
}
}
}
}
}
else
{
//使用自定義的結構元素進行膨脹
for(j = 1; j <lHeight-1; j++)
{
for(i = 0;i <lWidth; i++)
{
//由于使用3×3的結構元素，為防止越界，所以不處理最左邊和最右邊的兩列像素
//和最上邊和最下邊的兩列像素
// 指向源圖像倒數第j行，第i個象素的指針
lpSrc = (char *)lpDIBBits + lWidth * j + i;
// 指向目標圖像倒數第j行，第i個象素的指針
lpDst = (char *)lpNewDIBBits + lWidth * j + i;
//取得當前指針處的像素值，注意要轉換為unsigned char型
pixel = (unsigned char)*lpSrc;
//目標圖像中含有0和255外的其它灰度值
if(pixel != 255 && *lpSrc != 0)
return FALSE;
//目標圖像中的當前點先賦成白色
*lpDst = (unsigned char)255;
//原圖像中對應結構元素中為黑色的那些點中只要有一個是黑色，
//則將目標圖像中的當前點賦成黑色
//注意在DIB圖像中內容是上下倒置的
for (m = 0;m < 3;m++ )
{
for (n = 0;n < 3;n++)
{
if( structure[m][n] == -1)
continue;
pixel = *(lpSrc + ((2-m)-1)*lWidth + (n-1));
if (pixel == 0 )
{
*lpDst = (unsigned char)0;
break;
}
}
}
}
}
}
// 復制膨脹后的圖像
memcpy(lpDIBBits, lpNewDIBBits, lWidth * lHeight);
// 釋放內存
LocalUnlock(hNewDIBBits);
LocalFree(hNewDIBBits);
// 返回
return TRUE;
}
/*************************************************************************
*
* 函數名稱?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -