利用Arcpy批量对影像唯一值进行重分类(二分类)

在一些情况下,我们需要将唯一值影像在一定范围内的分为一类,在范围外的分为另外一类,Acrgis Pro也提供了唯一值的重分类的操作,但是对于数量很多的唯一值,一个个手动更改显得费时费力;利用Arcpy可以批量对影像的唯一值进行重分类,这里以二分类为例。


本文的分类要求:将影像值为0的分类为1,将值不为0的分类为-1。

Arcgis Pro唯一值重分类界面:

 显然手动更改费时费力!


利用Arcpy进行重分类代码如下:

import arcpy

# 设置环境参数
arcpy.env.workspace = "./workspace"
arcpy.env.overwriteOutput = True

# 输入栅格影像
# 这里使用的栅格影像值范围是-712~715
input_raster = "./data/raw.tif"

# 设置重分类规则
# 在重分类规则中,定义了像素值为0的为一类,值不为0的为另一类
reclass_field = "Value"
remap = arcpy.sa.RemapValue([[0, 1], [-712, -1, -1], [1, 715, -1]]) 

# 进行重分类
reclassified_raster = arcpy.sa.Reclassify(input_raster, reclass_field, remap)

# 保存重分类后的栅格影像
output_raster = "./reclassed.tif"
reclassified_raster.save(output_raster)

print("Reclassification complete.")

  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
使用arcpy进行批量采样的方法如下所示:首先,需要导入arcpy模块并设置工作空间和采样方法。然后,使用ListRasters函数获取要采样的栅格文件列表。接下来,使用ExtractByMask_sa函数将每个栅格文件按照指定的掩膜进行采样,并将结果保存到指定的输出路径中。最后,通过循环遍历每个栅格文件,输出采样完成的提示信息。以下是一个示例代码: ```python import arcpy arcpy.CheckOutExtension("spatial") arcpy.env.workspace = "E:\\ANUSSPLIN\\1000mPRE" # 栅格文件路径 rasters = arcpy.ListRasters("*", "tif") # 获取栅格文件列表 mask = "E:\\China_map\\长株潭垃圾站点信息\\长沙.shp" # 掩膜文件路径 for raster in rasters: out = "E:\\ANUSSPLIN\\PRE_output\\" + raster arcpy.gp.ExtractByMask_sa(raster, mask, out) print("ma_" + raster + " has done!") print("ok!!!") ``` 这段代码使用了arcpy的ExtractByMask_sa函数来进行采样操作,其中raster是要采样的栅格文件,mask是用来进行掩膜的矢量文件。通过循环遍历每个栅格文件,将其按照指定的掩膜进行采样,并将结果保存到指定的输出路径中。最后,输出采样完成的提示信息。\[1\] 希望对你有帮助! #### 引用[.reference_title] - *1* *3* [Python地理数据处理 十五:基于arcpy批量操作](https://blog.csdn.net/amyniez/article/details/127537354)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [Python遥感开发之arcpy批量采样](https://blog.csdn.net/qq_32306361/article/details/128088402)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^insertT0,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值