contour.cpp

基于GRID格式DEM的等值线追踪算法（已删除编译目录）

// Contour.cpp: implementation of the CContour class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "ContourGenerator.h"
#include "ContourTracer.h"
#include "ContourValues.h"
#include "Contour.h"

#include "EditContourLevelDlg.h"


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

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

CContour::CContour(CGeoData* pOwner)
{
	m_pDataOwner = pOwner;
	m_pContourDrawer = NULL;
	m_Contours.RemoveAll();
	m_bReGenerated = FALSE;

}


CContour::~CContour()
{
	RemoveAllContours();
	
	if( m_pContourDrawer != NULL )
		delete m_pContourDrawer;
}

void CContour::GenerateContours()
{//生成等值线

	int numContours = GetNumberOfValues();
	if( numContours <= 0 )
		return;


	if( this->IsHavingContours() )
	{
		this->RemoveAllContours();
	}

	CContourTracer Tracer;
	
	//给等值线追踪器设置输入
	Tracer.SetGridDataInfo(m_pDataOwner->GetDataInfo());
	Tracer.SetInput(m_pDataOwner->GetGridData());

	CCurveList *pCurveList;
	for(int i=0; i<numContours; i++)
	{
		pCurveList = new CCurveList;
		if(pCurveList==NULL)
		{
			AfxMessageBox("内存分配失败!");
			return;
		}

		Tracer.SetOutput(pCurveList);
		if( !Tracer.ExecuteTracing(m_contourValues.GetValue(i)) )
		{
			delete pCurveList;
		}
		else
		{
			m_Contours.Add(pCurveList);
		}		

	}

	m_bReGenerated = TRUE; //重新生成了等值线

}

void CContour::RemoveAllContours()
{
	int num = m_Contours.GetSize();
	if( num > 0 )
	{
		CCurveList* pCurveList;
		for(int i=0; i<num; i++)
		{
			pCurveList = this->m_Contours.GetAt(i);	
			ASSERT( pCurveList!=NULL );
			
			POSITION pos = pCurveList->GetHeadPosition();						
			while (pos!=NULL) 
			{
				delete pCurveList->GetNext(pos);
			}
			delete pCurveList;
		}

		m_Contours.RemoveAll();
	}
}

void CContour::EditContourLevels()
{
	CEditContourLevelDlg levelDlg(m_contourValues);


	ASSERT(m_pDataOwner!=NULL);
	levelDlg.SetZMinMax(m_pDataOwner->GetDataInfo().zMin,m_pDataOwner->GetDataInfo().zMax);

	if( levelDlg.DoModal() == IDOK )
	{
		m_contourValues = levelDlg.m_Levels;

		GenerateContours();

	}

}

void CContour::GenerateDefaultContourValues()
{//产生缺省值(10等份)

	ASSERT(m_pDataOwner!=NULL);
	float zMin = m_pDataOwner->GetDataInfo().zMin;
	float zMax = m_pDataOwner->GetDataInfo().zMax;

	float inc = (zMax-zMin)/9; //10等份

	float val;

	for(val= zMin; val<=zMax; val+=inc)
	{
		m_contourValues.AddValue(val);
	}
	

}

void CContour::Draw(CDC* pDC,const CRect& drawRect)
{//绘制等值线

    if( m_pContourDrawer == NULL)
	{
		m_pContourDrawer = new CContourDrawer(this);
	}

	ASSERT(m_pContourDrawer!=NULL);

	m_pContourDrawer->DoDraw(pDC,drawRect);

	m_bReGenerated = FALSE;

}