使用R光栅化带有孔的shapefile?
我正试图使用R来栅格化一些下载的形状文件,其中包括有洞的多边形,而R无法识别这些洞 这些孔似乎被指定为hole=T,因此我认为这个问题不是问题中解决的问题: 还可能发生什么?如何解决此问题使用R光栅化带有孔的shapefile?,r,shapefile,rasterize,R,Shapefile,Rasterize,我正试图使用R来栅格化一些下载的形状文件,其中包括有洞的多边形,而R无法识别这些洞 这些孔似乎被指定为hole=T,因此我认为这个问题不是问题中解决的问题: 还可能发生什么?如何解决此问题 我一直在使用rasterize功能,没有对默认设置进行任何更改,并且我所处理的数据都是从这里的shapefile中重新投影的:,例如hemhwdpine.shp您是否知道该文件中的20598个多边形中的哪一个会给您带来问题?这样做有助于我们帮助您(如果您在缩小问题范围的同时没有碰巧发现问题)。FGID 19
我一直在使用
rasterize
功能,没有对默认设置进行任何更改,并且我所处理的数据都是从这里的shapefile中重新投影的:,例如hemhwdpine.shp
您是否知道该文件中的20598个多边形中的哪一个会给您带来问题?这样做有助于我们帮助您(如果您在缩小问题范围的同时没有碰巧发现问题)。FGID 19334的多边形是一个违法者,但我确信它不是唯一的一个。@Josh O'Brien,我刚刚阅读了您对这个问题的回答:,并且意识到我已经使用maptools而不是rgdal阅读了我的ShapeFile。这可能是问题的一部分吗?为什么rgdal是更好的选择?一个更好的原因是它读入并保留层的CRS信息。要亲自查看,请将proj4string()
返回的值与maptools::readShapePoly()
和rgdal::readOGR()读入的对象进行比较。实际上,我在使用raster::shapefile()
中阅读了该文件,它是readOGR()
的一个方便的包装器,并且使用FGID==19334对复杂的多边形进行光栅化时没有遇到任何问题,因此看起来使用readOGR()
可能足以解决您的问题。尝试了切换,但仍然没有找到漏洞。现在想知道这是否与临时文件的磁盘空间有限有关,或者可能是某种内存问题。