arcpy之Cursor访问矢量-进行批量裁剪

最新推荐文章于 2024-10-17 16:38:46 发布

云朵不吃雨

最新推荐文章于 2024-10-17 16:38:46 发布

阅读量211

点赞数 1

CC 4.0 BY-SA版权

分类专栏： arcgis 文章标签： arcgis

本文链接：https://blog.youkuaiyun.com/d8dongdong/article/details/134488791

arcgis 专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一个使用Python和ArcGIS库进行矢量裁剪的详细教程，通过SearchCursor、MakeFeatureLayer_management和ExtractByMask函数实现按形状文件裁剪影像。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

系列文章目录

前言

Cursor 是一种数据访问对象，可用于在表中迭代一组行或者向表中插入新行。游标有三种形式：搜索、插入或更新。游标通常用于读取和更新属性。
本章内容主要使用的函数包括：
arcpy.SearchCursor：创建查找Cursor
SearchCursor 函数用于在要素类或表上建立只读游标。SearchCursor 可用于遍历 Row 对象并提取字段值。可以使用 where 子句或字段限制搜索，并对结果排序。

arcpy.MakeFeatureLayer_management：创建图层文件
arcpy.SelectLayerByAttribute_management：按图层属性进行选择
arcpy.sa.ExtractByMask:按掩膜进行裁剪

一、完整矢量裁剪代码

代码如下：

# -*- coding: utf-8 -*-

import arcpy
from arcpy.sa import *
def batchclipimgbyshp(img_dir, shp_dir,save_dir):
    """

    :param img_dir: 待裁剪影像路径
    :param shp_dir: 裁剪矢量路径
    :param save_dir: 裁剪后影像输出路径
    :return: None
    """
    rows = arcpy.SearchCursor(r"E:\changdong\HENANSample\cloudandshadow.shp")
    # 创建矢量图层
    layer = arcpy.MakeFeatureLayer_management(shp_dir, "lyr_shp")

    for row in rows:
        arcpy.SelectLayerByAttribute_management("lyr_shp", "NEW_SELECTION", ' "FID" = {} '.format(row.getValue("FID")))
        print(row.getValue("FID"))

        # 导出并查看矢量

        print(row.getValue("clsid"))
        outExtractByMask = ExtractByMask(img_dir, layer)
        outExtractByMask.save(save_dir.format(row.getValue("FID"), row.getValue("clsid")))
        arcpy.SelectLayerByAttribute_management("lyr_shp", "CLEAR_SELECTION")
        del outExtractByMask


if __name__ == '__main__':
    # :param img_dir: 待裁剪影像路径
    # :param shp_dir: 裁剪矢量路径
    # :param save_dir: 裁剪后影像输出路径
    # :return: None
    arcpy.env.overwriteOutput = True
    img_dir = r'D:\temp\oCB04A_20230709_L1A0000482611_8bit.tif'
    shp_dir = r"E:*\cloudandshadow.shp"
    save_dir = r"E:*\fid{}_clsid{}.tif"

    batchclipimgbyshp(img_dir, shp_dir,save_dir)