?投影

Python+ArcPy 對柵格數據進行批量重投影_糯米丸子TAT的博客-CSDN博客

轉tif同時投影?

批量提取HDF文件中某子數據集中特定波段_疾風勁草心兒_新浪博客?

需要注意的就是?arcpy.ProjectRaster_management這個函數的參數里面，第三個是目標數據的坐標系統，這個不能直接寫字符串“CGCS2000 3 Degree GK Zone 38”，而是一個“坐標系統”對象，需要用arcpy.SpatialReference獲取。?

基于Python批量將HDF文件轉為GeoTiff格式并進行拼接、投影轉換和矢量裁剪 - 走看看?

但是我直接用字符串也沒有報錯?

# -*- coding:utf-8 -*- ##======================== ##批量柵格數據投影轉換import arcpy,os,os.path,initdef projectRaster(rootPath):try:##arcpy工作目錄root_path = rootPatharcpy.env.workspace = root_path##待處理文件所在目錄(相對于根目錄)input_path = "tif-source"output_path = "tif-target"##源坐標系 "CGCS2000_3_Degree_GK_CM_123E" sourceSR = arcpy.SpatialReference("CGCS2000 3 Degree GK CM 123E")##目標坐標系(WGS 1984 Web Mercator Auxiliary Sphere)targetSR = arcpy.SpatialReference("WGS 1984 Web Mercator (auxiliary sphere)")##遍歷目錄，查找柵格數據files = os.listdir(root_path+os.sep+input_path)for f in files:if os.path.splitext(f)[1].upper() == ".TIF":fileName = os.path.splitext(f)[0] + ".tif"in_dataset = input_path + os.sep + fileNameout_dataset = output_path + os.sep + fileNameprint "begin project "+in_dataset+" from: " +sourceSR.name+" to: "+targetSR.namearcpy.ProjectRaster_management(in_dataset, out_dataset, targetSR, "NEAREST",\"#", "#", "#",sourceSR)print "project success!"except arcpy.ExecuteError:print "Project Raster example failed."print arcpy.GetMessages()################################################ if __name__ == '__main__':#指定處理文件根目錄root_path = r"E:\gis"projectRaster(root_path)

?

?

總結

以上是生活随笔為你收集整理的arcpy.ProjectRaster_management的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。