void RGISAnNetworkAccumulateCBK (Widget widget,RGISWorkspace *workspace,XmAnyCallbackStruct *callData) { DBDataset *dataset = UIDataset (); DBObjData *netData = dataset->Data (); DBObjData *inGridData, *outGridData; if ((inGridData = netData->LinkedData ()) == (DBObjData *) NULL) { CMmsgPrint (CMmsgAppError,"Null Linked Data in: %s %d",__FILE__,__LINE__); return; } outGridData = DBNetworkToGrid (netData,DBTypeGridContinuous); outGridData->Document (DBDocSubject,inGridData->Document (DBDocSubject)); if (UIDataHeaderForm (outGridData)) { UIPauseDialogOpen ((char *) "Accumulating Grid"); if (RGlibNetworkAccumulate (netData,inGridData,outGridData) == DBSuccess) workspace->CurrentData (outGridData); else delete outGridData; UIPauseDialogClose (); } else delete outGridData; }
int main (int argc,char *argv []) { int argPos, argNum = argc, ret, verbose = false; float coeff = 0.000001; char *title = (char *) NULL, *subject = (char *) NULL; char *domain = (char *) NULL, *version = (char *) NULL; char *netName = (char *) NULL, *stnName = (char *) NULL; char *tsName = (char *) NULL; char *fields [5] = { (char *) NULL, (char *) NULL, (char *) NULL, (char *) NULL, (char *) NULL }; DBInt shadeSet = DBFault; bool areaMult = true, correction = false, allowNegative = true; DBObjData *data, *netData, *grdData, *stnData = (DBObjData *) NULL, *tsData = (DBObjData *) NULL; for (argPos = 1;argPos < argNum; ) { if (CMargTest (argv [argPos],"-n","--network")) { if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing network!"); return (CMfailed); } netName = argv [argPos]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-a","--stations")) { if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing monitoring stations!"); return (CMfailed); } stnName = argv [argPos]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-i","--timeseries")) { if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing time series!"); return (CMfailed); } tsName = argv [argPos]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-r","--relate")) { if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing relate field!"); return (CMfailed); } fields [0] = argv [argPos]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-x","--nextstation")) { if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing next station field!"); return (CMfailed); } fields [1] = argv [argPos]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-j","--join")) { if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing join field!"); return (CMfailed); } fields [2] = argv [argPos]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-e","--date")) { if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing date field!"); return (CMfailed); } fields [3] = argv [argPos]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-g","--discharge")) { if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing discharge field!"); return (CMfailed); } fields [4] = argv [argPos]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-c","--correction")) { int corCodes [] = { true, false }; const char *corModes [] = { "on", "off", (char *) NULL }; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing correcition mode!"); return (CMfailed); } if ((correction = CMoptLookup (corModes,argv [argPos],true)) == DBFault) { CMmsgPrint (CMmsgUsrError,"Invalid correction mode!"); return (CMfailed); } correction = corCodes [correction]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-m","--mode")) { int modeCodes [] = { true, false }; const char *modes [] = { "rate", "volume", (char *) NULL }; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing accumulation mode!"); return (CMfailed); } if ((areaMult = CMoptLookup (modes,argv [argPos],true)) == DBFault) { CMmsgPrint (CMmsgUsrError,"Invalid accumulation mode!"); return (CMfailed); } areaMult = modeCodes [areaMult]; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-f","--coefficient")) { if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing coefficient!"); return (CMfailed); } if (sscanf (argv [argPos],"%f",&coeff) != 1) { CMmsgPrint (CMmsgUsrError,"Invalid coefficient"); return (CMfailed); } if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) break; continue; } if (CMargTest (argv [argPos],"-l","--negative")) { int negCodes [] = { true, false }; const char *modes [] = { "allow", "eliminate", (char *) NULL }; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing negative mode!"); return (CMfailed); } if ((allowNegative = CMoptLookup (modes,argv [argPos],true)) == DBFault) { CMmsgPrint (CMmsgUsrError,"Invalid negative mode!"); return (CMfailed); } allowNegative = negCodes [allowNegative]; 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],"-s","--shadeset")) { int shadeCodes [] = { DBDataFlagDispModeContStandard, DBDataFlagDispModeContGreyScale, DBDataFlagDispModeContBlueScale, DBDataFlagDispModeContBlueRed, DBDataFlagDispModeContElevation }; const char *shadeSets [] = { "standard","grey","blue","blue-to-red","elevation", (char *) NULL }; if ((argNum = CMargShiftLeft (argPos,argv,argNum)) <= argPos) { CMmsgPrint (CMmsgUsrError,"Missing shadeset!"); return (CMfailed); } if ((shadeSet = CMoptLookup (shadeSets,argv [argPos],true)) == DBFault) { CMmsgPrint (CMmsgUsrError,"Invalid shadeset!"); return (CMfailed); } shadeSet = shadeCodes [shadeSet]; 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 grid>",CMprgName(argv[0])); CMmsgPrint (CMmsgInfo," -n,--network [network coverage]"); CMmsgPrint (CMmsgInfo," -a,--stations [station coverage]"); CMmsgPrint (CMmsgInfo," -i,--timeseries [time series data]"); CMmsgPrint (CMmsgInfo," -r,--relate [relate field]"); CMmsgPrint (CMmsgInfo," -j,--join [join field]"); CMmsgPrint (CMmsgInfo," -e,--date [date field]"); CMmsgPrint (CMmsgInfo," -x,--nextstation [next station field]"); CMmsgPrint (CMmsgInfo," -g,--discharge [discharge field]"); CMmsgPrint (CMmsgInfo," -c,--correction [on|off (default is off)]"); CMmsgPrint (CMmsgInfo," -m,--mode [rate|volume]"); CMmsgPrint (CMmsgInfo," -f,--coefficient [conversion coefficient]"); CMmsgPrint (CMmsgInfo," -l,--negative [allow|eliminate]"); CMmsgPrint (CMmsgInfo," -t,--title [dataset title]"); CMmsgPrint (CMmsgInfo," -u,--subject [subject]"); CMmsgPrint (CMmsgInfo," -d,--domain [domain]"); CMmsgPrint (CMmsgInfo," -v,--version [version]"); CMmsgPrint (CMmsgInfo," -s,--shadeset [standard|grey|blue|blue-to-red|elevation]"); 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 (netName == (char *) NULL) { CMmsgPrint (CMmsgUsrError,"Network is not specified"); return (CMfailed); } if ((stnName == (char *) NULL) && (tsName != (char *) NULL)) { CMmsgPrint (CMmsgUsrError,"Missing station coverage!"); return (CMfailed); } if ((stnName != (char *) NULL) && (tsName == (char *) NULL)) { CMmsgPrint (CMmsgUsrError,"Missing time series data!"); return (CMfailed); } netData = new DBObjData (); if ((netData->Read (netName) == DBFault) || (netData->Type () != DBTypeNetwork)) { delete netData; return (CMfailed); } if (stnName != (char *) NULL) { stnData = new DBObjData (); if (stnData->Read (stnName) == DBFault) { delete netData; delete stnData; return (CMfailed); } } if (tsName != (char *) NULL) { tsData = new DBObjData (); if (tsData->Read (tsName) == DBFault) { delete netData; delete stnData; delete tsData; 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 netData; delete grdData; return (CMfailed); } if (title == (char *) NULL) title = (char *) "Accumulated Grid"; 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 = DBNetworkToGrid (netData,DBTypeGridContinuous); data->Name (title); data->Document (DBDocSubject,subject); data->Document (DBDocGeoDomain,domain); data->Document (DBDocVersion,version); if (shadeSet != DBFault) { data->Flags (DBDataFlagDispModeContShadeSets,DBClear); data->Flags (shadeSet,DBSet); } ret = RGlibNetworkAccumulate (netData,grdData,stnData,tsData, fields,(DBFloat) coeff, areaMult, correction, allowNegative,data); if (ret == DBSuccess) ret = (argNum > 2) && (strcmp (argv [2],"-") != 0) ? data->Write (argv [2]) : data->Write (stdout); delete data; delete netData; delete grdData; if (verbose) RGlibPauseClose (); return (ret); }