GIS经纬度转桩号工具

GIS经纬度转桩号工具

需求

现有一条公路，有公里桩和百米桩对应的经纬度数据；需要根据这个数据得到，输入一个经纬度得到对应的桩号，或者输入一个桩号得到对应的经纬度

实现思路

1. 根据原始数据的公里桩和百米桩对应的经纬度数据，打断获取到一米对应一个经纬度，即每米的桩号对应的经纬度
2. 根据第一步得到的结果值，编写查询服务接口，输入经纬度得到桩号、输入桩号得到经纬度。

具体实现

var stepNum int32 = 1func main() {fileName := "./zh_data_02.csv"// 读取csv文件zhDataList := readCSV(fileName)// 按照一米打断resultDataList := insertPoint(zhDataList)//保存成csv文件saveToFile(resultDataList)}

• 读取原桩号CSV文件

func readCSV(fileName string) []bean.CsvZHData {zhDataList := make([]bean.CsvZHData, 0)fs, err := os.Open(fileName)if err != nil {log.Fatalf("can not open the file, err is %+v", err)}defer fs.Close()r := csv.NewReader(fs)//针对大文件，一行一行的读取文件for {row, err := r.Read()if err != nil && err != io.EOF {log.Fatalf("can not read, err is %+v", err)}if err == io.EOF {break}zhStr := row[0]lonStr := row[1]latStr := row[2]eleStr := row[3]zh, err := strconv.ParseFloat(zhStr, 64)lon, err := strconv.ParseFloat(lonStr, 64)lat, err := strconv.ParseFloat(latStr, 64)ele, err := strconv.ParseFloat(eleStr, 64)zhDataList = append(zhDataList, bean.CsvZHData{ZHStr: zhStr, LonStr: lonStr, LatStr: latStr, EleStr: eleStr,ZH: zh, Lon: lon, Lat: lat, Ele: ele})}return zhDataList
}

• 按1米打断

func insertPoint(zhDataList []bean.CsvZHData) []bean.CsvZHData {resultZHList := make([]bean.CsvZHData, 0)size := len(zhDataList)for i := 0; i < size; i++ {zhPoint := zhDataList[i]endPoint := bean.CsvZHData{}if i == (size - 1) {endPoint = zhDataList[i]} else {resultZHList = append(resultZHList, zhPoint)endPoint = zhDataList[i+1]}startLon := zhPoint.LonstartLat := zhPoint.LatstartZH := zhPoint.ZHendLon := endPoint.LonendLat := endPoint.LatendZH := endPoint.ZHdistance := endZH - startZHazimuth := mytools.ComputeAzimuth(startLon, startLat, endLon, endLat)if distance <= 1000 { // 只有百米桩在计算，其他超过1百米情况不计算for addLen := 0.0; addLen < distance-1; {addLen = addLen + 1stepLenKM := addLen / 1000lonlat := mytools.ConvertDistanceToLogLat(startLon, startLat, stepLenKM, azimuth)lonlatArr := strings.Split(lonlat, ",")lon, _ := strconv.ParseFloat(lonlatArr[0], 64)lat, _ := strconv.ParseFloat(lonlatArr[1], 64)zh := startZH + addLen//startLon = lon//startLat = latresultZHList = append(resultZHList, bean.CsvZHData{Lon: lon, Lat: lat, Ele: 0, ZH: zh})}} else {}//resultZHList = append(resultZHList, endPoint)}resultZHList = append(resultZHList, zhDataList[len(zhDataList)-1])return resultZHList
}

• 保存成CSV

func saveToFile(zhList []bean.CsvZHData) {OUT_FILE_NAME := "over_zh_data_02.csv"fileCSV, err := os.OpenFile(OUT_FILE_NAME, os.O_CREATE|os.O_TRUNC|os.O_RDWR, 0666)if err != nil {log.Fatal(err)}// 关闭文件defer fileCSV.Close()for _, zhBean := range zhList {//fileCSV.WriteString(string(result) + "," + "rn")zhStr := strconv.FormatFloat(zhBean.ZH, 'f', -1, 32)lonStr := strconv.FormatFloat(zhBean.Lon, 'f', -1, 64)latStr := strconv.FormatFloat(zhBean.Lat, 'f', -1, 64)eleStr := strconv.FormatFloat(zhBean.Ele, 'f', -1, 64)fileCSV.WriteString(zhStr + "," + lonStr + "," + latStr + "," + eleStr + "rn")}
}

• 经纬度转桩号服务

router.GET("/tozh", func(c *gin.Context) {lonStr := c.Query("lon")latStr := c.Query("lat")precisionStr := c.Query("precision")fmt.Println("lon=", lonStr, "  lat=", latStr)resultObj := make(map[string]interface{})if lonStr == "" || latStr == "" {resultObj["msg"] = "请传入正确的经纬度!"resultObj["code"] = 401c.JSON(200, resultObj)} else {lon, err := strconv.ParseFloat(lonStr, 64)lat, err := strconv.ParseFloat(latStr, 64)if err != nil {fmt.Println(err)}var precision float64if precisionStr == "" {precision = 50}zhGJD02 := myConfig.ZH_GJD_02zhList, _ := myUtils.ReadGJDZH(zhGJD02)zhBean := myUtils.LonlatToZH(bean.LonLatToZHParamBean{Lon: lon, Lat: lat, Precision: precision}, zhList)c.JSON(200, zhBean)}})

• 桩号转经纬度服务

router.GET("/zhto", func(c *gin.Context) {zhStr := c.Query("zh")roadCode := c.Query("roadCode")resultObj := make(map[string]interface{})if zhStr == "" {resultObj["msg"] = "请传入正确的经纬度!"resultObj["code"] = 401c.JSON(200, resultObj)} else {if roadCode == "" {roadCode = "G320"}zhGJD02 := myConfig.ZH_GJD_02_, zhMap := myUtils.ReadGJDZH(zhGJD02)zhint, err := strconv.ParseInt(zhStr, 10, 64)if err != nil {fmt.Println(err)}lonlatBean := myUtils.ZhToLonLat(bean.ZHToLonLatParamBean{RoadCode: roadCode,ZHInt:    zhint,}, zhMap)c.JSON(200, lonlatBean)}})

最终编译成的工具

工具下载