遙感IDL二次開發(大氣校正)

2021-10-18 16:41:15 字數 2118 閱讀 8011

1.程式功能:

用envi開啟tm資料,將資料傳到idl工作空間,在idl環境下進行大氣校正,並將結果回傳給envi。

2.程式源**:

function atmocorr32,test ;*

****

****

****

*** 讀入資料 ***

****

****

****

****

**;size_image=

size

(test)

samples=size_image[1]

;tm列數

lines=size_image[2]

;tm行數

bands=size_image[3]

;tm波段數

;讀取 tm 的 gain 和 offset 值到 data 中

fn_calib=

'輻射定標引數txt檔案儲存路徑'

; openr,lun,fn_calib,

/get_lun

data=

fltarr(2

,6) data1=

fltarr(3

,6) readf,lun,data0

readf,lun,data1

free_lun,lun

gain=data0[0,

0:5]

;增益值

offset=data0[1,

0:5]

;偏移值

;用定標係數將原始dn值轉換為大氣層頂太陽輻亮度l

l=fltarr

(size_image[1]

,size_image[2]

,size_image[3]

);定義輻亮度

for i=

0,bands-

1do begin

l[*,

*, i]

=gain[i]

*test[*,

*,i]

+offset[i]

; endfor

;輸出輻射定標後檔案

envi_write_envi_file,out_name=

'大氣校正引數txt檔案路徑'

,l

;讀取大氣校正引數

a=data1[0,

*] b=data1[1,

*] c=data1[2,

*];設定乙個陣列準備存放大氣校正的結果

e=fltarr

(size_image[1]

,size_image[2]

,size_image[3]

);大氣校正

for i=

0,bands-

1do begin

y=a[i]

*l[*,*

,i]-offset[i]

; e[*,

*,i]

=y/(

1+c[i]

*y);

endfor

;輸出大氣校正檔案

envi_write_envi_file,out_name=

'輸出大氣校正結果txt檔案路徑'

,e

return

,e end

**注意:**首先要開啟 envi+idl聯合程式, 然後用envi載入影像,並將需要使用的影像匯入到 idl 中,最後在idl命令視窗中輸入函式:如 e=atmocorr32(test)

3.結果

postek二次開發 POSTEK二次開發介面

例項簡介 例項截圖 核心 using system using system.collections.generic using system.componentmodel using system.data using system.drawing using system.text using ...

mcms二次開發 YzmCMS二次開發手冊

摘要 yzmcms是採用mvc設計模式開發,基於模組和操作的方式進行訪問,採用單一入口模式進行專案部署和訪問,無論訪問任何乙個模組或者功能,只有一.yzmcms是採用mvc設計模式開發,基於模組和操作的方式進行訪問,採用單一入口模式進行專案部署和訪問,無論訪問任何乙個模組或者功能,只有乙個統一的入口...

mcms二次開發 YzmCMS二次開發手冊

yzmcms是採用mvc設計模式開發,基於模組和操作的方式進行訪問,採用單一入口模式進行專案部署和訪問,無論訪問任何乙個模組或者功能,只有乙個統一的入口。一 基本目錄結構.根目錄 cache 快取目錄 必須可寫入 common 全站公共目錄 必須可寫入 uploads 預設上傳目錄 必須可寫入 yz...