3
将复杂的shapefile划分为网格
我有一个具有多边形/多多边形特征的相当详细的shapefile（文件约为500mb）。它实际上是整个世界的shapefile，其特征代表海岸线。我需要使用网格划分此数据。明确地说，我不想对数据进行“排序”，而是将多边形实际切成小块。我知道之前曾有人问过这个问题，但是我发现的解决方案对我不起作用。 我试过了： 使用QGIS并将我的shapefile内容与矢量网格相交-结果非常糟糕。大多数主要陆地都神奇地消失了，尽管有时似乎有较小的一块土地。我应该注意，这种方法对简单得多的数据（例如，较少的点）非常有效 使用OGR的相交工具。我通过ogr2ogr甚至通过滚动自己的C ++工具都进行了尝试。两者都具有与QGIS相同的问题。对于简单的文件，它们也不会出现此问题，但是对于较复杂的文件，则不会通过。作为参考，我使用的是澳大利亚和新西兰的shapefile，大小不超过20mb，并且QGIS和OGR都无法对其进行“网格化”。 有人建议一次使用PostGIS，因为它具有交集功能-但PostGIS的ST_Intersect使用与OGR相同的GEOS后端。实际上，据我所知，它们都调用相同的函数，因此我认为PostGIS不会产生不同的结果。 我一直在寻找关于其他尝试方法的建议。我需要一个功能强大的应用程序或工具包，可以将高度详细的shapefile划分为图块。 编辑：添加一些更多的信息 针对Simbamangu： shapefile基本上是来自OpenStreetMap的海岸线数据。这是“ processed_p”文件的合并版本（因此不会拆分为图块），我通过给他们的开发人员列表发送了电子邮件。请注意，他们将瓦片分割成100km x 100km重叠的块并不一定是我想要的-我不希望重叠，我想自由选择网格大小，或者我只想使用默认处理的_p。 默认情况下，海岸线数据具有QGIS报告的几何误差。我使用一个小的工具将这些错误汇总在一起，使用发现的一些代码专门解决了这个问题（修复海岸线数据中的几何错误：https : //github.com/tudelft-gist/prepair）。使用此工具在文件上运行几乎可以修复QGIS发现的所有错误。我仅在清理文件后尝试进行相交。 我使用QGIS所做的正是：打开数据以确保它在QGIS中看起来不错。尝试通过使用具有指定间距的Vector Grid创建一层瓷砖，然后将这两层相交，将其分成多个瓷砖-不行。尝试使用较小的数据集-选择大洋洲（Aus，NZ）中的要素以尝试使用较小的数据集-此形状文件的大小小于20mb。再次尝试将其分割，不起作用。 我对OGR所做的事情：ogr2ogr直接使用带有spat_extent的'-spat'和'-clipsrc'选项。还写了一个可以在WKT上运行的小型C ++工具，因此我使用ogr2ogr将shapefile转换为WKT，然后将文本文件提供给我的应用程序。它遍历该文件并调用此处记录的Intersection（）方法：http : //www.gdal.org/ogr/classOGRGeometry.html。我认为最终结果与直接使用ogr2ogr完全一样。 针对布伦特： 是的 一切都在WGS84纬度/经度中 我以为是相反的事实-对于给定的一组网格图块，与一个巨型多面体相交将花费更长的时间，而不是一堆可以在空间上更局部地定位于每个图块的零碎要素，但这是一个有趣的建议-我会尝试并报告。 在此过程中不会保留任何属性字段，我只对几何感兴趣。 我不确定，但是我想你是说我应该选择与给定网格图块重叠的多边形，然后执行相交。使用QGIS手动操作太麻烦了。我的工具已经通过边界框检查在一定程度上做到了这一点。速度有所提高，但是最终结果仍然很差，并且没有明显的不同。 这不是一个选择。现在，我正在尝试将数据划分为1度经度x 1度经度，并且我正在寻找一种适用于所有情况的通用/稳健方法。我尝试增加网格尺寸（即10x10）以查看是否可以获得更好的结果，并且看不到网格尺寸与输出质量之间的任何关联。 编辑＃2： 我尝试了更多，总的来说，使用GEOS和QGIS（使用fTools，我不知道是否又使用GEOS）时，结果似乎都不可靠。我错误地指出网格的大小与结果无关—网格越大，结果越好（这是很好的认识，但仍然不是解决方案）。这是一个实际间隔开的网格的屏幕截图，该网格大部分可用，但在一个磁贴中部分失败： 几何图形很干净-QGIS使用“检查有效性”工具显示0错误。我不是要逐步解决这个问题。验证某些特征是否在视觉上不明显（且不会使用较小的图块）的情况下，无法通过如此大的数据集的交集不可行。