图像高提升滤波的vc实现
void CSDIELSView::OnSharpeningEnhance(){//程序编制：李立宗 lilizong@gmail.com//2012-8-12if(myImage1.IsNull())OnOpenResourceFile();if(!myImage2.IsNull())myImage2.Destroy();if(myImage2.IsNull()){myImage2.Create(myImage1.GetWidth(),myImage1.GetHeight(),24,0);}//COLORREF pixel; int maxY = myImage1.GetHeight();int maxX=myImage1.GetWidth();byte* pRealData;byte* pRealData2;pRealData=(byte*)myImage1.GetBits();pRealData2=(byte*)myImage2.GetBits();int pit=myImage1.GetPitch();int pit2=myImage2.GetPitch();//需要注意,pit和pit2的值并不一样,所以如果使用一个值,会导致不同的结果出现//CString str;//str.Format(TEXT("%d"),pit);//MessageBox(str);//str.Format(TEXT("%d"),pit2);//MessageBox(str);int bitCount=myImage1.GetBPP()/8;int bitCount2=myImage2.GetBPP()/8;int tempR,tempG,tempB;int tempX,tempY;int temp;float tempA;//int M[3][3]={{1,2,1},{2,4,2},{1,2,1}};int t=100; //门限int f=1; //参数，原图像的比例值float maxR,minR,maxG,minG,maxB,minB;float spanR,spanG,spanB;static int tempArray[1000][1000][3]={0}; //必须定义为static，否则无法支持大数组//tempR=tempG=tempG=0;//说明：将生产的图像作为24位图处理。for (int y=1; y<maxY-1; y++) {for (int x=1; x<maxX-1; x++) {temp=*(pRealData+pit*(y-1)+(x-1)*bitCount)+*(pRealData+pit*(y-1)+(x)*bitCount)+*(pRealData+pit*(y-1)+(x+1)*bitCount)+*(pRealData+pit*(y)+(x-1)*bitCount)-8*(*(pRealData+pit*(y)+(x)*bitCount))+*(pRealData+pit*(y)+(x+1)*bitCount)+*(pRealData+pit*(y+1)+(x-1)*bitCount)+*(pRealData+pit*(y+1)+(x)*bitCount)+*(pRealData+pit*(y+1)+(x+1)*bitCount);if(temp>255)temp=255;if(temp<-255)temp=-255;tempArray[y][x][0]=f*(*(pRealData+pit*(y)+(x)*bitCount))+temp+0.5;if(bitCount==1){/*tempG=tempR;tempB=tempR;*/tempArray[y][x][1]=tempArray[y][x][0];tempArray[y][x][2]=tempArray[y][x][0];}else{temp=*(pRealData+pit*(y-1)+(x-1)*bitCount+1)+*(pRealData+pit*(y-1)+(x)*bitCount+1)+*(pRealData+pit*(y-1)+(x+1)*bitCount+1)+*(pRealData+pit*(y)+(x-1)*bitCount+1)-8*(*(pRealData+pit*(y)+(x)*bitCount+1))+*(pRealData+pit*(y)+(x+1)*bitCount+1)+*(pRealData+pit*(y+1)+(x-1)*bitCount+1)+*(pRealData+pit*(y+1)+(x)*bitCount+1)+*(pRealData+pit*(y+1)+(x+1)*bitCount+1);if(temp>255)temp=255;if(temp<-255)temp=-255;tempArray[y][x][1]=f*(*(pRealData+pit*(y)+(x)*bitCount+1))+temp+0.5;temp=*(pRealData+pit*(y-1)+(x-1)*bitCount+2)+*(pRealData+pit*(y-1)+(x)*bitCount+2)+*(pRealData+pit*(y-1)+(x+1)*bitCount+2)+*(pRealData+pit*(y)+(x-1)*bitCount+2)-8*(*(pRealData+pit*(y)+(x)*bitCount+2))+*(pRealData+pit*(y)+(x+1)*bitCount+2)+*(pRealData+pit*(y+1)+(x-1)*bitCount+2)+*(pRealData+pit*(y+1)+(x)*bitCount+2)+*(pRealData+pit*(y+1)+(x+1)*bitCount+2);if(temp>255)temp=255;if(temp<-255)temp=-255;tempArray[y][x][2]=f*(*(pRealData+pit*(y)+(x)*bitCount+2))+temp+0.5;}}}maxR=maxG=maxB=0;minR=minG=minB=65535;for (int y=1; y<maxY-1; y++) {for (int x=1; x<maxX-1; x++) {if(maxR<(tempArray[y][x][0]))maxR=tempArray[y][x][0];if(maxG<(tempArray[y][x][1]))maxG=tempArray[y][x][1];if(maxB<(tempArray[y][x][2]))maxB=tempArray[y][x][2];if(minR>(tempArray[y][x][0]))minR=tempArray[y][x][0];if(minG>(tempArray[y][x][1]))minG=tempArray[y][x][1];if(minB>(tempArray[y][x][2]))minB=tempArray[y][x][2];}}spanR=maxR-minR;spanG=maxG-minG;spanB=maxB-minB;for (int y=1; y<maxY-1; y++) {for (int x=1; x<maxX-1; x++) {if(spanR>0)tempR=(int)(tempArray[y][x][0]-minR)*255/spanR;else if(tempArray[y][x][0]<=255)tempR=tempArray[y][x][0];elsetempR=255;if(spanG>0)tempG=(int)(tempArray[y][x][1]-minR)*255/spanR;else if(tempArray[y][x][1]<=255)tempG=tempArray[y][x][1];elsetempG=255;if(spanB>0)tempB=(int)(tempArray[y][x][2]-minR)*255/spanR;else if(tempArray[y][x][2]<=255)tempB=tempArray[y][x][2];elsetempB=255;*(pRealData2+pit2*y+x*bitCount2)=tempR;*(pRealData2+pit2*y+x*bitCount2+1)=tempG;*(pRealData2+pit2*y+x*bitCount2+2)=tempB;}}Invalidate();}