CoverageFunctions.java

/*
 * Copyright (c) 2022 Martin Davis.
 *
 * All rights reserved. This program and the accompanying materials
 * are made available under the terms of the Eclipse Public License 2.0
 * and Eclipse Distribution License v. 1.0 which accompanies this distribution.
 * The Eclipse Public License is available at http://www.eclipse.org/legal/epl-v20.html
 * and the Eclipse Distribution License is available at
 *
 * http://www.eclipse.org/org/documents/edl-v10.php.
 */
package org.locationtech.jtstest.function;

import java.util.Arrays;
import java.util.List;

import org.locationtech.jts.coverage.CoverageGapFinder;
import org.locationtech.jts.coverage.CoveragePolygonValidator;
import org.locationtech.jts.coverage.CoverageSimplifier;
import org.locationtech.jts.coverage.CoverageUnion;
import org.locationtech.jts.coverage.CoverageValidator;
import org.locationtech.jts.geom.Geometry;
import org.locationtech.jts.geom.util.PolygonExtracter;
import org.locationtech.jtstest.geomfunction.Metadata;

public class CoverageFunctions {
  
  public static Geometry validatePolygon(Geometry geom, Geometry adjacentPolys) {
    return CoveragePolygonValidator.validate(geom, toGeometryArray(adjacentPolys));
  }
  
  public static Geometry validatePolygonWithGaps(Geometry geom, Geometry adjacentPolys, 
      @Metadata(title="Gap width")
      double gapWidth) {
    return CoveragePolygonValidator.validate(geom, toGeometryArray(adjacentPolys), gapWidth);
  }
  
  public static Geometry validate(Geometry geom) {
    Geometry[] invalid = CoverageValidator.validate(toGeometryArray(geom));
    return FunctionsUtil.buildGeometryCollection(invalid, geom.getFactory().createLineString());
  }

  public static Geometry validateWithGaps(Geometry geom, 
      @Metadata(title="Gap width")
      double gapWidth) {
    Geometry[] invalid = CoverageValidator.validate(toGeometryArray(geom), gapWidth);
    return FunctionsUtil.buildGeometryCollection(invalid, geom.getFactory().createLineString());
  }

  public static Geometry findGaps(Geometry geom, 
      @Metadata(title="Gap width")
      double gapWidth) {
    return CoverageGapFinder.findGaps(toGeometryArray(geom),gapWidth);
  }

  @Metadata(description="Fast Union of a coverage")
  public static Geometry union(Geometry coverage) {
    Geometry[] cov = toGeometryArray(coverage);
    return CoverageUnion.union(cov);
  }
  
  @Metadata(description="Simplify a coverage")
  public static Geometry simplify(Geometry coverage, double tolerance) {
    Geometry[] cov = toGeometryArray(coverage);
    Geometry[] result =  CoverageSimplifier.simplify(cov, tolerance);
    return coverage.getFactory().createGeometryCollection(result);
  }
  
  @Metadata(description="Simplify a coverage with a smoothness weight")
  public static Geometry simplifySharp(Geometry coverage, 
      @Metadata(title="Distance tol")
      double tolerance, 
      @Metadata(title="Weight")
      double weight) {
    Geometry[] cov = toGeometryArray(coverage);
    CoverageSimplifier simplifier = new CoverageSimplifier(cov);
    simplifier.setSmoothWeight(weight);
    Geometry[] result = simplifier.simplify(tolerance);
    return coverage.getFactory().createGeometryCollection(result);
  }
  
  @Metadata(description="Simplify a coverage with a ring removal size factor")
  public static Geometry simplifyRemoveRings(Geometry coverage, 
      @Metadata(title="Distance tol")
      double tolerance, 
      @Metadata(title="Removal Size Factor")
      double factor) {
    Geometry[] cov = toGeometryArray(coverage);
    CoverageSimplifier simplifier = new CoverageSimplifier(cov);
    simplifier.setRemovableRingSizeFactor(factor);
    Geometry[] result = simplifier.simplify(tolerance);
    return coverage.getFactory().createGeometryCollection(result);
  }
  
  @Metadata(description="Simplify inner edges of a coverage")
  public static Geometry simplifyInner(Geometry coverage, double tolerance) {
    Geometry[] cov = toGeometryArray(coverage);
    Geometry[] result =  CoverageSimplifier.simplifyInner(cov, tolerance);
    return coverage.getFactory().createGeometryCollection(result);
  }
  
  @Metadata(description="Simplify outer edges of a coverage")
  public static Geometry simplifyOuter(Geometry coverage, double tolerance) {
    Geometry[] cov = toGeometryArray(coverage);
    Geometry[] result =  CoverageSimplifier.simplifyOuter(cov, tolerance);
    return coverage.getFactory().createGeometryCollection(result);
  }
  
  @Metadata(description="Simplify inner and outer edges of a coverage differently")
  public static Geometry simplifyInOut(Geometry coverage, 
      @Metadata(title="Inner Distance tol")
      double toleranceInner, 
      @Metadata(title="Outer Distance tol")
      double toleranceOuter) {
    Geometry[] cov = toGeometryArray(coverage);
    CoverageSimplifier simplifier = new CoverageSimplifier(cov);
    Geometry[] result = simplifier.simplify(toleranceInner, toleranceOuter);
    return coverage.getFactory().createGeometryCollection(result);
  }
  
  @Metadata(description="Simplify a coverage with per-geometry tolerances")
  public static Geometry simplifyTolerances(Geometry coverage, 
      @Metadata(title="Tolerances (comma-sep)")
      String tolerancesCSV) {
    Geometry[] cov = toGeometryArray(coverage);
    double[] tolerances = tolerances(tolerancesCSV, cov.length);
    Geometry[] result =  CoverageSimplifier.simplify(cov, tolerances);
    return FunctionsUtil.buildGeometry(result);
  }
  
  private static double[] tolerances(String csvList, int len) {
    Double[] tolsDouble = toDoubleArray(csvList);
    double[] tols = new double[len];
    for (int i = 0; i < tolsDouble.length; i++) {
      tols[i] = tolsDouble[i];
    }
    return tols;
  }  
  
  private static Double[] toDoubleArray(String csvList) {
    return Arrays.stream(csvList.split(",")).map(Double::parseDouble).toArray(Double[]::new);
  }
  
  private static Geometry[] toGeometryArray(Geometry geom) {
    Geometry[] geoms = new Geometry[geom.getNumGeometries()];
    for (int i = 0; i < geom.getNumGeometries(); i++) {
      geoms[i]= geom.getGeometryN(i);
    }
    return geoms;
  }
  
}