如何使用idl做GLT重投影（调用envi接口）

下面是调用envi接口来实现的GLT，效果比较好，就是稍微来说速度有点慢。（课上学的）

function get_hdf4_data,hd_name,sds_namehd_id=hdf_sd_start(hd_name)sd_id=hdf_sd_nametoindex(hd_id,sds_name)sds_id=hdf_sd_select(hd_id,sd_id)hdf_sd_getdata,sds_id,datahdf_sd_endaccess,sds_idhdf_sd_end,hd_idreturn,data
end
pro build_GLT_envicompile_opt idl2scale_factor=0.001dir='E:\modis_swath\'out_dir='E:\out\'out_dir1='E:\out_img\'file_list=file_search(dir,'*.hdf')print,file_listfile_n=n_elements(file_list)for file_i=0,file_n-1 do beginname=file_list[file_i]need_lon_max=118need_lon_min=115need_lat_max=42need_lat_min=39AOD_data=get_hdf4_data(name,'Image_Optical_Depth_Land_And_Ocean')lon=get_hdf4_data(name,'Longitude')lat=get_hdf4_data(name,'Latitude')lat_min=min(lat)lat_max=max(lat)lon_min=min(lon)lon_max=max(lon)if ((lat_min le need_lat_min) and (lat_max ge need_lat_max)and(lon_min le need_lon_min)and(lon_max ge need_lon_max)) then beginlon_name=out_dir+file_basename(name,'.hdf')+'lon'+'.tiff'lat_name=out_dir+file_basename(name,'.hdf')+'lat'+'.tiff'aod_name=out_dir+file_basename(name,'.hdf')+'aod'+'.tiff'write_tiff,lon_name,lon,/float;open函数不能打开hdf函数，故要转换为tiff文件write_tiff,lat_name,lat,/floatwrite_tiff,aod_name,AOD_data,/floatenvi,/restore_base_save_filesenvi_batch_initENVI_OPEN_FILE,lon_name,R_FID=fid_xENVI_OPEN_FILE,lat_name,R_FID=fid_yENVI_OPEN_FILE,aod_name,R_FID=fid_aodi_proj=envi_proj_create(/geographic)o_proj=envi_proj_create(/geographic)envi_file_query,fid_aod,sname=snameout_glt_img=out_dir+sname+'_glt_'+'.img'out_glt_hdr=out_dir+sname+'_glt_'+'.hdr'pixel_size=0.03rotation=0.0x_pos=0y_pos=0pos=0envi_glt_doit,i_proj=i_proj,$o_proj=o_proj,out_name=out_glt_img,$pixel_size=pixel_size,r_fid=glt_fid,rotation=rotation,$x_fid=fid_x,y_fid=fid_y,x_pos=x_pos,$y_pos=y_pos;glt_fid为输出参数out_name=out_dir1+sname+'_georeproj.img'out_name_hdr=out_dir1+sname+'_georeproj.hdr'envi_doit,'envi_georef_from_glt_doit',fid=fid_aod,$glt_fid=glt_fid,out_name=out_name,pos=pos,$r_fid=r_fid;其中r_fid为输出参数，可以不管。endif else begincontinueendelse   endfor
end

如何使用idl做GLT重投影（调用envi接口）相关推荐

Python实现VIIRS气溶胶产品重投影-类GLT实现
Python实现VIIRS气溶胶产品重投影-类GLT实现前言代码说明用到的外部python包代码实现注意事项后记前言最近用IDL写了个VIIRS气溶胶产品重投影的代码,之后心想着用py ...
ENVI_IDL:批量重投影ModisSwath产品(调用二次开发接口)+解析
目录 1. 课堂内容 1. 获取Modis Swath数据(这里只获取Lat.Lon.Aod(气溶胶厚度)三个数据集以及aod数据集的两个属性),并对aod数据进行简单的处理 2. 调用二次开发接口以 ...
IDL实现MODIS Grid（正弦投影）产品的重投影及拼接处理
IDL实现MODIS Grid(正弦投影)产品的重投影及拼接处理前言 map_proj_image函数使用关键单个文件的重投影示例多个文件的重投影+拼接后记前言关于MODIS正弦投影产品的 ...
IDL对HDF5图像做GLT校正
最近需要对HY-1C的HDF5图像做几何校正,并简单的显示,但是查找了大量的资料,并没有找到相应的代码可供参考,只找到了一位朋友的博客,是对HDF影像做处理,在此感谢,地址在这儿http://blog ...
IDL实现遥感数据的快速重投影（几何校正）- 以MODIS Swath产品为例
IDL实现遥感数据的快速重投影(几何校正)- 以MODIS Swath产品为例网上关于遥感数据的重投影资料其实已经很多了,比如基于GDAL的实现.基于ENVI GLT的实现.基于一些专用处理工具的实 ...
ENVI_IDL:批量重投影Modis Swath产品并指定范围输出为Geotiff格式+解析
目录 1. 实验内容 2. 知识储备 3. ENVI实操对对应DL代码部分 4. 编程 5. 题外话 5.1 n_element()与一些函数自带的count参数返回的区别 5.2 发现一个难以发现 ...
Kalibr源码学习(一): 重投影误差
Kalibr源码学习(一): 重投影误差给自己挖一个大坑, 从标定结果来学习Kalibr的标定源码, 这里基本以KB模型为例, 也就是标定时, kalibr的模型设定为 --model pinhol ...
在 VSLAM 的后端优化中的重投影误差的雅可比计算详细推导
对于相机位姿的变换可以通过旋转矩阵或者四元数进行表示,对于旋转矩阵的定义满足: R{∣R∣=1RRT=IR \begin{cases} |R| = 1 \\ RR^T = I\\ \end{cases ...
利用MODIS TOOL批量重投影MODIS
MODIS TOOL介绍 MODIS TOOL全称是MODIS Reprojection Tools,里面主要是重投影和镶嵌拼接和重采样,如果不嫌麻烦的话你可以用GUI一个一个操作,但是效率非常的低, ...

如何使用idl做GLT重投影（调用envi接口）

如何使用idl做GLT重投影（调用envi接口）相关推荐

最新文章

热门文章