可视域分析R3算法代码

数字地形分析中的R3可视域分析算法代码共享，提供数据检测

//DEM.CPP
/****************************************
**
**      REFERENCE PLANE ALGORITHM
**
*****************************************/
//DEM.CPP

#include“DEM.h“
#include“StdAfx.h“
#include“INPUT.h“
#include“OUTPUT.h“
using namespace std;
DEM::DEM（）{}

DEM::~DEM（）{}

	// Intial function receive the information of input object.
	// Then intialize the member of class DEM.
void DEM::Intial（Input& input）{

	m_ncols=input.GetNcols（）;
	m_nrows=input.GetNrows（）;
	m_xllcorner=input.GetXCorner（）;
	m_yllcorner=input.GetYCorner（）;
	m_cellsize=input.GetCellsize（）;
	m_NODATA_VALUE=input.GetNODATA_VALUE（）;
	counter=0;

	Viewshed=new int*[m_nrows];Aux=new AUXILIARY*[m_nrows];
	for（int i=0;i		Viewshed[i]=new int[m_ncols];
		Aux[i]=new AUXILIARY[m_ncols];
	}

	for（int i=0;i		for（int j=0;j			（*（Aux+i）+j）->m_row=i;
			（*（Aux+i）+j）->m_col=j;
			（*（Aux+i）+j）->m_altd=input.GetAltitude（ij）;
			//infile>>s;
			Viewshed[i][j]=0;
		}
		// Caution:The xllcorner and yllcorner coordinates are in southwestbut the altitude array begins in northwest.

}

int DEM::GetNrows（）{
	return m_nrows;
}

int DEM::GetNcols（）{
	return m_ncols;
}

int DEM::GetCellsize（）{
	return m_cellsize;
}

int DEM::GetNODATA_VALUE（）{
	return m_NODATA_VALUE;
}
	// delieve row and column if need get altitude
status DEM::GetAltitude（int iint j）{
	return （*（Aux+i）+j）->m_altd;
}

float DEM::GetXCorner（）{
	return m_xllcorner;
}

float DEM::GetYCorner（）{
	return m_yllcorner;
}

float DEM::IChangeY（int i）{
	return m_yllcorner+float（（m_nrows-i）*m_cellsize）;
}

float DEM::JChangeX（int j）{
	return m_xllcorner+float（（j-1）*m_cellsize）;
}

AUXILIARY *DEM::GetAltitudePnt（int iint j）{
	return *（Aux+i）+j;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    1401028  2012-11-14 20:23  R3 algrothm\Debug\m.txt

     文件     118272  2012-12-24 16:10  R3 algrothm\Debug\R3 algrothm.exe

     文件     860984  2012-12-24 16:10  R3 algrothm\Debug\R3 algrothm.ilk

     文件    1018880  2012-12-24 16:10  R3 algrothm\Debug\R3 algrothm.pdb

     文件     389802  2012-12-12 19:06  R3 algrothm\Debug\R3_t1.txt

     文件       1056  2012-11-05 18:55  R3 algrothm\Debug\t.txt

     文件   46530560  2012-12-24 16:07  R3 algrothm\ipch\r3 algrothm-3f32fe11\r3 algrothm-3bd3734f.ipch

     文件       4186  2012-12-24 16:09  R3 algrothm\R3 algrothm\Debug\cl.command.1.tlog

     文件     121078  2012-12-24 16:09  R3 algrothm\R3 algrothm\Debug\CL.read.1.tlog

     文件       4522  2012-12-24 16:09  R3 algrothm\R3 algrothm\Debug\CL.write.1.tlog

     文件      54913  2012-12-10 19:25  R3 algrothm\R3 algrothm\Debug\DEM.obj

     文件     235615  2012-12-08 10:19  R3 algrothm\R3 algrothm\Debug\INPUT.obj

     文件          2  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\link-cvtres.read.1.tlog

     文件          2  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\link-cvtres.write.1.tlog

     文件          2  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\link.1896-cvtres.read.1.tlog

     文件          2  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\link.1896-cvtres.write.1.tlog

     文件          2  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\link.1896.read.1.tlog

     文件          2  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\link.1896.write.1.tlog

     文件       2474  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\link.command.1.tlog

     文件       4352  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\link.read.1.tlog

     文件       1584  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\link.write.1.tlog

     文件     117314  2012-12-24 16:09  R3 algrothm\R3 algrothm\Debug\main.obj

     文件        434  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\mt.command.1.tlog

     文件        378  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\mt.read.1.tlog

     文件        378  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\mt.write.1.tlog

     文件     219999  2012-12-10 19:25  R3 algrothm\R3 algrothm\Debug\OUTPUT.obj

     文件       2524  2012-12-08 10:19  R3 algrothm\R3 algrothm\Debug\R3 algrothm.Build.CppClean.log

     文件        406  2012-12-08 10:19  R3 algrothm\R3 algrothm\Debug\R3 algrothm.exe.embed.manifest

     文件        472  2012-12-08 10:19  R3 algrothm\R3 algrothm\Debug\R3 algrothm.exe.embed.manifest.res

     文件        381  2012-12-24 16:10  R3 algrothm\R3 algrothm\Debug\R3 algrothm.exe.intermediate.manifest

............此处省略40个文件信息