rgdal

4熱度

1回答

我想鞏固我加載到rgdal爲r的形狀文件之間的非nonded交集，但我得到了以下錯誤： "Error: TopologyException: found non-noded intersection between LINESTRING (34.7279 1.59723, 34.7278 1.59729) and LINESTRING (34.7278 1.59723, 3

0熱度

1回答

爲什麼gdal給了我R的變量的不同最大值和最小值？

這三條線將NetCDF文件打開和讀取一個變量： f=open.ncdf("C:\\BR_Ji1.nc") A = get.var.ncdf(nc=f,varid="date",verbose=TRUE) B = get.var.ncdf(nc=f,varid="GPP",verbose=TRUE) [1] "get.var.ncdf: entering. Here is varid:" [

0熱度

1回答

是否有可能將（二進制）柵格數據直接傳遞給空間對象？

我使用curl包的getBinaryURL從Web接收柵格數據（png，tiff，...）並將數據保存到磁盤。此柵格數據稍後用於空間分析，其中我當前使用柵格（和rgdal）數據包加載柵格數據並創建柵格對象。 library(RCurl) url<-"http://demo.mapserver.org/cgi-bin/wms?version=1.1.1&service=WMS&request=Ge

10熱度

1回答

在shapefile背後繪製柵格

如何繪製shapefile對象後面的「柵格」對象？對自己的兩個情節不錯，但點不超過光柵情節： require(rgdal) require(maptools) require(raster) myproj = "+proj=utm +zone=12 +north +ellps=WGS84 +units=m" shp = readShapeSpatial(fn.shp, proj4stri

7熱度

1回答

如何使用rgdal包將NAD 83座標轉換爲經度和緯度？

我有座標，所有這些應該位於DC，但我無法弄清楚如何將它們從NAD 83轉換爲R中的經度和緯度。我在rgdal軟件包中使用spTransform()函數，並獲取有關不符合的數據。 library(rgdal) nad83_coords <- data.frame(x=c(396842.6, 397886.9, 398315.5, 398154.3, 398010.3), y=c(140887.1,