simstadt_openlayers.js 13.3 KB
Newer Older
eric.duminil's avatar
eric.duminil committed
1
2
3
proj4.defs("EPSG:3068", "+proj=cass +lat_0=52.41864827777778 +lon_0=13.62720366666667 +x_0=40000 +y_0=10000 +ellps=bessel +datum=potsdam +units=m +no_defs"); // http://spatialreference.org/ref/epsg/3068/proj4js/
proj4.defs("EPSG:32632", "+proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"); // http://spatialreference.org/ref/epsg/32632/proj4js/
proj4.defs("EPSG:31463", "+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/31463/proj4js/
4
proj4.defs("EPSG:31466", "+proj=tmerc +lat_0=0 +lon_0=6 +k=1 +x_0=2500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs "); // http://spatialreference.org/ref/epsg/31466/proj4js/
eric.duminil's avatar
eric.duminil committed
5
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/
6
proj4.defs("EPSG:31468",  "+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 +y_0=0 +ellps=bessel +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs"); 
eric.duminil's avatar
eric.duminil committed
7
8
9
10
proj4.defs("EPSG:32118", "+proj=lcc +lat_1=41.03333333333333 +lat_2=40.66666666666666 +lat_0=40.16666666666666 +lon_0=-74 +x_0=300000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"); // http://spatialreference.org/ref/epsg/32118/proj4js/
proj4.defs("EPSG:2263", "+proj=lcc +lat_1=41.03333333333333 +lat_2=40.66666666666666 +lat_0=40.16666666666666 +lon_0=-74 +x_0=300000.0000000001 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs "); //  http://www.spatialreference.org/ref/epsg/nad83-new-york-long-island-ftus/proj4/
//NOTE: Proj4 string for 28992 is wrong at http://spatialreference.org/ref/epsg/amersfoort-rd-new/
//NOTE: Corrected version from https://oegeo.wordpress.com/2008/05/20/note-to-self-the-one-and-only-rd-projection-string/
11
proj4.defs("EPSG:25832", "+proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs "); // http://spatialreference.org/ref/epsg/etrs89-utm-zone-32n/proj4/
eric.duminil's avatar
eric.duminil committed
12
13
proj4.defs("EPSG:28992", "+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.999908 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812 +no_defs <>"); // 

14
var regionChooser = (function(){
15
	var publicScope = {};
16
	var fromJavaFX = navigator.userAgent.indexOf('JavaFX') !== -1;
17
18
	var dataPanel = $('#dataPanel');
	var wgs84Sphere = new ol.Sphere(6378137);
19
20
21
22
	
	if (fromJavaFX){
		$("html").addClass("wait");
	}
23
24
25

	var osm_layer = new ol.layer.Tile({
		source: new ol.source.OSM()
eric.duminil's avatar
eric.duminil committed
26
27
	});

28
29
30
31
32
33
34
35
36
	function read_kml(url){
		return new ol.source.KML({
			projection : ol.proj.get('EPSG:3857'),
			url : url,
			extractAttributes : false,
			extractStyles : false
		});
	}
	
37
	var kml_source = read_kml(fromJavaFX ? undefined : 'data/citygml_hulls.kml');
eric.duminil's avatar
eric.duminil committed
38

39
40
41
42
43
44
45
46
47
48
49
50
	function polygon_style(color, alpha) {
		return new ol.style.Style({
			fill : new ol.style.Fill({
				color : 'rgba(255, 255, 255,' + alpha + ')'
			}),
			stroke : new ol.style.Stroke({
				color : color,
				width : 2,
				lineDash : [ 5, 10 ]
			}),
		});
	}
eric.duminil's avatar
eric.duminil committed
51

52
53
	var kml_layer = new ol.layer.Vector({
		source : kml_source,
54
		style : polygon_style('#447744', 0.2)
55
56
57
58
59
60
61
62
63
64
	});

	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.2)'
			})
eric.duminil's avatar
eric.duminil committed
65
66
		})
	});
67
68
69
70
71

	novafactory_vectors = new ol.source.Vector({
		features : []
	});

72
	publicScope.addCityGmlHull = function(kmlString) {
73
		options = {featureProjection: ol.proj.get('EPSG:3857')};
74
		feature = kmlFormat.readFeature(kmlString, options);
75
76
		kml_source.addFeature(feature);
		dataPanel.append('.');
77
78
79
80
		srsName = feature.get("srsName") || "EPSG:31467";
		if (proj4.defs(srsName) === undefined){
			console.warning(srsName + " isn't defined by Proj4js!")
		}
81
82
	};
	
83
84
85
86
87
88
89
90
	publicScope.addNovaFactoryProduct = function(xmin, ymin, xmax, ymax, name, epsgId) {
		var box = new ol.geom.Polygon(
				[ [ [ xmin, ymin ], [ xmin, ymax ], [ xmax, ymax ], [ xmax, ymin ], [ xmin, ymin ] ] ]);
		box.transform('EPSG:' + epsgId, 'EPSG:3857');
		var feature = new ol.Feature({
			geometry : box,
			name : name,
		});
91
		feature["geoJSON"] = geoJsonFormat.writeFeatureObject(feature);
eric.duminil's avatar
eric.duminil committed
92
		feature["area"] = feature.getGeometry().getArea();
93
94
		feature["description"] = "novaFACTORY>" + name;
		feature["source"] = "NovaFACTORY";
95
		
96
97
		novafactory_vectors.addFeature(feature);
	};
98
	
99
100
101
	var novafactory_layer = new ol.layer.Vector({
		source : novafactory_vectors,
		style : polygon_style('#ff7700', 0.1)
eric.duminil's avatar
eric.duminil committed
102
	});
103
104
105

	var map = new ol.Map({
		target : 'map',
106
		layers : [ osm_layer, kml_layer, novafactory_layer, intersections_layer ],
107
108
109
110
111
		interactions : ol.interaction.defaults({
			keyboard : true
		})
	});

112
	var geoJsonFormat = new ol.format.GeoJSON();
113
	var kmlFormat = new ol.format.KML({extractStyles: false});
114
115
116
117
118
119
120

	kml_layer.addEventListener("change", function() {
		map.getView().fitExtent(kml_source.getExtent(), (map.getSize()));
	});

	function updateGMLPolygons() {
		kml_source.forEachFeature(function(feature) {
121
			feature["geoJSON"] = geoJsonFormat.writeFeatureObject(feature);
122
123
124
125
126
127
128
129
130
131
132
133
134
135
			feature["area"] = feature.getGeometry().getArea();
			var project = feature.get("project");
			var name = feature.get("name");
			feature["description"] = project + ">" + name;
			feature["source"] = "CityGML";
		});
	}

	// 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({
eric.duminil's avatar
eric.duminil committed
136
			fill : new ol.style.Fill({
137
138
139
140
141
142
143
144
145
146
147
				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'
				})
eric.duminil's avatar
eric.duminil committed
148
149
150
			})
		})
	});
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
	featureOverlay.setMap(map);

	var selected_features = featureOverlay.getFeatures();
	selected_features.on('add', function(event) {
		var feature = event.element;
		feature.on("change", function() {
			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);

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

	var sketch;
eric.duminil's avatar
eric.duminil committed
179

180
181
182
183
184
185
186
187
	draw.on('drawstart', function(evt) {
		sketch = evt.feature;
		updateGMLPolygons();
	});
	var sourceProj = map.getView().getProjection();

	function findIntersections() {
		var sketch_area = sketch.getGeometry().getArea();
188
		var poly1 = geoJsonFormat.writeFeatureObject(sketch);
189
190
191
192
193
194
195
196
197
198
		var intersection_found = false;
		intersections.clear();
		function findIntersection(feature) {
			try {
				var jsonIntersection = turf.intersect(poly1, feature["geoJSON"]);
				if (undefined !== jsonIntersection) {
					if (!intersection_found) {
						dataPanel.append("Intersection found with :<br/>\n");
						intersection_found = true;
					}
199
					var intersection = geoJsonFormat.readFeature(jsonIntersection);
200
201
202
203
					var intersectionArea = intersection.getGeometry().getArea();
					var citygml_percentage = Math.round(intersectionArea / feature["area"] * 100);
					var sketch_percentage = Math.round(intersectionArea / sketch_area * 100);
					intersections.addFeature(intersection);
204
					var link = '<li>'
205
206

					var description;
207
					if (fromJavaFX) {
208
						link += "<a href=\"#\" onclick=\"regionChooser.downloadRegionFrom" + feature["source"] + "(" + i
209
210
								+ ");return false;\">" + feature["description"] + "</a>";
					} else {
211
						link += feature['description'];
212
213
					}

214
					link += " (" + citygml_percentage + "%";
215
					if (sketch_percentage == 100) {
216
						link += ", all inside";
217
					}
218
					dataPanel.append(link + ")\n");
eric.duminil's avatar
eric.duminil committed
219
220
				}

221
222
223
224
225
226
227
228
229
230
231
			} catch (err) {
				console.log(feature.get('description') + " - " + err);
			}
			i++;
		}
		var i = 0;
		novafactory_vectors.forEachFeature(findIntersection);
		i = 0;
		kml_source.forEachFeature(findIntersection);
		if (!intersection_found) {
			dataPanel.append("No intersection found with any CityGML or NovaFactory product<br/>\n");
eric.duminil's avatar
eric.duminil committed
232
233
		}
	}
234

235
236
237
238
	publicScope.display = function(text){
		dataPanel.append(text + "<br/>\n");
	}
	
239
240
241
242
243
244
245
246
247
248
249
	publicScope.downloadRegionFromCityGML = function(i) {
		// TODO: Disable all links
		// TODO: DRY
		var feature = kml_source.getFeatures()[i];
		// Waiting 100ms in order to let the cursor change
		setTimeout(function() {
			var start = new Date().getTime();
			var srsName = feature.get("srsName") || "EPSG:31467";
			if (proj4.defs(srsName)){
				$("html").addClass("wait");
				console.log("Selected region is written in " + srsName + " coordinate system.");
250
				fxapp.downloadRegionFromCityGML(sketchAsWKT(srsName), feature.get("project"), feature.get("name"), srsName);
251
252
253
254
255
256
257
258
259
260
261
262
263
				var end = new Date().getTime();
				var time = end - start;
				console.log('DL Execution time: ' + time);
				setTimeout(function() {
					$("html").removeClass("wait");
					dataPanel.append("Done<br/>\n");
				}, 100);
			} else {
				var msg = "ERROR : Unknown coordinate system : \"" + srsName + "\". Cannot extract any region";
				console.log(msg);
				dataPanel.append(msg + "<br/>\n");
			}
		}, 100);
eric.duminil's avatar
eric.duminil committed
264
265
	}

266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
	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();
		if (!fromJavaFX) {
			var wgs84_coords = "";
			var n = coords.length;
			for (var i = 0; i < n; i++) {
				var wgs84_coord = coords[i];
				// wgs84_coords += "regionPolygon.add(new Coord(" + wgs84_coord[1] +
				// "," + wgs84_coord[0] + "));<br/>";
				wgs84_coords += "(" + wgs84_coord[1] + "," + wgs84_coord[0] + ")<br/>";

			}
			dataPanel.append("WGS84 Coordinates<br/>");
			dataPanel.append(wgs84_coords + "<br/>\n");
eric.duminil's avatar
eric.duminil committed
286
		}
287
288
289
290
291
292
		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);
eric.duminil's avatar
eric.duminil committed
293
	}
294
295
296
297
298

	draw.on('drawend', function() {
		displayInfo();
		draw.setActive(false);
	});
299
300
301
302
303
304
305
306
	
	// Pressing ESCAPE or DELETE resets the drawing.
	// With OpenLayers 3.9, draw_interaction.removeLastPoint(); might be better.
	document.addEventListener('keydown', function(e) {
		if (e.which == 27 || e.which == 46){
			resetDrawing();
		}
	});
307

308
309
	function resetDrawing(){
		console.log("Reset drawing");
310
311
312
		try {
			draw.finishDrawing();
		} finally {
Eric Duminil's avatar
Eric Duminil committed
313
			displayHelp();
314
315
316
317
318
319
			$("html").removeClass("wait");
			draw.setActive(true);
			featureOverlay.getFeatures().clear();
			intersections.clear();
			focusOnMap();
		}
320
	}
321
	
322
323
324
	novafactory_layer.downloadFinished = function() {
		// FIXME: Weird <br>s are inserted between lines
		// FIXME: Doesn't stop waiting cursor
eric.duminil's avatar
eric.duminil committed
325
		$("html").removeClass("wait");
326
327
328
329
330
331
332
333
334
335
		setTimeout(function() {
			dataPanel.append("NovaFactory : DONE <br/>\n");
		}, 100);
	};

	novafactory_layer.updateStatus = function(status) {
		dataPanel.append("NovaFactory : " + status + "<br/>\n");
	};

	novafactory_layer.selectSaveFile = function(zipFilename) {
336
		fxapp.extractZIPtoGML(zipFilename);
337
338
339
340
341
342
343
	};

	publicScope.downloadRegionFromNovaFACTORY = function(i) {
		$("html").addClass("wait");
		var feature = novafactory_vectors.getFeatures()[i];
		// Waiting 100ms in order to let the cursor change
		setTimeout(function() {
344
			fxapp.downloadRegion(sketchAsWKT(), feature.get('name'), novafactory_layer);
345
		}, 100);
eric.duminil's avatar
eric.duminil committed
346
	}
347
348
349
350
351
352
353
354
355

	function sketchAsWKT(srsName) {
		srsName = (typeof srsName === 'undefined') ? 'EPSG:4326' : srsName;
		var wktFormat = new ol.format.WKT();
		return wktFormat.writeFeature(sketch, {
			dataProjection : ol.proj.get(srsName),
			featureProjection : ol.proj.get('EPSG:3857')
		});
	}
356
	
357
358
359
360
361
	function focusOnMap() {
		$('#map').focus();
		// $('#map').scrollIntoView();
	}
	
362
363
364
365
	var fxapp = undefined;
	
	publicScope.setFxApp = function(app){
		fxapp = app;
366
367
368
		console.log = function(message){
			fxapp.log(message);
		}
369
370
371
372
		
		console.warning = function(message){
			fxapp.warning(message);
		}
373
	}
Eric Duminil's avatar
Eric Duminil committed
374
375
376
377
378
379
380
381
382
383
	
	function displayHelp(){
		dataPanel.empty();
		dataPanel.append("Welcome to Region Chooser!<br><br>\n");
		dataPanel.append("You can draw a polygon on the map by clicking.<br>\n");
		dataPanel.append("You can add a new point to an existing edge by clicking and dragging.<br>\n");
		dataPanel.append("You can remove a point with SHIFT + clicking.<br>\n");
		dataPanel.append("You can cancel drawing with ESC or DEL.<br><br>\n");
		dataPanel.append("After drawing a polygon which intersects with a GML file, you can download the corresponding part by clicking on the filename.<br>\n");
	}
384

385
386
	// Executed by JavaFX when whole page is loaded.
	publicScope.ready = function() {
387
		updateGMLPolygons();
Eric Duminil's avatar
Eric Duminil committed
388
		displayHelp();
389
		$("html").removeClass("wait");
390
		console.log("READY!");
391
392
	}

393
	focusOnMap();
394
395
396
	//var regionChooser = publicScope; //NOTE: In order to open closure. For debugging
	return publicScope;
})();