學習筆記，僅供參考
學習書目：《R語言空間數據處理與分析實踐教程》–盧賓賓;

---

準備工作

設置工作路徑,并導包

workL = "F:/MyStudio/Rstudio/RSpaceMetrology/myRdoc/C3" setwd(workL) getwd()

讀取ESRI shapefile(.shp)空間格式數據

注意，如果這里報錯誤，就需要檢查一下工作目錄中是不是缺少其他文件(比如.dbf, .shx),如果缺失,就將這些文件放到工作目錄下(),再次導入.

#設置EPSG:27700 LNHP <- readShapePoints("LNHP", verbose=T,proj4string = CRS("+init=epsg:27700")) #summary(LNHP) plot(LNHP) LNNT <- readShapeLines("LNNT", verbose=T,proj4string = CRS("+init=epsg:27700")) #summary(LNNT) plot(LNNT)LN.bou <- readShapePoly("LondonBorough", verbose=T,proj4string = CRS("+init=epsg:27700")) #summary(LN.bou) plot(LN.bou)

關于EPSP的小知識

EPSP的英文全稱是European Petroleum Survey Group,中文名稱為歐洲石油調查組織。這個組織成立于1986年，2005年并入IOGP(International Association of Oil & Gas Producers)，中文名稱為國際油氣生產者協會。它負責維護并發布坐標參照系統的數據集參數，以及坐標轉換描述，該數據集被廣泛接受并使用，通過一個Web發布平臺進行分發，同時提供了微軟Acess數據庫的存儲文件，通過SQL 腳本文件，mySQL, Oracle 和PostgreSQL等數據庫也可使用。
目前已有的橢球體，投影坐標系等不同組合都對應著不同的ID號，這個號在EPSG中被稱為EPSG code，它代表特定的橢球體、單位、地理坐標系或投影坐標系等信息。EPSG對世界的每一個地方都制定了地圖，但是由于座標系不同，所以地圖也各不相同。

中國地圖比如對于中國來講，以地球的幾何球心為中心的地圖就是EPSG:4479，以地球的橢球焦點為中心就是EPSG:4480，此外還有EPSG:4490，因為選擇不同的座標系對于油氣勘探的成本至關重要，所以有不同的座標系。

世界地圖EPSG:4326在世界地圖方面，EPSG:4326是比較著名的一個，因為由美國主導的GPS系統就是在用它，它還有一個名氣更大的別名叫作WGS84，WGS(World Geodetic System)是世界大地測量系統的意思，由于是1984年定義的，所以叫WGS84，之前的版本還有WGS72、WGS66、WGS60。EPSG:3857另一個比較知名的編碼是EPSG:3857，這也是一張世界地圖，目前主要是各大互聯網地圖公司以它為基準，例如Google地圖，Microsoft地圖都在用它。

空間數據導出

writePointsShape(LNHP, fn="LNHP_w") writeLinesShape(LNNT, fn="LNNT_w") writePolyShape(LN.bou, fn="LondonBorough_w")

總結

以上是生活随笔為你收集整理的R语言空间数据处理(part2)--空间数据读写的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。