查找两个栅格的交点

9

我有两个栅格图像，按布尔分类（1或2）。它们在给定区域重叠。

我需要只有重叠部分的栅格。

使用Arcinfo10。Answer只能使用栅格方法，不能转换为多边形，点等。仅供参考，以下是屏幕截图。

arcgis-10.0 raster

— 萨德
source

两个答案都是正确的，第一个答案比较优雅，但是我决定给予第二个荣誉，因为它对我最有帮助。

— 塔德

8

这似乎是一个非常简单的剪辑操作。使用剪辑（数据管理）工具。

从帮助文档中：

现有的栅格或矢量层可用作裁剪范围。如果将要素类用作输出范围，则可以选择按要素类的最小边界矩形或要素的多边形几何来裁剪栅格。如果使用裁剪几何，则可以提高输出的像素深度。因此，您需要确保输出格式可以支持正确的像素深度。

实际上，这应该比您的要求更简单，因为它使用了核心功能，甚至不像其他方法那样需要Spatial Analyst扩展。

--------------编辑

另一方面，如果您的目标是创建一个栅格，以显示两个输入均为null的位置，或者一个输入具有值，或者两个输入都具有值的栅格，则可以使用Raster Multiplication。本质上，这将创建一个输出栅格，该栅格栅格的值是输入栅格的数值的倍数。请注意，以下内容需要Spatial Analyst扩展。

以您的情况为例，使用布尔栅格。

Raster 1:  1,2  
Raster 2:  1,2

Output Raster:  Possible combinations:  1x1,1x2,2x1,2x2 with output values: 1,2,4  
Value 1:  Both Boolean False  
Value 2:  1 or the other, Boolean True  
Value 4:  Both Boolean True

如果要确定哪个输入具有真实值，则需要对输入栅格之一进行重新分类，以使相乘后的值完全唯一。两个输入都需要重新分类，以使值不为“ 1”，因为这将允许在相乘的栅格中有重复的值。

例：

Raster 1: 1,2 - Use [Raster Addition][4], with a constant of 1 - Output Raster 3  
Raster 2: 1,2 - Use [Multiply Raster][5], with constant of 2 - Output Raster 4  
Raster 3: 2,3  
Raster 4: 2,4

Now, multiply Raster 3 and Raster 4  
Output Raster:  Possible Combinations:  2x2,3x2,4x2,4x3 with output values: 4,6,8,12

Value 4:  Both Boolean False  
Value 6:  Raster 3 True, Raster 4 False  
Value 8:  Raster 4 True, Raster 3 False  
Value 12: both Boolean True

您可以将这些操作组合在一起，在其中将一个栅格与另一个栅格一起裁剪，以仅显示重叠区域。然后重新分类并乘以栅格，以显示两个都为False，一个或另一个为True或两个均为True的地方。

希望这可以帮助。

— 获取空间
source

我不断收到输出范围错误。或蒙版错误。太令人沮丧了...。

— 塔德

nvm。我想我有一个投影问题。也就是说，这些功能不支持即时投影...因此，它们在相关功能上不会重叠。

— 塔德

@Thad-然后我将尝试重新投影其中一个栅格以匹配另一个栅格。使用数据管理工具下的项目栅格工具。或者，在您的工作流程中使用python代码以编程方式执行此操作。

— 获得空间

是的，我做到了，而且效果很好。

— 塔德

9

我假设在您的栅格中，白色（以某种方式表示无数据区域）的值为1，蓝色/绿色的值为2。

因此，您需要获得两个栅格都具有值2的区域。

使用ArcGIS中的栅格计算器可以轻松地执行此操作。

表达式将如下所示：

("raster1" == 2) & ("raster2" == 2)

要么

Con(("raster1" == 2) & ("raster2" == 2), 1, 0)

PS此工具需要Spatial Analyst扩展

— 亚历克斯·马可夫
source

0

我认为这种方法可以帮助您：

arctoolbox >> Data management tools >> Raster >> Rasterprocessing >> Clip 将栅格拖放到输入和输出中（无论是哪一个）。现在，您有一个来自重叠区域的shapefile（我们称其为overlay.shp）。
对于两个栅格，均使用按遮罩提取 arctoolbox >> Spatial analyst tools >> Extraction >> Extract by mask：：输入：栅格1，遮罩数据：overlay.shp。现在，您有了一个重叠栅格1的栅格（我们称它为overlay1）。 arctoolbox >> Spatial analyst tools >> Extraction >> Extract by mask：输入：栅格2.遮罩数据：overlap.shp。现在，您有了重叠栅格1的栅格（我们称其为overlay2）。
最后一步：：arctoolbox >> Spatial analyst tools >> Map algebra拖放overlay2和overlay1。

现在您可以做任何事情。您需要的每个统计信息。

— 无意义的
source