map alignment experiment

src/commonTest/kotlin/com/jillesvangurp/geogeometry/MapAlignmentTest.kt

@@ -0,0 +1,98 @@
+package com.jillesvangurp.geogeometry
+
+import com.jillesvangurp.geo.GeoGeometry
+import com.jillesvangurp.geo.GeoGeometry.Companion.roundDecimals
+import com.jillesvangurp.geojson.FeatureCollection
+import com.jillesvangurp.geojson.Geometry
+import com.jillesvangurp.geojson.asFeature
+import com.jillesvangurp.geojson.geoJsonIOUrl
+import com.jillesvangurp.geojson.latitude
+import com.jillesvangurp.geojson.longitude
+import kotlin.test.Test
+import kotlinx.serialization.json.JsonPrimitive
+import kotlinx.serialization.json.buildJsonObject
+
+
+enum class MapProviders {
+    Here, Google, Osm, Apple
+}
+
+fun DoubleArray.roundAndReverse() = this.map { it.roundDecimals(8) }.toDoubleArray().also {it.reverse()}
+
+class MapAlignmentTest {
+
+    val referencePoints = mapOf(
+        "Ahoy Berlin NW Building Corner in Berlin" to mapOf(
+            MapProviders.Osm to doubleArrayOf(52.54150, 13.38982).roundAndReverse(),
+            MapProviders.Google to doubleArrayOf(52.541510535244676, 13.38980753280997).roundAndReverse(),
+//            MapProviders.Here to doubleArrayOf(52.54131,13.38965).roundAndReverse(),
+            MapProviders.Here to doubleArrayOf(52.54150630324809, 13.389818199084072).roundAndReverse(),
+            MapProviders.Apple to doubleArrayOf(52.541497, 13.389822).roundAndReverse()
+        ),
+        "Siege Säule Centroid" to mapOf(
+            MapProviders.Osm to doubleArrayOf(52.51450, 13.35012).roundAndReverse(),
+            MapProviders.Google to doubleArrayOf(52.51451224763399, 13.35010044553336).roundAndReverse(),
+            MapProviders.Here to doubleArrayOf(52.514514513886496, 13.35009726850834).roundAndReverse(),
+            MapProviders.Apple to doubleArrayOf(52.514524, 13.350096).roundAndReverse(),
+        ),
+        "Galeria Kaufhof SW corner" to mapOf(
+            MapProviders.Osm to doubleArrayOf(52.52185,13.41171).roundAndReverse(),
+            MapProviders.Google to doubleArrayOf(52.52184966204656, 13.411722336815041).roundAndReverse(),
+            MapProviders.Here to doubleArrayOf(52.52184025422203, 13.411664047844162).roundAndReverse(),
+            MapProviders.Apple to doubleArrayOf(52.521853, 13.411718).roundAndReverse(),
+        )
+    )
+
+    @Test
+    fun distanceForGoogleHereAndOsmForReferencePoints() {
+        referencePoints.forEach { (testCase, points) ->
+            println(
+                """
+                ### $testCase
+                
+                | Provider 1 | Point | Provider 2 | Point | Distance (meters) |
+                |----|----|----|----|----|
+            """.trimIndent()
+            )
+
+
+            val ps = points.toList()
+            val seen = mutableSetOf<List<MapProviders>>()
+
+            ps.indices.forEach { index ->
+                val (refProvider, refPoint) = ps[index]
+                ps.indices.forEach { i ->
+                    if (i != index) {
+                        val (otherProvider, otherPoint) = ps[i]
+                        val providerPair = listOf(refProvider, otherProvider).sorted()
+                        if (!seen.contains(providerPair)) {
+                            println(
+                                "| $refProvider | ${refPoint.latitude}, ${refPoint.longitude} | $otherProvider | ${otherPoint.latitude}, ${otherPoint.longitude} | ${
+                                    GeoGeometry.distance(
+                                        refPoint,
+                                        otherPoint
+                                    ).roundDecimals(2)
+                                } |"
+                            )
+                            seen.add(providerPair)
+                        }
+                    }
+                }
+            }
+            val url = ps.map { (provider, point) ->
+                Geometry.Point(point).asFeature(
+                    buildJsonObject {
+                        put("title", JsonPrimitive(provider.name))
+                    }
+                )
+            }.let {
+                FeatureCollection(it).geoJsonIOUrl
+            }
+            println("""
+                
+                [Review on geojson.io]($url)
+                
+            """.trimIndent())
+        }
+    }
+}