2. 程式人生 > 其它 >【IDL程式碼庫】如何用IDL獲取柵格範圍平均海拔

【IDL程式碼庫】如何用IDL獲取柵格範圍平均海拔

阿新 • • 發佈：2022-05-27

FLAASH已提供IDL介面（http://blog.sina.com.cn/s/blog_764b1e9d0102xxrk.html），但在FLAASH大氣校正中需要輸入研究區平均海拔，我們已經介紹過如何利用ENVI的統計功能獲取研究區高程（http://blog.sina.com.cn/s/blog_764b1e9d0101drqg.html），下面介紹如何一種使用IDL獲取高程方法，以下為功能原始碼及測試程式碼：

  ;+ ; :description: ; :獲取柵格範圍的平均海拔，單位為Km ; 失敗則返回-1,否則返回浮點型柵格範圍平均海拔 ; :author: Hanzt ; :2017-11-28
; - function Get_Raster_Msl,raster     ;編譯規則優化   compile_opt idl2   ;獲取envi物件   if isa(e,'envi') then e=envi(/current) else e=envi()     ;獲取ENVI自帶~900m解析度DEM資料絕對路徑（如果有更高解析度DEM也可替換）  root_Dir=e.root_dir  demFile=filepath('GMTED2010.jp2',$    root_Dir=root_Dir,subdirectory=['data'])     ;如果DEM不存在，則手動選DEM資料
   if ~file_test(demfile) then begin    title='Select a DEM file which covered the extent of the input raster.'    demFile=envi_pickfile(title=title)    if demFile eq '' then return,-1   endif     ;開啟DEM，並記錄錯誤資訊err0  demRaster=e.openraster(demFile,error=err0)     ;獲取柵格行列數  spatialref1=raster.spatialref
  ns=raster.ncolumns  nl=raster.nrows      ;獲取柵格四個角點投影座標，mapx,mapy  spatialref1.convertfiletomap,[0,ns,ns,0],[0,0,nl,nl],mapx,mapy  coordSys1=envicoordsys(coord_sys_str=spatialref1.coord_sys_str)     ;將mapx,mapy轉換到DEM柵格投影座標，用於DEM裁剪  spatialref2=demraster.spatialref  coordSys2=envicoordsys(coord_sys_str=spatialref2.coord_sys_str)  coordSys1.convertmaptomap,mapx,mapy, Map2X, Map2Y, coordSys2     ;獲取範圍最值  x_max=max(Map2X,min=x_min)  y_max=max(Map2Y,min=y_min)     ;獲取裁剪範圍  sub_rect=[x_min,y_min,x_max,y_max]     ;用柵格範裁剪DEM資料  rasterSub=envisubsetraster( $    demRaster,               $    spatialref=spatialref2,  $    sub_rect=sub_rect,       $    error=err1)     ;如果裁剪失敗，則返回-1   if rasterSub eq !null then begin    msg='No intersections between the input raster and the dem raster.'    !null=dialog_message(msg,/info)    return,-1   endif     ;裁剪後DEM統計  Statistics =ENVIRasterStatistics(rasterSub,error=err2)     ;如果失敗返回-1，否則返回平高程   if err0 ne '' || err1 ne '' || err2 ne '' then $    return,-1 else $    return,float(Statistics['MEAN']/1000) end     ;測試程式碼 pro test_Get_Raster_Msl   compile_opt idl2   if isa(e,'envi') then e=envi(/current) else e=envi()  file='C:\Program Files\Exelis\ENVI53\data\qb_boulder_msi'  raster=e.openraster(file)  msl=get_raster_msl(raster)     ;如果成功則控制檯列印   if msl.typecode eq 4 then  print,'Mean Sea Leavel:',msl   end

相關推薦

【IDL程式碼庫】如何用IDL獲取柵格範圍平均海拔

FLAASH已提供IDL介面（http://blog.sina.com.cn/s/blog_764b1e9d0102xxrk.html），但在FLAASH大氣校正中需要輸入研究區平均海拔，我們已經介紹過如何利用ENVI的統計功能獲取研究區高程（http://blog.sina.com.cn/s/b

【IDL程式碼庫】使用IDL視覺化引力波資料

新年新氣象，來點（看上去）高階的東西。 本文內容參考此文編寫（http://www.codingpy.com/article/gwpy-ligo-analyze-gravitational-waves-data/，作者：EarlGrey@程式設計派）。引力波科普知識，以及Python視覺化

【IDL程式碼庫】解決IDL生成多邊形向量在ArcMap中不閉合問題

偶然發現的問題。 比如寫出一個4個頂點的矩形向量，在ArcGIS中如果設定不填充，竟然只有三個邊。。。但是在ENVI中是顯示正常的。後來請教同事後發現，是由於不閉合的原因。

【IDL程式碼庫】我用IDL畫了個表

IDL的特長就是資料分析和視覺化，可以做很多有趣的事情。 比如本人在2014年就寫了個2048遊戲，當時確實很火，感興趣的可以去了解下。

【IDL程式碼庫】Colibri開源軟體包

Colibri介紹 Colibri是一個基於ENVI/IDL平臺實現的一個開源軟體包。裡邊包含了多種多樣的演算法和程式。

【IDL程式碼庫】IDL實現MODIS HDF檔案的輻射定標（氣溶膠系統）

一、介紹： MODIS 1KM資料的HDF包括250米和500米的兩個通道的資料，所有通道的解析度均為1000米。HDF的科學資料集包含的波段及波段型別見圖1。

【IDL程式碼庫】IDL呼叫6S.exe生成查詢表原始碼 （氣溶膠系統）

最新下載地址，包涵了原始碼和6s.exe程式：http://pan.baidu.com/s/1ppk1C ;+ ;:Description:

【IDL程式碼庫】IDL實現MODIS HDF檔案的幾何校正（氣溶膠系統）

一、介紹 MODIS HDF包含經緯度資料集，因此可以進行幾何校正。系統使用的方法是先讀取HDF中的經緯度資料集，然後建立51*51的GCP控制點（見圖1），通過對GCP點的投影轉換來對角度資料集和輻射定標後的科學資料集進行

【IDL程式碼庫】IDL實現MODIS資料雲檢測（氣溶膠系統）

一、介紹 系統採用的雲檢測演算法是多光譜綜合去雲法。該演算法從可見光反射率、紅外波段亮溫值以及亮溫差等方面綜合考慮，逐步建立一個雲檢測掩膜。對於預測的MODIS資料必須是經過了輻射定標的，熱輻射強度轉換成了

【IDL程式碼庫】利用IDLDrawWidget實現C#與IDL混合程式設計

在博文 http://blog.sina.com.cn/s/blog_764b1e9d01018v35.html 中較為詳細地介紹了IDL混合程式設計的基本手段和方法。本文以 IDLDrawWidget 為例，介紹如何進行C#與IDL之間的功能呼叫。

【IDL程式碼庫】環境衛星CCD資料氣溶膠反演工具原始碼分享

ENVI/IDL實現HJ衛星氣溶反演：http://blog.sina.com.cn/s/blog_764b1e9d01019hdw.html 這裡將環境衛星氣溶膠反演的三個工具和查詢表建立源程式分享給大家，不同於之前的modis氣溶膠反演程式，該程式做了查詢表插值，

【IDL程式碼庫】IDL中顯示分類影象的方法

IDL作為第四代視覺化語言，具有強大的視覺化能力。根據方法的不同，可以分為四類視覺化手段，分別為快速視覺化（IDL 8.0版本開始提供）、直接圖形法、物件圖形法和智慧化程式設計工具。

【IDL程式碼庫】大資料分塊寫入HDF5檔案示例程式碼

IDL提供了專門針對HDF5科學資料格式的讀寫函式庫。可以參考IDL幫助的Routines (alphabetical) > Routines: H > HDF5 Routines 這個章節檢視函式列表。

【IDL程式碼庫】陣列編輯器

利用此程式碼可以實現手動編輯一維和二維陣列。程式碼中使用了 Clipboard 物件，支援使用者可以從其他程式（如Excel或寫字板等）複製和貼上文字。

【IDL程式碼庫】IDL中實現bwareaopen功能（過濾斑塊）

Matlab中的bwareaopen函式功能可檢視百度百科： http://baike.baidu.com/link?url=zvIOHl35oJmAcA6pSIQh1INKyfppE8IkidfVEw_p0y2JtgFz1_WKfCB1B59ZU2_3GemYozbfzNUVAL2bLVaJNq

【IDL程式碼庫】一個完整的ENVI擴充套件工具原始碼

以TVDI VTCI擴充套件工具為例，為廣大遙感愛好者提供一個完整ENVI/IDL二次開發示例，包括演算法編寫、資料分塊處理、繪圖、IDL介面搭建、事件響應和自定義ENVITask等內容。

【IDL程式碼庫】基於6S模型生成MODIS氣溶膠反演查詢表

;+ ; :Author: Hanzt ; :Email: [email protected]，歡迎討論交流 ; 更新日誌 ;   2017-04-27 新增元資訊描述

【IDL程式碼庫】柵格轉向量內部介面

自ENVI 5.2版本開始提供了ENVIClassificationToShapefileTask，可以實現由分類影象轉換為shapefile向量檔案的功能。

【IDL程式碼庫】IDL實現氣泡排序演算法

本文使用IDL實現經典的氣泡排序演算法，雖然IDL內建了排序函式sort()，但是經典的氣泡排序演算法還是值得學習。氣泡排序類似於漢諾塔遊戲，關於氣泡排序wiki百科解釋如下：

【工具類庫】Guava 程式碼參考

https://github.com/google/guava/blob/master/guava/src/com/google/common/util/concurrent/RateLimiter.java