GIS之arcgis系列10:arcpy实现批量掩膜提取

按掩膜提取 (Spatial Analyst)

提取掩膜所定义区域内的相应栅格像元。

OutRas = ExtractByMask(InRas1, InMsk1, "INSIDE")

使用情况

  • 输入栅格中的其他属性(若有的话)将按照原样添加到输出栅格属性表。 根据所记录的属性,某些属性值可能需要重新计算。

  • 将多波段栅格指定为输入栅格(Python 中的 in_raster)值时,将使用所有波段。

    要处理一系列来自多波段栅格的波段,请首先使用波段合成工具创建由这些特定波段组成的栅格数据集。 然后将结果用作输入栅格(Python 中的 in_raster)值。

    默认输出格式为地理数据库栅格。 如果将 Esri 格网堆栈指定为输出格式,则堆栈的名称不能以数字开头、不能使用空格且长度不得超过九个字符。

  • 当为输入栅格掩膜指定多波段栅格时,将只在运算中使用第一个波段。

  • 当输入栅格(Python 中的 in_raster)值和输入栅格或要素掩膜数据数据(Python 中的 in_mask_data)具有相同的像元大小,且像元对齐,则可将其直接应用于工具中。 在工具操作期间,系统不会对其进行内部重采样。

    如果像元大小不同,则输出像元大小将为输入的最大值,而输入栅格值将在内部用作捕捉栅格。 如果像元大小相同,但是像元未对齐,则输入栅格值将在内部用作捕捉栅格。 在执行提取操作之前,上述任一情况都会触发内部重采样。

    有关详细信息,请参阅像元大小捕捉栅格环境主题。

  • 默认情况下,如果掩膜输入为一个要素,则可使用像元大小和输入栅格值中的像元对齐(捕捉栅格)从内部将其转换为栅格。

  • 运行按掩膜提取工具时,如果在环境设置中指定了掩膜,则输出栅格仅包含位于环境掩膜和输入掩膜数据交集内区域的像元值。

  • 您可以使用分析范围(Python 中为 analysis_extent)参数为独立工具操作明确指定输出分析区域,或将覆盖环境设置作为工作流的一部分。 您可以通过输入值、选择显示范围、选择图层或浏览至输入数据集的方式指定范围。

  • 默认分析范围值是从输入栅格值和输入栅格或要素掩膜数据值的交集计算得出的。

  • 如果分析范围没有明确指定为参数值,则其会从分析环境设置中派生。

  • 有关适用于此工具的地理处理环境的详细信息,请参阅分析环境和 Spatial Analyst


完整代码如下:

python 复制代码
import os
import arcpy

arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput = 1

input_dir = arcpy.GetParameterAsText(0)
output_dir = arcpy.GetParameterAsText(1)
mask = arcpy.GetParameterAsText(2)

tif_files = []
for dirpath, dirnames, filenames in os.walk(input_dir):
    for filename in filenames:
        if filename.endswith('.tif'):
            tif_files.append(os.path.join(dirpath, filename))

error_files = []

for i, tif_file in enumerate(tif_files, start=1):

    relative_path = os.path.relpath(tif_file, input_dir)

    base_name, ext = os.path.splitext(relative_path)
    out_file = os.path.join(output_dir, base_name +ext)

    if not os.path.exists(os.path.dirname(out_file)):
        os.makedirs(os.path.dirname(out_file))
        arcpy.gp.ExtractByMask_sa(tif_file, mask, out_file)
相关推荐
疯狂学习GIS2 小时前
ArcGIS填补面图层的细小空白并删除主体部分外的零散部分
arcgis·gis·学术工作效率·gis数据
sinat_384241095 小时前
在有网络连接的机器上打包 electron 及其依赖项,在没有网络连接的机器上安装这些离线包
javascript·arcgis·electron
GIS 数据栈6 小时前
每日一书 《基于ArcGIS的Python编程秘笈》
开发语言·python·arcgis
GIS思维21 小时前
ArcGIS定义投影与投影的区别(数据和底图不套合的原因和解决办法)
arcgis·gis·地理信息·arcgis坐标系·动态投影
阮少年、1 天前
ArcGIS Pro的arpx项目在ArcGIS Server中发布要素服务(FeatureServer)
arcgis·arcgis server·gis服务发布·要素服务
圆弧YH1 天前
Arcgis 地图制作
arcgis
sinat_384241091 天前
使用 cnpm 安装 Electron,才是正确快速的方法
arcgis
中科GIS地理信息培训2 天前
ArcGIS Pro 3.4新功能1:唯一值符号化增加复选框,可在内容窗格和布局视图中控制类别的可见性。
arcgis·arcgis pro
黄鹂绿柳2 天前
学了Arcgis的水文分析——捕捉倾泻点,河流提取与河网分级,3D图层转要素失败的解决方法,测量学综合实习网站存着
arcgis
科研online2 天前
ArcGIS计算水库库容量
arcgis·学习方法