simstadt_openlayers.js 8.67 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
//TODO: Add zoom to extent as control
//TODO: Don't allow multiple sketch
//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)'
		})
	})
});

49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
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

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

vectorSource.addMarker = function(x, y, name) {
	this.addFeature(new ol.Feature({
		geometry : new ol.geom.Point(ol.proj.transform([ x, y ], 'EPSG:31467', 'EPSG:3857')),
		name : name,
	}));
}

var createTextStyle = function(feature, resolution) {
	var align = "Normal";
	var baseline = "Middle";
	var size = "12px";
	var weight = "Normal";
	var font = weight + ' ' + size + ' ' + "Arial";
	var fillColor = '#aa3300';
	var outlineColor = "#ffffff";
	var outlineWidth = 3;

	return new ol.style.Text({
		textAlign : align,
		textBaseline : baseline,
		font : font,
		text : feature.get('name'),
		fill : new ol.style.Fill({
			color : fillColor
		}),
		stroke : new ol.style.Stroke({
			color : outlineColor,
			width : outlineWidth
		}),
	});
};

var createPointStyleFunction = function() {
	return function(feature, resolution) {
		var style = new ol.style.Style({
			image : new ol.style.Circle({
				radius : 10,
				fill : new ol.style.Fill({
					color : 'rgba(255, 0, 0, 0.1)'
				}),
				stroke : new ol.style.Stroke({
					color : 'red',
					width : 1
				})
			}),
			text : createTextStyle(feature, resolution)
		});
		return [ style ];
	};
};

var labels = new ol.layer.Vector({
	source : vectorSource,
	style : createPointStyleFunction()
});

110
111
var map = new ol.Map({
	target : 'map',
112
	layers : [ osm_layer, kml_layer, intersections_layer, labels ],
113
114
115
	interactions : ol.interaction.defaults({
		keyboard : true
	})
116
117
});

duminil's avatar
duminil committed
118
119
120
121
122
123
124
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
125
126
127
		var project = feature.get("project");
		var name = feature.get("name");
		feature["description"] = project + ">" + name;
duminil's avatar
duminil committed
128
129
		var citygmlHere;
		if (fromJavaFX) {
130
			citygmlHere = fxapp.checkIfCityGMLSAreAvailable(project, name);
duminil's avatar
duminil committed
131
132
133
134
135
		}
		feature["available"] = citygmlHere;
	});
});

136
137
138
139
140
141
142
143
144
145
146
147
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
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
// 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
196
197
	var i = 0;
	kml_source.forEachFeature(function(feature) {
198
		try {
duminil's avatar
duminil committed
199
			var jsonIntersection = turf.intersect(poly1, feature["geoJSON"]);
200
201
202
203
204
205
206
			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
207
				var citygml_percentage = Math.round(intersectionArea / feature["area"] * 100);
208
209
210
				var sketch_percentage = Math.round(intersectionArea / sketch_area * 100);
				intersections.addFeature(intersection);

duminil's avatar
duminil committed
211
212
				var description;
				if (feature["available"]) {
duminil's avatar
duminil committed
213
214
					description = "<a href=\"#\" onclick=\"downloadRegionFromCityGML(" + i + ");return false;\">"
							+ feature["description"] + "</a>";
duminil's avatar
duminil committed
215
				} else {
duminil's avatar
duminil committed
216
					description = feature['description'];
duminil's avatar
duminil committed
217
218
219
				}

				$('#dataPanel').append(description + " (" + citygml_percentage + "%");
220
221
222
223
224
225
226
227
228
				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
229
230
		i++;
	})
231
232
233
234
	if (!intersection_found) {
		$('#dataPanel').append("No intersection found with any CityGML<br/>\n");
	}
}
duminil's avatar
duminil committed
235
236

function downloadRegionFromCityGML(i) {
237
	// TODO: Disable all links
238
	$("html").addClass("wait");
duminil's avatar
duminil committed
239
	var feature = kml_source.getFeatures()[i];
240
241
242
243
	// 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
244
				.get("name"), vectorSource);
245
246
247
		var end = new Date().getTime();
		var time = end - start;
		console.log('DL Execution time: ' + time);
248
//		$('#dataPanel').append("Imported buildings : " + buildings_count);
249
250
		$("html").removeClass("wait");
	}, 100);
duminil's avatar
duminil committed
251
252
}

253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
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/>";
	}
270
271
272
273
//	$('#dataPanel').append("WGS84 Coordinates<br/>");
//	$('#dataPanel').append(wgs84_coords + "<br/>\n");
//	$('#dataPanel').append("GSK3 Coordinates<br/>");
//	$('#dataPanel').append(gsk3_coords + "<br/>\n");
274
275
276
277
278
279
280
281
282
283
284
285
286
287
	$('#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
288
289
290
291
	try {
		draw.finishDrawing();
	} finally {
		$('#dataPanel').empty();
292
		$("html").removeClass("wait");
duminil's avatar
duminil committed
293
294
295
296
		featureOverlay.getFeatures().clear();
		intersections.clear();
		reset_btn.disabled = true;
		send_btn.disabled = true;
297
		focusOnMap();
duminil's avatar
duminil committed
298
	}
299
300
301
});

$('#send').click(function() {
duminil's avatar
duminil committed
302
303
304
305
306
	focusOnMap();
	fxapp.downloadRegion(sketchAsWKT());
});

function sketchAsWKT() {
307
	var wktFormat = new ol.format.WKT();
duminil's avatar
duminil committed
308
	return wktFormat.writeFeature(sketch, {
309
		dataProjection : ol.proj.get('EPSG:31467'),
310
		featureProjection : ol.proj.get('EPSG:3857')
duminil's avatar
duminil committed
311
312
	})
}
313
314
315

function focusOnMap() {
	$('#map').focus();
duminil's avatar
duminil committed
316
	// $('#map').scrollIntoView();
317
318
}

duminil's avatar
duminil committed
319
focusOnMap();