Python地理数据处理 24：基于arcpy批量操作（六）

1. 同年tif数据作差

对于同一年成像的两景遥感影像分别进行做差处理。如将2001q.tif文件减去2001m.tif文件，将2002q.tif文件减去2002m.tif文件，以此类推。

import arcpy
import os
arcpy.CheckOutExtension('Spatial')

# 设置工作空间
workspace = r"D:\dataset"  # 文件夹路径
arcpy.env.workspace = workspace

# 获取文件夹中的所有tif文件
tif_files = arcpy.ListFiles("*.tif")

# 存储每个年份的eos和sos文件路径
year_files = {
    
    }

# 遍历tif文件，根据文件名将eos和sos文件分组
for tif_file in tif_files:
    basename = os.path.splitext(tif_file)[0]
    year = basename[:4]
    file_type = basename[4:]  # 文件类型，如eos、sos
    
    if year in year_files:
        year_files[year][file_type] = tif_file
    else:
        year_files[year] = {
    
    file_type: tif_file}

# 对每个年份进行相减操作
for year in year_files:
    if "m" in year_files[year] and "q" in year_files[year]:
        m_raster = year_files[year]["m"]
        q_raster = year_files[year]["q"]
        
        # 生成新的栅格图像名称
        output_raster = year + "m_q.tif"
        
        # 相减操作
        out_raster = arcpy.sa.Minus(m_raster, q_raster)
        
        # 保存新的栅格图像
        output_path = "D:\\dataset\\s\\" + output_raster
        out_raster.save(output_path)  
        print year

2. Calculate Statistics

有时候tif影像无法建立属性表，或者无法分类，这种情况可能是由于图像未进行统计计算而导致的，以下为arcgis中该功能的界面，但是无法进行批量处理：

批量处理代码：

import arcpy
import os

# set workspace
workspace = r"D:\dataset"
arcpy.env.workspace = workspace

# get all .tif image
tif_files = arcpy.ListFiles("*.tif")

# run the Calculate Statistics in arcgis
for tif_file in tif_files:
    tif_path = os.path.join(workspace, tif_file)
    arcpy.CalculateStatistics_management(tif_path)
    print "to {} run the Calculate Statistics function".format(tif_file)