【ArcGIS】批量对栅格图像按要素掩膜提取

news2024/9/22 7:34:27

要把一张大的栅格图裁成分省或者分县市的栅格集,一般是用ArcGIS里的按掩膜提取。
在这里插入图片描述
但是有的时候所要求的栅格集量非常大,所以用代码来做批量掩膜(按字段)会非常方便。

import arcpy , shutil , os
from arcpy import env
from arcpy.sa import *


#使用说明
print "开始使用前,请认真阅读使用说明"
print "\n"
print "使用说明 1 :在工作空间内存放所需要掩膜提取的栅格文件。"
print "使用说明 2 :按字段提取的要素一般是.shp文件,在输入路径时务必确保以.shp为结尾。"
print "使用说明 3 :字段名是 使用说明 2 中要素里拥有的字段,一般使用名称作为输入要素。输入之前,务必确保与要素文件中字段名一致。"
print "使用说明 5 :栅格须存放在工作空间内,在这里,只能输入一个栅格,所输入的栅格名不能加文件后缀。"
print "使用说明 4 :输出文件夹作为存储掩膜后的栅格文件,文件夹里是一个栅格集。"
print "\n"


WorkSpace = raw_input("请输入工作空间路径:")
FieldFeature = raw_input("请输入按字段提取的要素路径:")
Field = raw_input("请输入要素字段名:")
Raster = raw_input("请输入需要裁剪的栅格名:")
OutputFile = raw_input("请输入输出文件夹路径:")
Free = raw_input("----------- 按回车开始执行 -----------")

print "\n"




if not os.path.exists(WorkSpace):
    print "工作空间不存在"
else:
    print "工作空间读取成功"
        
    
if not os.path.exists(OutputFile):
    print "输出路径不存在"
else:
    print "输出路径读取成功"



if os.path.exists(WorkSpace) and os.path.exists(OutputFile):         #当工作空间与临时文件夹都存在方能执行
    print "\n"
    print "正在检查程序运行环境,请稍等"
    print "\n"



    arcpy.env.workspace = WorkSpace                         #设置工作空间

    


    TempFile = WorkSpace + "/" + "TempFile"                   #设置临时文件夹路径
    

    if os.path.exists(TempFile):
        Enquire1 = raw_input("已存在临时文件夹TempFile,是否重建临时文件夹:")
        if Enquire1 == "是":
            shutil.rmtree(TempFile)
            print "临时文件夹删除成功"

                
    if not os.path.exists(TempFile):                         
        os.makedirs(TempFile)                                  #创建临时文件夹,分割输出文件夹,掩膜提取裁剪要素文件夹
        print "临时文件夹创建成功"                   
        print TempFile                                           #输出临时文件夹位置
        print "\n"


        print "正在执行程序,请稍等"
        print "\n"



                             

        #分割
        arcpy.Split_analysis(FieldFeature , FieldFeature , Field , TempFile)        #分割要素



        #按掩膜提取
        dirs = os.listdir(TempFile)                  
        for i in dirs:                           
            if os.path.splitext(i)[1] == ".shp":         #提取裁剪要素文件夹中.shp结尾的文件   分离path里的文件名和扩展名,文件名=[0],扩展名=[1]
                masks = TempFile + "/" + str(i)       #裁剪要素mask路径
                print str(i) + " 掩膜要素读取完毕"
                name = i.strip(".shp")                                                #去掉要素名称后缀.shp
                outExtractByMask = ExtractByMask(Raster, masks)                   #掩膜提取
                outpath = OutputFile + "/" + str(name)                           #掩膜提取输出路径
                outExtractByMask.save(outpath)                                  #保存输出文件
                print str(name) + " 掩膜提取执行完毕"
                print "\n"



        print "正在删除临时文件夹TempFile"
        print "\n"
        shutil.rmtree(TempFile)                                #强制删除临时文件夹
        print "临时文件夹TempFile已删除"
        print "\n"


            
        print "-----------按要素字段掩膜提取执行完毕-----------"





    else:
        print "已存在临时文件夹TempFile,无法执行"

else:
    print "执行未能成功"

直接复制到IDLE里运行即可,使用说明写得很详细。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1229217.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【GEE】基于GEE进行非监督学习

1 简介与摘要 之前写了多季节叠加的监督学习,所以这次简单写一个非监督学习吧。。 这次为了简单明了,就不整那么多虚的了,在这里我不叠图层了,有需要的可以参考前一篇博客自己添加输入的图层。 2 制作输入影像 首先&#xff0c…

C++大神之路——环境篇

序 在我还在做后端的时候,当时程序员圈里就有个梗很火,说的是当时几种常用编程语言的鄙视链:做C的鄙视做Java的,做Java的鄙视做C#的,而我很不幸,当时在鄙视链最底层。一开始只是当个笑话听听就算了&#x…

IC卡操作软件支持PN532

IC卡操作软件,在知道卡片密码的情况下,可以对卡片修改数据,格式化清卡,修改UID卡和CUID卡的卡号,锁UFUID卡等 卡片dump文件拖进软件,即可打开文件,编辑修改文件,写卡,就…

一文读懂:为什么GPU比CPU更快?

大家好,我是老猫,猫头鹰的猫。 在过去几十年里,GPU变得越来越流行,尤其是最近ChatGPT大火,背后训练大模型的硬件设备GPU达到了一片难求的地步。 你有没有好奇:为什么必须要用GPU?CPU被淘汰了吗…

Polygon Miden VM中的哈希函数对比

1. 引言 在Polygon Miden VM中,使用了多个不同的哈希函数: 1)“传统”哈希函数,如BLAKE3:对STARK之外的性能进行了优化。2)algebraic哈希函数,如Rescue Prime:对STARK内部优化&…

idea项目中java类名出现带 j 小红点,如何解决?

目录 一、问题描述 二、问题解决方案 1、寻找异常问题 2、解决方案 2.1常规操作方法 2.2 快速操作方法 一、问题描述 一打开idea的java项目,发现所有的文件边上都有带J的大红点 虽然,在 git bash 中进行编译时无异常。 但是视觉上给人的感受就是…

非遗拓鱼丨以传统文化助力守护美丽长江

长江生物多样性极为丰富,是中华民族赖以生存的生命之源,值得人类更多的关注与保护。为传播非遗文化、宣传长江生态保护,积极响应野生动物保护宣传的号召,以非遗手工助力守护长江,11月18日,在重庆市陆海国际…

[解决] 问题:ImportError: cannot import name ‘Callable‘ from ‘collections‘

问题 我在运行yolov8的代码时,出现了ImportError: cannot import name Callable from collections的错误 原因 版本问题:以下collections的方法都在Python3.10版本后被取消了 ["Awaitable", "Coroutine", "AsyncIterable", "A…

数据资产入表,给企业带来的机遇和挑战

作为推动数字经济发展的核心要素,近年来,数据资源对于企业特别是相关数据企业的价值和作用正日益凸显。 数据资产入表之后,能够为企业经营带来实质性的收益。“随着数据资产的纳入,企业的资产也出现了新标的。在资产负债表中&…

mybatisPlus的简单使用

封装实体类 编写Mapper service层 controller层

图片地址GPS经纬度查询

先打开exif图片查询的网站: 改图宝的:https://www.gaitubao.com/exif图虫de的:EXIF信息查看器 (tuchong.com) 将这个地点:51 deg 30 51.90" N, 0 deg 5 38.73" W 修改为:5130 51.90" N, 05 38.73" W 到谷…

优步让一切人工智能化

优步(Uber)的商业模式建立在对数据的颠覆性使用上--通过将双方智能手机的位置数据关联起来,将出租车司机与乘客配对。这意味着,它可以比传统出租车公司更快地安排司机去接乘客,极大地冲击了传统出租车公司的业务。 优步自成立以来&#xff0…

MySQL 之多版本并发控制 MVCC

MySQL 之多版本并发控制 MVCC 1、MVCC 中的两种读取方式1.1、快照读1.2、当前读 2、MVCC实现原理之 ReadView2.1、隐藏字段2.2、ReadView2.3、读已提交和可重复读隔离级别下,产生 ReadView 时机的区别 3、MVCC 解决幻读4、总结 MVCC(多版本并发控制&…

国外客户发开发信怎么发?写外贸邮件方法?

如何给国外客户发开发信?给老外发外贸开发信的技巧? 与国外客户保持联系并积极开展业务合作变得至关重要。一种有效的方式是通过发展信件,这可以帮助建立关系、传递信息和促成商机。蜂邮将介绍如何向国外客户发开发信,以确保您的…

c语言:矩阵交换

题目&#xff1a; 代码和思路&#xff1a; #define _CRT_SECURE_NO_WARNINGS #include<stdio.h>int main() {int n 0;int m 0;int arr[10][10] { 0 }; // 输入行和列scanf("%d%d", &n, &m);int i 0;int j 0;//读取数组for (i 0; i < n; i)…

nginx学习(3)

Nginx 负载均衡 实战案例 实现效果 浏览器地址栏输入地址 http://172.31.0.99/oa/a.html&#xff0c;负载均衡效果&#xff0c;平均 8083 和 8084 端口中 一、配置 1、先创建2个文件夹&#xff0c;并将apache-tomcat-8.5.87解压到tomcat8083和tomcat8084中 &#xff08;或…

Linux中安装部署环境(JAVA)

目录 在Linux中安装jdk 包管理器yum安装jdk JDK安装过程中的问题 验证安装jdk 在Linux中安装tomcat 安装mysql 在Linux中安装jdk jdk在Linux中的安装方式有很多种, 这里介绍最简单的方法, 也就是包管理器方法: 包管理器yum安装jdk Linux中常见的包管理器有: yumaptp…

Nginx部署前端项目

Nginx部署前端项目 1.在nginx官网http://nginx.org/en/download.html &#xff0c;下载稳定版本&#xff1a; 2.解压后&#xff0c;点击根目录中的nginx.exe即可启动Nginx&#xff0c;或是在nginx安装目录中启动cmd并输入以下命令启动&#xff1a; nginx.exe 或 start nginx3…

本地部署 EmotiVoice易魔声 多音色提示控制TTS

本地部署 EmotiVoice易魔声 多音色提示控制TTS EmotiVoice易魔声 介绍ChatGLM3 Github 地址部署 EmotiVoice准备模型文件准备预训练模型推理 EmotiVoice易魔声 介绍 EmotiVoice是一个强大的开源TTS引擎&#xff0c;支持中英文双语&#xff0c;包含2000多种不同的音色&#xff…

深信服测开实习

感觉会有人对这种类型的博客感兴趣&#xff0c;所以想了想还是发上来了。 今天早上十点收到了通知&#xff0c;这周五报道。 大致说工作时长&#xff1a; 周一到周五 一二四 早上九点到中午十二点。两小时午休。下午两点到晚上六点半。一小时晚休。七点半到晚上八点半下班。三…