📄 filter.cpp

📁 一些自己做的关于图象处理的程序（如傅立叶变换
💻 CPP
📖 第 1 页 / 共 2 页
字号:
12 下一页
// Filter.cpp: implementation of the CFilter class.
//
//////////////////////////////////////////////////////////////////////
/**********************************************************************
* Copyright (c) 2003, Medical Image Processing Lab, Sichuan University 
* All rights reserved.
* 
* 文件名称：Equalize.cpp
* 摘    要：图像中值滤波CFilter的源文件
* 
* 当前版本：1.1 
* 作    者：jian wei zhang 
* 完成日期：2004年10月30日
************************************************************************/

#include "stdafx.h"
#include "photostar.h"
#include "Filter.h"

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////



CFilter::~CFilter()
{

}

CFilter::CFilter(CImageObject *pImageObject, int iWinSize)
{
	ASSERT(pImageObject!=NULL);

	m_pImageObject=pImageObject;
	m_iWinSize=iWinSize;

}

/**********************************************************************
* Copyright (c) 2003, Medical Image Processing Lab, Sichuan University 
* All rights reserved.
* 
* 文件名称：Equalize.cpp
* 摘    要：图像中值滤波CFilter的源文件
* 
* 当前版本：1.1 
* 作    者：jian wei zhang 
* 完成日期：2004年10月30日
************************************************************************/


BOOL CFilter::MedianFilter()
{
	if (m_pImageObject == NULL)
	{
		return FALSE;
	}
	
	int iWidth = m_pImageObject->GetWidth();
	int iHeight = m_pImageObject->GetHeight();
	
	/******************************the old image*****************************/
	unsigned char *pOldBuffer, *pNewBuffer, *pOldBits, *pNewBits;
	BITMAPFILEHEADER *pOldBFH, *pNewBFH;
	BITMAPINFOHEADER *pOldBIH, *pNewBIH;
	RGBQUAD *pOldPalette, *pNewPalette;
	int nWidthBytes, nNumColors;
	
	pOldBuffer = (unsigned char *) m_pImageObject->GetDIBPointer(&nWidthBytes, m_pImageObject->GetNumBits());
	if (pOldBuffer == NULL) 
	{
		return FALSE;
	}
	
	pOldBFH = (BITMAPFILEHEADER *) pOldBuffer;
	pOldBIH = (BITMAPINFOHEADER *) &pOldBuffer[sizeof(BITMAPFILEHEADER)];
	nNumColors = m_pImageObject->GetNumColors();
	pOldPalette = (RGBQUAD *)  &pOldBuffer[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)];
	pOldBits = (unsigned char *) &pOldBuffer[sizeof(BITMAPFILEHEADER)
		+sizeof(BITMAPINFOHEADER)+nNumColors*sizeof(RGBQUAD)];
	
	/*****************************the new image*****************************/
	DWORD dwNewSize;
	HGLOBAL hNewDib;
	
	//Allocate a global memory block for the new image
	dwNewSize = sizeof(BITMAPFILEHEADER) + sizeof(BITMAPINFOHEADER) 
		+ nNumColors * sizeof(RGBQUAD)
		+ nWidthBytes * (m_pImageObject->GetHeight());
	hNewDib = ::GlobalAlloc(GMEM_MOVEABLE | GMEM_ZEROINIT, dwNewSize);
	if (hNewDib == NULL)
	{
		m_pImageObject->m_nLastError = IMAGELIB_MEMORY_ALLOCATION_ERROR;
		::GlobalUnlock(m_pImageObject->GetDib());
		return FALSE;
	}
	
	//Get the pointer to the new memory block
	pNewBuffer = (unsigned char *) ::GlobalLock(hNewDib);
	if (pNewBuffer == NULL)
	{
		::GlobalFree( hNewDib );
		m_pImageObject->m_nLastError = IMAGELIB_MEMORY_LOCK_ERROR;
		::GlobalUnlock(m_pImageObject->GetDib());
		return FALSE;
	}
	
	pNewBFH = (BITMAPFILEHEADER *) pNewBuffer;
	pNewBIH = (BITMAPINFOHEADER *) &pNewBuffer[sizeof(BITMAPFILEHEADER)];
	pNewPalette = (RGBQUAD *) &pNewBuffer[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)];
	pNewBits = (unsigned char *) &pNewBuffer[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+nNumColors*sizeof(RGBQUAD)];
	
	//Make a copy of the old image
	memcpy(pNewBuffer,pOldBuffer,dwNewSize);
	
	/*******************************Median filter****************************/
	//define a template of the median filter
	int iAmount = m_iWinSize*m_iWinSize;
	int* piTemplate = new int[iAmount];
	for (int k=0; k<iAmount; k++)
	{
		piTemplate[k] = 0;
	}
	
	//the original processing point
	int iX = m_iWinSize / 2;
	int iY = m_iWinSize / 2;
	int iMedIndex = iAmount / 2;
	
	//begin to process
	int i, j, ii, jj;
	long a=m_pImageObject->GetNumBits();	
	switch (m_pImageObject->GetNumBits())
	{
	case 8:
		//image loop
		for (j=iY; j<iHeight-iY; j++)
		{
			for (i=iX; i<iWidth-iX; i++)
			{
				//template loop
				for (jj=0; jj<m_iWinSize; jj++)
				{
					for (ii=0; ii<m_iWinSize; ii++)
					{
						//Get gray value of the pixel in line (j-iY+jj) and column (i-iX+ii)
						piTemplate[jj*m_iWinSize+ii] = 
							pOldBits[(iHeight-1-(j-iY+jj))*nWidthBytes + (i-iX+ii)];
						
					}
				}
				
				//Sort the elements in the template
				BubbleSort(piTemplate, iAmount);
				//Get the median value in the template
				pNewBits[(iHeight-1-j)*nWidthBytes + i] = piTemplate[iMedIndex];
				 
			}
		}
		break;
	case 16:
		break;
	case 24:
		break;
	case 32:
		break;
	}
	
	//Free the temporary memory
	delete [] piTemplate;
	piTemplate = NULL;
	
	//Free the memory for the old image
	::GlobalUnlock(m_pImageObject->GetDib());
	::GlobalFree(m_pImageObject->GetDib());
	
	//Update the DIB in the object pointed by m_pImaeObject
	::GlobalUnlock(hNewDib); 
	m_pImageObject->SetDib(hNewDib);
	
	return TRUE;
	
}


/**********************************************************************
* Copyright (c) 2003, Medical Image Processing Lab, Sichuan University 
* All rights reserved.
* 
* 文件名称：Equalize.cpp
* 摘    要：图像最大值滤波CFilter的源文件
* 
* 当前版本：1.1 
* 作    者：jian wei zhang 
* 完成日期：2004年11月8日
************************************************************************/

BOOL CFilter::MaxFilter()
{
  	if (m_pImageObject == NULL)
	{
		return FALSE;
	}
	
	int iWidth = m_pImageObject->GetWidth();
	int iHeight = m_pImageObject->GetHeight();
	int iNumBits=m_pImageObject->GetNumBits();
	/******************************the old image*****************************/
	unsigned char *pOldBuffer, *pNewBuffer, *pOldBits, *pNewBits;
	BITMAPFILEHEADER *pOldBFH, *pNewBFH;
	BITMAPINFOHEADER *pOldBIH, *pNewBIH;
	RGBQUAD *pOldPalette, *pNewPalette;
	int nWidthBytes, nNumColors;
	
	pOldBuffer = (unsigned char *) m_pImageObject->GetDIBPointer(&nWidthBytes, m_pImageObject->GetNumBits());
	if (pOldBuffer == NULL) 
	{
		return FALSE;
	}
	
	pOldBFH = (BITMAPFILEHEADER *) pOldBuffer;
	pOldBIH = (BITMAPINFOHEADER *) &pOldBuffer[sizeof(BITMAPFILEHEADER)];
	nNumColors = m_pImageObject->GetNumColors();
	pOldPalette = (RGBQUAD *)  &pOldBuffer[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)];
	pOldBits = (unsigned char *) &pOldBuffer[sizeof(BITMAPFILEHEADER)
		+sizeof(BITMAPINFOHEADER)+nNumColors*sizeof(RGBQUAD)];
	
	/*****************************the new image*****************************/
	DWORD dwNewSize;
	HGLOBAL hNewDib;
	
	//Allocate a global memory block for the new image
	dwNewSize = sizeof(BITMAPFILEHEADER) + sizeof(BITMAPINFOHEADER) 
		+ nNumColors * sizeof(RGBQUAD)
		+ nWidthBytes * (m_pImageObject->GetHeight());
	hNewDib = ::GlobalAlloc(GMEM_MOVEABLE | GMEM_ZEROINIT, dwNewSize);
	if (hNewDib == NULL)
	{
		m_pImageObject->m_nLastError = IMAGELIB_MEMORY_ALLOCATION_ERROR;
		::GlobalUnlock(m_pImageObject->GetDib());
		return FALSE;
	}
	
	//Get the pointer to the new memory block
	pNewBuffer = (unsigned char *) ::GlobalLock(hNewDib);
	if (pNewBuffer == NULL)
	{
		::GlobalFree( hNewDib );
		m_pImageObject->m_nLastError = IMAGELIB_MEMORY_LOCK_ERROR;
		::GlobalUnlock(m_pImageObject->GetDib());
		return FALSE;
	}
	
	pNewBFH = (BITMAPFILEHEADER *) pNewBuffer;
	pNewBIH = (BITMAPINFOHEADER *) &pNewBuffer[sizeof(BITMAPFILEHEADER)];
	pNewPalette = (RGBQUAD *) &pNewBuffer[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)];
	pNewBits = (unsigned char *) &pNewBuffer[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+nNumColors*sizeof(RGBQUAD)];
	
	//Make a copy of the old image
	memcpy(pNewBuffer,pOldBuffer,dwNewSize);
	
	/*******************************Max filter****************************/
	//define a template of the median filter
	int iAmount = m_iWinSize*m_iWinSize;
	int* piTemplate = new int[iAmount];
	for (int k=0; k<iAmount; k++)
	{
		piTemplate[k] = 0;
	}
	
	//the original processing point
	int iX = m_iWinSize / 2;
	int iY = m_iWinSize / 2;
	
	
	//begin to process
	int i, j, ii, jj;
	long a=iNumBits;	
	switch (iNumBits)
	{
	case 8:
		//image loop
		for (j=iY; j<iHeight-iY; j++)
		{
			for (i=iX; i<iWidth-iX; i++)
			{
				//template loop
				for (jj=0; jj<m_iWinSize; jj++)
				{
					for (ii=0; ii<m_iWinSize; ii++)
					{
						//Get gray value of the pixel in line (j-iY+jj) and column (i-iX+ii)
						piTemplate[jj*m_iWinSize+ii] = 
							pOldBits[(iHeight-1-(j-iY+jj))*nWidthBytes + (i-iX+ii)];
						
					}
				}
				
				//Sort the elements in the template
				BubbleSort(piTemplate, iAmount);
				//Get the median value in the template
				pNewBits[(iHeight-1-j)*nWidthBytes + i] = piTemplate[iAmount-1];
				 
			}
		}
		break;
	case 16:
		break;
	case 24:
		break;
	case 32:
		break;
	}
	
	//Free the temporary memory
	delete [] piTemplate;
	piTemplate = NULL;
	
	//Free the memory for the old image
	::GlobalUnlock(m_pImageObject->GetDib());
	::GlobalFree(m_pImageObject->GetDib());
	
	//Update the DIB in the object pointed by m_pImaeObject
	::GlobalUnlock(hNewDib); 
	m_pImageObject->SetDib(hNewDib);
	
	return TRUE;
}
/**********************************************************************
* Copyright (c) 2003, Medical Image Processing Lab, Sichuan University 
* All rights reserved.
* 
* 文件名称：Equalize.cpp
* 摘    要：图像最小值滤波CFilter的源文件
* 
* 当前版本：1.1 
* 作    者：jian wei zhang 
* 完成日期：2004年11月8日
************************************************************************/
BOOL CFilter::MinFilter()
{
if (m_pImageObject == NULL)
	{
		return FALSE;
	}
	
	int iWidth = m_pImageObject->GetWidth();
	int iHeight = m_pImageObject->GetHeight();
	int iNumBits=m_pImageObject->GetNumBits();
	/******************************the old image*****************************/
	unsigned char *pOldBuffer, *pNewBuffer, *pOldBits, *pNewBits;
	BITMAPFILEHEADER *pOldBFH, *pNewBFH;
	BITMAPINFOHEADER *pOldBIH, *pNewBIH;
	RGBQUAD *pOldPalette, *pNewPalette;
	int nWidthBytes, nNumColors;
	
	pOldBuffer = (unsigned char *) m_pImageObject->GetDIBPointer(&nWidthBytes, m_pImageObject->GetNumBits());
	if (pOldBuffer == NULL) 
	{
		return FALSE;
	}
	
	pOldBFH = (BITMAPFILEHEADER *) pOldBuffer;
	pOldBIH = (BITMAPINFOHEADER *) &pOldBuffer[sizeof(BITMAPFILEHEADER)];
	nNumColors = m_pImageObject->GetNumColors();
	pOldPalette = (RGBQUAD *)  &pOldBuffer[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)];
	pOldBits = (unsigned char *) &pOldBuffer[sizeof(BITMAPFILEHEADER)
		+sizeof(BITMAPINFOHEADER)+nNumColors*sizeof(RGBQUAD)];
	
	/*****************************the new image*****************************/
	DWORD dwNewSize;
	HGLOBAL hNewDib;
	
	//Allocate a global memory block for the new image
	dwNewSize = sizeof(BITMAPFILEHEADER) + sizeof(BITMAPINFOHEADER) 
		+ nNumColors * sizeof(RGBQUAD)
		+ nWidthBytes * (m_pImageObject->GetHeight());
	hNewDib = ::GlobalAlloc(GMEM_MOVEABLE | GMEM_ZEROINIT, dwNewSize);
	if (hNewDib == NULL)
	{
		m_pImageObject->m_nLastError = IMAGELIB_MEMORY_ALLOCATION_ERROR;
		::GlobalUnlock(m_pImageObject->GetDib());
		return FALSE;
	}
	
	//Get the pointer to the new memory block
	pNewBuffer = (unsigned char *) ::GlobalLock(hNewDib);
	if (pNewBuffer == NULL)
	{
		::GlobalFree( hNewDib );
		m_pImageObject->m_nLastError = IMAGELIB_MEMORY_LOCK_ERROR;
		::GlobalUnlock(m_pImageObject->GetDib());
		return FALSE;
	}
	
	pNewBFH = (BITMAPFILEHEADER *) pNewBuffer;
	pNewBIH = (BITMAPINFOHEADER *) &pNewBuffer[sizeof(BITMAPFILEHEADER)];
	pNewPalette = (RGBQUAD *) &pNewBuffer[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)];
	pNewBits = (unsigned char *) &pNewBuffer[sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+nNumColors*sizeof(RGBQUAD)];
	
	//Make a copy of the old image
	memcpy(pNewBuffer,pOldBuffer,dwNewSize);
	
	/*******************************Max filter****************************/
	//define a template of the median filter
	int iAmount = m_iWinSize*m_iWinSize;
	int* piTemplate = new int[iAmount];
	for (int k=0; k<iAmount; k++)
	{
		piTemplate[k] = 0;
	}
	
	//the original processing point
	int iX = m_iWinSize / 2;
	int iY = m_iWinSize / 2;
	
	
	//begin to process
	int i, j, ii, jj;
	long a=iNumBits;	
	switch (iNumBits)
	{
12 下一页
💿 文件大小 507 K
👤 上传用户 wugu0
📂 所属分类数学计算
🏷️ 相关标签

#图象处理 #傅立叶变换 #程序
⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -