IDL实现NPP卫星数据地理校正

NPP（National Polar-orbiting Operational Environmental Satellite System Preparatory Project，国家极地轨道运行环境卫星系统预备工程卫）卫星于2011年10月28日，从加利福尼亚州范登堡空军基地发射成功,NPP卫星数据主要有22个波段，包括16个M波段数据，5个I波段，1个day and night波段，如下图

具体卫星的介绍，在我另外一篇博文中有介绍，大家可以参考。

NPP卫星数据自带经纬度信息，如图中的gdnbo,gimgo,gmodo文件就是包含经纬度信息的波段数据，分别对应day and night,M波段，I波段数据，具体介绍可以参考官网介绍，https://ladsweb.modaps.eosdis.nasa.gov/missions-and-measurements/level-0-level-1

因此可以利用自带的经纬度信息，建立GLT地理查找表文件，然后将GLT文件应用于NPP卫星数据，从而实现对NPP 卫星数据的几何校正，主要是利用ENVI_GLT_DOIT和ENVI_GEOREF_FROM_GLT_DOIT这两个函数

下面就以day and night 波段数据为例进行说明：

1;获取校正参数文件 root_path为存放文件的目录
gdnbo=file_search(root_path,'*gdnbo.h5')
2;读取不同类型文件的经纬度数据
lat_gdnbo=e.OpenRaster(gdnbo[0],DATASET_NAME='All_Data/VIIRS-DNB-GEO_All/Latitude')
lon_gdnbo=e.OpenRaster(gdnbo[0],DATASET_NAME='/All_Data/VIIRS-DNB-GEO_All/Longitude')
3；输入待校正的文件input_file为输入的待校正的文件
dnb_raster_ref=e.OpenRaster(input_file,DATASET_NAME='/All_Data/VIIRS-DNB-SDR_All/Radiance')
dbn_id=ENVIRASTERTOFID(svm_raster_ref)
envi_file_query,dbn_id,dims=dims_dbn
4；构建GLT地理文件，r_fid为返回的GLT文件的fid
  lat_id=ENVIRASTERTOFID(lat_gdnbo)
  lon_id=ENVIRASTERTOFID(lon_gdnbo)
  envi_file_query,lat_id,nb=nb_lat
  y_pos=lindgen(nb_lat)
  envi_file_query,lon_id,nb=nb_lon
  x_pos=lindgen(nb_lon)
  input_prj=ENVI_PROJ_CREATE(/geographic,datum='WGS-84')
  output_prj=ENVI_PROJ_CREATE(/UTM,zone=48,datum='North America 1927')
  ENVI_DOIT,'ENVI_GLT_DOIT',dims=dims,I_PROJ=input_prj,O_PROJ=output_prj,$
    OUT_NAME=out_name,r_fid=glt_id,$
    ROTATION=0,X_FID=lon_id, X_POS=x_pos[0], Y_FID=lat_id, Y_POS=y_pos[0]
5；应用GLT文件实现地理校正 out_name为输出的文件名
ENVI_FILE_QUERY,glt_fid,dims=dims
svi_id=ENVIRASTERTOFID(dnb_raster_ref)
ENVI_DOIT,'ENVI_GEOREF_FROM_GLT_DOIT',FID=svi_id,GLT_DIMS=dims,GLT_FID=glt_fid,OUT_NAME=out_name,pos=0,r_fid=out_fid