//---------------------------------------------------------
void CD8_Flow_Analysis::Get_Basins(void)
{
Process_Set_Text(_TL("Drainage Basins"));
for(int y=0; y<Get_NY() && Set_Progress(y); y++)
{
for(int x=0; x<Get_NX(); x++)
{
Get_Basin(x, y);
}
}
//-----------------------------------------------------
CSG_Shapes *pBasins = Parameters("BASINS")->asShapes();
if( pBasins )
{
bool bResult;
SG_RUN_MODULE(bResult, "shapes_grid", 6,
pModule->Get_Parameters()->Set_Parameter(SG_T("GRID") , m_pBasins)
&& pModule->Get_Parameters()->Set_Parameter(SG_T("POLYGONS"), pBasins)
)
pBasins->Set_Name(_TL("Drainage Basins"));
}
}
//---------------------------------------------------------
bool CCost_Accumulated::Get_Destinations(CPoints &Points)
{
Points.Clear();
m_pAccumulated->Set_NoData_Value(-1.0); m_pAccumulated->Assign(-1.0);
m_pAllocation ->Set_NoData_Value(-1.0); m_pAllocation ->Assign( 0.0);
if( Parameters("DEST_TYPE")->asInt() == 0 ) // Point
{
CSG_Shapes *pDestinations = Parameters("DEST_POINTS")->asShapes();
for(int i=0, x, y; i<pDestinations->Get_Count(); i++)
{
if( Get_System().Get_World_to_Grid(x, y, pDestinations->Get_Shape(i)->Get_Point(0)) && !m_pCost->is_NoData(x, y) )
{
Points.Add(x, y); m_pAllocation->Set_Value(x, y, Points.Get_Count()); m_pAccumulated->Set_Value(x, y, 0.0);
}
}
}
else // Grid
{
CSG_Grid *pDestinations = Parameters("DEST_GRID")->asGrid();
for(int y=0; y<Get_NY(); y++) for(int x=0; x<Get_NX(); x++)
{
if( !pDestinations->is_NoData(x, y) && !m_pCost->is_NoData(x, y) )
{
Points.Add(x, y); m_pAllocation->Set_Value(x, y, Points.Get_Count()); m_pAccumulated->Set_Value(x, y, 0.0);
}
}
}
return( Points.Get_Count() > 0 );
}
bool CCreatePointGrid::On_Execute(void){
CSG_Shape *pShape;
double dXMin = ((CSG_Parameter_Range *) Parameters("X_EXTENT")->Get_Data())->Get_LoVal();
double dYMin = ((CSG_Parameter_Range *) Parameters("Y_EXTENT")->Get_Data())->Get_LoVal();
double dXMax = ((CSG_Parameter_Range *) Parameters("X_EXTENT")->Get_Data())->Get_HiVal();
double dYMax = ((CSG_Parameter_Range *) Parameters("Y_EXTENT")->Get_Data())->Get_HiVal();
double dDistance = Parameters("DIST")->asDouble();
if (dDistance<=0){
return false;
}//if
CSG_Shapes *pShapes = Parameters("POINTS")->asShapes();
pShapes->Create(SHAPE_TYPE_Point, _TL("Point Grid"));
pShapes->Add_Field("X", SG_DATATYPE_Double);
pShapes->Add_Field("Y", SG_DATATYPE_Double);
for (double x=dXMin; x<dXMax; x=x+dDistance){
for (double y=dYMin; y<dYMax; y=y+dDistance){
pShape = pShapes->Add_Shape();
pShape->Add_Point(x,y);
pShape->Set_Value(0, x);
pShape->Set_Value(1, y);
}//for
}//for
return true;
}//method
//---------------------------------------------------------
bool C_Kriging_Base::_Get_Points(void)
{
int iShape, iPart, iPoint;
CSG_Shape *pShape , *pPoint;
CSG_Shapes *pPoints;
m_pShapes = Parameters("SHAPES") ->asShapes();
m_zField = Parameters("FIELD") ->asInt();
if( m_pShapes->Get_Type() != SHAPE_TYPE_Point )
{
pPoints = SG_Create_Shapes(SHAPE_TYPE_Point, SG_T(""), m_pShapes);
for(iShape=0; iShape<m_pShapes->Get_Count() && Set_Progress(iShape, m_pShapes->Get_Count()); iShape++)
{
pShape = m_pShapes->Get_Shape(iShape);
if( !pShape->is_NoData(m_zField) )
{
for(iPart=0; iPart<pShape->Get_Part_Count(); iPart++)
{
for(iPoint=0; iPoint<pShape->Get_Point_Count(iPart); iPoint++)
{
pPoint = pPoints->Add_Shape(pShape, SHAPE_COPY_ATTR);
pPoint->Add_Point(pShape->Get_Point(iPoint, iPart));
}
}
}
}
m_pShapes = pPoints;
}
return( m_pShapes->Get_Count() > 1 );
}
//---------------------------------------------------------
bool CGridding_Spline_Base::_Get_Grid(void)
{
CSG_Grid *pGrid = m_bGridPoints ? Parameters("GRIDPOINTS")->asGrid () : NULL;
CSG_Shapes *pShapes = !m_bGridPoints ? Parameters("SHAPES") ->asShapes() : NULL;
//-----------------------------------------------------
m_pGrid = NULL;
switch( Parameters("TARGET")->asInt() )
{
case 0: // user defined...
if( m_Grid_Target.Init_User(m_bGridPoints ? pGrid->Get_Extent() : pShapes->Get_Extent()) && Dlg_Parameters("USER") )
{
m_pGrid = m_Grid_Target.Get_User();
}
break;
case 1: // grid...
if( Dlg_Parameters("GRID") )
{
m_pGrid = m_Grid_Target.Get_Grid();
}
break;
}
//-------------------------------------------------
if( m_pGrid )
{
m_pGrid->Set_Name(CSG_String::Format(SG_T("%s (%s)"), m_bGridPoints ? pGrid->Get_Name() : pShapes->Get_Name(), Get_Name()));
m_pGrid->Assign_NoData();
}
//-----------------------------------------------------
return( m_pGrid != NULL );
}
bool CPoints_From_Table::On_Execute(void)
{
int iRecord, iXField, iYField;
double dX, dY;
CSG_Table *pTable;
CSG_Table_Record *pRecord;
CSG_Shapes *pShapes;
CSG_Shape *pShape;
pTable = Parameters("TABLE") ->asTable();
pShapes = Parameters("POINTS") ->asShapes();
iXField = Parameters("X") ->asInt();
iYField = Parameters("Y") ->asInt();
pShapes->Create(SHAPE_TYPE_Point, pTable->Get_Name(), pTable);
if( pTable->Get_Field_Count() > 1 && pTable->Get_Record_Count() > 0 )
{
for(iRecord=0; iRecord<pTable->Get_Record_Count() && Set_Progress(iRecord, pTable->Get_Record_Count()); iRecord++)
{
pRecord = pTable->Get_Record(iRecord);
dX = pRecord->asDouble(iXField);
dY = pRecord->asDouble(iYField);
pShape = pShapes->Add_Shape(pRecord, SHAPE_COPY_ATTR);
pShape->Add_Point(dX, dY);
}//for
return( true );
}//if
return( false );
}//method
//---------------------------------------------------------
bool CShapes_Create_Empty::On_Execute(void)
{
TSG_Vertex_Type Vertex;
switch( Parameters("VERTEX")->asInt() )
{
default:
case 0: Vertex = SG_VERTEX_TYPE_XY; break;
case 1: Vertex = SG_VERTEX_TYPE_XYZ; break;
case 2: Vertex = SG_VERTEX_TYPE_XYZM; break;
}
//-----------------------------------------------------
CSG_Shapes *pShapes;
switch( Parameters("TYPE")->asInt() )
{
default: return( false );
case 0: pShapes = SG_Create_Shapes(SHAPE_TYPE_Point , Parameters("NAME")->asString(), NULL, Vertex); break;
case 1: pShapes = SG_Create_Shapes(SHAPE_TYPE_Points , Parameters("NAME")->asString(), NULL, Vertex); break;
case 2: pShapes = SG_Create_Shapes(SHAPE_TYPE_Line , Parameters("NAME")->asString(), NULL, Vertex); break;
case 3: pShapes = SG_Create_Shapes(SHAPE_TYPE_Polygon, Parameters("NAME")->asString(), NULL, Vertex); break;
}
//-----------------------------------------------------
int i, n;
TSG_Data_Type Type;
CSG_Parameters *pAttributes;
pAttributes = Parameters("FIELDS")->asParameters();
n = pAttributes->Get_Count() / 3;
for(i=0; i<n; i++)
{
switch( pAttributes->Get_Parameter(GET_TYPE(i))->asInt() )
{
default:
case 0: Type = SG_DATATYPE_String; break;
case 1: Type = SG_DATATYPE_Date; break;
case 2: Type = SG_DATATYPE_Color; break;
case 3: Type = SG_DATATYPE_Byte; break;
case 4: Type = SG_DATATYPE_Char; break;
case 5: Type = SG_DATATYPE_Word; break;
case 6: Type = SG_DATATYPE_Short; break;
case 7: Type = SG_DATATYPE_DWord; break;
case 8: Type = SG_DATATYPE_Int; break;
case 9: Type = SG_DATATYPE_ULong; break;
case 10: Type = SG_DATATYPE_Long; break;
case 11: Type = SG_DATATYPE_Float; break;
case 12: Type = SG_DATATYPE_Double; break;
case 13: Type = SG_DATATYPE_Binary; break;
}
pShapes->Add_Field(pAttributes->Get_Parameter(GET_NAME(i))->asString(), Type);
}
Parameters("SHAPES")->Set_Value(pShapes);
return( true );
}
//---------------------------------------------------------
bool CGrid_Fill::On_Execute(void)
{
CSG_Shapes *pPoints = Parameters("POINTS")->asShapes();
if( !pPoints->is_Valid() || !pPoints->Get_Extent().Intersects(Get_System().Get_Extent()) || !Parameters_Set(Parameters) )
{
return( false );
}
//-----------------------------------------------------
int nReplaced = 0;
#define GET_NPOINTS (bSelection ? pPoints->Get_Selection_Count() : pPoints->Get_Count())
#define GET_POINT(i) (bSelection ? pPoints->Get_Selection(i)->Get_Point(0) : pPoints->Get_Shape(i)->Get_Point(0))
bool bSelection = pPoints->Get_Selection_Count() > 0;
for(int i=0; i<GET_NPOINTS && Process_Get_Okay(); i++)
{
nReplaced += Fill(GET_POINT(i));
}
Message_Fmt("\n%d %s\n", nReplaced, _TL("replacements"));
return( true );
}
//---------------------------------------------------------
CSG_Shapes * COGR_DataSource::Read_Shapes(int iLayer)
{
OGRLayer *pLayer = Get_Layer(iLayer);
//-----------------------------------------------------
if( pLayer && Get_Type(iLayer) != SHAPE_TYPE_Undefined )
{
int iField;
OGRFeature *pFeature;
OGRFeatureDefn *pDef = pLayer->GetLayerDefn();
CSG_Shapes *pShapes = SG_Create_Shapes(Get_Type(iLayer), CSG_String(pDef->GetName()));
for(iField=0; iField<pDef->GetFieldCount(); iField++)
{
OGRFieldDefn *pDefField = pDef->GetFieldDefn(iField);
pShapes->Add_Field(pDefField->GetNameRef(), COGR_Driver::Get_Type(pDefField->GetType()));
}
pLayer->ResetReading();
//-------------------------------------------------
while( (pFeature = pLayer->GetNextFeature()) != NULL && SG_UI_Process_Get_Okay(false) )
{
OGRGeometry *pGeometry = pFeature->GetGeometryRef();
if( pGeometry != NULL )
{
CSG_Shape *pShape = pShapes->Add_Shape();
for(iField=0; iField<pDef->GetFieldCount(); iField++)
{
OGRFieldDefn *pDefField = pDef->GetFieldDefn(iField);
switch( pDefField->GetType() )
{
default: pShape->Set_Value(iField, SG_STR_MBTOSG(pFeature->GetFieldAsString (iField))); break;
case OFTString: pShape->Set_Value(iField, SG_STR_MBTOSG(pFeature->GetFieldAsString (iField))); break;
case OFTInteger: pShape->Set_Value(iField, pFeature->GetFieldAsInteger(iField)); break;
case OFTReal: pShape->Set_Value(iField, pFeature->GetFieldAsDouble (iField)); break;
}
}
//-----------------------------------------
if( _Read_Geometry(pShape, pGeometry) == false )
{
pShapes->Del_Shape(pShape);
}
}
OGRFeature::DestroyFeature(pFeature);
}
return( pShapes );
}
//-----------------------------------------------------
return( NULL );
}
//---------------------------------------------------------
void CGrid_Cross_Profiles::Make_Report(const SG_Char *FileName, CSG_Grid *pDEM, CSG_Shapes *pLines, CSG_Shapes *pProfiles, double Distance)
{
if( FileName )
{
int iProfile, iPoint, nSamples, iBox;
CSG_Rect r;
CSG_Shape *pProfile, *pLine;
CSG_Shapes Profile;
CSG_Doc_PDF pdf;
pdf.Open(PDF_PAGE_SIZE_A4, PDF_PAGE_ORIENTATION_PORTRAIT, _TL("Cross Profiles"));
pdf.Layout_Add_Box(5, 5, 95, 20);
pdf.Layout_Add_Box(5, 25, 95, 45);
pdf.Layout_Add_Box(5, 50, 95, 70);
pdf.Layout_Add_Box(5, 75, 95, 90);
nSamples = pProfiles->Get_Field_Count() - OFFSET;
Distance = Distance / (nSamples - 1);
iBox = 0;
for(iProfile=0; iProfile<pProfiles->Get_Count() && Set_Progress(iProfile, pProfiles->Get_Count()); iProfile++)
{
pProfile = pProfiles->Get_Shape(iProfile);
Profile.Create(SHAPE_TYPE_Line);
pLine = Profile.Add_Shape();
for(iPoint=0; iPoint<nSamples; iPoint++)
{
if( !pProfile->is_NoData(OFFSET + iPoint) )
{
pLine->Add_Point(
(iPoint - nSamples / 2) * Distance,
pProfile->asDouble(OFFSET + iPoint)
);
}
}
if( pLine->Get_Point_Count(0) > 1 )
{
if( iBox >= NBOXES )
{
pdf.Add_Page();
iBox = 0;
}
// pdf.Draw_Graticule (pdf.Layout_Get_Box(iBox), Profile.Get_Extent(), 20);
// pdf.Draw_Shapes (pdf.Layout_Get_Box(iBox), &Profile);
r = pLine->Get_Extent(); r.Inflate(10);
pdf.Draw_Graticule (pdf.Layout_Get_Box(iBox), r, 15);
pdf.Draw_Shape (pdf.Layout_Get_Box(iBox), pLine, PDF_STYLE_POLYGON_FILLSTROKE, SG_COLOR_GREEN, SG_COLOR_BLACK, 1, &r);
iBox++;
}
}
pdf.Save(Parameters("DOCUMENT")->asString());
}
}
//---------------------------------------------------------
bool C_Kriging_Base::_Get_Grid(void)
{
CSG_Shapes *pShapes = Parameters("SHAPES")->asShapes();
m_pGrid = NULL;
m_pVariance = NULL;
//-------------------------------------------------
switch( Parameters("TARGET")->asInt() )
{
case 0: // user defined...
if( Dlg_Parameters("USER") )
{
m_pGrid = _Get_Grid(pShapes->Get_Extent());
}
break;
case 1: // grid system...
if( Dlg_Parameters("SYSTEM") )
{
m_pGrid = SG_Create_Grid(*Get_Parameters("SYSTEM")->Get_Parameter("SYSTEM")->asGrid_System(), SG_DATATYPE_Float);
}
break;
case 2: // grid...
if( Dlg_Parameters("GRID") )
{
m_pGrid = Get_Parameters("GRID")->Get_Parameter("GRID") ->asGrid();
m_pVariance = Get_Parameters("GRID")->Get_Parameter("VARIANCE") ->asGrid();
}
break;
}
//-------------------------------------------------
if( m_pGrid )
{
if( !m_pVariance && Parameters("BVARIANCE")->asBool() )
{
m_pVariance = SG_Create_Grid(m_pGrid, SG_DATATYPE_Float);
}
m_pGrid->Set_Name(CSG_String::Format(SG_T("%s (%s)"), pShapes->Get_Name(), Get_Name()));
Parameters("GRID")->Set_Value(m_pGrid);
if( m_pVariance )
{
m_pVariance->Set_Name(CSG_String::Format(SG_T("%s (%s - %s)"), pShapes->Get_Name(), Get_Name(), _TL("Variance")));
Parameters("VARIANCE")->Set_Value(m_pVariance);
}
if( Parameters("TARGET")->asInt() == 2 )
{
Get_Parameters("GRID")->Get_Parameter("VARIANCE")->Set_Value(m_pVariance);
}
}
//-----------------------------------------------------
return( m_pGrid != NULL );
}
//---------------------------------------------------------
bool CLine_Properties::On_Execute(void)
{
CSG_Shapes *pLines = Parameters("LINES")->asShapes();
if( pLines->is_Valid() && pLines->Get_Count() > 0 )
{
if( Parameters("OUTPUT")->asShapes() && Parameters("OUTPUT")->asShapes() != pLines )
{
pLines = Parameters("OUTPUT")->asShapes();
pLines->Create(*Parameters("LINES")->asShapes());
}
//-------------------------------------------------
int iOffset = pLines->Get_Field_Count();
pLines->Add_Field(SG_T("N_VERTICES"), SG_DATATYPE_Int);
pLines->Add_Field(SG_T("LENGTH") , SG_DATATYPE_Double);
//-------------------------------------------------
for(int iLine=0; iLine<pLines->Get_Count() && Set_Progress(iLine, pLines->Get_Count()); iLine++)
{
CSG_Shape *pLine = pLines->Get_Shape(iLine);
pLine->Set_Value(iOffset + 0, pLine->Get_Point_Count());
pLine->Set_Value(iOffset + 1, ((CSG_Shape_Line *)pLine)->Get_Length());
}
return( true );
}
return( false );
}
//---------------------------------------------------------
bool CGPX_Import::On_Execute(void)
{
CSG_Shapes *pWay;
CSG_MetaData GPX;
//-----------------------------------------------------
m_Name = Parameters("FILE") ->asString();
m_pShapes = Parameters("SHAPES") ->asShapesList();
m_bTime = Parameters("TIME") ->asBool();
//-----------------------------------------------------
if( !GPX.Create(m_Name) || GPX.Get_Name().CmpNoCase(SG_T("gpx")) )
{
return( false );
}
//-----------------------------------------------------
pWay = SG_Create_Shapes(SHAPE_TYPE_Point, m_Name);
m_Name = SG_File_Get_Name(m_Name, false);
m_pShapes->Del_Items();
//-----------------------------------------------------
for(int i=0; i<GPX.Get_Children_Count(); i++)
{
CSG_MetaData *pChild = GPX.Get_Child(i);
if( pChild->Get_Name().CmpNoCase(SG_T("wpt")) == 0 )
{
Add_Point(pChild, pWay);
}
else if( pChild->Get_Name().CmpNoCase(SG_T("rte")) == 0 )
{
Add_Route(pChild);
}
else if( pChild->Get_Name().CmpNoCase(SG_T("trk")) == 0 )
{
Add_Track(pChild);
}
}
//-----------------------------------------------------
if( pWay->Get_Count() > 0 )
{
m_pShapes->Add_Item(pWay);
}
else
{
delete(pWay);
}
return( m_pShapes->Get_Count() > 0 );
}
//---------------------------------------------------------
void CSG_Network::_On_Construction(void)
{
m_Nodes.Create(sizeof(CSG_Network_Node **), 0, SG_ARRAY_GROWTH_1);
m_Edges.Create(SHAPE_TYPE_Line , SG_T("EDGES"));
m_Edges.Add_Field(SG_T("ID") , SG_DATATYPE_Int);
m_Edges.Add_Field(SG_T("NODE_A") , SG_DATATYPE_Int);
m_Edges.Add_Field(SG_T("NODE_B") , SG_DATATYPE_Int);
m_Edges.Add_Field(SG_T("SHAPE_TYPE"), SG_DATATYPE_Int);
m_Edges.Add_Field(SG_T("PROCESSED") , SG_DATATYPE_Int);
}
//---------------------------------------------------------
bool CSelection_Invert::On_Execute(void)
{
CSG_Shapes *pInput;
pInput = Parameters("INPUT") ->asShapes();
pInput->Inv_Selection();
DataObject_Update(pInput);
return( true );
}
//---------------------------------------------------------
bool CGridding_Spline_Base::_Get_Points(CSG_Points_Z &Points, bool bInGridOnly)
{
Points.Clear();
if( m_bGridPoints )
{
int x, y;
TSG_Point p;
CSG_Grid *pGrid = Parameters("GRIDPOINTS") ->asGrid();
for(y=0, p.y=pGrid->Get_YMin(); y<pGrid->Get_NY() && Set_Progress(y, pGrid->Get_NY()); y++, p.y+=pGrid->Get_Cellsize())
{
for(x=0, p.x=pGrid->Get_XMin(); x<pGrid->Get_NX(); x++, p.x+=pGrid->Get_Cellsize())
{
if( !pGrid->is_NoData(x, y) && (!bInGridOnly || m_pGrid->is_InGrid_byPos(p)) )
{
Points.Add(p.x, p.y, pGrid->asDouble(x, y));
}
}
}
}
else
{
CSG_Shapes *pShapes = Parameters("SHAPES") ->asShapes();
int zField = Parameters("FIELD") ->asInt();
for(int iShape=0; iShape<pShapes->Get_Count() && Set_Progress(iShape, pShapes->Get_Count()); iShape++)
{
CSG_Shape *pShape = pShapes->Get_Shape(iShape);
if( !pShape->is_NoData(zField) )
{
double zValue = pShape->asDouble(zField);
for(int iPart=0; iPart<pShape->Get_Part_Count(); iPart++)
{
for(int iPoint=0; iPoint<pShape->Get_Point_Count(iPart); iPoint++)
{
TSG_Point p = pShape->Get_Point(iPoint, iPart);
if( !bInGridOnly || m_pGrid->is_InGrid_byPos(p) )
{
Points.Add(p.x, p.y, zValue);
}
}
}
}
}
}
return( Points.Get_Count() >= 3 );
}
//---------------------------------------------------------
bool CAddCoordinates::On_Execute(void)
{
CSG_Shapes *pShapes = Parameters("OUTPUT")->asShapes();
if( pShapes )
{
pShapes->Assign(Parameters("INPUT")->asShapes());
}
else
{
pShapes = Parameters("INPUT")->asShapes();
}
//-----------------------------------------------------
int xField = pShapes->Get_Field_Count();
pShapes->Add_Field("X", SG_DATATYPE_Double);
int yField = pShapes->Get_Field_Count();
pShapes->Add_Field("Y", SG_DATATYPE_Double);
//-----------------------------------------------------
for(int i=0; i<pShapes->Get_Count(); i++)
{
CSG_Shape *pShape = pShapes->Get_Shape(i);
pShape->Set_Value(xField, pShape->Get_Point(0).x);
pShape->Set_Value(yField, pShape->Get_Point(0).y);
}
return( true );
}
//---------------------------------------------------------
bool CPolygon_Vertex_Check::On_Execute(void)
{
//-----------------------------------------------------
CSG_Shapes *pPolygons = Parameters("POLYGONS")->asShapes();
if( Parameters("CHECKED")->asShapes() && Parameters("CHECKED")->asShapes() != pPolygons )
{
CSG_Shapes *pCopy = Parameters("CHECKED")->asShapes();
pCopy->Create(*pPolygons);
pCopy->Set_Name(CSG_String::Format(SG_T("%s [%s]"), pPolygons->Get_Name(), _TL("checked")));
pPolygons = pCopy;
}
if( (m_pAdded = Parameters("ADDED")->asShapes()) != NULL )
{
m_pAdded->Create(SHAPE_TYPE_Point, _TL("Added"));
}
double Epsilon = Parameters("EPSILON")->asDouble();
//-----------------------------------------------------
for(int iPolygon=0; iPolygon<pPolygons->Get_Count()-1 && Set_Progress(iPolygon, pPolygons->Get_Count()-1); iPolygon++)
{
CSG_Shape_Polygon *pA = (CSG_Shape_Polygon *)pPolygons->Get_Shape(iPolygon);
for(int jPolygon=iPolygon+1; jPolygon<pPolygons->Get_Count() && Process_Get_Okay(); jPolygon++)
{
CSG_Shape_Polygon *pB = (CSG_Shape_Polygon *)pPolygons->Get_Shape(jPolygon);
for(int iPart=0; iPart<pA->Get_Part_Count() && Process_Get_Okay(); iPart++)
{
for(int jPart=0; jPart<pB->Get_Part_Count() && Process_Get_Okay(); jPart++)
{
if( pA->Get_Part(iPart)->Get_Extent().Intersects(pB->Get_Part(jPart)->Get_Extent()) )
{
Check_Vertices(pA->Get_Part(iPart), pB->Get_Part(jPart), Epsilon);
Check_Vertices(pB->Get_Part(jPart), pA->Get_Part(iPart), Epsilon);
}
}
}
}
}
//-----------------------------------------------------
return( true );
}
//---------------------------------------------------------
bool CGSPoints_Distances::On_Execute(void)
{
//-----------------------------------------------------
CSG_Shapes *pPoints = Parameters("POINTS") ->asShapes();
CSG_Table *pTable = Parameters("TABLE") ->asTable();
//-----------------------------------------------------
CSG_PRQuadTree QT(pPoints, 0);
CSG_Simple_Statistics s;
double x, y, z;
for(int iPoint=0; iPoint<pPoints->Get_Count() && Set_Progress(iPoint, pPoints->Get_Count()); iPoint++)
{
TSG_Point p = pPoints->Get_Shape(iPoint)->Get_Point(0);
if( QT.Select_Nearest_Points(p.x, p.y, 2) && QT.Get_Selected_Point(1, x, y, z) && (x != p.x || y != p.y) )
{
s.Add_Value(SG_Get_Distance(x, y, p.x, p.y));
}
}
//-----------------------------------------------------
if( s.Get_Count() > 0 )
{
CSG_Table_Record *pRecord;
pTable->Destroy();
pTable->Set_Name(CSG_String::Format(SG_T("%s [%s]"), _TL("Minimum Distance Analysis"), pPoints->Get_Name()));
pTable->Add_Field(SG_T("NAME") , SG_DATATYPE_String);
pTable->Add_Field(SG_T("VALUE") , SG_DATATYPE_Double);
SET_VALUE(_TL("Mean Average") , s.Get_Mean());
SET_VALUE(_TL("Minimum") , s.Get_Minimum());
SET_VALUE(_TL("Maximum") , s.Get_Maximum());
SET_VALUE(_TL("Standard Deviation") , s.Get_StdDev());
SET_VALUE(_TL("Duplicates") , pPoints->Get_Count() - s.Get_Count());
DataObject_Update(pTable, SG_UI_DATAOBJECT_SHOW);
return( true );
}
Message_Dlg(_TL("not enough observations"));
return( false );
}
bool CCreateChartLayer::On_Execute(void){
CSG_Shapes *pInput;
int i=0;
int iType;
int iSizeField;
if (GetExtraParameters()){
iSizeField = Parameters("SIZE")->asInt();
m_fMaxSize = (float) Parameters("MAXSIZE")->asDouble();
m_fMinSize = (float) Parameters("MINSIZE")->asDouble();
if (m_fMinSize > m_fMaxSize){
m_fMinSize = m_fMaxSize;
}//if
iType = Parameters("TYPE")->asInt();
pInput = Parameters("INPUT")->asShapes();
m_fMaxValue = pInput->Get_Maximum(iSizeField);
m_fMinValue = pInput->Get_Minimum(iSizeField);
if (iType == TYPE_PIE){
m_pOutput = SG_Create_Shapes(SHAPE_TYPE_Polygon, _TL("Chart (sectors):"));
}//if
else{
m_pOutput = SG_Create_Shapes(SHAPE_TYPE_Polygon, _TL("Chart (bars):"));
}//else
m_pOutput->Add_Field(_TL("Field (ID)"), SG_DATATYPE_Int);
m_pOutput->Add_Field(_TL("Field (Name)"), SG_DATATYPE_String);
for (i = 0; i < pInput->Get_Count(); i++){
if (iType == TYPE_PIE){
// AddPieChart(pInput->Get_Shape(i),m_pOutput->Get_Type());
AddPieChart(pInput->Get_Shape(i),pInput->Get_Type());
}//if
else{
// AddBarChart(pInput->Get_Shape(i),m_pOutput->Get_Type());
AddBarChart(pInput->Get_Shape(i),pInput->Get_Type());
}//else
}//for
DataObject_Add(m_pOutput, false);
delete [] m_bIncludeParam;
return true;
}//if
delete [] m_bIncludeParam;
return false;
}//method
//---------------------------------------------------------
int CShapes_Save::On_Parameter_Changed(CSG_Parameters *pParameters, CSG_Parameter *pParameter)
{
if( !SG_STR_CMP(pParameter->Get_Identifier(), "SHAPES") )
{
CSG_Shapes *pShapes = pParameter->asShapes() ? pParameter->asShapes() : NULL;
if( pShapes && *pShapes->Get_Name() )
{
pParameters->Get_Parameter("NAME")->Set_Value(pShapes->Get_Name());
}
Set_SRID(pParameters, pShapes ? pShapes->Get_Projection().Get_EPSG() : -1);
}
return( CSG_PG_Module::On_Parameter_Changed(pParameters, pParameter) );
}
//---------------------------------------------------------
bool CXYZ_Import::On_Execute(void)
{
CSG_File Stream;
if( !Stream.Open(Parameters("FILENAME")->asString(), SG_FILE_R) )
{
Error_Set(_TL("file could not be opened"));
return( false );
}
//-----------------------------------------------------
CSG_Shapes *pPoints = Parameters("POINTS")->asShapes();
pPoints->Create(SHAPE_TYPE_Point, SG_File_Get_Name(Parameters("FILENAME")->asString(), false));
pPoints->Add_Field("Z", SG_DATATYPE_Double);
//-----------------------------------------------------
if( Parameters("HEADLINE")->asBool() )
{
CSG_String sLine;
if( !Stream.Read_Line(sLine) )
{
Error_Set(_TL("could not read headline"));
return( false );
}
}
//-----------------------------------------------------
sLong Length = Stream.Length();
double x, y, z;
while( Stream.Scan(x) && Stream.Scan(y) && Stream.Scan(z) && Set_Progress((double)Stream.Tell(), (double)Length) )
{
CSG_Shape *pPoint = pPoints->Add_Shape();
pPoint->Add_Point(x, y);
pPoint->Set_Value(0, z);
}
return( pPoints->Get_Count() > 0 );
}
//---------------------------------------------------------
bool CSG_Shapes::Assign(CSG_Data_Object *pObject)
{
int iShape;
CSG_Shape *pShape;
CSG_Shapes *pShapes;
//-----------------------------------------------------
if( pObject && pObject->is_Valid() && (pObject->Get_ObjectType() == DATAOBJECT_TYPE_Shapes || pObject->Get_ObjectType() == DATAOBJECT_TYPE_PointCloud) )
{
pShapes = (CSG_Shapes *)pObject;
Create(pShapes->Get_Type(), pShapes->Get_Name(), pShapes, pShapes->Get_Vertex_Type());
for(iShape=0; iShape<pShapes->Get_Count() && SG_UI_Process_Set_Progress(iShape, pShapes->Get_Count()); iShape++)
{
pShape = Add_Shape();
pShape->Assign(pShapes->Get_Shape(iShape));
}
SG_UI_Process_Set_Ready();
Update();
Get_History() = pObject->Get_History();
return( true );
}
return( false );
}
//---------------------------------------------------------
int CGW_Multi_Regression::On_Parameter_Changed(CSG_Parameters *pParameters, CSG_Parameter *pParameter)
{
if( m_Grid_Target.On_User_Changed(pParameters, pParameter) )
{
return( true );
}
if( !SG_STR_CMP(pParameter->Get_Identifier(), SG_T("POINTS")) )
{
CSG_Shapes *pPoints = pParameters->Get_Parameter("POINTS") ->asShapes();
CSG_Parameters *pAttributes = pParameters->Get_Parameter("PREDICTORS") ->asParameters();
pAttributes->Destroy();
pAttributes->Set_Name(_TL("Predictors"));
for(int i=0; pPoints && i<pPoints->Get_Field_Count(); i++)
{
switch( pPoints->Get_Field_Type(i) )
{
default: // not numeric
break;
case SG_DATATYPE_Byte:
case SG_DATATYPE_Char:
case SG_DATATYPE_Word:
case SG_DATATYPE_Short:
case SG_DATATYPE_DWord:
case SG_DATATYPE_Int:
case SG_DATATYPE_ULong:
case SG_DATATYPE_Long:
case SG_DATATYPE_Float:
case SG_DATATYPE_Double:
pAttributes->Add_Value(
NULL, SG_Get_String(i, 0), pPoints->Get_Field_Name(i), _TL(""), PARAMETER_TYPE_Bool, false
);
break;
}
}
return( true );
}
return( false );
}
//---------------------------------------------------------
bool CCRS_Transform_Shapes::On_Execute_Transformation(void)
{
if( m_bList )
{
CSG_Parameter_Shapes_List *pSources, *pTargets;
pSources = Parameters("SOURCE")->asShapesList();
pTargets = Parameters("TARGET")->asShapesList();
pTargets->Del_Items();
for(int i=0; i<pSources->Get_Count() && Process_Get_Okay(false); i++)
{
CSG_Shapes *pSource = pSources->asShapes(i);
CSG_Shapes *pTarget = SG_Create_Shapes(pSource);
if( Transform(pSource, pTarget) )
{
pTargets->Add_Item(pTarget);
}
else
{
delete(pTarget);
}
}
return( pTargets->Get_Count() > 0 );
}
else
{
CSG_Shapes *pSource = Parameters("SOURCE")->asShapes();
CSG_Shapes *pTarget = Parameters("TARGET")->asShapes();
if( pSource == pTarget )
{
pTarget = SG_Create_Shapes(pSource);
if( Transform(pSource, pTarget) )
{
pSource->Assign(pTarget);
return( true );
}
else
{
delete(pTarget);
return( false );
}
}
else
{
pTarget->Create(pSource->Get_Type(), pSource->Get_Name(), pSource, pSource->Get_Vertex_Type());
}
return( Transform(pSource, pTarget) );
}
}
//---------------------------------------------------------
bool CGW_Multi_Regression::Get_Predictors(void)
{
int i;
CSG_Shapes *pPoints = Parameters("POINTS") ->asShapes();
CSG_Parameters *pAttributes = Parameters("PREDICTORS") ->asParameters();
m_nPredictors = 0;
m_iPredictor = new int[pPoints->Get_Field_Count()];
for(i=0; i<pAttributes->Get_Count(); i++)
{
if( pAttributes->Get_Parameter(i)->asBool() )
{
m_iPredictor[m_nPredictors++] = CSG_String(pAttributes->Get_Parameter(i)->Get_Identifier()).asInt();
}
}
CSG_Parameters *pGrids = Get_Parameters("GRID"), Tmp;
Tmp.Assign(pGrids);
pGrids->Create(this, Tmp.Get_Name(), Tmp.Get_Description(), Tmp.Get_Identifier(), false);
m_Grid_Target.Add_Grid_Parameter(SG_T("QUALITY") , _TL("Quality") , false);
m_Grid_Target.Add_Grid_Parameter(SG_T("INTERCEPT") , _TL("Intercept"), false);
pGrids->Get_Parameter("QUALITY")->Get_Parent()->asGrid_System()->Assign(*Tmp("QUALITY")->Get_Parent()->asGrid_System());
pGrids->Get_Parameter("QUALITY") ->Set_Value(Tmp("QUALITY") ->asGrid());
pGrids->Get_Parameter("INTERCEPT")->Set_Value(Tmp("INTERCEPT")->asGrid());
for(i=0; i<m_nPredictors; i++)
{
m_Grid_Target.Add_Grid_Parameter(SG_Get_String(i, 0),
CSG_String::Format(SG_T("%s [%s]"), _TL("Slope"), pPoints->Get_Field_Name(m_iPredictor[i])), false
);
if( Tmp(SG_Get_String(i, 0)) )
{
pGrids->Get_Parameter(SG_Get_String(i, 0))->Set_Value(Tmp(SG_Get_String(i, 0))->asGrid());
}
}
return( m_nPredictors > 0 );
}
//---------------------------------------------------------
bool CSelection_Delete::On_Execute(void)
{
CSG_Shapes *pInput;
pInput = Parameters("INPUT") ->asShapes();
if( pInput->Get_Selection_Count() <= 0 )
{
Error_Set(_TL("no shapes in selection"));
return( false );
}
pInput->Del_Selection();
DataObject_Update(pInput);
return( true );
}
//---------------------------------------------------------
void CGrid_3D_Image::_Set_Shapes(CSG_Shapes *pInput)
{
int iShape, iPart, iPoint;
double x, y, z, dx, dy;
T3DPoint p;
TSG_Point Point;
CSG_Shape *pShape;
CSG_Shapes *pOutput;
if( pInput && pInput->is_Valid() )
{
Process_Set_Text("%s \"%s\"", _TL("Project"), pInput->Get_Name());
pOutput = SG_Create_Shapes(*pInput);
dx = (double)Get_NX() / Get_System().Get_XRange();
dy = (double)Get_NY() / Get_System().Get_YRange();
for(iShape=0; iShape<pOutput->Get_Count() && Set_Progress(iShape, pOutput->Get_Count()); iShape++)
{
pShape = pOutput->Get_Shape(iShape);
for(iPart=0; iPart<pShape->Get_Part_Count(); iPart++)
{
for(iPoint=0; iPoint<pShape->Get_Point_Count(iPart); iPoint++)
{
Point = pShape->Get_Point(iPoint, iPart);
x = dx * (Point.x - Get_XMin());
y = dy * (Point.y - Get_YMin());
z = m_pDEM->is_InGrid((int)x, (int)y, true) ? m_pDEM->asDouble((int)x, (int)y) : 0.0;
_Get_Position(x, y, z, p);
pShape->Set_Point(p.x, p.y, iPoint, iPart);
}
}
}
DataObject_Add(pOutput);
}
}
bool CCreateChartLayer::GetExtraParameters(){
int i;
CSG_Shapes *pInput;
CSG_Parameter *pParam;
CSG_String sName;
bool bIsValidSelection = false;
pInput = Parameters("INPUT")->asShapes();
m_pExtraParameters->Create(this, _TL("Fields for diagram"), _TL(""), SG_T("EXTRA"));
m_bIncludeParam = new bool [pInput->Get_Field_Count() ];
for (i = 0; i < pInput->Get_Field_Count(); i++)
{
switch( pInput->Get_Field_Type(i) )
{
default:
break;
case SG_DATATYPE_Byte:
case SG_DATATYPE_Char:
case SG_DATATYPE_Word:
case SG_DATATYPE_Short:
case SG_DATATYPE_DWord:
case SG_DATATYPE_Int:
case SG_DATATYPE_ULong:
case SG_DATATYPE_Long:
case SG_DATATYPE_Float:
case SG_DATATYPE_Double: // is numeric field
m_pExtraParameters->Add_Value(
NULL, SG_Get_String(i,0).c_str(), pInput->Get_Field_Name(i), _TL(""), PARAMETER_TYPE_Bool, false
);
break;
}
}//for
if(Dlg_Parameters("EXTRA")){
for (i = 0; i < pInput->Get_Field_Count(); i++){
sName = SG_Get_String(i,0);
if (pParam = Get_Parameters("EXTRA")->Get_Parameter(sName.c_str())){
m_bIncludeParam[i] = pParam->asBool();
bIsValidSelection = true;
}//try
else{
m_bIncludeParam[i] = false;
}//else
}//for
m_pExtraParameters->Destroy();
return bIsValidSelection;
}//if
m_pExtraParameters->Destroy();
return false;
}//method
//---------------------------------------------------------
bool CCollect_Points::On_Execute_Finish(void)
{
CSG_Shapes *pTarget = Parameters("REF_TARGET")->asShapes();
if( pTarget != NULL )
{
pTarget->Create(SHAPE_TYPE_Point, _TL("Reference Points (Projection)"));
pTarget->Add_Field("X_SRC", SG_DATATYPE_Double);
pTarget->Add_Field("Y_SRC", SG_DATATYPE_Double);
pTarget->Add_Field("X_MAP", SG_DATATYPE_Double);
pTarget->Add_Field("Y_MAP", SG_DATATYPE_Double);
pTarget->Add_Field("RESID", SG_DATATYPE_Double);
for(int iPoint=0; iPoint<m_pPoints->Get_Count(); iPoint++)
{
CSG_Shape *pPoint = pTarget->Add_Shape(m_pPoints->Get_Shape(iPoint), SHAPE_COPY_ATTR);
pPoint->Add_Point(
pPoint->asDouble(2),
pPoint->asDouble(3)
);
}
}
m_Engine.Destroy();
return( true );
}
