基于idl的遙感影像波段運(yùn)算

1、基于IDL的遙感影像波段運(yùn)算1.IDL語言特征在IDL中，矩陣是按照列的方式處理的，即矩陣是以n列、m行的方式表示的，這一點(diǎn)與Ftran語言一樣，但與C語言（按m行、n列的維度來標(biāo)記）表示矩陣的方式不同。只有熟悉IDL的數(shù)據(jù)存儲(chǔ)、處理的特點(diǎn)，才能方便地進(jìn)行分析。以300列、200行的遙感數(shù)據(jù)b1為例，其數(shù)據(jù)的記錄特征如下：IDL圖形坐標(biāo)b1DN值坐標(biāo)圖1.IDL中像元坐標(biāo)與亮度矩陣的對(duì)應(yīng)關(guān)系（1）數(shù)據(jù)按列存儲(chǔ)；最先存儲(chǔ)的是第一列的數(shù)據(jù)（

2、圖象第一列的像元亮度值），接著是第二列，直至最后一列。（show:printdata[0]）（2）采用順序顯示數(shù)據(jù)的方式。將第一行的亮度值按順序從左到右以行的方式顯示或存儲(chǔ)，不管b1矩陣有多少行、多少列，總是先顯示、存儲(chǔ)完第一行，然后是下一行。由于IDL的默認(rèn)顯示寬度是80列，所以在輸出窗口中要用幾行才能顯示完矩陣的第一行亮度值。然后另起一行顯示亮度矩陣的下一行值，直到最后一行。這一點(diǎn)對(duì)于圖象顯示和分析相當(dāng)重要。復(fù)習(xí)：復(fù)習(xí)：IDL矩陣的

3、顯示方式，矩陣的顯示方式，data（show:printdata[0]）（3）默認(rèn)情況下，圖象數(shù)組的第一行和第一個(gè)元素顯示在屏幕的左下角左下角。后面的行將從下往上顯示。在圖形窗口顯示圖象時(shí)，圖象的每一個(gè)元素在屏幕中顯示為一個(gè)像素[1]。在IDL中圖象是從下往上畫的。即把b1亮度矩陣的第一行值顯示在圖象最下端的一行像元中，把b1亮度矩陣的最后一行值顯示在圖象最上端的一行像元中。這是圖象顯示最重要的特征。區(qū)別兩幅圖象的顯示效果：（區(qū)別兩幅圖

4、象的顯示效果：（tvscldata）與（）與（tvscldata300，0，der=1）（4）矩陣轉(zhuǎn)置與圖象的關(guān)系像元坐標(biāo)亮度矩陣1121……30019294……901222……30029692……91…………………………12002200……30020098101……90從左到有一一對(duì)應(yīng)，亮度矩陣中的第二行值，與圖象中的第二行像元從左到有一一對(duì)應(yīng)，依次類推，直到最后一行。在IDL中打開b1的亮度矩陣并且顯示在IDL繪圖環(huán)境中的時(shí)候，要正

5、確地顯示圖象（見圖2），必須要確保亮度矩陣與圖象坐標(biāo)系統(tǒng)的正確關(guān)系，如亮度矩陣中的左下三角（或右上三角）一定要嚴(yán)格對(duì)應(yīng)圖象像元矩陣的左下三角（或右上三角），否則就會(huì)發(fā)生圖象翻轉(zhuǎn)、倒置（見圖3），其實(shí)質(zhì)就是將波段亮度矩陣進(jìn)行了轉(zhuǎn)置或行列倒置。練習(xí)：程序調(diào)用。調(diào)用b4的II文件的代碼如下：file=filepath(b1_matrix.txtsubdir=examplesdata)openrlunfileget_lundata=lonarr

6、(300200)事先知道圖象大小readflundataprintsize(datan_elements)；data就是b4圖象的亮度矩陣free_lunlun3.波段運(yùn)算（BMath）BMathTM功能允許你處理導(dǎo)致單個(gè)波段輸出的復(fù)雜表達(dá)式。這些數(shù)學(xué)表達(dá)式也可以應(yīng)用于一個(gè)多波段文件中的所有波段，providing“FileMath”。關(guān)于使用波段運(yùn)算的更多信息，請(qǐng)參閱ENVIProgrammer’sGuide第29頁的“BMathBa

7、sics”。1可利用的波段運(yùn)算功能（AvailableBMathFunctions）BMath功能為用戶提供一個(gè)靈活的圖像處理工具，其中許多功能是無法在任何其它的圖像處理系統(tǒng)中獲得的。該功能的能力與IDL語言的能力直接相關(guān)?？捎玫暮瘮?shù)包括但不僅限于表42中列出的數(shù)學(xué)表達(dá)式。表42:一些可用的波段運(yùn)算函數(shù)。SeriesScalar數(shù)學(xué)數(shù)學(xué)三角函數(shù)三角函數(shù)其它波段運(yùn)算選項(xiàng)其它波段運(yùn)算選項(xiàng)加（）正弦（sin（x））關(guān)系運(yùn)算符（EQ、NE、LE