Package org.noise_planet.noisemodelling.pathfinder.utils.geometry

Class JTSUtility

java.lang.Object
org.noise_planet.noisemodelling.pathfinder.utils.geometry.JTSUtility
public class JTSUtility extends Object
Author:
Nicolas Fortin

  • Field Summary

    Fields
    Modifier and Type
    Field
    Description
    static final double
    TRIANGLE_INTERSECTION_EPSILON
    Intersection test for topology triangle navigation

  • Method Summary

    Modifier and Type
    Method
    Description
    static Double
    dist2D(org.locationtech.jts.geom.Coordinate c0, org.locationtech.jts.geom.Coordinate c1)
    compute the distance between two points in two dimensions
    static double
    dist3D(org.locationtech.jts.geom.Coordinate c0, org.locationtech.jts.geom.Coordinate c1)
    compute the distance between two points of dimension three
    static boolean
    dotInTri(org.locationtech.jts.geom.Coordinate p, org.locationtech.jts.geom.Coordinate a, org.locationtech.jts.geom.Coordinate b, org.locationtech.jts.geom.Coordinate c)
     
    static boolean
    dotInTri(org.locationtech.jts.geom.Coordinate p, org.locationtech.jts.geom.Coordinate a, org.locationtech.jts.geom.Coordinate b, org.locationtech.jts.geom.Coordinate c, AtomicReference<Double> error)
    Fast dot in triangle test http://www.blackpawn.com/texts/pointinpoly/default.html
    static double[]
    getLinearFunction(org.locationtech.jts.geom.Coordinate p1, org.locationtech.jts.geom.Coordinate p2)
    Compute a and b linear function of the line p1 p2
    static double[]
    getLinearRegressionPolyline(List<org.locationtech.jts.geom.Coordinate> xzList)
    NFS 31-133 P.69 Annex E
    static double[]
    getMeanPlaneCoefficients(org.locationtech.jts.geom.Coordinate[] profile)
    calculate the mean plane y = A.x + B for a sequence of terrain points projected on the unfolded propagation plane ; using (x,y) coordinates as in section VI.2.2.c of the JRC-2012 reference report.
    static org.locationtech.jts.geom.Coordinate
    getNearestPoint(org.locationtech.jts.geom.Coordinate from, org.locationtech.jts.geom.LineString to)
     
    static org.locationtech.jts.geom.Coordinate
    getNearestPoint(org.locationtech.jts.geom.LineSegment segment, org.locationtech.jts.geom.Coordinate p)
     
    static List<org.locationtech.jts.geom.Coordinate>
    getNewCoordinateSystem(List<org.locationtech.jts.geom.Coordinate> listPoints)
    ChangeCoordinateSystem, use original coordinate in 3D to change into a new markland in 2D with new x' computed by algorithm and y' is original height of point.
    static List<org.locationtech.jts.geom.Coordinate>
    getNewCoordinateSystem(List<org.locationtech.jts.geom.Coordinate> listPoints, double tolerance)
    ChangeCoordinateSystem, use original coordinate in 3D to change into a new markland in 2D with new x' computed by algorithm and y' is original height of point.
    static org.locationtech.jts.geom.Coordinate
    getOldCoordinateSystem(org.locationtech.jts.geom.Coordinate Point, double angle)
    ChangeCoordinateSystem, use original coordinate in 3D to change into a new markland in 2D with new x' computed by algorithm and y' is original height of point.
    static List<org.locationtech.jts.geom.Coordinate>
    getOldCoordinateSystemList(List<org.locationtech.jts.geom.Coordinate> listPoints, double angle)
    ChangeCoordinateSystem, use original coordinate in 3D to change into a new markland in 2D with new x' computed by algorithm and y' is original height of point.
    static double
    getSlope(org.locationtech.jts.geom.Coordinate p0, org.locationtech.jts.geom.Coordinate p1)
     
    static org.locationtech.jts.math.Vector3D
    getTriangleNormal(org.locationtech.jts.geom.Coordinate p1, org.locationtech.jts.geom.Coordinate p2, org.locationtech.jts.geom.Coordinate p3)
     
    static List<org.locationtech.jts.geom.Coordinate>
    getXAscendingHullPoints(org.locationtech.jts.geom.Coordinate[] coordinates)
     
    static org.locationtech.jts.geom.Coordinate
    makePointImage(double a, double b, org.locationtech.jts.geom.Coordinate point)
    Make image (like line y=ax+b would be a mirror) of point
    static org.locationtech.jts.geom.Coordinate
    makeProjectedPoint(double a, double b, org.locationtech.jts.geom.Coordinate point)
    Project point on y=ax+b plane

    Methods inherited from class java.lang.Object

    clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

  • Field Details

    • TRIANGLE_INTERSECTION_EPSILON

      public static final double TRIANGLE_INTERSECTION_EPSILON
      Intersection test for topology triangle navigation
      See Also:

  • Method Details

    • getLinearFunction

      public static double[] getLinearFunction(org.locationtech.jts.geom.Coordinate p1, org.locationtech.jts.geom.Coordinate p2)
      Compute a and b linear function of the line p1 p2
      Parameters:
      p1 - p1
      p2 - p2 with p2.x != p1.x
      Returns:
      [a,b] linear function parameters.

    • getNearestPoint

      public static org.locationtech.jts.geom.Coordinate getNearestPoint(org.locationtech.jts.geom.LineSegment segment, org.locationtech.jts.geom.Coordinate p)
      Parameters:
      segment -
      p -
      Returns:

    • getNearestPoint

      public static org.locationtech.jts.geom.Coordinate getNearestPoint(org.locationtech.jts.geom.Coordinate from, org.locationtech.jts.geom.LineString to)
      Parameters:
      from -
      to -
      Returns:

    • getLinearRegressionPolyline

      public static double[] getLinearRegressionPolyline(List<org.locationtech.jts.geom.Coordinate> xzList)
      NFS 31-133 P.69 Annex E
      Parameters:
      xzList - Line coordinates in the same plan of the line formed by the first and the last point. X must be incremental.
      Returns:
      [a,b] Linear function parameters produced by least square regression of provided points.

    • makeProjectedPoint

      public static org.locationtech.jts.geom.Coordinate makeProjectedPoint(double a, double b, org.locationtech.jts.geom.Coordinate point)
      Project point on y=ax+b plane
      Parameters:
      a - a plane parameter
      b - b plane parameter
      point - Point to project
      Returns:
      Projected point

    • makePointImage

      public static org.locationtech.jts.geom.Coordinate makePointImage(double a, double b, org.locationtech.jts.geom.Coordinate point)
      Make image (like line y=ax+b would be a mirror) of point
      Parameters:
      a - a linear parameter
      b - b linear parameter
      point - Point to project
      Returns:
      Mirrored point position

    • dotInTri

      public static boolean dotInTri(org.locationtech.jts.geom.Coordinate p, org.locationtech.jts.geom.Coordinate a, org.locationtech.jts.geom.Coordinate b, org.locationtech.jts.geom.Coordinate c)

    • getTriangleNormal

      public static org.locationtech.jts.math.Vector3D getTriangleNormal(org.locationtech.jts.geom.Coordinate p1, org.locationtech.jts.geom.Coordinate p2, org.locationtech.jts.geom.Coordinate p3)

    • dotInTri

      public static boolean dotInTri(org.locationtech.jts.geom.Coordinate p, org.locationtech.jts.geom.Coordinate a, org.locationtech.jts.geom.Coordinate b, org.locationtech.jts.geom.Coordinate c, AtomicReference<Double> error)
      Fast dot in triangle test http://www.blackpawn.com/texts/pointinpoly/default.html
      Parameters:
      p - coordinate of the point
      a - coordinate of the A vertex of triangle
      b - coordinate of the B vertex of triangle
      c - coordinate of the C vertex of triangle
      Returns:
      True if dot is in triangle

    • getNewCoordinateSystem

      public static List<org.locationtech.jts.geom.Coordinate> getNewCoordinateSystem(List<org.locationtech.jts.geom.Coordinate> listPoints)
      ChangeCoordinateSystem, use original coordinate in 3D to change into a new markland in 2D with new x' computed by algorithm and y' is original height of point.
      Parameters:
      listPoints - X Y Z points, all should be on the same plane as first and last points.
      Returns:
      X Z projected points

    • getNewCoordinateSystem

      public static List<org.locationtech.jts.geom.Coordinate> getNewCoordinateSystem(List<org.locationtech.jts.geom.Coordinate> listPoints, double tolerance)
      ChangeCoordinateSystem, use original coordinate in 3D to change into a new markland in 2D with new x' computed by algorithm and y' is original height of point.
      Parameters:
      listPoints - X Y Z points, all should be on the same plane as first and last points.
      tolerance - Simplify the point list by not adding points where the distance from the line segments formed from the previous and the next point is inferior to this tolerance (remove intermediate collinear points)
      Returns:
      X Z projected points

    • getOldCoordinateSystemList

      public static List<org.locationtech.jts.geom.Coordinate> getOldCoordinateSystemList(List<org.locationtech.jts.geom.Coordinate> listPoints, double angle)
      ChangeCoordinateSystem, use original coordinate in 3D to change into a new markland in 2D with new x' computed by algorithm and y' is original height of point. Attention this function can just be used when the points in the same plane. "http://en.wikipedia.org/wiki/Rotation_matrix"
      Parameters:
      listPoints - X Y Z points, all should be on the same plane as first and last points.
      Returns:
      X Z projected points

    • getOldCoordinateSystem

      public static org.locationtech.jts.geom.Coordinate getOldCoordinateSystem(org.locationtech.jts.geom.Coordinate Point, double angle)
      ChangeCoordinateSystem, use original coordinate in 3D to change into a new markland in 2D with new x' computed by algorithm and y' is original height of point. Attention this function can just be used when the points in the same plane. "http://en.wikipedia.org/wiki/Rotation_matrix"
      Parameters:
      Point - X Y Z points, all should be on the same plane as first and last points.
      Returns:
      X Z projected points

    • getMeanPlaneCoefficients

      public static double[] getMeanPlaneCoefficients(org.locationtech.jts.geom.Coordinate[] profile)
      calculate the mean plane y = A.x + B for a sequence of terrain points projected on the unfolded propagation plane ; using (x,y) coordinates as in section VI.2.2.c of the JRC-2012 reference report.
      Parameters:
      profile - u v coordinates @see {getNewCoordinateSystem(List)}
      Returns:
      Coefficient A and B

    • getXAscendingHullPoints

      public static List<org.locationtech.jts.geom.Coordinate> getXAscendingHullPoints(org.locationtech.jts.geom.Coordinate[] coordinates)
      Parameters:
      coordinates - Coordinates
      Returns:
      Parts of the clock-wise ConvexHull where x value are increasing from the minimum X value

    • dist3D

      public static double dist3D(org.locationtech.jts.geom.Coordinate c0, org.locationtech.jts.geom.Coordinate c1)
      compute the distance between two points of dimension three
      Parameters:
      c0 -
      c1 -
      Returns:
      the distance in double

    • dist2D

      public static Double dist2D(org.locationtech.jts.geom.Coordinate c0, org.locationtech.jts.geom.Coordinate c1)
      compute the distance between two points in two dimensions
      Parameters:
      c0 -
      c1 -
      Returns:
      the distance in double

    • getSlope

      public static double getSlope(org.locationtech.jts.geom.Coordinate p0, org.locationtech.jts.geom.Coordinate p1)
      Parameters:
      p0 -
      p1 -
      Returns: