栅格投影的意义:将分辨率单位为度的变为米
定义投影的意义:将某个栅格的坐标系变为想要的坐标系
批量栅格投影的代码如下:
# -*- 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!")
希望以上对大家有所帮助!
长路漫漫…
唯有坚持…