add gis

Author: zing-dev

modules/gis/distance/distance.go

@@ -0,0 +1,42 @@
+package main
+
+import (
+	"fmt"
+	"github.com/asmarques/geodist"
+	geo "github.com/kellydunn/golang-geo"
+	"gitlab.com/psns/geocalc"
+)
+
+type Dis struct {
+	Start    geocalc.Point
+	End      geocalc.Point
+	Distance float64
+}
+
+// Distance calculation
+func Distance() {
+	for _, dis := range []Dis{
+		{
+			//41.53786, 121.5665	40.99264, 123.1513	145.661公里 / 90.509英里
+			Start:    geocalc.Point{Lat: 41.53786, Lon: 121.5665},
+			End:      geocalc.Point{Lat: 40.99264, Lon: 123.1513},
+			Distance: 145661,
+		},
+		{
+			Start:    geocalc.Point{Lat: 30.551589, Lon: 114.202234},
+			End:      geocalc.Point{Lat: 30.557181, Lon: 114.503553},
+			Distance: 28841,
+		},
+	} {
+		d1, _ := geodist.VincentyDistance(geodist.Point{Lat: dis.Start.Lat, Long: dis.Start.Lon}, geodist.Point{Lat: dis.End.Lat, Long: dis.End.Lon})
+		d2 := geodist.HaversineDistance(geodist.Point{Lat: dis.Start.Lat, Long: dis.Start.Lon}, geodist.Point{Lat: dis.End.Lat, Long: dis.End.Lon})
+		d3 := D2(dis.Start.Lat, dis.Start.Lon, dis.End.Lat, dis.End.Lon)
+		d4 := geocalc.Distance(dis.Start, dis.End)
+		d5 := geo.NewPoint(dis.Start.Lat, dis.Start.Lon).GreatCircleDistance(geo.NewPoint(dis.End.Lat, dis.End.Lon))
+		fmt.Printf("d1:%.4f,d2:%.4f,d3:%.4f,d4:%.4f,d5:%.4f\n", d1, d2, d3/1000, d4/1000, d5)
+	}
+}
+
+func main() {
+	Distance()
+}

modules/gis/distance/f.go

@@ -0,0 +1,21 @@
+package main
+
+import "math"
+
+func D2(lat1, lon1, lat2, lon2 float64) float64 {
+	const r = 6378137 // Earth radius in METERS https://en.wikipedia.org/wiki/Earth_radius
+
+	// Distance between lat and longs
+	dLat := (lat2 - lat1) * math.Pi / 180.0
+	dLon := (lon2 - lon1) * math.Pi / 180.0
+
+	// convert to radians
+	lat1 = lat1 * math.Pi / 180.0
+	lat2 = lat2 * math.Pi / 180.0
+
+	// apply formulae
+	a := (math.Pow(math.Sin(dLat/2), 2)) + math.Pow(math.Sin(dLon/2), 2)*math.Cos(lat1)*math.Cos(lat2)
+	c := 2 * math.Asin(math.Sqrt(a))
+
+	return r * c
+}

modules/gis/go.mod

@@ -0,0 +1,16 @@
+module psns.geocalc
+
+go 1.17
+
+require (
+	github.com/asmarques/geodist v1.0.1
+	github.com/kellydunn/golang-geo v0.7.0
+	gitlab.com/psns/geocalc v1.3.1
+)
+
+require (
+	github.com/erikstmartin/go-testdb v0.0.0-20160219214506-8d10e4a1bae5 // indirect
+	github.com/kylelemons/go-gypsy v1.0.0 // indirect
+	github.com/lib/pq v1.10.3 // indirect
+	github.com/ziutek/mymysql v1.5.4 // indirect
+)