鉴于前面两篇中存在背景值干扰问题,这里采用CopyRaster方法赋值栅格数据并去除背景值
代码如下:
# -*- coding: utf-8 -*-
import os,arcpy,datetime
#读取和存放文件夹的位置
inputPath1 = r'H:\input1'#镶嵌数据文件夹1
inputPath2 =r'H:\input2'#镶嵌数据文件夹2
path=r'H:\outputMocaic'#输出结果存放文件夹
# start calculate time
startTime = datetime.datetime.now()
print startTime
# set the intermediate data folder #这段根据需要使用
intermediateDataPath = path+"\\"+"IntermediateData" #如果有中间文件就存放在这里
# set result data folder
resultDataPath = path+"\\"+"Result"
# determine if the folder exists
if os.path.exists(intermediateDataPath):
print "IntermediateData floder exists"
else:
# create a intermediate data floder
arcpy.CreateFolder_management(path,"IntermediateData")
if os.path.exists(resultDataPath):
print("Result floder exists")
else:
本文介绍使用Python的arcpy模块进行遥感图像处理,通过CopyRaster方法实现栅格数据的背景值去除。在实际操作中发现,尽管结果显示正常,但在图层列表中仍检测到-3000的背景值,暗示可能存在异常值问题。
订阅专栏 解锁全文
1459





