📄 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 + -