?? tiduruihuadib.cpp
字號:
void TiDuRuiHuaDib::Leijianzuida(int *tongji)
{
int i,j;
int height,wide;
LPBYTE p_data;
LPBYTE p_temp;
height = this->GetHeight();
wide = this->GetDibWidthBytes();
p_data=this->GetData(); //獲得圖像數據
p_temp = new BYTE[height*wide];
memcpy(p_temp,p_data,height*wide);
//域值
int Th; //使相對距離最大的域值
int Th0; //域值初始值；
//兩類的均值；
int u1,u2;
//相對距離
int s;
int s0; //初始域值所對應的相對距離
//臨時變量
int Temp0, Temp1, Temp2, Temp3;
Temp0 = Temp1 = Temp2 = Temp3 = 0;
//初始化Th,Th0
Th=Th0=127;
//計算u1
for (i = 0; i < Th0 + 1; i++)
{
Temp0 += tongji[i] * i;
Temp1 += tongji[i];
}
u1=Temp0/Temp1;
//計算u2
for (i = Th0 + 1; i < 256; i++)
{
Temp2 += tongji[i] * i;
Temp3 += tongji[i];
}
u2=Temp2/Temp3;
//計算相對距離
s0=(u2-Th0)*(Th0-u1)/((u2-u1)*(u2-u1));
for (Th0=1;Th0<255;Th0++)
{
//計算u1
Temp0 = Temp1 = Temp2 = Temp3 = 0;
for (i = 0; i < Th0 + 1; i++)
{
Temp0 += tongji[i] * i;
Temp1 += tongji[i];
}
if (Temp1==0)
{
Temp1=1;
}
u1=Temp0/Temp1;
//計算u2
for (i = Th0 + 1; i < 256; i++)
{
Temp2 += tongji[i] * i;
Temp3 += tongji[i];
}
if (Temp3==0)
{
Temp3=1;
}
u2=Temp2/Temp3;
//計算相對距離
s=(u2-Th0)*(Th0-u1)/((u2-u1)*(u2-u1));
//與s0比較大小
if (s>s0)
{
s0=s;
Th=Th0;
}
}
//根據域值Th對圖像二值化
for (j=0;j<height;j++)
{
for (i=0;i<wide;i++)
{
if(*(p_temp+j*wide+i)<Th)
*(p_temp+j*wide+i)=0;
else
*(p_temp+j*wide+i)=255;
}
}
memcpy(p_data,p_temp,height*wide);
delete p_temp;
}
void TiDuRuiHuaDib::Yuzhifenge(int Yuzhi)
{
// 指向源圖像的指針
LPBYTE p_data;
LPBYTE lpSrc;
// 指向緩存圖像的指針
LPBYTE lpDst;
// 指向緩存DIB圖像的指針
LPBYTE temp;
//循環變量
long i;
long j;
//圖像的高和寬
long wide;
long height;
p_data=GetData();
if(m_pBitmapInfoHeader->biBitCount<9) //灰度圖像
wide=GetWidth();
else //24位彩色
wide=GetDibWidthBytes();
height=GetHeight();
temp=new BYTE[wide*height];
memset(temp,255,wide*height);
for (j=0;j<height;j++)
{
for(i=0;i<wide;i++)
{
lpSrc=p_data+wide*j+i;
lpDst=temp+wide*j+i;
if(abs(*lpSrc-Yuzhi)<30)
//*lpDst=Yuzhi;
*lpDst=0;
}
}
memcpy(p_data,temp,wide*height);
delete temp;
}
void TiDuRuiHuaDib::nnzhong(int n)
{
DWORD size;
size=GetSize(); //獲得圖像打小
BYTE* p_temp=new BYTE [size]; //開辟內存
memset(p_temp,255,size);
int yy,xx,n2,nn,chuo,chg,m,medi,madom,mado[1000];
BYTE *p_data; //原圖數據區指針
int wide,height; //原圖長、寬
if(n<3||n%2!=1)//檢查取值是否為3.5.7等的奇數
AfxMessageBox("請輸入一個大于等于3的奇數");
if(n>=3&&n%2==1)
n2=(n-1)/2;
nn=n*n;
chuo=(nn-1)/2;
p_data=this->GetData();//取得原圖的數據區指針
wide=this->GetWidth(); //取得原圖的數據區寬
height=this->GetHeight(); //取得原圖的數據區高
//n*n中值濾波
for(int j=n2;j<height-n2;j++)
for(int i=n2;i<wide-n2;i++)
{
//把n*n屏蔽窗口部分的所有像素值放入mado[m]
m=0;
for(yy=j-n2;yy<=j+n2;yy++)
for(xx=i-n2;xx<=i+n2;xx++)
{
mado[m]=p_data[yy*wide+xx];
m++;
}
//把mado[m]中的值按下降順序用冒泡法排序
do{
chg=0;
for(m=0;m<nn-1;m++)
{
if(mado[m]<mado[m+1])
{
madom=mado[m];
mado[m]=mado[m+1];
mado[m+1]=madom;
chg=1;
}
}
}while(chg==1);
//求中值medi
medi=mado[chuo];
//把中值代入顯示圖像中
p_temp[j*wide+i]=medi;
}
memcpy(p_data,p_temp,size);
delete p_temp;
}
void TiDuRuiHuaDib::nnjunzhi(int n)
{
DWORD size;
size=GetSize();
BYTE *p_data;
int xx,yy,n2,sum;
int wide,height; //原圖長、寬
BYTE* p_temp=new BYTE [size];
memset(p_temp,255,size);
if(n<3||n%2!=1)//確認n為奇數
AfxMessageBox("請輸入一個大于等于3的奇數");
if(n>=3&&n%2==1) //如果n是大于等于3的奇數
n2=(n-1)/2;
p_data=this->GetData();//取得原圖的數據區指針
wide=this->GetWidth(); //取得原圖的數據區寬
height=this->GetHeight(); //取得原圖的數據區高
//用N*N屏蔽窗口平均化
for(int j=n2;j<height-n2;j++)
{
for(int i=n2;i<wide-n2;i++)
{ //求N*N屏蔽窗口內的總和sum
sum=0;
for(yy=j-n2;yy<=j+n2;yy++)
for(xx=i-n2;xx<=i+n2;xx++)
sum+=p_data[yy*wide+xx];
//把n*n屏蔽窗口內的平均值四舍五入后作為顯示圖像像素值
p_temp[j*wide+i]=(int)((float)sum/(n*n)+0.5);
}
}
memcpy(p_data,p_temp,size);
delete p_temp;
}
void TiDuRuiHuaDib::Jubupingjun()
{
DWORD size;
size=GetSize();
BYTE *p_data;//原圖數據區指針
BYTE* p_temp=new BYTE[size]; //緩存圖像指針
int wide,height;
int n,ji[9],nmin;
float mean[9],bunsan[9],bmin;
p_data=this->GetData();//取得原圖的數據區指針
wide=this->GetWidth(); //取得原圖的數據區寬
height=this->GetHeight(); //取得原圖的數據區高
memset(p_temp,255,size); //初始化緩存區域
for(int j=2;j<=height-3;j++)
for(int i=2;i<=wide-3;i++)
{
//求9種近鄰區域的均值及其方差
//第1近鄰區域
ji[0]=p_data[(j-1)*wide+(i-1)];
ji[1]=p_data[(j-1)*wide+i];
ji[2]=p_data[(j-1)*wide+(i+1)];
ji[3]=p_data[j*wide+(i-1)];
ji[4]=p_data[j*wide+i];
ji[5]=p_data[j*wide+(i+1)];
ji[6]=p_data[(j+1)*wide+(i-1)];
ji[7]=p_data[(j+1)*wide+i];
ji[8]=p_data[(j+1)*wide+(i+1)];
mean[0]=(float)(ji[0]+ji[1]+ji[2]+ji[3]+ji[4]+ji[5]+ji[6]+ji[7]+ji[8])/9;
bunsan[0]=0;
for(n=0;n<=8;n++)
bunsan[0]+=ji[n]*ji[n]-mean[0]*mean[0];
//第2近鄰區域
ji[0]=p_data[(j-2)*wide+(i-1)];
ji[1]=p_data[(j-2)*wide+i];
ji[2]=p_data[(j-2)*wide+(i+1)];
ji[3]=p_data[(j-1)*wide+(i-1)];
ji[4]=p_data[(j-1)*wide+i];
ji[5]=p_data[(j-1)*wide+(i+1)];
ji[6]=p_data[j*wide+i];
mean[1]=(float)(ji[0]+ji[1]+ji[2]+ji[3]+ji[4]+ji[5]+ji[6])/7;
bunsan[1]=0;
for(n=0;n<=6;n++)
bunsan[1]+=ji[n]*ji[n]-mean[1]*mean[1];
//第3近鄰區域
ji[0]=p_data[(j-1)*wide+(i-2)];
ji[1]=p_data[(j-1)*wide+(i-1)];
ji[2]=p_data[j*wide+(i-2)];
ji[3]=p_data[j*wide+(i-1)];
ji[4]=p_data[j*wide+i];
ji[5]=p_data[(j+1)*wide+(i-2)];
ji[6]=p_data[(j+1)*wide+(i-1)];
mean[2]=(float)(ji[0]+ji[1]+ji[2]+ji[3]+ji[4]+ji[5]+ji[6])/7;
bunsan[2]=0;
for(n=0;n<=6;n++)
bunsan[2]+=ji[n]*ji[n]-mean[2]*mean[2];
//第4近鄰區域
ji[0]=p_data[j*wide+i];
ji[1]=p_data[(j+1)*wide+(i-1)];
ji[2]=p_data[(j+1)*wide+i];
ji[3]=p_data[(j+1)*wide+(i+1)];
ji[4]=p_data[(j+2)*wide+(i-1)];
ji[5]=p_data[(j+2)*wide+i];
ji[6]=p_data[(j+2)*wide+(i+1)];
mean[3]=(float)(ji[0]+ji[1]+ji[2]+ji[3]+ji[4]+ji[5]+ji[6])/7;
bunsan[3]=0;
for(n=0;n<=6;n++)
bunsan[3]+=ji[n]*ji[n]-mean[3]*mean[3];
//第5近鄰區域
ji[0]=p_data[(j-1)*wide+(i+1)];
ji[1]=p_data[(j-1)*wide+(i+2)];
ji[2]=p_data[j*wide+i];
ji[3]=p_data[j*wide+(i+1)];
ji[4]=p_data[j*wide+(i+2)];
ji[5]=p_data[(j+1)*wide+(i+1)];
ji[6]=p_data[(j+1)*wide+(i+2)];
mean[4]=(float)(ji[0]+ji[1]+ji[2]+ji[3]+ji[4]+ji[5]+ji[6])/7;
bunsan[4]=0;
for(n=0;n<=6;n++)
bunsan[4]+=ji[n]*ji[n]-mean[4]*mean[4];
//第6近鄰區域
ji[0]=p_data[(j-2)*wide+(i+1)];
ji[1]=p_data[(j-2)*wide+(i+2)];
ji[2]=p_data[(j-1)*wide+i];
ji[3]=p_data[(j-1)*wide+(i+1)];
ji[4]=p_data[(j-1)*wide+(i+2)];
ji[5]=p_data[j*wide+i];
ji[6]=p_data[j*wide+(i+1)];
mean[5]=(float)(ji[0]+ji[1]+ji[2]+ji[3]+ji[4]+ji[5]+ji[6])/7;
bunsan[5]=0;
for(n=0;n<=6;n++)
bunsan[5]+=ji[n]*ji[n]-mean[5]*mean[5];
//第7近鄰區域
ji[0]=m_pData[(j-2)*wide+(i-2)];
ji[1]=p_data[(j-2)*wide+(i-1)];
ji[2]=p_data[(j-1)*wide+(i-2)];
ji[3]=p_data[(j-1)*wide+(i-1)];
ji[4]=p_data[(j-1)*wide+i];
ji[5]=p_data[j*wide+(i-1)];
ji[6]=p_data[j*wide+i];
mean[6]=(float)(ji[0]+ji[1]+ji[2]+ji[3]+ji[4]+ji[5]+ji[6])/7;
bunsan[6]=0;
for(n=0;n<=6;n++)
bunsan[6]+=ji[n]*ji[n]-mean[6]*mean[6];
//第8近鄰區域
ji[0]=m_pData[j*wide+(i-1)];
ji[1]=p_data[j*wide+i];
ji[2]=p_data[(j+1)*wide+(i-2)];
ji[3]=p_data[(j+1)*wide+(i-1)];
ji[4]=p_data[(j+1)*wide+i];
ji[5]=p_data[(j+2)*wide+(i-2)];
ji[6]=p_data[(j+2)*wide+(i-1)];
mean[7]=(float)(ji[0]+ji[1]+ji[2]+ji[3]+ji[4]+ji[5]+ji[6])/7;
bunsan[7]=0;
for(n=0;n<=6;n++)
bunsan[7]+=ji[n]*ji[n]-mean[7]*mean[7];
//第9近鄰區域
ji[0]=p_data[j*wide+i];
ji[1]=p_data[j*wide+(i+1)];
ji[2]=p_data[(j+1)*wide+i];
ji[3]=p_data[(j+1)*wide+(i+1)];
ji[4]=p_data[(j+1)*wide+(i+2)];
ji[5]=p_data[(j+2)*wide+(i+1)];
ji[6]=p_data[(j+2)*wide+(i+2)];
mean[8]=(float)(ji[0]+ji[1]+ji[2]+ji[3]+ji[4]+ji[5]+ji[6])/7;
bunsan[8]=0;
for(n=0;n<=6;n++)
bunsan[8]+=ji[n]*ji[n]-mean[8]*mean[8];
//求方差最小的近鄰區域nmin
bmin=bunsan[0];
nmin=0;
for(n=0;n<=8;n++)
{
if(bmin>bunsan[n])
{
bmin=bunsan[n];
nmin=n;
}
//把nmin的值四舍五入后作為顯示圖像的值
p_temp[j*wide+i]=(int)(mean[nmin]+0.5);
}
}
memcpy(p_data,p_temp,size);
delete p_temp;
}
void TiDuRuiHuaDib::Templat(BYTE *m_pdata, int wide, int height, int tempH, int tempW, int tempMX, int tempMY, float *fpArray, float fCoef)
{
int i,j,k,l; //循環變量
BYTE* temp=new BYTE[wide*height]; //新圖像緩沖區
//初始化新圖像為原始圖像
memcpy( temp,m_pdata,wide*height);
float fResult; //像素值計算結果
for(j=tempMY;j<height-tempH+tempMY+1;j++)
for(i=tempMX;i<wide-tempW+tempMX+1;i++)
{
//計算像素值
fResult=0;
for(k=0;k<tempH;k++)
for(l=0;l<tempW;l++)
fResult=fResult+m_pdata[(j-tempMY+k)*wide+(i-tempMX+l)]*fpArray[k*tempW+l];
//乘上系數
fResult*=fCoef;
//取絕對值
fResult=(float)fabs(fResult);
//判斷是否超過255
if(fResult>255)
//若超過255，直接賦值為255
temp[j*wide+i]=255;
else
//未超過255，賦值為計算結果
temp[j*wide+i]=(int)(fResult+0.5);
}
memcpy(m_pdata, temp,wide*height); //復制處理后的圖像
delete temp;
}
void TiDuRuiHuaDib::Laplacian()
{
int tempH; //模板高度
int tempW; //模板寬度
float tempC; //模板系數
int tempMY; //模板中心元素Y坐標
int tempMX; //模板中心元素X坐標
float Template[9]; //模板數組
LPBYTE p_data; //原圖數據區指針
int wide,height; //原圖長、寬
p_data=this->GetData ();
wide=this->GetWidth ();
height=this->GetHeight ();
LPBYTE temp1=new BYTE[wide*height]; //新圖像緩沖區
//拷貝原圖像到緩存圖像
memcpy( temp1,p_data,wide*height);
//設置Laplacian模板 參數
tempW=3;
tempH=3;
tempC=1.0;
tempMY=1;
tempMX=1;
Template[0]=-1.0;
Template[1]=-1.0;
Template[2]=-1.0;
Template[3]=-1.0;
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -