经纬度坐标和笛卡尔坐标的相互转换--r语言

写论文想找GPS坐标转笛卡尔坐标的方法，搜了一圈，大多是通过编程、arcgis等方式，转换方式很复杂，小叶子现将发现的一种简单的方法分享给大家。

第一步，将你的经纬度坐标按图1所示整理，lon代表经度，lat代表纬度。ps：我一般是将excel表格保存成csv格式，习惯啦。后面需要用到“SoDA"包，先提前下载好噢。

图1

第二步 当然是直接上手啦。

```

library("SoDA")

#导入经纬度坐标数据

spa<-read.csv("behticspa.csv",header = T)

xy <- geoXY(spa$lon, spa$lat, unit = 1)

xy

#输出到文件夹spa

write.csv(xy,file = "spa.csv")

```

图2

好啦，如图2所示，转换完成啦。

```

plot(xy[,1], xy[,2], asp = 1)

```

图3就是笛卡尔坐标在平面上的图形啦

图3