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)
相关推荐
新中地GIS开发老师14 小时前
WebGIS和WebGL的基本概念介绍和差异对比
学习·arcgis·webgl
王解14 小时前
Jest项目实战(2): 项目开发与测试
前端·javascript·react.js·arcgis·typescript·单元测试
GISerQ.2 天前
ArcGIS/QGIS按掩膜提取或栅格裁剪后栅格数据的值为什么变了?
arcgis·栅格数据·按掩膜提取
坫鹏加油站2 天前
【ArcGIS】绘制各省碳排放分布的中国地图
arcgis
前端 贾公子2 天前
koa项目实战 == 实现注册登录鉴权
arcgis
jyan_敬言4 天前
【Linux】centos7安装配置DHCP、FTP服务
linux·运维·服务器·c语言·开发语言·arcgis·青少年编程
yngsqq5 天前
arcgis坐标系问题
arcgis
GIS思维5 天前
ArcGIS影像调色(三原色)三原色调整
arcgis
新中地GIS开发老师6 天前
地理信息科学专业想搞GIS开发:学前端还是后端?
前端·javascript·arcgis·前端框架·cesium
Z_W_H_7 天前
【ArcGISPro】宣布推出适用于 ArcGIS 的 AI 助手
人工智能·arcgis·ai