forked from MapComplete/MapComplete
		
	
		
			
				
	
	
		
			364 lines
		
	
	
	
		
			15 KiB
		
	
	
	
		
			TypeScript
		
	
	
	
	
	
			
		
		
	
	
			364 lines
		
	
	
	
		
			15 KiB
		
	
	
	
		
			TypeScript
		
	
	
	
	
	
| import * as Assert from "assert";
 | |
| import {equal} from "assert";
 | |
| import T from "./TestHelper";
 | |
| import {GeoOperations} from "../Logic/GeoOperations";
 | |
| import {BBox} from "../Logic/BBox";
 | |
| 
 | |
| export default class GeoOperationsSpec extends T {
 | |
| 
 | |
|     private static polygon = {
 | |
|         "type": "Feature",
 | |
|         "properties": {},
 | |
|         "geometry": {
 | |
|             "type": "Polygon",
 | |
|             "coordinates": [
 | |
|                 [
 | |
|                     [
 | |
|                         1.8017578124999998,
 | |
|                         50.401515322782366
 | |
|                     ],
 | |
|                     [
 | |
|                         -3.1640625,
 | |
|                         46.255846818480315
 | |
|                     ],
 | |
|                     [
 | |
|                         5.185546875,
 | |
|                         44.74673324024678
 | |
|                     ],
 | |
|                     [
 | |
|                         1.8017578124999998,
 | |
|                         50.401515322782366
 | |
|                     ]
 | |
|                 ]
 | |
|             ]
 | |
|         }
 | |
|     };
 | |
|     private static multiPolygon = {
 | |
|         "type": "Feature",
 | |
|         "properties": {},
 | |
|         "geometry": {
 | |
|             "type": "MultiPolygon",
 | |
|             "coordinates": [[
 | |
|                 [
 | |
|                     [
 | |
|                         1.8017578124999998,
 | |
|                         50.401515322782366
 | |
|                     ],
 | |
|                     [
 | |
|                         -3.1640625,
 | |
|                         46.255846818480315
 | |
|                     ],
 | |
|                     [
 | |
|                         5.185546875,
 | |
|                         44.74673324024678
 | |
|                     ],
 | |
|                     [
 | |
|                         1.8017578124999998,
 | |
|                         50.401515322782366
 | |
|                     ]
 | |
|                 ],
 | |
|                 [
 | |
|                     [
 | |
|                         1.0107421875,
 | |
|                         48.821332549646634
 | |
|                     ],
 | |
|                     [
 | |
|                         1.329345703125,
 | |
|                         48.25394114463431
 | |
|                     ],
 | |
|                     [
 | |
|                         1.988525390625,
 | |
|                         48.71271258145237
 | |
|                     ],
 | |
|                     [
 | |
|                         0.999755859375,
 | |
|                         48.86471476180277
 | |
|                     ],
 | |
|                     [
 | |
|                         1.0107421875,
 | |
|                         48.821332549646634
 | |
|                     ]
 | |
|                 ]
 | |
|             ]]
 | |
|         }
 | |
|     };
 | |
| 
 | |
|     private static inHole = [1.42822265625, 48.61838518688487]
 | |
|     private static inMultiPolygon = [2.515869140625, 47.37603463349758]
 | |
|     private static outsidePolygon = [4.02099609375, 47.81315451752768]
 | |
| 
 | |
|     constructor() {
 | |
|         super(
 | |
|             "GeoOperationsSpec", [
 | |
|                 ["Point out of polygon", () => {
 | |
|                     GeoOperationsSpec.isFalse(GeoOperations.inside([
 | |
|                         3.779296875,
 | |
|                         48.777912755501845
 | |
|                     ], GeoOperationsSpec.polygon), "Point is outside of the polygon");
 | |
|                 }
 | |
| 
 | |
|                 ],
 | |
|                 ["Point inside of polygon", () => {
 | |
| 
 | |
|                     GeoOperationsSpec.isTrue(GeoOperations.inside([
 | |
|                         1.23046875,
 | |
|                         47.60616304386874
 | |
|                     ], GeoOperationsSpec.polygon), "Point is inside of the polygon");
 | |
|                 }
 | |
|                 ],
 | |
|                 ["MultiPolygonTest", () => {
 | |
| 
 | |
|                     const isTrue = GeoOperationsSpec.isTrue;
 | |
|                     const isFalse = GeoOperationsSpec.isFalse;
 | |
| 
 | |
|                     isFalse(GeoOperations.inside(GeoOperationsSpec.inHole, GeoOperationsSpec.multiPolygon), "InHole was detected as true");
 | |
|                     isTrue(GeoOperations.inside(GeoOperationsSpec.inMultiPolygon, GeoOperationsSpec.multiPolygon), "InMultiPolgyon was not detected as true");
 | |
|                     isFalse(GeoOperations.inside(GeoOperationsSpec.outsidePolygon, GeoOperationsSpec.multiPolygon), "OutsideOfPolygon was detected as true");
 | |
| 
 | |
|                 }],
 | |
|                 ["Intersection between a line and a polygon", () => {
 | |
|                     const line = {
 | |
|                         "type": "Feature",
 | |
|                         "properties": {},
 | |
|                         "geometry": {
 | |
|                             "type": "LineString",
 | |
|                             "coordinates": [
 | |
|                                 [
 | |
|                                     3.779296875,
 | |
|                                     48.777912755501845
 | |
|                                 ],
 | |
|                                 [
 | |
|                                     1.23046875,
 | |
|                                     47.60616304386874
 | |
|                                 ]
 | |
|                             ]
 | |
|                         }
 | |
|                     };
 | |
| 
 | |
|                     const overlap = GeoOperations.calculateOverlap(line, [GeoOperationsSpec.polygon]);
 | |
|                     Assert.equal(1, overlap.length)
 | |
|                 }],
 | |
|                 ["Fully enclosed", () => {
 | |
|                     const line = {
 | |
|                         "type": "Feature",
 | |
|                         "properties": {},
 | |
|                         "geometry": {
 | |
|                             "type": "LineString",
 | |
|                             "coordinates": [
 | |
|                                 [
 | |
|                                     0.0439453125,
 | |
|                                     47.31648293428332
 | |
|                                 ],
 | |
|                                 [
 | |
|                                     0.6591796875,
 | |
|                                     46.77749276376827
 | |
|                                 ]
 | |
|                             ]
 | |
|                         }
 | |
|                     };
 | |
| 
 | |
|                     const overlap = GeoOperations.calculateOverlap(line, [GeoOperationsSpec.polygon]);
 | |
|                     Assert.equal(1, overlap.length)
 | |
|                 }],
 | |
|                 ["overlapWith matches points too", () => {
 | |
|                     const point = {
 | |
|                         "type": "Feature",
 | |
|                         "properties": {},
 | |
|                         "geometry": {
 | |
|                             "type": "Point",
 | |
|                             "coordinates": [
 | |
|                                 2.274169921875,
 | |
|                                 46.76244305208004
 | |
|                             ]
 | |
|                         }
 | |
|                     };
 | |
| 
 | |
|                     const overlap = GeoOperations.calculateOverlap(point, [GeoOperationsSpec.polygon]);
 | |
|                     Assert.equal(1, overlap.length)
 | |
|                 }],
 | |
|                 ["bbox bounds test",
 | |
|                     () => {
 | |
|                         const bbox = BBox.fromTile(16, 32754, 21785)
 | |
|                         equal(-0.076904296875, bbox.minLon)
 | |
|                         equal(-0.0714111328125, bbox.maxLon)
 | |
|                         equal(51.5292513551899, bbox.minLat)
 | |
|                         equal(51.53266860674158, bbox.maxLat)
 | |
|                     }
 | |
|                 ],
 | |
|                 ["Regression test: intersection/overlap", () => {
 | |
| 
 | |
|                     const polyGrb ={
 | |
|                         "type": "Feature",
 | |
|                         "properties": {
 | |
|                             "osm_id": "25189153",
 | |
|                             "size_grb_building": "217.14",
 | |
|                             "addr:housenumber": "173",
 | |
|                             "addr:street": "Kortrijksestraat",
 | |
|                             "building": "house",
 | |
|                             "source:geometry:entity": "Gbg",
 | |
|                             "source:geometry:date": "2015/02/27",
 | |
|                             "source:geometry:oidn": "1729460",
 | |
|                             "source:geometry:uidn": "8713648",
 | |
|                             "H_DTM_MIN": "17.28",
 | |
|                             "H_DTM_GEM": "17.59",
 | |
|                             "H_DSM_MAX": "29.04",
 | |
|                             "H_DSM_P99": "28.63",
 | |
|                             "HN_MAX": "11.45",
 | |
|                             "HN_P99": "11.04",
 | |
|                             "detection_method": "from existing OSM building source: house ,hits (3)",
 | |
|                             "auto_building": "house",
 | |
|                             "size_shared": "210.68",
 | |
|                             "size_source_building": "212.63",
 | |
|                             "id": "https://betadata.grbosm.site/grb?bbox=360935.6475626023,6592540.815539878,361088.52161917265,6592693.689596449/37",
 | |
|                             "_lat": "50.83736194999996",
 | |
|                             "_lon": "3.2432137000000116",
 | |
|                             "_layer": "GRB",
 | |
|                             "_length": "48.51529464293261",
 | |
|                             "_length:km": "0.0",
 | |
|                             "_now:date": "2021-12-05",
 | |
|                             "_now:datetime": "2021-12-05 21:51:40",
 | |
|                             "_loaded:date": "2021-12-05",
 | |
|                             "_loaded:datetime": "2021-12-05 21:51:40"
 | |
|                         },
 | |
|                         "geometry": {
 | |
|                             "type": "Polygon",
 | |
|                             "coordinates": [
 | |
|                                 [
 | |
|                                     [
 | |
|                                         3.2431059999999974,
 | |
|                                         50.83730270000021
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.243174299999987,
 | |
|                                         50.83728850000007
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.2432116000000173,
 | |
|                                         50.83736910000003
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.2433214000000254,
 | |
|                                         50.83740350000011
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.24329779999996,
 | |
|                                         50.837435399999855
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.2431881000000504,
 | |
|                                         50.83740090000025
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.243152699999997,
 | |
|                                         50.83738980000017
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.2431059999999974,
 | |
|                                         50.83730270000021
 | |
|                                     ]
 | |
|                                 ]
 | |
|                             ]
 | |
|                         },
 | |
|                         "id": "https://betadata.grbosm.site/grb?bbox=360935.6475626023,6592540.815539878,361088.52161917265,6592693.689596449/37",
 | |
|                         "_lon": 3.2432137000000116,
 | |
|                         "_lat": 50.83736194999996,
 | |
|                         "bbox": {
 | |
|                             "maxLat": 50.837435399999855,
 | |
|                             "maxLon": 3.2433214000000254,
 | |
|                             "minLat": 50.83728850000007,
 | |
|                             "minLon": 3.2431059999999974
 | |
|                         }
 | |
|                     }
 | |
|                     const polyHouse = {
 | |
|                         "type": "Feature",
 | |
|                         "id": "way/594963177",
 | |
|                         "properties": {
 | |
|                             "timestamp": "2021-12-05T04:04:55Z",
 | |
|                             "version": 3,
 | |
|                             "changeset": 114571409,
 | |
|                             "user": "Pieter Vander Vennet",
 | |
|                             "uid": 3818858,
 | |
|                             "addr:housenumber": "171",
 | |
|                             "addr:street": "Kortrijksestraat",
 | |
|                             "building": "house",
 | |
|                             "source:geometry:date": "2018-10-22",
 | |
|                             "source:geometry:ref": "Gbg/5096537",
 | |
|                             "_last_edit:contributor": "Pieter Vander Vennet",
 | |
|                             "_last_edit:contributor:uid": 3818858,
 | |
|                             "_last_edit:changeset": 114571409,
 | |
|                             "_last_edit:timestamp": "2021-12-05T04:04:55Z",
 | |
|                             "_version_number": 3,
 | |
|                             "id": "way/594963177",
 | |
|                             "_backend": "https://www.openstreetmap.org",
 | |
|                             "_lat": "50.83736395",
 | |
|                             "_lon": "3.2430937",
 | |
|                             "_layer": "OSM-buildings",
 | |
|                             "_length": "43.561938680928506",
 | |
|                             "_length:km": "0.0",
 | |
|                             "_now:date": "2021-12-05",
 | |
|                             "_now:datetime": "2021-12-05 21:51:40",
 | |
|                             "_loaded:date": "2021-12-05",
 | |
|                             "_loaded:datetime": "2021-12-05 21:51:39",
 | |
|                             "_surface": "93.32785810484549",
 | |
|                             "_surface:ha": "0"
 | |
|                         },
 | |
|                         "geometry": {
 | |
|                             "type": "Polygon",
 | |
|                             "coordinates": [
 | |
|                                 [
 | |
|                                     [
 | |
|                                         3.2429993,
 | |
|                                         50.8373243
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.243106,
 | |
|                                         50.8373027
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.2431527,
 | |
|                                         50.8373898
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.2431881,
 | |
|                                         50.8374009
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.2431691,
 | |
|                                         50.8374252
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.2430936,
 | |
|                                         50.837401
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.243046,
 | |
|                                         50.8374112
 | |
|                                     ],
 | |
|                                     [
 | |
|                                         3.2429993,
 | |
|                                         50.8373243
 | |
|                                     ]
 | |
|                                 ]
 | |
|                             ]
 | |
|                         },
 | |
|                         "_lon": 3.2430937,
 | |
|                         "_lat": 50.83736395,
 | |
|                         "bbox": {
 | |
|                             "maxLat": 50.8374252,
 | |
|                             "maxLon": 3.2431881,
 | |
|                             "minLat": 50.8373027,
 | |
|                             "minLon": 3.2429993
 | |
|                         }
 | |
|                     }
 | |
| 
 | |
|                     const overlaps = GeoOperations.calculateOverlap(polyGrb, [polyHouse])
 | |
|                     Assert.equal(overlaps.length, 1)
 | |
|                     const overlapsRev= GeoOperations.calculateOverlap(polyHouse, [polyGrb])
 | |
|                     Assert.equal(overlaps.length, 1)
 | |
| 
 | |
|                 }]
 | |
|             ]
 | |
|         )
 | |
| 
 | |
|     }
 | |
| }
 |