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回答

    这三条线将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对象后面的“栅格”对象?对自己的两个情节不错,但点不超过光栅情节: require(rgdal) require(maptools) require(raster) myproj = "+proj=utm +zone=12 +north +ellps=WGS84 +units=m" shp = readShapeSpatial(fn.shp, proj4stri

    7热度

    1回答

    我有坐标,所有这些应该位于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,

    1热度

    1回答

    我试图将地理坐标(度)转换为UTM坐标(米),并不断收到错误消息,指出“检测到非有限变换”。你知道我该如何解决这个问题吗?下面是我使用的代码: > GPS.Points <- Gomer.Data[, c('Longitude', 'Latitude')] > head(GPS.Points) Longitude Latitude 1 23.85474 -19.52211 2 2