大型shapefile到栅格

我有一个很大的shapefile（1 gb），需要对其进行栅格化。我已经尝试了以下方法。

1.）使用v.in.ogr此失败将其导入GRASS，并显示以下错误消息：错误：G_realloc：无法在break_polygons.c：188上分配498240036字节

2.）我的第二个想法是使用PostGIS。导入shapefile，在网格的x，y位置对其重新采样，然后导出这些点并从xyz创建网格。我成功导入了shapefile（多边形和点），但是将300万个点与100万个多边形相交似乎很慢。我使用了以下PostGIS条件，也许还有改进的余地。

select polygons.land_id,grid.geom from grid,polygons where grid.geom && polygons.geom and within(grid.geom,polygons.geom)

3.）我也曾尝试simplify()在PostGIS中使用。但是我迷失了许多小多边形（即某些仅由小多边形覆盖的区域变为空）。

任何想法将不胜感激。

您可以尝试gdal_rasterize，尽管我还没有使用过这么大的shapefile，所以您可能会遇到与GRASS相同的问题。我认为类似以下内容的东西应该可以工作（GDAL> = 1.8.0）：

gdal_rasterize -a AN_ATTRIB -l THE_LAYER -a_nodata -9999 -a_srs EPSG:27700 -co TILED=YES -tr 10 10 -ot Float32 src.shp dest.tif

当然，您将不得不根据源shapefile尝试一些选项。最重要的参数是-tr，它指定像素的分辨率；没有它，您可能会发现自己的栅格很大。

如果要坚持使用GRASS，请尝试为栅格化设置较小的范围，然后将过程分成可管理的块，然后将栅格拼接成一个。

关于您的初衷，您可以尝试在具有更多RAM或交换磁盘空间的计算机上运行v.in.ogr命令吗？

如果没有，您可以将其拆分为多个文件，然后在重新合并之前对其栅格化。

如果在进行栅格化时需要更多控制，请查看perrygeo的poly_density.py脚本，该脚本在后台使用GDAL，但可用于重叠要素或添加条件评估，而超出了gdal_rasterize单独评估的范围。

就像开膛手杰克一样，分部分完成。从postgis导出切片的数据导入草并进行转换。

这是有关如何在PostGIS中栅格化矢量表的指南 ...