批量转换遥感影像：WGS84至CGCS2000 3度带坐标

pro step3
  compile_opt IDL2
  envi,/restore_base_save_files
  envi_batch_init
  ;打开要转换投影的影像
  image_dir='D:\test\input\' ;根据文件存放的目录进行相应修改
  image_files=file_search(image_dir,'*.img',count=numfiles)  ;根据相应的文件格式修改过滤条件
      ;****************添加功能代码******************
       for i=0,numfiles-1 do begin
        e = ENVI(/headless)
        image_file=image_files[i]
        envi_open_file,image_file,r_fid=fid
        envi_file_query,fid,ns=ns,nl=nl,nb=nb,dims=dims
         raster1 = e.OpenRaster(image_file)
         NCOLUMNS=raster1.NCOLUMNS                ;影像的行列数
         NROWS=raster1.NROWS
         spatialRef = raster1.SPATIALREF          ;空间参考
         File1X = round(NCOLUMNS/2)               ;影像中心坐标
         File1Y = round(NROWS/2)
         spatialRef.ConvertFileToMap, File1X, File1Y, MapX, MapY
         spatialRef.ConvertMapToLonLat, MapX, MapY, Lon, Lat        ;地图坐标转经纬度
         a=round(Lon/3)                           ;影像中心经度进行3度取整
         Lonint=a*3
         print,Lonint
         print,file_baseName(image_file,'.img')+'.img'
;*******************CGCS2000-57E**************************         
         if Lonint EQ 57 then begin
          ;定义CGCS2000-57E投影
           name='CGCS_2000_3_Degree_GK_CM_57E'
           datum='D_China_2000'

           params=[6378137.0,6356752.3,0.000000, 57.000000,500000.0, 0.0, 1.000000]
           o_proj=envi_proj_create(type=3,name=name,params=params,datum=datum)
           ;投影转换
           pos=lindgen(nb)
           o_ps=[250,250]
           grid=[25,25]
           outfiledir='D:\test\true'
           out_name = outfiledir +'\' +file_baseName(image_file,'.img')+'.img'
           envi_convert_file_map_projection,fid=fid,r_fid=fid_out,o_proj=o_proj,dims=dims,pos=pos,out_name=out_name,background=0,o_pixel_size=o_ps,grid=grid,warp_method=2,resampling=0
         endif
 ;*******************CGCS2000-60E**************************          
          if Lonint EQ 60 then begin
           ;定义CGCS2000-60E投影
           name='CGCS_2000_3_Degree_GK_CM_60E'
           datum='D_China_2000'
           params=[6378137.0,6356752.3,0.000000, 60.000000,500000.0, 0.0, 1.000000]
           o_proj=envi_proj_create(type=3,name=name,params=params,datum=datum)
           ;投影转换
           pos=lindgen(nb)
           o_ps=[250,250]
           grid=[25,25]
           outfiledir='D:\test\true'
           out_name = outfiledir +'\' +file_baseName(image_file,'.img')+'.img'
           envi_convert_file_map_projection,fid=fid,r_fid=fid_out,o_proj=o_proj,dims=dims,pos=pos,out_name=out_name,background=0,o_pixel_size=o_ps,grid=grid,warp_method=2,resampling=0
         endif
  ;*******************CGCS2000-63E**************************  
     if Lonint EQ 63 then begin
    ;定义CGCS2000-63E投影
    name='CGCS_2000_3_Degree_GK_CM_63E'
    datum='D_China_2000'