void RGISAnalysePointSubbasinStatsCBK (Widget widget, RGISWorkspace *workspace,XmAnyCallbackStruct *callData)
{
char *selection;
DBDataset *dataset = UIDataset ();
DBObjData *pntData = dataset->Data ();
DBObjData *grdData = pntData->LinkedData ();
DBObjData *netData;
DBObjData *tblData;
static Widget selectWidget = (Widget) NULL;
if (selectWidget == (Widget) NULL) selectWidget = UISelectionCreate ((char *) "Select Network");
if ((selection = UISelectObject (selectWidget,(DBObjectLIST<DBObject> *) dataset->DataList (),_RGISAnalyseTestNetwork)) == (char *) NULL) return;
if ((netData = dataset->Data (selection)) == (DBObjData *) NULL) return;
tblData = new DBObjData ("",DBTypeTable);
tblData->Document (DBDocGeoDomain,pntData->Document (DBDocGeoDomain));
tblData->Document (DBDocSubject,(char *) "Subbasin Statistics");
if (UIDataHeaderForm (tblData))
{
UIPauseDialogOpen ((char *) "Calculating Statistics");
if (RGlibPointSubbasinStats (pntData,netData,grdData,tblData) == DBSuccess)
workspace->CurrentData (tblData);
else delete tblData;
UIPauseDialogClose ();
}
else delete tblData;
return;
}
int main (int argc,char *argv [])
{
int argPos, argNum = argc, ret, verbose = false;
char *title = (char *) NULL, *subject = (char *) NULL;
char *domain = (char *) NULL, *version = (char *) NULL;
char *pointName = (char *) NULL, *networkName = (char *) NULL;
DBObjData *data, *pntData, *netData, *grdData;
for (argPos = 1;argPos < argNum; )
{
if (CMargTest (argv [argPos],"-p","--point"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing point coverage!"); return (CMfailed); }
pointName = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-n","--network"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing network coverage!"); return (CMfailed); }
networkName = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-t","--title"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing title!"); return (CMfailed); }
title = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-u","--subject"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing subject!"); return (CMfailed); }
subject = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-d","--domain"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing domain!"); return (CMfailed); }
domain = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-v","--version"))
{
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos)
{ CMmsgPrint (CMmsgUsrError,"Missing version!"); return (CMfailed); }
version = argv [argPos];
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-V","--verbose"))
{
verbose = true;
if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break;
continue;
}
if (CMargTest (argv [argPos],"-h","--help"))
{
CMmsgPrint (CMmsgInfo,"%s [options] <input grid> <output table>",CMfileName(argv[0]));
CMmsgPrint (CMmsgInfo," -p,--point [point coverage]");
CMmsgPrint (CMmsgInfo," -n,--network [network coverage]");
CMmsgPrint (CMmsgInfo," -t,--title [dataset title]");
CMmsgPrint (CMmsgInfo," -u,--subject [subject]");
CMmsgPrint (CMmsgInfo," -d,--domain [domain]");
CMmsgPrint (CMmsgInfo," -v,--version [version]");
CMmsgPrint (CMmsgInfo," -V,--verbose");
CMmsgPrint (CMmsgInfo," -h,--help");
return (DBSuccess);
}
if ((argv [argPos][0] == '-') && (strlen (argv [argPos]) > 1))
{ CMmsgPrint (CMmsgUsrError,"Unknown option: %s!",argv [argPos]); return (CMfailed); }
argPos++;
}
if (argNum > 3) { CMmsgPrint (CMmsgUsrError,"Extra arguments!"); return (CMfailed); }
if (verbose) RGlibPauseOpen (argv[0]);
if (pointName == (char *) NULL)
{ CMmsgPrint (CMmsgUsrError,"Point coverage is not specified"); return (CMfailed); }
if (networkName == (char *) NULL)
{ CMmsgPrint (CMmsgUsrError,"Network coverage is not specified"); return (CMfailed); }
pntData = new DBObjData ();
if ((pntData->Read (pointName) == DBFault) || (pntData->Type () != DBTypeVectorPoint))
{ delete pntData; return (CMfailed); }
netData = new DBObjData ();
if ((netData->Read (networkName) == DBFault) || (netData->Type () != DBTypeNetwork))
{ delete pntData; delete netData; return (CMfailed); }
grdData = new DBObjData ();
ret = (argNum > 1) && (strcmp (argv [1],"-") != 0) ? grdData->Read (argv [1]) : grdData->Read (stdin);
if ((ret == DBFault) || (grdData->Type () != DBTypeGridContinuous))
{ delete pntData, delete netData; delete grdData; return (CMfailed); }
if (title == (char *) NULL) title = (char *) "Subbasin Histogram";
if (subject == (char *) NULL) subject = grdData->Document (DBDocSubject);
if (domain == (char *) NULL) domain = netData->Document (DBDocGeoDomain);
if (version == (char *) NULL) version = (char *) "0.01pre";
data = new DBObjData (title,DBTypeTable);
data->Document (DBDocSubject,subject);
data->Document (DBDocGeoDomain, domain);
data->Document (DBDocVersion, version);
if ((ret = RGlibPointSubbasinStats (pntData,netData,grdData,data)) == DBSuccess)
ret = (argNum > 2) && (strcmp (argv [2],"-") != 0) ? data->Write (argv [2]) : data->Write (stdout);
delete netData; delete pntData; delete grdData; delete data;
if (verbose) RGlibPauseClose ();
return (ret);
}
