35得票8回答
readOGR()无法打开文件。

wmap <- readOGR(dsn="~/R/funwithR/data/ne_110m_land", layer="ne_110m_land") 这段代码没有加载形状文件,会生成错误Error in ogrInfo(dsn = dsn, layer = layer, encodin...

26得票4回答
安装rgdal遇到问题

我想为"R version 3.2.3 (2015-12-10)"安装rgdal。我已经下载并安装了: GDAL 1.11完整版 PROJ框架v4.9.2-2 GEOS框架v3.5.0-1 从KyngChaos下载。 然后在RStudio中输入install.packages("rg...

20得票3回答
PROJ4升级到PROJ6和“废弃数据”警告

Context 这些问题与从PROJ4升级到PROJ6所引起的更改及其对各种R空间包(sp,sf,raster)的影响有关。 我们现在收到了很多关于“已放弃基准”的警告,看起来有点令人担忧,我有点困惑我应该做什么。我可以看到这可能会在某些情况下产生可怕的后果,并且在其他情况下可以忽略它们。 看...

20得票4回答
如何在使用fortify()函数后保留shapefile中的信息?

如何在 shapefile 后保留多边形信息?让我试着解释一下: 我有一个包含以下数据的 shapefile:> head(mapa@data) ID CD_GEOCODI TIPO CD_GEOCODB NM_BAIRRO CD_GEOCOD...

16得票1回答
空间多边形 - 在R中从坐标创建一组多边形

我正在尝试从保存为X,Y格式的顶点位置创建一组多边形。 这是我的数据示例 - 每行表示一个多边形的顶点。这些多边形是正方形。square <- rbind(c(255842.4, 4111578, 255862.4, 4111578, 255862.4, 4111558, ...

15得票5回答
将纬度和经度点转换为UTM

我找到了一个相当简单的例子来演示如何做这件事,但我无法让它为我工作。我对R比较陌生。library(rgdal) xy <- cbind(c(118, 119), c(10, 50)) project(xy, "+proj=utm +zone=51 ellps=WGS84") ...

14得票3回答
使用readOGR时如何编码德语umlauts

我试图使用rgdal包提供的readOGR函数读取OGR矢量地图,但是在处理德语umlauts时遇到了一些麻烦。下面是数据示例,像ö这样的umlauts被替换为\303\266。map <-readOGR("/path/to/data.gdb", layer = "layer") map...

14得票4回答
使用readOGR和readShapePoly读取shape文件的区别

我使用maptools包中的readShapePoly函数读取了一个shapefile文件,但无法使用readOGR函数读取同一文件。我希望有人能帮助我使用readOGR函数读取这个shapefile文件。 我从这里下载了文件orcounty.shp:http://geography.uor...

12得票1回答
在一个矢量文件的背后绘制一个栅格图像

我该如何在一个"raster"对象后面绘制一个shapefile对象?它们各自单独绘制得很好,但点无法覆盖在栅格上:require(rgdal) require(maptools) require(raster) myproj = "+proj=utm +zone=12 +north +el...

11得票1回答
在Linux上安装rgeos包时出现错误 [R]

我正在尝试在Linux上安装rgeos软件包。但是我遇到了以下错误:system("sudo apt-get update") system("sudo apt install libgdal-dev -y") install.packages("rgeos") collect2: error...