借助Go语言开发地理信息系统,涵盖坐标系转换、空间关系判断、地图服务发布等核心功能

# 安装Go语言（以Ubuntu为例）
sudo apt install golang-go

# 验证GIS核心库
go get -u github.com/twpayne/go-geom       # 几何对象处理
go get -u github.com/paulmach/orb         # 地理坐标系转换

package main

import (
	"fmt"
	"github.com/go-spatial/proj"
)

func main() {
	// 创建坐标转换器（WGS84转Web墨卡托）
	transformer, _ := proj.NewTransformer(proj.MustNewCRS("EPSG:4326"), proj.MustNewCRS("EPSG:3857"))

	// 故宫博物院坐标（经度, 纬度）
	wgs84Point := proj.Point{X: 116.397, Y: 39.916}
	
	// 执行坐标转换
	webMercatorPoint, _ := transformer.TransformPoint(wgs84Point)
	
	fmt.Printf("转换结果: X=%.2f, Y=%.2f\n", webMercatorPoint.X, webMercatorPoint.Y)
	// 输出: 转换结果: X=12958234.62, Y=4852834.51
}

package main

import (
	"fmt"
	"github.com/twpayne/go-geom"
	"github.com/twpayne/go-geom/xy"
)

func main() {
	// 创建多边形（模拟北京市五环范围）
	polygon := geom.NewPolygon(geom.XY).MustSetCoords([][]geom.Coord{
		{{116.2, 39.8}, {116.5, 39.8}, {116.5, 40.0}, {116.2, 40.0}, {116.2, 39.8}},
	})

	// 待检测点（天安门坐标）
	point := geom.NewPoint(geom.XY).MustSetCoords(geom.Coord{116.397, 39.903})

	// 判断点是否在多边形内
	isInside := xy.IsPointInRing(
		point.Coord(),
		polygon.Coord(0), // 取多边形外环
	)

	fmt.Printf("天安门是否在五环内: %v\n", isInside) // 输出: true
}

package main

import (
	"github.com/gin-gonic/gin"
	"github.com/twpayne/go-geojson"
)

func main() {
	r := gin.Default()

	// 定义GeoJSON端点
	r.GET("/buildings", func(c *gin.Context) {
		featureCollection := geojson.NewFeatureCollection([]*geojson.Feature{
			{
				Geometry: geojson.NewPoint([]float64{116.397, 39.903}),
				Properties: map[string]interface{}{
					"name": "人民大会堂",
					"height": 46.5,
				},
			},
		})

		c.JSON(200, featureCollection)
	})

	r.Run(":8080")
}

func queryBuildings(db *sql.DB) {
	rows, _ := db.Query(`
		SELECT ST_AsGeoJSON(geom), name 
		FROM buildings 
		WHERE ST_Within(geom, ST_MakeEnvelope($1,$2,$3,$4, 4326))`,
		116.3, 39.8, 116.5, 40.0)

	for rows.Next() {
		var geojsonStr string
		var name string
		rows.Scan(&geojsonStr, &name)
		// 解析GeoJSON并处理...
	}
}

// 使用gdal.org/go包读取遥感影像
dataset, _ := gdal.Open("beijing.tif", gdal.ReadOnly)
band := dataset.RasterBand(1)

// 读取指定区域像素值
buffer := make([]float32, 100*100)
band.IO(gdal.Read, 0, 0, 100, 100, buffer, 100, 100, 0, 0)