R语言空间数据处理(part2)--空间数据读写
學(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)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 幽灵行动荒野怎么隐身
- 下一篇: R语言观察日志(part10)--fil