simstadt_openlayers.js 12 KB
Newer Older
Eric Duminil's avatar
Eric Duminil committed
1
var regionChooser = (function(){
2
	var publicScope = {};
3
	var fromJavaFX = navigator.userAgent.indexOf('JavaFX') !== -1;
4
	//NOTE: Could do without jQuery
5
6
	var dataPanel = $('#dataPanel');
	var wgs84Sphere = new ol.Sphere(6378137);
7
	var gmlId = 0;
8
9
10
11
	
	if (fromJavaFX){
		$("html").addClass("wait");
	}
12

13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
    // Hash function. For testing purposes in browser (no Feature ID is available)
	const cyrb53 = (str, seed = 0) => {
	  let h1 = 0xdeadbeef ^ seed,
		h2 = 0x41c6ce57 ^ seed;
	  for (let i = 0, ch; i < str.length; i++) {
		ch = str.charCodeAt(i);
		h1 = Math.imul(h1 ^ ch, 2654435761);
		h2 = Math.imul(h2 ^ ch, 1597334677);
	  }
	  
	  h1 = Math.imul(h1 ^ (h1 >>> 16), 2246822507) ^ Math.imul(h2 ^ (h2 >>> 13), 3266489909);
	  h2 = Math.imul(h2 ^ (h2 >>> 16), 2246822507) ^ Math.imul(h1 ^ (h1 >>> 13), 3266489909);
	  
	  return 4294967296 * (2097151 & h2) + (h1 >>> 0);
	};
28

29
30
	var osm_layer = new ol.layer.Tile({
		source: new ol.source.OSM()
eric.duminil's avatar
eric.duminil committed
31
32
	});

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

44
45
46
47
48
49
50
51
52
53
54
55
	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
56

57
58
	var kml_layer = new ol.layer.Vector({
		source : kml_source,
59
		style : polygon_style('#447744', 0.2)
60
61
62
63
64
65
66
67
68
69
	});

	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
70
71
		})
	});
72

73
	publicScope.addCityGmlHull = function(kmlString) {
74
		options = {featureProjection: ol.proj.get('EPSG:3857')};
75
		feature = kmlFormat.readFeature(kmlString, options);
76
		feature.setId(gmlId++);
77
78
		kml_source.addFeature(feature);
		dataPanel.append('.');
79
80
81
82
		srsName = feature.get("srsName") || "EPSG:31467";
		if (proj4.defs(srsName) === undefined){
			console.warning(srsName + " isn't defined by Proj4js!")
		}
83
84
	};
	
85
86
	var map = new ol.Map({
		target : 'map',
87
		layers : [ osm_layer, kml_layer, intersections_layer ],
88
89
90
91
92
		interactions : ol.interaction.defaults({
			keyboard : true
		})
	});

93
	var geoJsonFormat = new ol.format.GeoJSON();
94
	var kmlFormat = new ol.format.KML({extractStyles: false});
95
96
97
98
99
100
101

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

	function updateGMLPolygons() {
		kml_source.forEachFeature(function(feature) {
102
			feature["geoJSON"] = geoJsonFormat.writeFeatureObject(feature);
103
			feature["area"] = feature.getGeometry().getArea();
Eric Duminil's avatar
Eric Duminil committed
104
105
			feature["project"] = feature.get("project");
			feature["name"] = feature.get("name");
106
107
108
109
110
111
112
113
114
115
			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
116
			fill : new ol.style.Fill({
117
118
119
120
121
122
123
124
125
126
127
				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
128
129
130
			})
		})
	});
131
132
	featureOverlay.setMap(map);

Eric Duminil's avatar
Eric Duminil committed
133
	//TODO: Rename to Javascript naming convention (CamelCase).
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
	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
160

161
162
163
164
165
	draw.on('drawstart', function(evt) {
		sketch = evt.feature;
		updateGMLPolygons();
	});
	var sourceProj = map.getView().getProjection();
166
167
168
169
170
171
172
173
	
	function showLinkToDownload(feature, jsonIntersection, polygonArea){
		var intersection = geoJsonFormat.readFeature(jsonIntersection);
		var intersectionArea = intersection.getGeometry().getArea();
		var citygml_percentage = Math.round(intersectionArea / feature["area"] * 100);
		var sketch_percentage = Math.round(intersectionArea / polygonArea * 100);
		intersections.addFeature(intersection);
		var link = '<li>'
Eric Duminil's avatar
todo    
Eric Duminil committed
174
175
		// TODO: Add checkbox + label + 
		// TODO: Add submit.
Eric Duminil's avatar
Notes    
Eric Duminil committed
176
        // TODO: If possible, highlight the corresponding polygon when hovering above a name.
177
178
        console.log("Feature name : "+ feature["name"]);
        console.log("Feature ID : "+ feature.getId());
179
		if (fromJavaFX) {
180
			link += '<input type="checkbox" id="citygml_' + feature.getId() + '" class="select_citygml"><label for="citygml_' + feature.getId() + '">' +  feature['name'] + '</label>';
181
		} else {
182
            h = cyrb53(feature['name']);
183
			link += '<input type="checkbox" id="citygml_' + h+ '" class="select_citygml"><label for="citygml_' + h + '">' +  feature['name'] + '</label>';
184
185
186
187
188
189
190
191
192
		}

		link += " (" + citygml_percentage + "%";
		if (sketch_percentage == 100) {
			link += ", all inside";
		}
		dataPanel.append(link + ")\n");
	}
	
193
	function findIntersection(feature, polygon) {
194
		try {
195
			return turf.intersect(polygon, feature["geoJSON"]);
196
197
198
199
		} catch (err) {
			console.log(feature.get('name') + " - " + err);
		}
	}
200
201

	function findIntersections() {
202
203
		var polygon = geoJsonFormat.writeFeatureObject(sketch);
		var polygonArea = sketch.getGeometry().getArea();
204
205
		var intersection_found = false;
		intersections.clear();
206
		//NOTE: getFeatures seems to not be sorted anymore. :-/
207
		features_by_project = groupBy(kml_source.getFeatures(), "project");
208
		
209
		Object.keys(features_by_project).forEach(function(project) {
210
		    features = features_by_project[project];
211
212
213
		    features_and_intersections = features.map(f=> [f, findIntersection(f,polygon)]).filter(l => l[1] !== undefined);
		    if (features_and_intersections.length > 0){
		    	intersection_found = true;
214
				dataPanel.append("<h2 class='info'>" + project);
215
216
		    	features_and_intersections.forEach(l => showLinkToDownload(l[0], l[1], polygonArea));
		    }
217
		});
218
		
219
220
221
222
223
		if (intersection_found) {
			document.getElementById("download").style.visibility = 'visible';
		}
		else {
			document.getElementById("download").style.visibility = 'hidden';
224
			dataPanel.append("No intersection found with any CityGML file.<br/>\n");
eric.duminil's avatar
eric.duminil committed
225
226
		}
	}
227

228
229
230
231
	publicScope.display = function(text){
		dataPanel.append(text + "<br/>\n");
	}
	
232
	publicScope.downloadRegionFromCityGML = function(checkbox_ids) {
233
234
		// TODO: Disable all links
		// TODO: DRY
235
236
237
238
239
240
241
242
		var features = checkbox_ids.map(checkbox_id => {
			var i = Number(checkbox_id.replace("citygml_", ""));
			return kml_source.getFeatureById(i);
		})
		
		var project = features[0].get("project"); //TODO: Check all the same
		var srsName = features[0].get("srsName") || "EPSG:31467"; //TODO: Check all the same
		var citygmlNames = features.map(f => f.get("name"));
243
244
245
246
247
248
		// Waiting 100ms in order to let the cursor change
		setTimeout(function() {
			var start = new Date().getTime();
			if (proj4.defs(srsName)){
				$("html").addClass("wait");
				console.log("Selected region is written in " + srsName + " coordinate system.");
249
				try {
250
251
252
					console.log("Before JAVA");
					fxapp.downloadRegionFromCityGML(sketchAsWKT(srsName), project, citygmlNames.join(";"), srsName);
					console.log("After JAVA");
253
254
					dataPanel.append("<h2 class='ok'>Done!</h2><br/>\n");
				} catch (e) {
255
					console.log("ERROR " + e);
256
257
					dataPanel.append("<h2 class='error'>Some problem occured!</h2><br/>\n");
				}
258
259
				var end = new Date().getTime();
				var time = end - start;
Eric Duminil's avatar
Eric Duminil committed
260
				console.log('Download Execution time: ' + (time / 1000).toFixed(3) + 's');
261
262
263
264
265
266
267
268
269
				setTimeout(function() {
					$("html").removeClass("wait");
				}, 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
270
271
	}

272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
	function displayInfo() {
		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 += "(" + 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
289
		}
290
		dataPanel.append("<h3 class='clean'>Area : " + (area / 10000).toFixed(1) + " ha\n");
291
292
		dataPanel.append('<button type="button" onclick="regionChooser.clickety_click()" id="download" style="visibility:hidden">Download Region</button>');
		dataPanel.append('<br/>\n');
293
		findIntersections();
eric.duminil's avatar
eric.duminil committed
294
	}
295
296
297
298
299

	draw.on('drawend', function() {
		displayInfo();
		draw.setActive(false);
	});
300
301
302
303
304
305
306
307
	
	// 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();
		}
	});
308

309
310
	function resetDrawing(){
		console.log("Reset drawing");
311
312
313
		try {
			draw.finishDrawing();
		} finally {
Eric Duminil's avatar
Eric Duminil committed
314
			displayHelp();
315
316
317
318
319
320
			$("html").removeClass("wait");
			draw.setActive(true);
			featureOverlay.getFeatures().clear();
			intersections.clear();
			focusOnMap();
		}
321
	}
322
	
323
324
325
326
327
328
329
330
	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')
		});
	}
331
	
332
333
334
335
	function focusOnMap() {
		$('#map').focus();
	}
	
336
337
338
339
	var fxapp = undefined;
	
	publicScope.setFxApp = function(app){
		fxapp = app;
340
341
342
		console.log = function(message){
			fxapp.log(message);
		}
343
344
345
346
		
		console.warning = function(message){
			fxapp.warning(message);
		}
347
	}
Eric Duminil's avatar
Eric Duminil committed
348
	
349
350
351
352
353
354
355
	groupBy = function(xs, key) {
		  return xs.reduce(function(rv, x) {
			(rv[x[key]] = rv[x[key]] || []).push(x);
			return rv;
		  }, {});
		};
	
Eric Duminil's avatar
Eric Duminil committed
356
357
	function displayHelp(){
		dataPanel.empty();
358
		dataPanel.append("<h2 class='info'>Welcome to Region Chooser!<br><br>\n");
Eric Duminil's avatar
Eric Duminil committed
359
360
361
362
363
364
		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");
	}
365

366
367
	// Executed by JavaFX when whole page is loaded.
	publicScope.ready = function() {
368
		updateGMLPolygons();
Eric Duminil's avatar
Eric Duminil committed
369
		displayHelp();
370
		$("html").removeClass("wait");
371
		console.log("READY!");
372
	}
373
374
375
	
	publicScope.clickety_click = function() {
		console.log("You clicked pretty well");
376
		var checkedBoxes = Array.from(document.querySelectorAll("input.select_citygml")).filter(c => c.checked);
Eric Duminil's avatar
Eric Duminil committed
377
		if (checkedBoxes.length === 0){
Eric Duminil's avatar
Eric Duminil committed
378
379
			console.log("You should select at least one citygml, though.");
		} else{
380
381
			checkedBoxes.forEach(c => {
				console.log("Nice! You checked Citygml " + c.id)
Eric Duminil's avatar
Eric Duminil committed
382
			});
383
				publicScope.downloadRegionFromCityGML(checkedBoxes.map(c => c.id));
Eric Duminil's avatar
Eric Duminil committed
384
		}
385
	}
386

387
	focusOnMap();
Eric Duminil's avatar
Eric Duminil committed
388
389
	//var regionChooser = publicScope; //NOTE: In order to open closure. For debugging
	return publicScope;
390
})();