templat.cpp

《医学图象的远程传输系统》

#include "stdafx.h"
#include "Templat.h"
/*Para:
lpBMP 指向图像的指针
Width 图像的宽度
Height 图像的高度
TempW 模板的宽度
TempH 模板的高度
TempMX 模板的中心元素x
TempMY 模板的中心元素y
fpTemp 指向模板的数组指针
Coef 模板系数

*/
#define BOUND(x,y,z) ((x)<(y))?(y):(((x)>(z))?(z):(x))//y=<x=<z 

BOOL WINAPI Templat(BYTE* lpBMP,long Width,long Height,int TempW,int TempH,int TempMX,
					int TempMY,float* fpTemp,float Coef)
{	
	float ResultR,ResultG,ResultB;
	ULONG nDibSize=Width*Height*3;
	BYTE* pDibTmp=(BYTE*)::GlobalAlloc(GMEM_FIXED,nDibSize);
	if (!pDibTmp){
		AfxMessageBox("Failed to allocate memory!");
		return FALSE;
	}
	memcpy (pDibTmp,lpBMP,nDibSize);
	for (int i=TempMY;i<Height-TempH+TempMY+1;i++){
		for (int j=TempMX;j<Width-TempW+TempMX+1;j++){
			ResultR=ResultG=ResultB=0;
			for (int k=0;k<TempH;k++){
				for (int l=0;l<TempW;l++){
					ResultR+=pDibTmp[(i-TempMY+k)*Width+(j-TempMX+l)]*fpTemp[k*TempW+l];
					ResultG+=pDibTmp[(i-TempMY+k)*Width+(j-TempMX+l)+1]*fpTemp[k*TempW+l];
					ResultB+=pDibTmp[(i-TempMY+k)*Width+(j-TempMX+l)+2]*fpTemp[k*TempW+l];
				}
			}
			pDibTmp[(i+TempMY)*Width+(j-TempMX)]=BOUND(int(Coef*ResultR),0,255);
			pDibTmp[(i+TempMY)*Width+(j-TempMX)+1]=BOUND(int(Coef*ResultG),0,255);
			pDibTmp[(i+TempMY)*Width+(j-TempMX)+2]=BOUND(int(Coef*ResultB),0,255);
			
		}
	}
	memcpy (lpBMP,pDibTmp,nDibSize);
	GlobalFree (pDibTmp);
	return TRUE;
}