學(xué)習(xí)筆記，僅供參考
學(xué)習(xí)書目：《R語言空間數(shù)據(jù)處理與分析實(shí)踐教程》–盧賓賓;

---

準(zhǔn)備工作

設(shè)置工作路徑,并導(dǎo)包

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

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

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

#設(shè)置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)

關(guān)于EPSP的小知識(shí)

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

中國地圖比如對(duì)于中國來講，以地球的幾何球心為中心的地圖就是EPSG:4479，以地球的橢球焦點(diǎn)為中心就是EPSG:4480，此外還有EPSG:4490，因?yàn)檫x擇不同的座標(biāo)系對(duì)于油氣勘探的成本至關(guān)重要，所以有不同的座標(biāo)系。

世界地圖EPSG:4326在世界地圖方面，EPSG:4326是比較著名的一個(gè)，因?yàn)橛擅绹鲗?dǎo)的GPS系統(tǒng)就是在用它，它還有一個(gè)名氣更大的別名叫作WGS84，WGS(World Geodetic System)是世界大地測(cè)量系統(tǒng)的意思，由于是1984年定義的，所以叫WGS84，之前的版本還有WGS72、WGS66、WGS60。EPSG:3857另一個(gè)比較知名的編碼是EPSG:3857，這也是一張世界地圖，目前主要是各大互聯(lián)網(wǎng)地圖公司以它為基準(zhǔn)，例如Google地圖，Microsoft地圖都在用它。

空間數(shù)據(jù)導(dǎo)出

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

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。