写论文想找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
