//---------------------------------------------------------
bool CGSGrid_Residuals::Get_Statistics(int x, int y, bool bCenter)
{
if( m_pGrid->is_InGrid(x, y) )
{
int i, ix, iy, nLower;
double z, iz, id, iw;
CSG_Simple_Statistics Statistics;
for(i=0, nLower=0, z=m_pGrid->asDouble(x, y); i<m_Cells.Get_Count(); i++)
{
if( m_Cells.Get_Values(i, ix = x, iy = y, id, iw, true) && (bCenter || id > 0.0) && m_pGrid->is_InGrid(ix, iy) )
{
Statistics.Add_Value(iz = m_pGrid->asDouble(ix, iy), iw);
if( z > iz )
{
nLower++;
}
}
}
//-------------------------------------------------
if( Statistics.Get_Weights() > 0.0 )
{
m_pMean ->Set_Value(x, y, Statistics.Get_Mean());
m_pDiff ->Set_Value(x, y, z - Statistics.Get_Mean());
m_pStdDev ->Set_Value(x, y, Statistics.Get_StdDev());
m_pRange ->Set_Value(x, y, Statistics.Get_Range());
m_pMin ->Set_Value(x, y, Statistics.Get_Minimum());
m_pMax ->Set_Value(x, y, Statistics.Get_Maximum());
m_pDevMean ->Set_Value(x, y, Statistics.Get_StdDev() > 0.0 ? ((z - Statistics.Get_Mean()) / Statistics.Get_StdDev()) : 0.0);
m_pPercent ->Set_Value(x, y, 100.0 * nLower / (double)Statistics.Get_Count());
return( true );
}
}
//-----------------------------------------------------
m_pMean ->Set_NoData(x, y);
m_pDiff ->Set_NoData(x, y);
m_pStdDev ->Set_NoData(x, y);
m_pRange ->Set_NoData(x, y);
m_pMin ->Set_NoData(x, y);
m_pMax ->Set_NoData(x, y);
m_pDevMean ->Set_NoData(x, y);
m_pPercent ->Set_NoData(x, y);
return( false );
}
//---------------------------------------------------------
bool CGPS_Track_Aggregation::Set_Statistic(CSG_Table_Record *pAggregate, CSG_Simple_Statistics &Statistic, CSG_Simple_Statistics &Time, int nDropped, bool bVerbose)
{
if( pAggregate )
{
pAggregate ->Set_Value(AGG_PARM , Statistic.Get_Mean());
pAggregate ->Set_Value(AGG_TIME , Time .Get_Mean());
if( bVerbose )
{
pAggregate ->Set_Value(AGG_MIN , Statistic.Get_Minimum());
pAggregate ->Set_Value(AGG_MAX , Statistic.Get_Maximum());
pAggregate ->Set_Value(AGG_RANGE , Statistic.Get_Range ());
pAggregate ->Set_Value(AGG_STDDEV , Statistic.Get_StdDev ());
pAggregate ->Set_Value(AGG_COUNT , Statistic.Get_Count ());
pAggregate ->Set_Value(AGG_DTIME , Time .Get_Range ());
pAggregate ->Set_Value(AGG_DROPPED, nDropped);
}
return( true );
}
return( false );
}
//---------------------------------------------------------
bool CGSGrid_Statistics::On_Execute(void)
{
//-----------------------------------------------------
CSG_Parameter_Grid_List *pGrids = Parameters("GRIDS")->asGridList();
if( pGrids->Get_Count() <= 1 )
{
Error_Set(_TL("no grids in selection"));
return( false );
}
//-----------------------------------------------------
CSG_Grid *pMean = Parameters("MEAN" )->asGrid();
CSG_Grid *pMin = Parameters("MIN" )->asGrid();
CSG_Grid *pMax = Parameters("MAX" )->asGrid();
CSG_Grid *pRange = Parameters("RANGE" )->asGrid();
CSG_Grid *pSum = Parameters("SUM" )->asGrid();
CSG_Grid *pVar = Parameters("VAR" )->asGrid();
CSG_Grid *pStdDev = Parameters("STDDEV" )->asGrid();
CSG_Grid *pStdDevLo = Parameters("STDDEVLO")->asGrid();
CSG_Grid *pStdDevHi = Parameters("STDDEVHI")->asGrid();
CSG_Grid *pPercentile = Parameters("PCTL" )->asGrid();
if( !pMean && !pMin && !pMax && !pRange && !pSum && !pVar && !pStdDev && !pStdDevLo && !pStdDevHi && !pPercentile )
{
Error_Set(_TL("no parameter output specified"));
return( false );
}
double dRank = Parameters("PCTL_VAL")->asDouble() / 100.0;
//-----------------------------------------------------
for(int y=0; y<Get_NY() && Set_Progress(y); y++)
{
#pragma omp parallel for
for(int x=0; x<Get_NX(); x++)
{
CSG_Table Values;
CSG_Simple_Statistics s;
for(int i=0; i<pGrids->Get_Count(); i++)
{
if( !pGrids->asGrid(i)->is_NoData(x, y) )
{
double z = pGrids->asGrid(i)->asDouble(x, y);
s.Add_Value(z);
if( pPercentile )
{
if( Values.Get_Field_Count() == 0 )
{
Values.Add_Field("Z", SG_DATATYPE_Double);
}
Values.Add_Record()->Set_Value(0, z);
}
}
}
//-----------------------------------------
if( s.Get_Count() <= 0 )
{
if( pMean ) pMean ->Set_NoData(x, y);
if( pMin ) pMin ->Set_NoData(x, y);
if( pMax ) pMax ->Set_NoData(x, y);
if( pRange ) pRange ->Set_NoData(x, y);
if( pSum ) pSum ->Set_NoData(x, y);
if( pVar ) pVar ->Set_NoData(x, y);
if( pStdDev ) pStdDev ->Set_NoData(x, y);
if( pStdDevLo ) pStdDevLo ->Set_NoData(x, y);
if( pStdDevHi ) pStdDevHi ->Set_NoData(x, y);
if( pPercentile ) pPercentile ->Set_NoData(x, y);
}
else
{
if( pMean ) pMean ->Set_Value(x, y, s.Get_Mean());
if( pMin ) pMin ->Set_Value(x, y, s.Get_Minimum());
if( pMax ) pMax ->Set_Value(x, y, s.Get_Maximum());
if( pRange ) pRange ->Set_Value(x, y, s.Get_Range());
if( pSum ) pSum ->Set_Value(x, y, s.Get_Sum());
if( pVar ) pVar ->Set_Value(x, y, s.Get_Variance());
if( pStdDev ) pStdDev ->Set_Value(x, y, s.Get_StdDev());
if( pStdDevLo ) pStdDevLo ->Set_Value(x, y, s.Get_Mean() - s.Get_StdDev());
if( pStdDevHi ) pStdDevHi ->Set_Value(x, y, s.Get_Mean() + s.Get_StdDev());
if( pPercentile )
{
Values.Set_Index(0, TABLE_INDEX_Ascending);
pPercentile->Set_Value(x, y, Values.Get_Record_byIndex((int)(dRank * s.Get_Count()))->asDouble(0));
}
}
}
}
//-----------------------------------------------------
return( true );
}
