simstadt_openlayers.js 7.38 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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
//TODO: Add zoom to extent as control
//TODO: Don't allow multiple sketch
//TODO: Add text to citygml vector
//TODO: Try to leave everything in 4326 

var reset_btn = $('#reset')[0];
var send_btn = $('#send')[0];

proj4.defs("EPSG:31467", "+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0"
		+ " +ellps=bessel +datum=potsdam +units=m +no_defs"); // http://spatialreference.org/ref/epsg/31467/proj4js/

var wgs84Sphere = new ol.Sphere(6378137);

var osm_layer = new ol.layer.Tile({
	source : new ol.source.OSM()
});

var kml_source = new ol.source.KML({
	projection : ol.proj.get('EPSG:3857'),
	url : 'data/citygml_hulls.kml',
	extractAttributes : false,
	extractStyles : false
})

var kml_layer = new ol.layer.Vector({
	source : kml_source,
	style : new ol.style.Style({
		fill : new ol.style.Fill({
			color : 'rgba(255, 255, 255, 0.2)'
		}),
		stroke : new ol.style.Stroke({
			color : '#777777',
			width : 2,
			lineDash : [ 5, 10 ]
		}),
	})
});

var intersections = new ol.source.Vector();

var intersections_layer = new ol.layer.Vector({
	source : intersections,
	style : new ol.style.Style({
		fill : new ol.style.Fill({
			color : 'rgba(255, 155, 51, 0.5)'
		})
	})
});

var map = new ol.Map({
	target : 'map',
	layers : [ osm_layer, kml_layer, intersections_layer ],
53
54
55
	interactions : ol.interaction.defaults({
		keyboard : true
	})
56
57
});

duminil's avatar
duminil committed
58
59
60
61
62
63
64
var geoJSONformat = new ol.format.GeoJSON();
kml_layer.addEventListener("change", function(event) {
	var fromJavaFX = (typeof fxapp !== 'undefined');
	map.getView().fitExtent(kml_source.getExtent(), (map.getSize()));
	kml_source.forEachFeature(function(feature) {
		feature["geoJSON"] = geoJSONformat.writeFeatureObject(feature);
		feature["area"] = feature.getGeometry().getArea();
duminil's avatar
duminil committed
65
66
67
		var project = feature.get("project");
		var name = feature.get("name");
		feature["description"] = project + ">" + name;
duminil's avatar
duminil committed
68
69
		var citygmlHere;
		if (fromJavaFX) {
70
			citygmlHere = fxapp.checkIfCityGMLSAreAvailable(project, name);
duminil's avatar
duminil committed
71
72
73
74
75
		}
		feature["available"] = citygmlHere;
	});
});

76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
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
// The features are not added to a regular vector layer/source,
// but to a feature overlay which holds a collection of features.
// This collection is passed to the modify and also the draw
// interaction, so that both can add or modify features.
var featureOverlay = new ol.FeatureOverlay({
	style : new ol.style.Style({
		fill : new ol.style.Fill({
			color : 'rgba(255, 155, 51, 0.5)'
		}),
		stroke : new ol.style.Stroke({
			color : '#ffcc33',
			width : 4
		}),
		image : new ol.style.Circle({
			radius : 5,
			fill : new ol.style.Fill({
				color : '#ffcc33'
			})
		})
	})
});
featureOverlay.setMap(map);

var selected_features = featureOverlay.getFeatures();
selected_features.on('add', function(event) {
	var feature = event.element;
	feature.on("change", function(event) {
		displayInfo();
	});
});

var modify = new ol.interaction.Modify({
	features : featureOverlay.getFeatures(),
	// the SHIFT key must be pressed to delete vertices, so
	// that new vertices can be drawn at the same position
	// of existing vertices
	deleteCondition : function(event) {
		return ol.events.condition.shiftKeyOnly(event) && ol.events.condition.singleClick(event);
	}
});
map.addInteraction(modify);

draw = new ol.interaction.Draw({
	features : featureOverlay.getFeatures(),
	type : 'Polygon'
});
map.addInteraction(draw);

var sketch;
draw.on('drawstart', function(evt) {
	sketch = evt.feature;
	reset_btn.disabled = false;
});
var sourceProj = map.getView().getProjection();

function findIntersections() {
	var sketch_area = sketch.getGeometry().getArea();
	var poly1 = geoJSONformat.writeFeatureObject(sketch);
	var intersection_found = false;
	intersections.clear();
duminil's avatar
duminil committed
136
137
	var i = 0;
	kml_source.forEachFeature(function(feature) {
138
		try {
duminil's avatar
duminil committed
139
			var jsonIntersection = turf.intersect(poly1, feature["geoJSON"]);
140
141
142
143
144
145
146
			if (undefined != jsonIntersection) {
				if (!intersection_found) {
					$('#dataPanel').append("Intersection found with :<br/>\n");
					intersection_found = true;
				}
				var intersection = geoJSONformat.readFeature(jsonIntersection);
				var intersectionArea = intersection.getGeometry().getArea();
duminil's avatar
duminil committed
147
				var citygml_percentage = Math.round(intersectionArea / feature["area"] * 100);
148
149
150
				var sketch_percentage = Math.round(intersectionArea / sketch_area * 100);
				intersections.addFeature(intersection);

duminil's avatar
duminil committed
151
152
				var description;
				if (feature["available"]) {
duminil's avatar
duminil committed
153
154
					description = "<a href=\"#\" onclick=\"downloadRegionFromCityGML(" + i + ");return false;\">"
							+ feature["description"] + "</a>";
duminil's avatar
duminil committed
155
				} else {
duminil's avatar
duminil committed
156
					description = feature['description'];
duminil's avatar
duminil committed
157
158
159
				}

				$('#dataPanel').append(description + " (" + citygml_percentage + "%");
160
161
162
163
164
165
166
167
168
				if (sketch_percentage == 100) {
					$('#dataPanel').append(", all inside");
				}
				$('#dataPanel').append(")<br/>\n");
			}

		} catch (err) {
			console.log(feature.get('description') + " - " + err);
		}
duminil's avatar
duminil committed
169
170
		i++;
	})
171
172
173
174
	if (!intersection_found) {
		$('#dataPanel').append("No intersection found with any CityGML<br/>\n");
	}
}
duminil's avatar
duminil committed
175
176

function downloadRegionFromCityGML(i) {
177
	//TODO: Disable all links
178
	$("html").addClass("wait");
duminil's avatar
duminil committed
179
	var feature = kml_source.getFeatures()[i];
180
181
182
183
184
185
186
187
188
189
190
	// Waiting 100ms in order to let the cursor change
	setTimeout(function() {
		var start = new Date().getTime();
		var buildings_count = fxapp.downloadRegionFromCityGML(sketchAsWKT(), feature.get("project"), feature
				.get("name"));
		var end = new Date().getTime();
		var time = end - start;
		console.log('DL Execution time: ' + time);
		$('#dataPanel').append("Imported buildings : " + buildings_count);
		$("html").removeClass("wait");
	}, 100);
duminil's avatar
duminil committed
191
192
}

193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
function displayInfo() {
	var start = new Date().getTime();
	$('#dataPanel').empty();
	var geom = /** @type {ol.geom.Polygon} */
	(sketch.getGeometry().clone().transform(sourceProj, 'EPSG:4326'));
	var coordinates = geom.getLinearRing(0).getCoordinates();
	var area = Math.abs(wgs84Sphere.geodesicArea(coordinates));
	var coords = geom.getLinearRing(0).getCoordinates();
	var gsk3_coords = "";
	var wgs84_coords = "";
	var n = coords.length - 1;
	for (var i = 0; i < n; i++) {
		var wgs84_coord = coords[i];
		wgs84_coords += "(" + wgs84_coord[1] + "," + wgs84_coord[0] + ")<br/>";
		var gsk3_coord = ol.proj.transform(coords[i], ol.proj.get('EPSG:4326'), ol.proj.get('EPSG:31467'))
		gsk3_coords += "(" + gsk3_coord[0] + "," + gsk3_coord[1] + ")<br/>";
	}
	$('#dataPanel').append("WGS84 Coordinates<br/>");
	$('#dataPanel').append(wgs84_coords + "<br/>\n");
	$('#dataPanel').append("GSK3 Coordinates<br/>");
	$('#dataPanel').append(gsk3_coords + "<br/>\n");
	$('#dataPanel').append("Area" + "<br/>\n");
	$('#dataPanel').append((Math.round(area / 1000) / 10).toString() + " ha<br/><br/>\n");
	findIntersections();
	var end = new Date().getTime();
	var time = end - start;
	console.log('Execution time: ' + time);
}

draw.on('drawend', function(e) {
	displayInfo();
	send_btn.disabled = false;
});

$('#reset').click(function() {
duminil's avatar
duminil committed
228
229
230
231
	try {
		draw.finishDrawing();
	} finally {
		$('#dataPanel').empty();
232
		$("html").removeClass("wait");
duminil's avatar
duminil committed
233
234
235
236
		featureOverlay.getFeatures().clear();
		intersections.clear();
		reset_btn.disabled = true;
		send_btn.disabled = true;
237
		focusOnMap();
duminil's avatar
duminil committed
238
	}
239
240
241
});

$('#send').click(function() {
duminil's avatar
duminil committed
242
243
244
245
246
	focusOnMap();
	fxapp.downloadRegion(sketchAsWKT());
});

function sketchAsWKT() {
247
	var wktFormat = new ol.format.WKT();
duminil's avatar
duminil committed
248
	return wktFormat.writeFeature(sketch, {
249
		dataProjection : ol.proj.get('EPSG:31467'),
250
		featureProjection : ol.proj.get('EPSG:3857')
duminil's avatar
duminil committed
251
252
	})
}
253
254
255

function focusOnMap() {
	$('#map').focus();
duminil's avatar
duminil committed
256
	// $('#map').scrollIntoView();
257
258
}

duminil's avatar
duminil committed
259
focusOnMap();