?投影

Python+ArcPy 對(duì)柵格數(shù)據(jù)進(jìn)行批量重投影_糯米丸子TAT的博客-CSDN博客

轉(zhuǎn)tif同時(shí)投影?

批量提取HDF文件中某子數(shù)據(jù)集中特定波段_疾風(fēng)勁草心兒_新浪博客?

需要注意的就是?arcpy.ProjectRaster_management這個(gè)函數(shù)的參數(shù)里面，第三個(gè)是目標(biāo)數(shù)據(jù)的坐標(biāo)系統(tǒng)，這個(gè)不能直接寫字符串“CGCS2000 3 Degree GK Zone 38”，而是一個(gè)“坐標(biāo)系統(tǒng)”對(duì)象，需要用arcpy.SpatialReference獲取。?

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

但是我直接用字符串也沒有報(bào)錯(cuò)?

# -*- coding:utf-8 -*- ##======================== ##批量柵格數(shù)據(jù)投影轉(zhuǎn)換import arcpy,os,os.path,initdef projectRaster(rootPath):try:##arcpy工作目錄root_path = rootPatharcpy.env.workspace = root_path##待處理文件所在目錄(相對(duì)于根目錄)input_path = "tif-source"output_path = "tif-target"##源坐標(biāo)系 "CGCS2000_3_Degree_GK_CM_123E" sourceSR = arcpy.SpatialReference("CGCS2000 3 Degree GK CM 123E")##目標(biāo)坐標(biāo)系(WGS 1984 Web Mercator Auxiliary Sphere)targetSR = arcpy.SpatialReference("WGS 1984 Web Mercator (auxiliary sphere)")##遍歷目錄，查找柵格數(shù)據(jù)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)

?

?

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。