rgdal

    0熱度

    2回答

    我有一個來自UTM區域33北部的大緯度經度信息。 我嘗試下面的命令來這​​個地理信息轉換爲UTM座標(我的數據集對象最初稱爲S3km): library(rgdal) UTM33N<-"+proj=utm+zone=33+north" UTM33N<-paste(UTM33N,"+ellps=WGS84",sep="") UTM33N<-paste(UTM33N,"+datum=WGS84"

    3熱度

    2回答

    我想在R中的光柵文件上繪製shapefile,但我無法使它們完美重疊:光柵似乎逆時針旋轉幾度。這是投影問題嗎? 請考慮以下MWE library(raster) library(rgdal) # Download from http://biogeo.ucdavis.edu/data/gadm2/shp/ITA_adm.zip shape_file = "ITA_adm1.shp" #

    1熱度

    2回答

    我想從shapefile中獲取緯度和經度。直到現在,我只知道如何閱讀shapefile。 library(rgdal) centroids.mp <- readOGR(".","35DSE250GC_SIR") 但我怎麼能從centroids.mp中提取經緯度?

    0熱度

    2回答

    我從1983年開始在國家平面座標系統中具有以下數據集。所有座標位於長島區(3104)。 dput(example) structure(c(1008031L, 1000852L, 1001869L, 1005306L, 986887L, 998031L, 1018703L, 1014319L, 1016186L, 1006977L, 1006891L, 1000883L, 1001403L,

    4熱度

    1回答

    我試圖在Rstudio(0.97.332)中使用Ubuntu 12.10(quantal)和R 2.15.1(烤棉花糖)安裝rgdal庫。 我在Ubuntu的安裝以下關於@中: $ apt-cache search gdal dans-gdal-scripts - GDAL contributed tools by Geographic Information Network of Alaska

    7熱度

    1回答

    我試圖做同樣的事情在這個問題Cartogram + choropleth map in R問,但從一個SpatialPolygonsDataFrame開始,並希望結束了相同類型的對象。 我可以將對象保存爲shape文件,使用scapetoad,重新打開它並將其轉換回來,但我寧願將它全部放在R中,以便過程完全可重現,並且可以自動編寫幾十個變體。 我已經在Github上分叉了Rcartogram代碼,

    2熱度

    1回答

    我成功導入了一個帶有readOGR的多邊形shapefile文件,但未導入.prj文件中的信息。 測試 類:SpatialPolygonsDataFrame n功能:19407 程度:35551.4,1585917,6318047,9408727(XMIN,XMAX,YMIN,YMAX) 座標。 REF。 :NA 摘要(測試)類SpatialPolygonsDataFrame 的 物體座標: 最小

    1熱度

    1回答

    我想知道如何根據特定顏色的矩陣編寫一個GTiff?這可能是顏色不代表數量變量,而是定性空間現象,如人口統計集羣。如果我有例如一個1000x1000的顏色值矩陣(像'#FF0000'這樣的字符串),我可以很容易地用grid.raster(mtrx, interpolate=FALSE)來繪製這些矩陣,但這很難與其他形狀文件疊加。我也可以將它們轉換爲一個SpatialPixelsDataFrame,字

    1熱度

    1回答

    我試圖將地理座標(度)轉換爲UTM座標(米),並不斷收到錯誤消息,指出「檢測到非有限變換」。你知道我該如何解決這個問題嗎?下面是我使用的代碼: > GPS.Points <- Gomer.Data[, c('Longitude', 'Latitude')] > head(GPS.Points) Longitude Latitude 1 23.85474 -19.52211 2 2

    11熱度

    4回答

    我發現了一個相當簡單的例子,說明如何做到這一點,但我不能讓它爲我工作。我很新的R library(rgdal) xy <- cbind(c(118, 119), c(10, 50)) project(xy, "+proj=utm +zone=51 ellps=WGS84") [,1] [,2] [1,] -48636.65 1109577 [2,] 213372.05 554