批量栅格投影（arcpy）

栅格投影的意义：将分辨率单位为度的变为米
定义投影的意义：将某个栅格的坐标系变为想要的坐标系

批量栅格投影的代码如下：

# -*- coding: UTF-8 -*-
#需要自行更改的如下：
#inws：输入路径
#outws：输出路径
#Coordinate_System：目标投影坐标系
import glob
import os
import arcpy

# 输入路径  应该注意，中文路径，会导致读不出文件
inws = r"------"

# 参考文件路径    使用栅格数据集（从其导入方形像元大小）的路径
Coordinate_System="PROJCS['UTM_Zone_46N',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['false_easting',500000.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',93.0],PARAMETER['scale_factor',0.9996],PARAMETER['latitude_of_origin',0.0],UNIT['Meter',1.0]]"
# 输出路径
outws = r"----"

# 利用glob包，将inws下的所有tif文件读存放到rasters中
rasters = glob.glob(os.path.join(inws, "*.tif"))

# 循环rasters中的所有影像，进行按掩模提取操作
for raster in rasters:
    print str(raster)
    nameT = os.path.basename(raster).split(".")[0] + ".tif"	
    outname = os.path.join(outws, nameT)  # 合并输出文件名+输出路径

    arcpy.ProjectRaster_management(raster, outname, Coordinate_System, "NEAREST","#","#","#","#")

print("OK!")

希望以上对大家有所帮助！

长路漫漫…
唯有坚持…

标签：rasters,os,arcpy,路径,投影,inws,栅格,PARAMETER
来源： https://blog.csdn.net/qq_41828818/article/details/120977814