IDL實現NPP衛星資料地理校正

2021-09-26 16:28:19 字數 1934 閱讀 1262

npp(national polar-orbiting operational environmental satellite system preparatory project,國家極地軌道執行環境衛星系統預備工程衛)衛星於2023年10月28日,從加利福尼亞州范登堡空軍基地發射成功,npp衛星資料主要有22個波段,包括16個m波段資料,5個i波段,1個day and night波段,如下圖

npp衛星資料自帶經緯度資訊,如圖中的gdnbo,gimgo,gmodo檔案就是包含經緯度資訊的波段資料,分別對應day and night,m波段,i波段資料,具體介紹可以參考官網介紹,

因此可以利用自帶的經緯度資訊,建立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

白鴿 使用IDL批量預處理國產高分光學衛星資料

乙個idl程式,用於批量預處理高分資料。衛星 載荷 wfvpms gf1?gf1b c d gf2 gf6?版本發布日期 v2.2 2021 01 08 v2.1 2020 11 14 v2.0 2020 10 31 v1.3 alpha 2020 06 16 v1.2 alpha 2020 05 ...

IDL 實現求算 DEM 坡度坡向

關於坡度坡向的定義,請google之。idl 原始碼 pro aspect slope,dem,aspect aspect,slope slope,pixelsize pixelsize ref to doi 10.1016 j.cageo.2003.07.005 aspect 1 slope 1 ...

IDL實現向量 shp 裁剪柵格TASK(一

隨著envi idl版本的更新,idl對向量和柵格資料的處理也變得越來越簡單化。其提供了很多方便的介面,使得使用者呼叫和學習練習便捷成為了可能。最近接觸idl,發現好多網上的 都是延後的,新的介面 理解和編寫起來都比較方便,尤其是在做大量資料研究和應用時,使用批處理的方式顯得尤其重要。新的介面還在摸...