資源簡(jiǎn)介
只要改變代碼的數(shù)據(jù)文件路徑，無(wú)論上千上萬(wàn)個(gè)數(shù)據(jù)文件，一次運(yùn)行即可機(jī)器處理。
Python批量讀取modis的hdf文件轉(zhuǎn)asc文件,修改文件夾即可，自動(dòng)創(chuàng)建子文件夾生產(chǎn)結(jié)果到子文件夾。
代碼片段和文件信息
#?Import?system?modules
import?arcpytimeos
from?arcpy.sa?import?*
arcpy.CheckOutExtension（“spatial“）
arcpy.env.overwriteOutput?=?True
year=2005?#modify
arcpy.env.workspace?=?r‘E:\Data\RSdata\MODIS\MOD09CMG_Reflectance_5600m_daily_Global\%d‘%year?#你的數(shù)據(jù)路徑，自己修改
file=“output%d“%year??#change?your?name
if?not?arcpy.Exists（arcpy.env.workspace+os.sep+file）:
????????arcpy.CreateFolder_management（arcpy.env.workspacefile）
????????
rasters?=?arcpy.ListRasters（“*“?“HDF“）
print?‘共有%s個(gè)文件處理‘%len（rasters）
for?iinput?in?enumerate（rasters[40:]）:
????t1=time.time（）
????outASCII?=?arcpy.env.workspace+os.sep+file+‘\\‘+input.replace（‘hdf‘‘a(chǎn)sc‘）
????#print?‘Execute?%d?RasterToASCII...‘%i
????arcpy.RasterToASCII_conversion（input?outASCII）???
????t2=time.time（）
????print?‘第%d個(gè)用時(shí)%s評(píng)論
共有 條評(píng)論