java - Geotools:GridCoverage 具有值 x 的总面积

Question

我正在使用 Java 进行 Geotools 工作。到目前为止，我有一个 GridCoverage2D 和一个几何图形列表。GridCoverage2D 是一种数字高程模型，源自 geotiff。到这里为止一切正常。

现在我想获取高程具有特定值的每个多边形的面积。例如，对于这个几何图形，我想知道海拔为 27 m 的总面积。我怎样才能实现这个目标？

我不知道如何开始:(

我想到的两个选择:

• 将几何图形分割成小部分(我该怎么做)，对于每个点获取中心点(我可以做到这一点)，然后评估这个 GridCoverage2D。这样我就有了一个包含非常小的几何图形的列表，以及与每个几何图形相对应的高程。进一步需要一些阵列魔法。这是一个好的计划吗？如何快速将几何图形分割成小部分？
• 使用过滤器/查询。但我不知道这是如何工作的，教程也没有帮助。是否有可能通过过滤器实现我想要的效果？

Answer 1

再见，另一种选择是使用区域统计操作或什至更好的过程。

该操作是覆盖模块的一部分:

https://github.com/geotools/geotools/blob/12.x/modules/library/coverage/src/main/java/org/geotools/coverage/processing/operation/ZonalStats.java https://github.com/geotools/geotools/blob/12.x/modules/library/coverage/src/test/java/org/geotools/coverage/processing/ZonalStasTest.java

该进程位于不受支持的进程栅格中: https://github.com/geotools/geotools/blob/12.x/modules/unsupported/process-raster/src/main/java/org/geotools/process/raster/RasterZonalStatistics.java https://github.com/geotools/geotools/blob/12.x/modules/unsupported/process-raster/src/test/java/org/geotools/process/raster/ZonalStatsProcessTest.java

我希望这会有所帮助。