RingSelfIntCheck.java 4.76 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
/*-
 *  Copyright 2020 Beuth Hochschule für Technik Berlin, Hochschule für Technik Stuttgart
 * 
 *  This file is part of CityDoctor2.
 *
 *  CityDoctor2 is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU Lesser General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  CityDoctor2 is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with CityDoctor2.  If not, see <https://www.gnu.org/licenses/>.
 */
package de.hft.stuttgart.citydoctor2.checks.geometry;

import java.util.ArrayList;
import java.util.Collections;
import java.util.List;
import java.util.Map;

import de.hft.stuttgart.citydoctor2.check.Check;
import de.hft.stuttgart.citydoctor2.check.CheckError;
import de.hft.stuttgart.citydoctor2.check.CheckId;
import de.hft.stuttgart.citydoctor2.check.CheckResult;
import de.hft.stuttgart.citydoctor2.check.CheckType;
import de.hft.stuttgart.citydoctor2.check.ResultStatus;
import de.hft.stuttgart.citydoctor2.check.error.PointTouchesEdgeError;
Matthias Betz's avatar
Matthias Betz committed
33
import de.hft.stuttgart.citydoctor2.check.error.RingEdgeIntersectionError;
Matthias Betz's avatar
Matthias Betz committed
34
import de.hft.stuttgart.citydoctor2.checks.Checks;
35
36
37
38
39
40
41
42
43
import de.hft.stuttgart.citydoctor2.datastructure.Edge;
import de.hft.stuttgart.citydoctor2.datastructure.Geometry;
import de.hft.stuttgart.citydoctor2.datastructure.LinearRing;
import de.hft.stuttgart.citydoctor2.datastructure.Vertex;
import de.hft.stuttgart.citydoctor2.math.DistanceResult;
import de.hft.stuttgart.citydoctor2.math.Segment3d;
import de.hft.stuttgart.citydoctor2.parser.ParserConfiguration;

public class RingSelfIntCheck extends Check {
Matthias Betz's avatar
Matthias Betz committed
44

45
46
47
48
49
50
51
52
53
54
55
56
	private static final String EPSILON_NAME = "minVertexDistance";

	private static final List<CheckId> dependencies;

	static {
		ArrayList<CheckId> deps = new ArrayList<>();
		deps.add(CheckId.C_GE_R_TOO_FEW_POINTS);
		deps.add(CheckId.C_GE_R_NOT_CLOSED);
		deps.add(CheckId.C_GE_R_DUPLICATE_POINT);
		dependencies = Collections.unmodifiableList(deps);
	}

Matthias Betz's avatar
Matthias Betz committed
57
	private double epsilon = Checks.MIN_VERTEX_DISTANCE_DEFAULT;
58
59
60
61

	@Override
	public void init(Map<String, String> parameters, ParserConfiguration config) {
		String epsilonString = parameters.get(EPSILON_NAME);
Matthias Betz's avatar
Matthias Betz committed
62
		if (epsilonString != null) {
63
64
65
66
67
68
69
70
71
72
73
			epsilon = Double.parseDouble(epsilonString);
		}
	}

	@Override
	public void check(LinearRing lr) {
		checkRingJava(lr);
	}

	private void checkRingJava(LinearRing lr) {
		List<Edge> edges = getEdgesForRing(lr);
Matthias Betz's avatar
Matthias Betz committed
74

75
76
77
		for (Edge e : edges) {
			if (checkForPointsTouchingEdge(lr, e)) {
				return;
Matthias Betz's avatar
Matthias Betz committed
78
			}
79
		}
Matthias Betz's avatar
Matthias Betz committed
80

81
82
		for (int i = 0; i < edges.size() - 1; i++) {
			Edge e1 = edges.get(i);
Matthias Betz's avatar
Matthias Betz committed
83

84
85
86
87
88
89
90
91
92
93
			for (int j = i + 1; j < edges.size(); j++) {
				Edge e2 = edges.get(j);
				if (e1.getConnectionPoint(e2) != null) {
					continue;
				}
				Segment3d s1 = new Segment3d(e1.getFrom(), e1.getTo());
				Segment3d s2 = new Segment3d(e2.getFrom(), e2.getTo());
				DistanceResult dr = s1.getDistanceResult(s2);
				if (dr.getDistance() < epsilon) {
					// intersection
Matthias Betz's avatar
Matthias Betz committed
94
					CheckError err = new RingEdgeIntersectionError(lr, e1, e2, dr.getPoint1());
95
96
97
98
99
100
					CheckResult cr = new CheckResult(this, ResultStatus.ERROR, err);
					lr.addCheckResult(cr);
					return;
				}
			}
		}
Matthias Betz's avatar
Matthias Betz committed
101

102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
		// no errors detected
		CheckResult cr = new CheckResult(this, ResultStatus.OK, null);
		lr.addCheckResult(cr);
	}

	private boolean checkForPointsTouchingEdge(LinearRing lr, Edge e1) {
		Segment3d seg = new Segment3d(e1.getFrom(), e1.getTo());
		for (Vertex v : lr.getVertices()) {
			if (v != e1.getFrom() && v != e1.getTo() && seg.getDistance(v) < epsilon) {
				CheckError err = new PointTouchesEdgeError(lr, e1, v);
				CheckResult cr = new CheckResult(this, ResultStatus.ERROR, err);
				lr.addCheckResult(cr);
				return true;
			}
		}
		return false;
	}

	private List<Edge> getEdgesForRing(LinearRing lr) {
		List<Edge> edges = new ArrayList<>();
		Geometry geom = lr.getParent().getParent();
		for (int i = 0; i < lr.getVertices().size() - 1; i++) {
			Vertex v1 = lr.getVertices().get(i);
			Vertex v2 = lr.getVertices().get(i + 1);
			Edge e = geom.getEdge(v1, v2);
			if (e == null) {
				throw new IllegalStateException("Edge between v1=" + v1 + " to v2=" + v2 + " is missing");
			}
			edges.add(e);
		}
		return edges;
	}

	@Override
	public List<CheckId> getDependencies() {
		return dependencies;
	}

	@Override
	public CheckType getType() {
		return CheckType.GEOMETRY;
	}

	@Override
	public Check createNewInstance() {
		return new RingSelfIntCheck();
	}
Matthias Betz's avatar
Matthias Betz committed
149
150
151
152
153

	@Override
	public CheckId getCheckId() {
		return CheckId.C_GE_R_SELF_INTERSECTION;
	}
154
}