Examples with OmsExtractBasin org.hortonmachine.hmachine.modules.demmanipulation.wateroutlet.OmsExtractBasin used on opensource projects

Example 1 with OmsExtractBasin

use of org.hortonmachine.hmachine.modules.demmanipulation.wateroutlet.OmsExtractBasin in project hortonmachine by TheHortonMachine.

the class ExtractBasin method process.

@Execute
public void process() throws Exception {
    OmsExtractBasin extractbasin = new OmsExtractBasin();
    extractbasin.pNorth = pNorth;
    extractbasin.pEast = pEast;
    extractbasin.inFlow = getRaster(inFlow);
    extractbasin.inNetwork = getVector(inNetwork);
    extractbasin.pSnapbuffer = pSnapbuffer;
    boolean doVector = false;
    if (outVectorBasin != null) {
        doVector = true;
    }
    extractbasin.doVector = doVector;
    extractbasin.doSmoothing = doSmoothing;
    extractbasin.pm = pm;
    extractbasin.doProcess = doProcess;
    extractbasin.doReset = doReset;
    extractbasin.process();
    dumpRaster(extractbasin.outBasin, outBasin);
    dumpVector(extractbasin.outOutlet, outOutlet);
    dumpVector(extractbasin.outVectorBasin, outVectorBasin);
    outArea = extractbasin.outArea;
    pm.message("Basin Area = " + extractbasin.outArea);
}

Example 2 with OmsExtractBasin

use of org.hortonmachine.hmachine.modules.demmanipulation.wateroutlet.OmsExtractBasin in project hortonmachine by TheHortonMachine.

the class RunMaximumDischarge method main.

public static void main(String[] args) throws Exception {
    String dtm = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/dtm.tiff";
    String aspect = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/aspect.tiff";
    String pit = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/pit.tiff";
    String flow = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/flow.tiff";
    String drain = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/drain.tiff";
    String tca = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/tca.tiff";
    String net = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/net300.tiff";
    String basin = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/basin.tiff";
    String basinShp = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/basin_vect.shp";
    String basinOutletShp = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/basin_outlet.shp";
    String basinPit = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/basin_cutout.tiff";
    String slope = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/slope.tiff";
    String topindex = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/topindex.tiff";
    String curvaturesPlan = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/curvatures_plan.tiff";
    String curvaturesProf = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/curvatures_prof.tiff";
    String curvaturesTan = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/curvatures_tan.tiff";
    String alung = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/alung.tiff";
    String alungB = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/alungB.tiff";
    String rescaled = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/flanginec/rescaled.tiff";
    // String dtm = "/media/hydrologis/Samsung_T3/MAZONE/DTM/dtm_toblino/dtm_toblino.tiff";
    // String pit = "/media/hydrologis/Samsung_T3/MAZONE/DTM/dtm_toblino/pit.tiff";
    // String flow = "/media/hydrologis/Samsung_T3/MAZONE/DTM/dtm_toblino/flow.tiff";
    // String drain = "/media/hydrologis/Samsung_T3/MAZONE/DTM/dtm_toblino/drain.tiff";
    // String tca = "/media/hydrologis/Samsung_T3/MAZONE/DTM/dtm_toblino/tca.tiff";
    // String dtm =
    // "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/DTM_calvello/dtm_all_float.tiff";
    // String pit = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/DTM_calvello/pit.tiff";
    // String flow = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/DTM_calvello/flow.tiff";
    // String drain = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/DTM_calvello/drain.tiff";
    // String tca = "/media/hydrologis/Samsung_T3/MAZONE/PITFILLE/DTM_calvello/tca.tiff";
    // OmsAspect aspectMod = new OmsAspect();
    // aspectMod.inElev = OmsRasterReader.readRaster(dtm);
    // // aspectMod.doRadiants = doRadiants;
    // aspectMod.doRound = true;
    // aspectMod.process();
    // OmsRasterWriter.writeRaster(aspect, aspectMod.outAspect);
    OmsDePitter pitfiller = new OmsDePitter();
    pitfiller.inElev = OmsRasterReader.readRaster("/home/hydrologis/Dropbox/hydrologis/lavori/2017_06_mapzone/test/dtm_test2.tiff");
    pitfiller.process();
    OmsRasterWriter.writeRaster(pit, pitfiller.outPit);
    OmsRasterWriter.writeRaster(flow, pitfiller.outFlow);
    if (true) {
        System.exit(0);
    }
    OmsDrainDir draindir = new OmsDrainDir();
    draindir.inPit = OmsRasterReader.readRaster(pit);
    draindir.inFlow = OmsRasterReader.readRaster(flow);
    // draindir.inFlownet = OmsRasterReader.readRaster(inFlownet);
    draindir.pLambda = 1f;
    draindir.process();
    OmsRasterWriter.writeRaster(drain, draindir.outFlow);
    OmsRasterWriter.writeRaster(tca, draindir.outTca);
    OmsExtractNetwork extractnetwork = new OmsExtractNetwork();
    extractnetwork.inTca = OmsRasterReader.readRaster(tca);
    extractnetwork.inFlow = OmsRasterReader.readRaster(drain);
    // extractnetwork.inSlope = OmsRasterReader.readRaster(inSlope);
    // extractnetwork.inTc3 = OmsRasterReader.readRaster(inTc3);
    extractnetwork.pThres = 300;
    extractnetwork.pMode = Variables.TCA;
    // extractnetwork.pExp = pExp;
    // extractnetwork.pm = pm;
    extractnetwork.process();
    OmsRasterWriter.writeRaster(net, extractnetwork.outNet);
    OmsExtractBasin extractbasin = new OmsExtractBasin();
    extractbasin.pNorth = 5112705.635839384;
    extractbasin.pEast = 1638565.1058415675;
    extractbasin.inFlow = OmsRasterReader.readRaster(drain);
    // extractbasin.inNetwork = getVector(inNetwork);
    // extractbasin.pSnapbuffer = pSnapbuffer;
    extractbasin.doVector = true;
    // extractbasin.doSmoothing = doSmoothing;
    extractbasin.process();
    OmsRasterWriter.writeRaster(basin, extractbasin.outBasin);
    OmsVectorWriter.writeVector(basinOutletShp, extractbasin.outOutlet);
    OmsVectorWriter.writeVector(basinShp, extractbasin.outVectorBasin);
    System.out.println("Basin Area = " + extractbasin.outArea);
    OmsCutOut cutPit = new OmsCutOut();
    cutPit.inRaster = OmsRasterReader.readRaster(pit);
    cutPit.inMask = OmsRasterReader.readRaster(basin);
    // c.pMax = pMax;
    // c.pMin = pMin;
    // c.doInverse = doInverse;
    cutPit.process();
    OmsRasterWriter.writeRaster(basinPit, cutPit.outRaster);
    OmsGradient gradient = new OmsGradient();
    gradient.inElev = OmsRasterReader.readRaster(basinPit);
    gradient.pMode = Variables.FINITE_DIFFERENCES;
    gradient.doDegrees = true;
    gradient.process();
    OmsRasterWriter.writeRaster(slope, gradient.outSlope);
    OmsTopIndex topindexMod = new OmsTopIndex();
    topindexMod.inTca = OmsRasterReader.readRaster(tca);
    topindexMod.inSlope = OmsRasterReader.readRaster(slope);
    topindexMod.process();
    OmsRasterWriter.writeRaster(topindex, topindexMod.outTopindex);
    OmsCurvatures curv = new OmsCurvatures();
    curv.inElev = OmsRasterReader.readRaster(basinPit);
    curv.process();
    OmsRasterWriter.writeRaster(curvaturesProf, curv.outProf);
    OmsRasterWriter.writeRaster(curvaturesPlan, curv.outPlan);
    OmsRasterWriter.writeRaster(curvaturesTan, curv.outTang);
    OmsAb ab = new OmsAb();
    ab.inTca = OmsRasterReader.readRaster(tca);
    ab.inPlan = OmsRasterReader.readRaster(curvaturesPlan);
    ab.process();
    OmsRasterWriter.writeRaster(alung, ab.outAb);
    OmsRasterWriter.writeRaster(alungB, ab.outB);
    GridCoverage2D cutDrain = OmsCutOut.cut(OmsRasterReader.readRaster(drain), OmsRasterReader.readRaster(basin));
    GridCoverage2D cutNet = OmsCutOut.cut(OmsRasterReader.readRaster(net), OmsRasterReader.readRaster(basin));
    OmsRescaledDistance rescaleddistance = new OmsRescaledDistance();
    rescaleddistance.inFlow = cutDrain;
    rescaleddistance.inNet = cutNet;
    rescaleddistance.inElev = OmsRasterReader.readRaster(basinPit);
    rescaleddistance.pRatio = 0.3;
    rescaleddistance.process();
    OmsRasterWriter.writeRaster(rescaled, rescaleddistance.outRescaled);
}

Example 3 with OmsExtractBasin

use of org.hortonmachine.hmachine.modules.demmanipulation.wateroutlet.OmsExtractBasin in project hortonmachine by TheHortonMachine.

the class TestWateroutlet method testExtractBasin.

public void testExtractBasin() throws Exception {
    double[][] flowData = HMTestMaps.flowData;
    HashMap<String, Double> envelopeParams = HMTestMaps.getEnvelopeparams();
    CoordinateReferenceSystem crs = HMTestMaps.getCrs();
    GridCoverage2D flowCoverage = CoverageUtilities.buildCoverage("flow", flowData, envelopeParams, crs, true);
    OmsExtractBasin extractBasin = new OmsExtractBasin();
    extractBasin.pm = pm;
    extractBasin.inFlow = flowCoverage;
    extractBasin.pNorth = 5139885.0;
    extractBasin.pEast = 1640724.0;
    extractBasin.process();
    GridCoverage2D basinCoverage = extractBasin.outBasin;
    // System.out.println(extractBasin.outArea);
    // PrintUtilities.printCoverageData(basinCoverage);
    checkMatrixEqual(basinCoverage.getRenderedImage(), HMTestMaps.basinWateroutletData, 0);
}