📄 flow_areaupslope.cpp
字号:
///////////////////////////////////////////////////////////
// //
// SAGA //
// //
// System for Automated Geoscientific Analyses //
// //
// Module Library: //
// ta_hydrology //
// //
//-------------------------------------------------------//
// //
// Flow_AreaUpslope.cpp //
// //
// Copyright (C) 2003 by //
// Olaf Conrad //
// //
//-------------------------------------------------------//
// //
// This file is part of 'SAGA - System for Automated //
// Geoscientific Analyses'. SAGA is free software; you //
// can redistribute it and/or modify it under the terms //
// of the GNU General Public License as published by the //
// Free Software Foundation; version 2 of the License. //
// //
// SAGA is distributed in the hope that it will be //
// useful, but WITHOUT ANY WARRANTY; without even the //
// implied warranty of MERCHANTABILITY or FITNESS FOR A //
// PARTICULAR PURPOSE. See the GNU General Public //
// License for more details. //
// //
// You should have received a copy of the GNU General //
// Public License along with this program; if not, //
// write to the Free Software Foundation, Inc., //
// 59 Temple Place - Suite 330, Boston, MA 02111-1307, //
// USA. //
// //
//-------------------------------------------------------//
// //
// e-mail: oconrad@saga-gis.org //
// //
// contact: Olaf Conrad //
// Institute of Geography //
// University of Goettingen //
// Goldschmidtstr. 5 //
// 37077 Goettingen //
// Germany //
// //
///////////////////////////////////////////////////////////
//---------------------------------------------------------
///////////////////////////////////////////////////////////
// //
// //
// //
///////////////////////////////////////////////////////////
//---------------------------------------------------------
#include "Flow_AreaUpslope.h"
///////////////////////////////////////////////////////////
// //
// //
// //
///////////////////////////////////////////////////////////
//---------------------------------------------------------
CFlow_AreaUpslope::CFlow_AreaUpslope(void)
{
m_pDTM = NULL;
m_pRoute = NULL;
m_pFlow = NULL;
}
//---------------------------------------------------------
CFlow_AreaUpslope::~CFlow_AreaUpslope(void)
{
Finalise();
}
///////////////////////////////////////////////////////////
// //
// //
// //
///////////////////////////////////////////////////////////
//---------------------------------------------------------
CSG_String CFlow_AreaUpslope::Get_Description(void)
{
return(
_TW(
"This module allows you to specify target cells, "
"for which the upslope contributing area shall be identified. "
"The result will give "
"for each cell the percentage of its flow that reaches the target cell(s).\n\n"
"References:\n\n"
"Deterministic 8\n"
"- O'Callaghan, J.F. / Mark, D.M. (1984):\n"
" 'The extraction of drainage networks from digital elevation data',\n"
" Computer Vision, Graphics and Image Processing, 28:323-344\n\n"
"Deterministic Infinity:\n"
"- Tarboton, D.G. (1997):\n"
" 'A new method for the determination of flow directions and upslope areas in grid digital elevation models',\n"
" Water Ressources Research, Vol.33, No.2, p.309-319\n\n"
"Multiple Flow Direction:\n"
"- Freeman, G.T. (1991):\n"
" 'Calculating catchment area with divergent flow based on a regular grid',\n"
" Computers and Geosciences, 17:413-22\n\n"
"- Quinn, P.F. / Beven, K.J. / Chevallier, P. / Planchon, O. (1991):\n"
" 'The prediction of hillslope flow paths for distributed hydrological modelling using digital terrain models',\n"
" Hydrological Processes, 5:59-79\n\n"
)
);
}
//---------------------------------------------------------
CSG_String CFlow_AreaUpslope::Get_Methods(void)
{
return( CSG_String::Format(SG_T("%s|%s|%s|"),
_TL("Deterministic 8"),
_TL("Deterministic Infinity"),
_TL("Multiple Flow Direction")
));
}
///////////////////////////////////////////////////////////
// //
// //
// //
///////////////////////////////////////////////////////////
//---------------------------------------------------------
bool CFlow_AreaUpslope::Initialise(int Method, CSG_Grid *pDTM, CSG_Grid *pRoute, CSG_Grid *pFlow, double MFD_Converge)
{
Finalise();
if( pDTM && pDTM->is_Valid() && pFlow && pFlow->is_Valid() && pFlow->Get_System() == pDTM->Get_System() )
{
m_Method = Method;
m_pDTM = pDTM;
m_pFlow = pFlow;
m_MFD_Converge = MFD_Converge;
if( pRoute && pRoute->is_Valid() && pRoute->Get_System() == pDTM->Get_System() )
{
m_pRoute = pRoute;
}
return( true );
}
return( false );
}
//---------------------------------------------------------
bool CFlow_AreaUpslope::Finalise(void)
{
m_pDTM = NULL;
m_pRoute = NULL;
m_pFlow = NULL;
return( true );
}
//---------------------------------------------------------
bool CFlow_AreaUpslope::Add_Target(int x, int y)
{
if( m_pFlow && m_pFlow->is_InGrid(x, y, false) )
{
m_pFlow->Set_Value(x, y, 100.0);
return( true );
}
return( false );
}
//---------------------------------------------------------
bool CFlow_AreaUpslope::Clr_Target(void)
{
if( m_pFlow )
{
m_pFlow->Assign(0.0);
return( true );
}
return( false );
}
//---------------------------------------------------------
bool CFlow_AreaUpslope::Get_Area(int x, int y)
{
return( Clr_Target() && Add_Target(x, y) && Get_Area() );
}
//---------------------------------------------------------
bool CFlow_AreaUpslope::Get_Area(void)
{
int i, x, y;
if( m_pDTM && m_pFlow )
{
for(i=0; i<m_pDTM->Get_NCells() && SG_UI_Process_Set_Progress(i, m_pDTM->Get_NCells()); i++)
{
m_pDTM->Get_Sorted(i, x, y, false);
if( m_pFlow->asDouble(x, y) > 0.0 )
{
break;
}
}
for(i++; i<m_pDTM->Get_NCells() && SG_UI_Process_Set_Progress(i, m_pDTM->Get_NCells()); i++)
{
m_pDTM->Get_Sorted(i, x, y, false);
Set_Value(x, y);
}
return( true );
}
return( false );
}
///////////////////////////////////////////////////////////
// //
// //
// //
///////////////////////////////////////////////////////////
//---------------------------------------------------------
void CFlow_AreaUpslope::Set_Value(int x, int y)
{
int i;
if( m_pRoute && (i = m_pRoute->asChar(x, y)) >= 0 )
{
int ix, iy;
double Flow;
ix = m_pDTM->Get_System().Get_xTo(i, x);
iy = m_pDTM->Get_System().Get_yTo(i, y);
if( m_pDTM->is_InGrid(ix, iy, true) && (Flow = m_pFlow->asDouble(ix, iy)) > 0.0 )
{
m_pFlow->Set_Value(x, y, Flow);
}
}
else if( !m_pDTM->is_NoData(x, y) )
{
switch( m_Method )
{
case 0: Set_D8 (x, y); break;
case 1: Set_DInf (x, y); break;
case 2: Set_MFD (x, y); break;
}
}
}
//---------------------------------------------------------
void CFlow_AreaUpslope::Set_D8(int x, int y)
{
int i;
if( (i = m_pDTM->Get_Gradient_NeighborDir(x, y, true)) >= 0 )
{
int ix, iy;
double Flow;
ix = m_pDTM->Get_System().Get_xTo(i, x);
iy = m_pDTM->Get_System().Get_yTo(i, y);
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -