灰度共生矩阵VC++实现

1、/ GrayShow.cpp: implementation of the CGrayShow class./ 灰灰灰灰灰灰灰灰灰灰 /#include "stdafx.h"#include "GrayShow.h"#include "math.h"#ifdef _DEBUG#undef THIS_FILEstatic char THIS_FILE=_FILE_;#define new DEBUG_NEW#endif/ Construction/Destruction/CGrayShow:CGrayShow()PMatrixRD =

2、NULL; /45灰度度度度灰灰灰灰灰灰 PMatrixLD = NULL; /135灰度度度度灰灰灰灰灰灰 PMatrixV = NULL; /90灰度度度度灰灰灰灰灰灰 PMatrixH = NULL; /0灰度度度度灰灰灰灰灰灰/ ImageArray = NULL; /原原度原原原原原原/m_tOriPixelArray = NULL;/原原原原原原 /m_tResPixelArray = NULL;/参参原原原原/ m_dib = new CDib();distance = 5;FilterWindowWidth = 16;GrayLayerNum = 16; /初原初初初8个灰灰个

3、，可可可可 /这这灰灰灰灰这初GrayLayerNum×GrayLayerNumPMatrixH = new intGrayLayerNum*GrayLayerNum;PMatrixLD= new intGrayLayerNum*GrayLayerNum; PMatrixRD= new intGrayLayerNum*GrayLayerNum;PMatrixV = new intGrayLayerNum*GrayLayerNum;/ for(i=0; i<GrayLayerNum; i+)/ / PMatrixHi = new intGrayLayerNum;/ PMatrix

4、LDi= new intGrayLayerNum;/ PMatrixRDi= new intGrayLayerNum;/ PMatrixVi = new intGrayLayerNum;/ CGrayShow:CGrayShow()if(PMatrixRD)delete PMatrixRD ;if(PMatrixLD)delete PMatrixLD ; /135灰度度度度灰灰灰灰灰灰 if(PMatrixV)delete PMatrixV ; /90灰度度度度灰灰灰灰灰灰 if(PMatrixH)delete PMatrixH;/功功：装装原原 /参原：FilePathName原原图图图图图

5、 / void CGrayShow:LoadImage(CString &FilePathName)/ / / / int i,j;/ /检检检检初检，不初检不不不不不不 / if(m_tOriPixelArray !=NULL)/ / for(i=0; i<ImageHeight; i+)/ / delete m_tOriPixelArrayi;/ delete m_tResPixelArrayi;/ delete ImageArrayi;/ / / m_dib->Open(FilePathName);/ ImageWidth = m_dib->GetWidth()

6、;/ ImageHeight = m_dib->GetHeight();/ BYTE *colorTable;/ colorTable = (BYTE *)m_dib->m_pDibBits;/ int byteBitCount = m_dib->GetBiBitCount()/8;/ m_tOriPixelArray = new RGBQUAD*ImageHeight;/ m_tResPixelArray = new RGBQUAD*ImageHeight;/ ImageArray = new BYTE*ImageHeight;/ for(int l=0 ; l<Im

7、ageHeight; l+)/ / m_tOriPixelArrayl = new RGBQUADImageWidth;/ m_tResPixelArrayl = new RGBQUADImageWidth;/ ImageArrayl = new BYTEImageWidth;/ / int count = 0;/ for(i=ImageHeight-1; i>=0; i-)/ / for(j=0; j<ImageWidth; j+)/ / m_tOriPixelArrayij.rgbBlue =colorTablecount+;/ m_tOriPixelArrayij.rgbGr

8、een=colorTablecount+;/ m_tOriPixelArrayij.rgbRed =colorTablecount+;/ m_tOriPixelArrayij.rgbReserved = 0;/ m_tResPixelArrayij=m_tOriPixelArrayij;/ ImageArrayij = (BYTE)(unsigned int)m_tOriPixelArrayij.rgbBlue/ +(unsigned int)m_tOriPixelArrayij.rgbGreen/ +(unsigned int)m_tOriPixelArrayij.rgbRed)/3);/

9、count += byteBitCount-3;/ / count += (4-(ImageWidth*byteBitCount)%4)%4;/ /函原功功：计灰计计计计 /参原：FeatureEnergy功能 / FeatureEntropy熵 / FeatureInertiaQuadrature惯惯灰 / FeatureCorrelation相相 / FeatureLocalCalm局局局局 / pMatrix灰灰灰灰 / dim灰灰灰灰度共原 /void CGrayShow:ComputeFeature(double &FeatureEnergy, double &Fea

10、tureEntropy,double &FeatureInertiaQuadrature, double &FeatureCorrelation, double &FeatureLocalCalm, int* pMatrix, int dim)int i,j;double *pdMatrix;pdMatrix = new doubledim*dim;int total = 0;for(i=0; i<dim; i+)for(j=0; j<dim; j+)total += pMatrixi*dim+j;for(i=0; i<dim; i+)for(j=0;

11、 j<dim; j+)pdMatrixi*dim + j = (double)pMatrixi*dim + j/(double)total;FeatureEnergy = 0.0;FeatureEntropy = 0.0;FeatureInertiaQuadrature = 0.0;FeatureLocalCalm = 0.0; /计灰功能、熵、惯惯灰、局局局局 for(i=0; i<dim; i+)for(j=0; j<dim; j+) /功能 FeatureEnergy += pdMatrixi*dim+j*pdMatrixi*dim+j; /熵 if(pdMatrixi

12、*dim+j>1e-12)FeatureEntropy -= pdMatrixi*dim+j*log(pdMatrixi*dim+j); /惯惯灰 FeatureInertiaQuadrature += (double)(i-j)*(double)(i-j)*pdMatrixi*dim+j; /局局局局 FeatureLocalCalm += pdMatrixi*dim+j/(1+(double)(i-j)*(double)(i-j); /计灰ux double ux = 0.0;double localtotal = 0.0;for(i=0; i<dim; i+)localtot

13、al = 0.0;for(j=0; j<dim; j+)localtotal += pdMatrixi*dim+j;ux += (double)i * localtotal; /计灰uy double uy = 0.0;for(j=0; j<dim; j+)localtotal = 0.0;for(i=0; i<dim; i+)localtotal += pdMatrixi*dim+j;uy += (double)j * localtotal; /计灰sigmax double sigmax = 0.0;for(i=0; i<dim; i+)localtotal = 0

14、.0;for(j=0; j<dim; j+)localtotal += pdMatrixi*dim+j;sigmax += (double)(i-ux) * (double)(i-ux) * localtotal; /计灰sigmay double sigmay = 0.0;for(j=0; j<dim; j+) localtotal = 0.0;for(i=0; i<dim; i+)localtotal += pdMatrixi*dim+j;sigmay += (double)(j-uy) * (double)(j-uy) * localtotal;/计灰相相 Featur

15、eCorrelation = 0.0;for(i=0; i<dim; i+)for(j=0; j<dim; j+)FeatureCorrelation += (double)(i-ux) * (double)(j-uy) * pdMatrixi*dim+j;FeatureCorrelation /= sigmax;FeatureCorrelation /= sigmay;delete pdMatrix;pdMatrix = NULL;/功功：计灰灰灰灰灰 /参原：LocalImage用用计灰度局局计计用用原原 / window局局计计用用局灰 /void CGrayShow:Com

16、puteMatrix(LPBYTE LocalImage,const WORD width,/local image 度局 const WORD height,/local image 度的 int i_count, /rowint j_count, /col int window )memset(PMatrixH,0,GrayLayerNum*GrayLayerNum*sizeof(int);memset(PMatrixLD,0,GrayLayerNum*GrayLayerNum*sizeof(int); memset(PMatrixRD,0,GrayLayerNum*GrayLayerNu

17、m*sizeof(int);memset(PMatrixV,0,GrayLayerNum*GrayLayerNum*sizeof(int);int i,j;/计灰0灰度灰灰灰灰灰 for( i=i_count; i<i_count+window; i+)for( j=j_count; j<j_count+window-distance; j+)LPBYTE lptemp = LocalImage + i*width*3 + j*3; BYTE gray_I_J =BYTE(*(lptemp)*0.114 + *(lptemp+1)*0.587+*(lptemp+2)*.299)/(

18、256/GrayLayerNum); LPBYTE lpdistance = lptemp +distance*3;BYTE gray_I_Distance = BYTE(*(lpdistance)*0.114+*(lpdistance+1)*0.587+*(lpdistance+2)*0.299)/(256/GrayLayerNum);PMatrixHgray_I_J*GrayLayerNum + gray_I_Distance += 1; PMatrixHgray_I_Distance*GrayLayerNum + gray_I_J += 1;/计灰90灰度灰灰灰灰灰 for(i=i_co

19、unt; i<i_count+window-distance; i+)for(j=j_count; j<j_count+window; j+) LPBYTE lptemp = LocalImage + i*width*3 + j*3;BYTE gray_I_J =BYTE(*(lptemp)*0.114 + *(lptemp+1)*0.587+*(lptemp+2)*.299)/(256/GrayLayerNum);LPBYTE lpdistance = /*lptemp +distance*3*/LocalImage + (i+distance)*width*3 + j*3;BY

20、TE gray_I_Distance = BYTE(*(lpdistance)*0.114+*(lpdistance+1)*0.587+*(lpdistance+2)*0.299)/(256/GrayLayerNum); PMatrixVgray_I_J*GrayLayerNum + gray_I_Distance += 1;PMatrixVgray_I_Distance*gray_I_J + gray_I_J += 1;/计灰135灰度灰灰灰灰灰 for(i=i_count; i<i_count+window-distance; i+)for(j=j_count; j<j_cou

21、nt+window-distance; j+)int newi, newj;newi = i+distance;newj = j+distance;LPBYTE lptemp = LocalImage + i*width*3 + j*3;BYTE gray_I_J =BYTE(*(lptemp)*0.114 + *(lptemp+1)*0.587+*(lptemp+2)*.299)/(256/GrayLayerNum);LPBYTE lpdistance = /*lptemp +distance*3*/LocalImage + (newi)*width*3 + newj*3;BYTE gray

22、_I_Distance = BYTE(*(lpdistance)*0.114+*(lpdistance+1)*0.587+*(lpdistance+2)*0.299)/(256/GrayLayerNum); PMatrixLDgray_I_J*GrayLayerNum + gray_I_Distance += 1;PMatrixLDgray_I_Distance*gray_I_J + gray_I_J += 1; /计灰45灰度灰灰灰灰灰 for(i=i_count+distance; i<i_count+window; i+) for(j=j_count; j<j_count+w

23、indow-distance; j+) int newi, newj;newi = i-distance;newj = j+distance;LPBYTE lptemp = LocalImage + i*width*3 + j*3;BYTE gray_I_J =BYTE(*(lptemp)*0.114 + *(lptemp+1)*0.587+*(lptemp+2)*.299)/(256/GrayLayerNum);LPBYTE lpdistance = /*lptemp +distance*3*/LocalImage + (newi)*width*3 + newj*3;BYTE gray_I_Distance = BYTE(*(lpdistance)*0.114+*(lpdistance+1)*0.587+*(lpdistance+2)*0.299)/(256/GrayLayerNum);PMatrixRDgr