批量栅格投影(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!")

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

长路漫漫…
唯有坚持…

上一篇:Verilog之参数传递实现条件编译


下一篇:[2021.10.21]<呆头熊的开发日记>怪物AI之有限状态机(1)