dlgintensity.cpp

该软件是我读硕士的时候写的

// DlgIntensity.cpp : implementation file
//

#include "stdafx.h"
#include "process.h"
#include "DlgIntensity.h"

#include <windows.h>
#include <math.h>
#include <malloc.h>
#include "resource.h"
#include <direct.h>

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

/////////////////////////////////////////////////////////////////////////////
// CDlgIntensity dialog


CDlgIntensity::CDlgIntensity(CWnd* pParent /*=NULL*/)
	: CDialog(CDlgIntensity::IDD, pParent)
{
	//{{AFX_DATA_INIT(CDlgIntensity)
	m_nDataFrom = 0.0;
	m_nDataTo = 0.0;
	//}}AFX_DATA_INIT
}


void CDlgIntensity::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CDlgIntensity)
	DDX_Text(pDX, IDC_EDIT2, m_nDataFrom);
	DDX_Text(pDX, IDC_EDIT3, m_nDataTo);
	//}}AFX_DATA_MAP
}


BEGIN_MESSAGE_MAP(CDlgIntensity, CDialog)
	//{{AFX_MSG_MAP(CDlgIntensity)
	ON_WM_PAINT()
	ON_BN_CLICKED(IDC_BUTTON1, OnShowIntnsity)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CDlgIntensity message handlers

BOOL CDlgIntensity::OnInitDialog()
{
 CDialog::OnInitDialog();

	// 返回TRUE
	return TRUE;

}

void CDlgIntensity::OnPaint() 
{
	m_nDataFrom=nDataFrom;
	m_nDataTo=nDataTo;
	UpdateData(false);

	// 字符串
	CString str;
	
	// 循环变量
	LONG i;
	
	// 最大计数
	int lMaxCount = 0;
	
	// 设备上下文
	CPaintDC dc(this);
	
	// 获取绘制坐标的文本框
	CWnd* pWnd = GetDlgItem(IDC_COORD);
	
	// 指针
	CDC* pDC = pWnd->GetDC();
	pWnd->Invalidate();
	pWnd->UpdateWindow();
	
	pDC->Rectangle(0,0,720,300);
	
	// 创建画笔对象
	CPen* pPenRed = new CPen;
	
	// 红色画笔
	pPenRed->CreatePen(PS_SOLID,2,RGB(255,0,0));
	
	// 创建画笔对象
	CPen* pPenBlue = new CPen;
	
	// 蓝色画笔
	pPenBlue->CreatePen(PS_SOLID,1,RGB(0,0, 255));
	
	// 创建画笔对象
	CPen* pPenGreen = new CPen;
	
	// 绿色画笔
	pPenGreen->CreatePen(PS_DOT,1,RGB(0,255,0));
	
	// 选中当前红色画笔，并保存以前的画笔
	CGdiObject* pOldPen = pDC->SelectObject(pPenRed);
	
	// 绘制坐标轴
	pDC->MoveTo(10,10);
	
	// 垂直轴
	pDC->LineTo(10,280);
	
	// 水平轴
	pDC->LineTo(710,280);
	
	// 写X轴刻度值
	str.Format("0");
	pDC->TextOut(10, 283, str);
	str.Format("50");
	pDC->TextOut(60, 283, str);
	str.Format("100");
	pDC->TextOut(110, 283, str);
	str.Format("150");
	pDC->TextOut(160, 283, str);
	str.Format("200");
	pDC->TextOut(210, 283, str);
	str.Format("250");
	pDC->TextOut(260, 283, str);

	str.Format("300");
	pDC->TextOut(310, 283, str);
	str.Format("350");
	pDC->TextOut(360, 283, str);
	str.Format("400");
	pDC->TextOut(410, 283, str);
	str.Format("450");
	pDC->TextOut(460, 283, str);
	str.Format("500");
	pDC->TextOut(510, 283, str);
	str.Format("550");
	pDC->TextOut(560, 283, str);
	str.Format("600");
	pDC->TextOut(610, 283, str);
	str.Format("650");
	pDC->TextOut(660, 283, str);
	// 绘制X轴刻度
	for (i = 0; i < 660; i += 5)
	{
		if ((i & 1) == 0)
		{
			// 10的倍数
			pDC->MoveTo(i + 10, 280);
			pDC->LineTo(i + 10, 284);
		}
		else
		{
			// 10的倍数
			pDC->MoveTo(i + 10, 280);
			pDC->LineTo(i + 10, 282);
		}
	}
	
	// 绘制X轴箭头
	pDC->MoveTo(705,275);
	pDC->LineTo(710,280);
	pDC->LineTo(705,285);
	
	// 绘制y轴箭头
	pDC->MoveTo(10,10);
	pDC->LineTo(5,15);
	pDC->MoveTo(10,10);
	pDC->LineTo(15,15);
	
    // 计算最大计数值
	for (i = 0; i < MaxData; i ++)
	{
		// 判断是否大于当前最大值
		if (m_lCount[i] > lMaxCount)
		{
			// 更新最大值
			lMaxCount = m_lCount[i];
		}
	}
	
	// 输出最大计数值
	pDC->MoveTo(10, 25);
	pDC->LineTo(14, 25);
	str.Format("%d", lMaxCount);
	pDC->TextOut(20, 26, str);
	
	
	// 更改成蓝色画笔
	pDC->SelectObject(pPenBlue);	
	
	// 判断是否有计数
	if (lMaxCount > 0)
	{
		// 绘制直方图
		for (i = 0; i < MaxData; i ++)
		{
			pDC->MoveTo(i + 10, 280);
			pDC->LineTo(i + 10, 281 - (int) (m_lCount[i] * 256 / lMaxCount));
		}
	}

	// 恢复以前的画笔
	pDC->SelectObject(pOldPen);	
	
	// 删除新的画笔
	delete pPenRed;
	delete pPenBlue;
	delete pPenGreen;
}

void CDlgIntensity::OnShowIntnsity() 
{    //显示局部直方图
	int i;
	int j;
	int temp;
	
    UpdateData(true);
	nDataFrom=m_nDataFrom;
	nDataTo=m_nDataTo;
    //计算所要显示的局部曲线的范围
	int temp2,temp1;
	temp1=0;
	temp2=DataNum-1;
    while(nDataFrom>DataDepth[temp1])
	{
		temp1++;
	}
	while(nDataTo<DataDepth[temp2])
	{
		temp2--;
	}

	MaxData=0;
	//测井曲线的最大值
	for (int ii = temp1; ii < temp2; ii++)
	{
        (int)temp=DataData[ii];
		if(temp>MaxData)  MaxData=temp;
	}
		
	
	for (i = 0; i < MaxData; i ++)
	{
		// 清零
		m_lCount[i] = 0;
	}
	// 计算各个灰度值的计数
	for (j =temp1; j <temp2; j ++)
	{
		// 计数加1
		m_lCount[(int)DataData[j]]++;
	}
	
   OnPaint(); 
}