ENVI中批量转TIFF，Arcgis批量地理校正
一、数据介绍
本次使用的数据为时空连续臭氧产品，从武汉大学“地学感知数据质量改善与融合应用研究室”网站（网址：http://sendimage.whu.edu.cn/send-resource-download/）下载，它是“利用发展的遥感信息重建理论与方法，基于Aura卫星上OMI传感器的臭氧总量三级产品OMTO3e，通过对该产品的缺失信息进行重建与修复，生成了全球（2004-2014）最高分辨率的时空连续臭氧制图产品。”——来自官网介绍。产品经解压后如下所示：
将上面的数据从ENVI里导出为TIFF，并在Arcgis中打开后如下所示，与左下方标准坐标数据（绿色坐标）在地理空间上不匹配，无法进行下一步处理：
由于数据量比较大，一个一个地georeferencing比较麻烦，因此需要采用批量处理的方法。
二、总体处理思路
首先是要在ENVI中，将图1所示的原始文件批量转换为TIFF，再在arcgis中批量地理校正。
三、ENVI中批量转TIFF
新建TXT文件，并重命名为“ENVIRaster_Study.pro”，用记事本或者notepad++或其他软件打开后，输入以下代码，其中第12行的输出文件路径根据实际情况调整。
PRO ENVIRaster_Study
e = ENVI() ; Launch the application
inputfilefolder=dialog_pickfile(/directory,title='Choose the input directory!');选择输入文件夹
inputfilefolderfiles=file_search(inputfilefolder,'*filled',count=num);读取文件夹下的文件
for j=0,num-1 do begin
inputfile=inputfilefolderfiles[j]
raster1 = e.OpenRaster(inputfile) ; 读入文件
inputfilesplit = strsplit(inputfile,'\',/extract)
;这里inputfile类似于"F:\Zhanglan\OMIO3\2004m1001-TFFSRC_filled"这样的字符串
;将字符串以"\"符号打断,<br /> ;以便取最后一个"\"后的字符串"2004m1001-TFFSRC_filled"经变换后做输出文件名
outputfilename=inputfilesplit[N_ELEMENTS(inputfilesplit)-1];取打断后最后一个"\"后的字符串做输出文件名
filepath_output = 'F:\Zhanglan\Postgraduate\Zhoumen\Tongji\O3Project\FN'+outputfilename+'.tif' ; 输出文件路径
raster1.Export, filepath_output, 'TIFF'; 输出为tiff格式
Print,'finished' ;打印finished，提示完成
endfor
end保存文件后，双击打开（前提是安装了IDL），这时IDL会自动打开pro文件，点击“运行”，选择要处理的数据所在的文件夹，即可。运行结果如下图所示：
四、ArcGIS中批量进行地理校正
经第三步转换后的TIFF文件，在ArcGIS中打开后，如第二张图所示，丢失了地理坐标信息，无法与正确的地理坐标匹配，所以，还需要进行地理校正。
传统的方法是使用georeferencing手动校正，这无法满足大量数据均需要校正的要求，所以，使用Python脚本做批量处理。
大致方法是在ArcMap中新建toolbox，右键新建的toolbox/add/script，添加写好的脚本。然后，再右键新建的toolbox/new/model，新建一个model，并把添加的脚本拖进model里，点击Model/Run Entire Model即可。具体自行百度。
那脚本如下，其中source_pnt为待校正图像的四至点，target_pnt为要校正到的四至点（目标四至点or正确的四至点）：
try:
import arcpy
arcpy.env.workspace = r"F:/Zhanglan/Postgraduate/Zhoumen/Tongji/O3Project"
files=arcpy.ListRasters()
##Warp a TIFF raster dataset with control points
##Define source control points
source_pnt = "'0 720';'1440 720';'0 0';'1440 0'"
##Define target control points
target_pnt = outpnts = "'-180 90';'180 90';'-180 -90';'180 -90'"
for file in files:
outputfilename="F:/Zhanglan/Postgraduate/Zhoumen/Tongji/O3Warp/"+file[:-4]+"WARP.tif"
arcpy.Warp_management(file, source_pnt, target_pnt, outputfilename, "POLYORDER1","NEAREST")
except:
print "Warp example failed."
print arcpy.GetMessages()等待结果即可。