//---------------------------------------------------------
bool CErosion_LS_Fields::Get_Statistics(void)
{
CSG_Shapes *pFields = Parameters("FIELDS" )->asShapes();
CSG_Shapes *pStatistics = Parameters("STATISTICS")->asShapes();
//-----------------------------------------------------
if( !pStatistics || !pFields || m_nFields <= 0 || pFields->Get_Count() != m_nFields )
{
return( false );
}
//-----------------------------------------------------
int i;
CSG_Simple_Statistics *Statistics = new CSG_Simple_Statistics[m_nFields];
//-----------------------------------------------------
for(int y=0; y<Get_NY() && Set_Progress(y); y++)
{
for(int x=0; x<Get_NX(); x++)
{
if( !m_pLS->is_NoData(x, y) && (i = m_Fields.asInt(x, y)) >= 0 && i < m_nFields )
{
Statistics[i] += m_pLS->asDouble(x, y);
}
}
}
//-----------------------------------------------------
pStatistics->Create(SHAPE_TYPE_Polygon, CSG_String::Format(SG_T("%s [%s]"), pFields->Get_Name(), _TL("LS")));
pStatistics->Add_Field("NCELLS", SG_DATATYPE_Int );
pStatistics->Add_Field("MEAN" , SG_DATATYPE_Double);
pStatistics->Add_Field("MIN" , SG_DATATYPE_Double);
pStatistics->Add_Field("MAX" , SG_DATATYPE_Double);
pStatistics->Add_Field("STDDEV", SG_DATATYPE_Double);
for(i=0; i<pFields->Get_Count() && Set_Progress(i, pFields->Get_Count()); i++)
{
CSG_Shape *pField = pStatistics->Add_Shape(pFields->Get_Shape(i));
if( Statistics[i].Get_Count() > 0 )
{
pField->Set_Value(0, Statistics[i].Get_Count ());
pField->Set_Value(1, Statistics[i].Get_Mean ());
pField->Set_Value(2, Statistics[i].Get_Minimum());
pField->Set_Value(3, Statistics[i].Get_Maximum());
pField->Set_Value(4, Statistics[i].Get_StdDev ());
}
else for(int j=0; j<pFields->Get_Field_Count(); j++)
{
pField->Set_NoData(j);
}
}
//-------------------------------------------------
delete[](Statistics);
return( true );
}
//---------------------------------------------------------
bool CCRS_Transform_Grid::Transform(CSG_Parameter_Grid_List *pGrids, CSG_Shapes *pPoints)
{
if( !pPoints || !pGrids || pGrids->Get_Count() < 1 )
{
return( false );
}
CSG_Grid *pGrid = pGrids->asGrid(0);
if( !m_Projector.Set_Source(pGrid->Get_Projection()) )
{
return( false );
}
int x, y, i;
TSG_Point Point;
pPoints->Create(SHAPE_TYPE_Point, _TL("Points"));
pPoints->Get_Projection() = m_Projector.Get_Target();
for(i=0; i<pGrids->Get_Count(); i++)
{
pPoints->Add_Field(pGrids->asGrid(i)->Get_Name(), pGrids->asGrid(i)->Get_Type());
}
for(y=0, Point.y=pGrid->Get_YMin(); y<pGrid->Get_NY() && Set_Progress(y, pGrid->Get_NY()); y++, Point.y+=pGrid->Get_Cellsize())
{
for(x=0, Point.x=pGrid->Get_XMin(); x<pGrid->Get_NX(); x++, Point.x+=pGrid->Get_Cellsize())
{
TSG_Point Point_Transformed = Point;
if( m_Projector.Get_Projection(Point_Transformed) )
{
CSG_Shape *pPoint = pPoints->Add_Shape();
pPoint->Add_Point(Point_Transformed);
for(i=0; i<pGrids->Get_Count(); i++)
{
if( !pGrids->asGrid(i)->is_NoData(x, y) )
{
pPoint->Set_Value(i, pGrids->asGrid(i)->asDouble(x, y));
}
else
{
pPoint->Set_NoData(i);
}
}
}
}
}
return( true );
}
//---------------------------------------------------------
bool CPROJ4_Grid::Set_Shapes(CSG_Parameter_Grid_List *pSources, CSG_Shapes *pTarget)
{
int x, y, i;
double z;
TSG_Point Pt_Source, Pt_Target;
CSG_Grid *pSource;
CSG_Shape *pShape;
if( pSources && pSources->Get_Count() > 0 && pTarget )
{
pSource = pSources->asGrid(0);
pTarget->Create(SHAPE_TYPE_Point, CSG_String::Format(SG_T("%s [%s]"), pSource->Get_Name(), Get_Proj_Name().c_str()));
for(i=0; i<pSources->Get_Count(); i++)
{
pTarget->Add_Field(pSources->asGrid(i)->Get_Name(), pSources->asGrid(i)->Get_Type());
}
for(y=0, Pt_Source.y=pSource->Get_YMin(); y<pSource->Get_NY() && Set_Progress(y, pSource->Get_NY()); y++, Pt_Source.y+=pSource->Get_Cellsize())
{
for(x=0, Pt_Source.x=pSource->Get_XMin(); x<pSource->Get_NX(); x++, Pt_Source.x+=pSource->Get_Cellsize())
{
if( !pSource->is_NoData(x, y) )
{
Pt_Target = Pt_Source;
if( Get_Converted(Pt_Target) )
{
pShape = pTarget->Add_Shape();
pShape->Add_Point(Pt_Target);
for(i=0; i<pSources->Get_Count(); i++)
{
if( pSources->asGrid(i)->Get_Value(Pt_Source, z, m_Interpolation) )
{
pShape->Set_Value(i, z);
}
else
{
pShape->Set_NoData(i);
}
}
}
}
}
}
return( true );
}
return( false );
}
//---------------------------------------------------------
bool CGrid_Values_AddTo_Points::On_Execute(void)
{
int iShape, iGrid, iField, Offset, Interpolation;
double Value;
CSG_Shapes *pShapes;
CSG_Parameter_Grid_List *pGrids;
//-----------------------------------------------------
pShapes = Parameters("RESULT") ->asShapes();
pGrids = Parameters("GRIDS" ) ->asGridList();
Interpolation = Parameters("INTERPOL") ->asInt();
//-----------------------------------------------------
if( pGrids->Get_Count() <= 0 )
{
return( false );
}
//-----------------------------------------------------
if( pShapes == NULL )
{
pShapes = Parameters("SHAPES")->asShapes();
}
else if( pShapes != Parameters("SHAPES")->asShapes() )
{
pShapes->Create(*Parameters("SHAPES")->asShapes());
}
Offset = pShapes->Get_Field_Count();
//-----------------------------------------------------
for(iGrid=0; iGrid<pGrids->Get_Count(); iGrid++)
{
pShapes->Add_Field(pGrids->asGrid(iGrid)->Get_Name(), SG_DATATYPE_Double);
}
//-----------------------------------------------------
for(iShape=0; iShape<pShapes->Get_Count() && Set_Progress(iShape, pShapes->Get_Count()); iShape++)
{
CSG_Shape *pShape = pShapes->Get_Shape(iShape);
for(iGrid=0, iField=Offset; iGrid<pGrids->Get_Count(); iGrid++, iField++)
{
CSG_Grid *pGrid = pGrids->asGrid(iGrid);
if( pGrid->Get_Value(pShape->Get_Point(0), Value, Interpolation) )
{
pShape->Set_Value(iField, Value);
}
else
{
pShape->Set_NoData(iField);
}
}
}
//-----------------------------------------------------
if( pShapes == Parameters("SHAPES")->asShapes() )
{
DataObject_Update(pShapes);
}
return( true );
}
//---------------------------------------------------------
bool CGrid_Values_AddTo_Shapes::On_Execute(void)
{
CSG_Parameter_Grid_List *pGrids;
CSG_Shapes *pShapes;
//-----------------------------------------------------
pShapes = Parameters("RESULT")->asShapes();
pGrids = Parameters("GRIDS" )->asGridList();
//-----------------------------------------------------
if( pGrids->Get_Count() <= 0 )
{
return( false );
}
//-----------------------------------------------------
if( pShapes == NULL )
{
pShapes = Parameters("SHAPES")->asShapes();
}
else if( pShapes != Parameters("SHAPES")->asShapes() )
{
pShapes->Create(*Parameters("SHAPES")->asShapes());
}
//-----------------------------------------------------
switch( Parameters("RESAMPLING")->asInt() )
{
default: m_Resampling = GRID_RESAMPLING_NearestNeighbour; break;
case 1: m_Resampling = GRID_RESAMPLING_Bilinear; break;
case 2: m_Resampling = GRID_RESAMPLING_BicubicSpline; break;
case 3: m_Resampling = GRID_RESAMPLING_BSpline; break;
}
//-----------------------------------------------------
for(int iGrid=0; iGrid<pGrids->Get_Count(); iGrid++)
{
CSG_Grid *pGrid = pGrids->asGrid(iGrid);
int Field = pShapes->Get_Field_Count();
pShapes->Add_Field(pGrid->Get_Name(), SG_DATATYPE_Double);
for(int iShape=0; iShape<pShapes->Get_Count() && Set_Progress(iShape, pShapes->Get_Count()); iShape++)
{
CSG_Simple_Statistics Statistics;
CSG_Shape *pShape = pShapes->Get_Shape(iShape);
if( pShape->Get_Extent().Intersects(pGrid->Get_Extent()) )
{
switch( pShapes->Get_Type() )
{
case SHAPE_TYPE_Point: default:
case SHAPE_TYPE_Points: Get_Data_Point (Statistics, pShape, pGrid); break;
case SHAPE_TYPE_Line: Get_Data_Line (Statistics, pShape, pGrid); break;
case SHAPE_TYPE_Polygon: Get_Data_Polygon(Statistics, pShape, pGrid); break;
}
}
if( Statistics.Get_Count() > 0 )
{
pShape->Set_Value(Field, Statistics.Get_Mean());
}
else
{
pShape->Set_NoData(Field);
}
}
}
//-----------------------------------------------------
if( pShapes == Parameters("SHAPES")->asShapes() )
{
DataObject_Update(pShapes);
}
return( true );
}
