在Repast Simphony中读取并显示一个或多个光栅文件

Question

我对在就餐方面使用GIS还不熟悉。代理必须从多个shapefile和两个光栅文件（基于lat和long）中获取像素值，以确定它们在每个滴答声中的移动路线。到目前为止，我首先创建了一个地理，然后将光栅覆盖添加到地理中。

我的问题如下：

>

我的代理将在这种地理环境中移动，所以我是否需要创建可写的网格覆盖（针对光栅和形状文件），然后将其添加到地理环境中？

是否可以向地理添加2个光栅文件和3-4个形状文件？

如何从光栅文件中读取数据？例如，植被光栅文件中的NDVI指数。（我正在看光栅图层类。）

上述所有文件都是按年创建的。所以，我需要根据适当的刻度重新创建这些层。这种方法正确吗？但是，我的主要问题是，如何在就餐时显示光栅数据？一些资源会很有用。

我用过这个组织。geotools软件包，用于读取geotif格式的光栅文件。

GeoTiffReader geoTiffReader = new GeoTiffReader(new File(filename));
GridCoverage2D coverage = (GridCoverage2D) geoTiffReader.read(null);
geography.addCoverage("coverage1", coverage);

现在，没有错误消息，所以我假设代码正确运行（显示打印stmts）并读取光栅文件。不过，显示光栅数据会很好。

Answer 1

作为参考，我建议先看一下就餐地理演示模型，即关于覆盖范围的GIS文档(https://repast.github.io/docs/RepastReference/RepastReference.html#gis-以及GeoTools GridCoverage API(http://docs.geotools.org/latest/userguide/library/coverage/grid.html).就餐地理投影和显示与GeoTools GridCoverage 2D实现完全兼容。您可以使用GeoTools阅读器创建示例中描述的覆盖范围，也可以使用RepastCoverageFactory创建可写GridCoverage2D，这只是Repast中GridCoverage2D的自定义实现，允许代理在内存中写入覆盖范围。

如果代理只需要读取覆盖范围，则可以使用GeoTools读卡器加载光栅文件，并通过“地理”指定给地理位置。addCoverage（“MyCoverage”，coverage），你可以随意添加。还可以根据需要从ShapeFile向地理位置添加任意多的代理层。就餐地理演示中也演示了这一点。

一个警告是覆盖范围的光栅数据和代理的形状文件应该使用相同的CRS和投影。Repast Geography假设地理中的所有GIS数据都在同一个投影中。此外，GIS3D显示要求数据在WGS84中。GeoTools提供了动态转换的实用程序，但如果需要，我建议在ArcGIS或QGIS等外部GIS工具中重新投影。

要引用加载的覆盖范围，代理可以调用geography。getCoverage（coverageName），您可以删除地理覆盖范围。移除覆盖（coverageName）。一旦代理拥有覆盖率实例，它就可以通过

double[] value = null;
DirectPosition pos = new DirectPosition2D(loc.x, loc.y);
coverage.evaluate(pos,value);

值的类型将取决于光栅文件中的数据，例如浮点数、int、字节等。值数组将包含与位置对应的每个波段的所有值，因此您需要知道光栅文件是如何索引的——通常规格将从您下载文件的位置提供。对于NDVI光栅文件，索引可能是双精度或浮点数，并且可能是数组中唯一的值。

您当然可以在不同的时间点添加和删除覆盖层，以模拟随时间的变化。只需确保代理行为的顺序考虑到这一点，例如，覆盖范围应在勾号开始时更新，以便代理在每个勾号使用新数据。可以创建负责更新的图层管理器代理。基于勾号删除和添加图层应该不会有任何问题，尽管我还没有用显示器测试过，所以我不确定这是否会对视觉产生影响。

要显示覆盖层，需要使用GIS3D显示，然后在Repast运行时使用显示向导简单地添加覆盖层。可以在一个显示器中包含多个覆盖层。覆盖层需要覆盖样式来确定光栅中的数据如何转换为图像。一些光栅文件在独立的特定颜色贴图上建立索引，您可以先尝试使用DefaultCoverageStyle，这是显示向导中的默认选项。或者，可以提供CoverageStyle类，并使用它指定光栅数据确定显示颜色的方式。CoverageStyle返回GeoTools RasterSymboler实例，BlueCoverageStyle类的地理演示中提供了一个示例。