3
* The JTS Topology Suite is a collection of Java classes that
4
* implement the fundamental operations required to validate a given
5
* geo-spatial data set to a known topological specification.
7
* Copyright (C) 2001 Vivid Solutions
9
* This library is free software; you can redistribute it and/or
10
* modify it under the terms of the GNU Lesser General Public
11
* License as published by the Free Software Foundation; either
12
* version 2.1 of the License, or (at your option) any later version.
14
* This library is distributed in the hope that it will be useful,
15
* but WITHOUT ANY WARRANTY; without even the implied warranty of
16
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
17
* Lesser General Public License for more details.
19
* You should have received a copy of the GNU Lesser General Public
20
* License along with this library; if not, write to the Free Software
21
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
23
* For more information, contact:
27
* 2328 Government Street
32
* www.vividsolutions.com
34
package com.vividsolutions.jtsexample.precision;
38
import com.vividsolutions.jts.io.WKTReader;
39
import com.vividsolutions.jts.geom.*;
40
import com.vividsolutions.jts.precision.EnhancedPrecisionOp;
43
* Example of using {@link EnhancedPrecisionOp} to avoid robustness problems
47
public class EnhancedPrecisionOpExample
49
public static void main(String[] args) throws Exception
51
EnhancedPrecisionOpExample example = new EnhancedPrecisionOpExample();
55
catch (Exception ex) {
60
private WKTReader reader = new WKTReader();
62
public EnhancedPrecisionOpExample() {
69
// two geometries which cause robustness problems
70
wkt1 = "POLYGON ((708653.498611049 2402311.54647056, 708708.895756966 2402203.47250014, 708280.326454234 2402089.6337791, 708247.896591321 2402252.48269854, 708367.379593851 2402324.00761653, 708248.882609455 2402253.07294874, 708249.523621829 2402244.3124463, 708261.854734465 2402182.39086576, 708262.818392579 2402183.35452387, 708653.498611049 2402311.54647056))";
71
wkt2 = "POLYGON ((708258.754920656 2402197.91172757, 708257.029447455 2402206.56901508, 708652.961095455 2402312.65463437, 708657.068786251 2402304.6356364, 708258.754920656 2402197.91172757))";
72
Geometry g1 = reader.read(wkt1);
73
Geometry g2 = reader.read(wkt2);
75
System.out.println("This call to intersection will throw a topology exception due to robustness problems:");
77
Geometry result = g1.intersection(g2);
79
catch (TopologyException ex) {
83
System.out.println("Using EnhancedPrecisionOp allows the intersection to be performed with no errors:");
84
Geometry result = EnhancedPrecisionOp.intersection(g1, g2);
85
System.out.println(result);