遙感課程設(shè)計

1、課 程 設(shè) 計 報 告課程名稱： 海洋遙感技術(shù) 實習(xí)名稱： 從TRMM/TMI衛(wèi)星數(shù)據(jù)反演全球11GHZ風(fēng)速 院 （系）： 海洋科學(xué)與工程學(xué)院 專業(yè)班級： 海洋技術(shù) 姓 名： 學(xué) 號： 指導(dǎo)教師： 2014年 7 月 5 日一、課程設(shè)計時間2014年6月30日至7月4日二、課程設(shè)計地點天津科技大學(xué)9-513海洋信息技術(shù)實驗室三、課程設(shè)計目的 理論與實驗課的綜合運用，提高課堂與實踐相結(jié)合的分析能力1、了解海洋衛(wèi)星反演*的物理原理以及兩種類型（熱紅外與微波輻射計）的衛(wèi)星傳感器；2、熟悉TRMM/TMI；3、掌握TMI反演*的算法；4、進一步熟悉IDL/matlab編程語言；5、了解TMI的數(shù)據(jù)結(jié)構(gòu)

2、以及數(shù)據(jù)的輸入與輸出；6、學(xué)會數(shù)據(jù)的圖形顯示，創(chuàng)建圖形輸出；7、掌握數(shù)據(jù)的圖象顯示；8、熟悉*數(shù)據(jù)的讀取，顯示，輸出，并注意colorbar的顯示；9、了解*對全球氣候的影響以及重要意義，著重說明*全球分布的特點。四、課程設(shè)計主要儀器設(shè)備，軟件1、硬件準(zhǔn)備：PC機；2、操作系統(tǒng)：Linux系統(tǒng)或Windows 2k以上系統(tǒng)；3、軟件開發(fā)工具：IDL/matlab語言4、2012-2013年日、周、月均的TRMM/TMI二進制衛(wèi)星數(shù)據(jù)五、微波輻射計反演的原理介紹吸收系數(shù)與電離層離子或大氣層分子的吸收有關(guān)，我們使用（z）描述在這個路徑位置z處的吸收系數(shù)。沿著電磁波輻射的傳播路徑對吸收的輻射通量積分

3、，可獲得輻射的在一個總衰減量。假設(shè)傳感器在一個吸收路徑的遠端測量一個具有表面輻亮度的光源。該光源輻射的電磁波在位置z處的輻亮度是L（z），經(jīng)介質(zhì)吸收而衰減，輻亮度L（z）在路徑微分元上為 （6-39）上式只考慮吸收引起的衰減而未考慮散射引起的衰減，是布格-郎伯-比爾定律在散射可以忽略條件下的特例。根據(jù)基爾霍夫定律，給出吸收氣體發(fā)射的輻亮度是：介質(zhì)吸收率a和吸收系數(shù)的關(guān)系因此，輻射傳輸方程給出的輻亮度變化的凈變化：衛(wèi)星輻射計探測的輻亮度有三個組成部分。第一部分是達到衛(wèi)星輻射計的海表面發(fā)射的輻亮度；第二部分是衛(wèi)星輻射計探測到的大氣自發(fā)輻射的輻亮度；第三部分是大氣向下發(fā)射的、達到海表面后又經(jīng)海表面反

4、射的輻亮度。我們得到輻射計探測到的輻亮度：微波頻率低于300GHZ，滿足瑞麗金斯定律的條件。因此，在微波輻射計對應(yīng)的輻射傳輸工程中，我們可以使用亮溫代替輻亮度。根據(jù)電磁輻射理論，輻射傳輸方程的微分形式是由于微波波長大于大氣層中各種粒子的尺度，散射衰減對于微波影響一般可以忽略，而僅僅需要考慮吸收的衰減。得到方程的解將瑞麗金斯定律帶入，我們得到考慮更多地輻射源，適合于微波輻射計的輻射傳輸方程可表示為平靜海面的微波發(fā)射率：根據(jù)適合兩介質(zhì)界面處的基爾霍夫定律，海面發(fā)射率e與菲涅爾發(fā)射率關(guān)系是：粗糙海面的微波發(fā)射率：輻射計接收到的海面輻射用輻亮度L表示。根據(jù)瑞麗金斯定律，輻亮度在微波波段與溫度成線性關(guān)系

5、。在不考慮大氣校正時，輻射計探測到的海面亮溫與海表面溫度有下列關(guān)系：式中，代表粗糙衛(wèi)星海面發(fā)射率，它是衛(wèi)星觀測角、微波頻率f、輻射計極化狀態(tài)、海表面溫度T、海表面鹽度、海面以上10m高出的風(fēng)速和風(fēng)向的函數(shù)。例如，基于小斜率近似六、TRMM/TMI衛(wèi)星及傳感器介紹1、寫出三顆微波輻射計及所在衛(wèi)星，并簡單介紹；衛(wèi)星TIROS-N 傳感器 MSU 上世紀(jì)70 年代末至90 年代末，以TIROS-N 發(fā)射為標(biāo)志，衛(wèi)星性能相對穩(wěn)定，形成了由大型計算機構(gòu)成的地面資料接收處理系統(tǒng)，真正進入了實際業(yè)務(wù)應(yīng)用階段。1978 年發(fā)射的第一次用于大氣探測業(yè)務(wù)的TIROS-N 衛(wèi)星，裝載有微波探測儀（MSU）。MSU

6、為4 通道微波輻射計（50.3, 53.7, 55.0, 57.9 GHz），采用雙旋轉(zhuǎn)鏡面天線，刈幅為2 300 km，空間分辨率為110 km，用于從較低對流層到較低平流層的大氣反演和大氣垂直溫度探測。從此，MSU 成為NOAA 系列衛(wèi)星的主要載荷，其通道主要在50 GHz 至60 GHz 之間，用以探測大氣溫度廓線。 衛(wèi)星 NOAA-15 傳感器 AMSU 1998 年，NOAA-15 衛(wèi)星上裝載新研制的20 通道先進微波探測載荷AMSU，取代了原較低分辨率微波探測裝置。自此，美國空間微波探測進入高速發(fā)展階段。AMSU 采用圓周掃描方式，觀測頻率從23.8 至183.0GHz ，主要用于

7、大氣垂直溫度和濕度的探測。 衛(wèi)星SZ-4 傳感器 多頻段微波輻射計（MFMR）2002 年12 月, 神舟4 號（SZ-4）飛船發(fā)射成功, 其主載荷為多模態(tài)微波遙感系統(tǒng), 包括了多頻段微波輻射計（MFMR）、雷達高度計和雷達散射計，實現(xiàn)了我國星載微波遙感器的突破。MFMR 微波輻射計主要用于探測土壤水分、降水、大氣水汽含量、積雪分布、海面溫度等，它是以功能性驗證和應(yīng)用試驗為主要目的的空間對地觀測微波遙感系統(tǒng)。MFMR 采用5 頻段（頻點包括6. 6, 13. 9, 19. 3, 23. 8 和37.0 GHz）、雙極化（其中6.6，19.3 和37.0 GHz 為水平或垂直極化）、全功率輻射計

8、體制，以及集中供電與統(tǒng)一數(shù)據(jù)管理設(shè)計方案，探測數(shù)據(jù)通過數(shù)據(jù)管理單元傳至衛(wèi)星公用分系統(tǒng)，并下傳至地面應(yīng)用系統(tǒng)。其中，13. 9 GHz 輻射計與雷達散射計公用一個接收機前端和天饋系統(tǒng)，采用圓錐掃描成像；而6. 6、19. 3、23. 8 和37.0 GHz 采取多頻共體饋源、公用偏置拋物面技術(shù)以及主反射器40°角斜入射和不掃描方案，頻段共用孔徑形成多頻段同軸輻射, 偏置拋物面將輻射流信號聚焦到饋源口徑上，再通過分波網(wǎng)絡(luò)將不同頻率、不同極化的信號分開, 最后饋入各自響應(yīng)頻段的接收機。2、 TRMM TRMM衛(wèi)星是從1997年11月27日發(fā)射成功的。它是由美國NASA和日本NASDA共同研

9、制的試驗衛(wèi)星。自發(fā)射成功以來，它為氣象工作者提供了大量熱帶海洋降水、云中液態(tài)水的含量、潛熱釋放等氣象數(shù)據(jù)。TRMM衛(wèi)星在2001年8月6日以前高度約為350km，軌道范圍為35S-35S。2001年8月以后，為減少功耗將軌道高度提升到距地面403km。 TRMM衛(wèi)星共搭載5種遙感儀器，分別為：可見光和紅外掃描儀VIRS、TRMM微波圖像儀TMI、降水雷達PR、閃電圖像儀LIS及云和地球輻射能量系統(tǒng)CERES，其中VIRS、TMI和PR為TRMM衛(wèi)星的基本降水測量儀器。74um光學(xué)探測儀，它被用來觀測閃電分布形式及其變化。CERES是工作波長分別為0.3um-0.5um和8.0-12.0um的可

10、見光和紅外探測器，它主要是用來探測地表及大氣發(fā)射和反射的輻射能量及其分布形式。3、 TMI 微波成像儀TMI（TMRR microwave imager）在五個波段設(shè)有10個通道，這10個通道是10.7vh、19.4vh、21.3vh、37.0vh和89.0vh等。七、TMI反演11GHZ風(fēng)速算法 表1.海洋產(chǎn)品的預(yù)計檢索精度地球物理參數(shù)均方根精度海表溫度Ts0.5 K近地表風(fēng)速W1.0 m/s垂直整合（即，柱狀）水蒸氣V1.0 mm垂直整合的云中液態(tài)水的L0.02 mm 表2. 過去和未來的衛(wèi)星輻射計系統(tǒng)的比較輻射計頻率/偏振角度SeaSat SMMR 6.6VH10.7VH 18.0VH

11、21.0VH 37.0VH 49°Nimbus-7 SMMR6.6VH10.7VH 18.0VH 21.0VH 37.0VH 51°SSM/I19.3VH22.2V 37.0VH 85.5VH 53°TRMM TMI10.7VH19.3VH21.0VH 37.0VH 85.5VH53°PM AMSR6.9VH10.7VH18.7VH23.8VH36.5VH89.0VH55°校準(zhǔn)輻射傳輸模型（RTM）使用SSM / I和SMMR觀測衛(wèi)星AMSR算法開發(fā)的第一個版本，通過RTM利用模擬TB的計算 版本1推出前的AMSR算法 利用TRMM結(jié)果指定10

12、.7 GH參數(shù) 版本2推出前的AMSR算法 最后的發(fā)射前的測試 在EOS PM最后發(fā)射前的AMSR算法 在6.9 GHz的指定的海表面發(fā)射率 使用實際的在軌AMSR觀測 版本1推出后的AMSR算法 圖1.海洋算法的開發(fā)步驟 環(huán)境場景42195個無線電探空儀5個云模型SST的隨機變化從0到30 C風(fēng)速隨機變化從0到20米/秒風(fēng)方向的隨機變化從0到360 真相：TS，W，V，L 完整的輻射傳輸模型 模擬 AMSR TB's 高斯噪聲 保留的數(shù)據(jù)集衍生物系數(shù)的多次線性回歸算法 算法的系數(shù) 運行算法 檢索值為TS，W，V，L Evalulate 算法性能 性能和Cross Talk統(tǒng)計 圖5.

13、求解和測試線性回歸算法如今伴隨著高精度的衛(wèi)星微波輻射計的出現(xiàn)，使得我們長期獲取一系列的地球物理參數(shù)變的可能，這對于研究全球水循環(huán)和地球輻射收支有著重要意義。環(huán)顧全球海洋，這些輻射計擁有空氣溫度的分布和三個階段的大氣水（蒸汽，液體，和冰）的能力。此外，如近地表風(fēng)速、海面溫度、海冰類型和濃度等地表參數(shù)也可以被糾正。各種各樣的水文和輻射過程中可以對這些測量結(jié)果進行研究，包括空氣和海水、空氣和冰的相互作用（即潛熱和感熱通量、淡水通量和表面應(yīng)力）和云對輻射通量的影響。微波輻射計是一個真正用于研究地球的獨特、有價值的工具。 在本文中，我們著重于獲取從特殊傳感器微波成像儀(SSM/I) Hollinger

14、et al. 1987上所觀測到的世界各地海洋中的地球物理參數(shù)的問題。SSM/I飛行于美國國防的極軌衛(wèi)星計劃（DMSP）的兩個運算極地軌道平臺間。7 SSM/I's系列中的第一個發(fā)射于1987年6月，并且SSM/ I系列可能將持續(xù)到2000年，到那時SSM/ I系列將會被一個稱為SSM/ IS的聯(lián)合成像儀/測深儀所替換。因此，我們將有機會獲取一個可以使用SSM/ IS觀察來進一步擴展的13年的全球時間系列的地球物理產(chǎn)品。 SSM / I的四個運行頻率為19.35、22.235、37和85.5 GHz。通過這些渠道，讓我們在海洋上獲取三個重要的地球物理參數(shù)變得可能:近地面風(fēng)速W（米/秒）

15、，柱狀水蒸氣V（毫米），以及柱狀云中液態(tài)水L（mm）。降雨同樣也可以被推算出，但是在本文中，我們的調(diào)查只考慮無雨的海洋場景。在無雨的情況下，由SSM / I測定的海洋亮度溫度（TB）和W、V和 L之間有一種相對簡單獨特的關(guān)系。雨的出現(xiàn)增加了相當(dāng)大的復(fù)雜性，就目前來看，我們要避免這個問題。 潛在的，由于TB和（W，V，L）之間的獨特關(guān)系，W、V和L可以被糾正成一個很高的精度。這個關(guān)系是由一個非降雨氣氛包圍在由海面下方的輻射傳遞方程（RTE）給出。它已經(jīng)表明， RTE可以通過一個相對簡單的閉形表達式（即，沒有積分），這就是所謂的TB模型函數(shù)Wentz, 1983。（W，V，L）的糾正是通過改變這些

16、參數(shù)，直到TB模型函數(shù)與SSM / I觀測值相匹配來實現(xiàn)的。因此，（W，V，L）的準(zhǔn)確度依賴于TB模型的準(zhǔn)確性。為了獲得盡可能高的糾正精度，TB模型必須包括RTE的所有相關(guān)參數(shù)的影響并且TB模型必須精確地校準(zhǔn)。對TB模型的功能和它的后續(xù)精密校準(zhǔn)的完全參數(shù)化是本文的主題。 本文用SSM / I傳感器的描述作為開始。然后，我們將討論TB模型的參數(shù)?？偣灿?個主要參數(shù)（W，V，L）和4次要參數(shù)：海面溫度TS（K），有效的大氣溫度TE（K），水蒸氣柱的有效大氣壓力P（MB），以及風(fēng)向f。TB對這些次要參數(shù)的依賴性相對于主要參數(shù)來說要較弱。然而，這些二次依賴關(guān)系仍然顯著并且必須要考慮到。第三節(jié)討論了從氣

17、候?qū)W得到的用于指定TS和TE和P的統(tǒng)計關(guān)系。 風(fēng)向太容易變以至于我們必須通過指定特定的氣候，因此我們將其包括為除了W，V，L之外的第四個糾正參數(shù)。將風(fēng)向帶入到TB模型中分為視線的風(fēng)分量WLS，其中風(fēng)矢量是沿著SSM / I觀測方向的分量。納入WLS作為第四個糾正有兩個好處，首先，對于其他參數(shù)（特別是W）它降低了糾正的錯誤，其次它還提供了海洋上空風(fēng)向的新信息。 第四節(jié)討論的是糾正算法。對于每個SSM / I象素，當(dāng)他們被取代成TB模型后，該算法所產(chǎn)生的（W，V，L，WLS）的值等于SSM / I觀測在19 GHz v-pol, 22 GHz v-pol, 37 GHz v-pol, 和 37 G

18、Hz h-pol的值。TB模型是準(zhǔn)線性的四個參數(shù)的計算，因此糾正中涉及解決四個未知數(shù)四個方程。第五節(jié)則完整的闡述了TB模型。 第6、7節(jié)描述了非常大的浮標(biāo)和用來校準(zhǔn)模型和糾正算法的無線電探空儀的數(shù)據(jù)集。全球浮標(biāo)和無線電探空儀觀測是在1987年至1990年4年間收集的?，F(xiàn)有66浮標(biāo)站點和55探空站點。這些現(xiàn)場測量并置與SSM / I立交橋。搭配方法一共有37,650 SSM / I浮標(biāo)天橋和35,108探空站點立交橋的產(chǎn)生。 第8節(jié)描述了TB模型中大氣系數(shù)的推導(dǎo)。理論上的亮度溫度由無線電探空儀觀測使用RTE的完整的積分公式計算得出。在TB模型中這些理論的TB用于計算大氣系數(shù)。由于TB模型的部分大

19、氣近似導(dǎo)致的糾正錯誤是由RTE TB的數(shù)據(jù)作為輸入輸入到糾正算法中模擬從而確定的。第8節(jié)還討論了糾正中空氣壓力變化的影響。 TB模型（及此后的糾正算法）的校準(zhǔn)經(jīng)由逆建模技術(shù)將會在第9節(jié)被討論。在TB模型中校準(zhǔn)是通過改變系數(shù)從而使得W和V匹配浮標(biāo)和無線電探空儀觀測。此外，L糾正的直方圖需要滿足一些統(tǒng)計學(xué)的約束。總之，W，V和L的糾正要求滿足19個統(tǒng)計條件。這種類型的校準(zhǔn)被稱為逆模型，這種類型的校準(zhǔn)被稱為逆模型，因為該模型的系數(shù)的推導(dǎo)是基于模型的逆（即糾正算法）的輸出。 第10節(jié)討論了由于風(fēng)向變化而導(dǎo)致的WLS和風(fēng)速錯誤的糾正。我們總結(jié)了一個完整的誤差分析。一個錯誤的預(yù)算表顯示了TB建模誤差、輻射

20、計的噪聲、和SSM/ I原位時空抽樣誤差，這些都有助于在糾正中觀察到的總RMS變化?；谶@個誤差分析，我們可以估計TB模型的準(zhǔn)確性和地球物理反演。 通過有競爭力的同行評審，NASA已經(jīng)選擇了本文描述用于生產(chǎn)的SMMR，SSM / I探路者數(shù)據(jù)集的算法。該數(shù)據(jù)集將是一個20年的地球物理參數(shù)的時間序列，這將成為一個廣泛分布的研究領(lǐng)域。在本文的分析是基于的1987-1990第一SSM / I的觀察，即飛赴強奎DMSP F08飛船Hollinger等。 1987。航天器F08的軌道是近圓形太陽同步和近極地，具有98.8°的傾角。海拔高度為860±25公里，軌道周期為102分鐘。在高

21、度的變化是由于軌道的偏心率和地球的扁率。本地時間的F08升序赤道橫穿是上午6:15am。 該SSM / I傳感器由7個獨立的總功率輻射計共享一個共同的嗽叭。這7個輻射計取雙偏振測量在19.35，37.0和85.5GHz，而僅僅一個垂直極化測量在22.235 GHz。該SSM / I使用尺寸61到66厘米的偏移拋物面反射鏡收集微波輻射。反射聚焦輻射包括瓦楞紙，寬帶，7端口饋源。反射器和喇叭天線旋轉(zhuǎn)在最低點軸使用一個單位。旋轉(zhuǎn)周期為1.9秒。將冷空間反射器和熱參考負載連接到旋轉(zhuǎn)軸，并且不旋轉(zhuǎn)。每一次掃描旋轉(zhuǎn)饋源觀察固定的冷反射和熱負荷。以這種方式，校準(zhǔn)化觀測采取每次掃描。地球觀測在102.4

22、76;旋轉(zhuǎn)段期間拍攝。在102.4°圓弧的中心位于飛船子軌跡，并對應(yīng)于地球的一個1400公里的大片表面。 1400公里的表面和98.8°軌道傾角為地球提供兩到三天的全覆蓋，除了以南北兩極2.4°為中心的兩個小圓形的部分。對地球反射觀看的最低點角為45°，從而導(dǎo)致地球入射角 在53.4°±0.25°。在頻率較低的通道（19，22，和37 千兆赫）進行采樣，使得像素間距為25公里，以一個12.5公里像素間距對85 GHz的通道進行采樣。 SSM / I測量來自于地球景物發(fā)射的輻射強度。瑞利 - 金斯近似Eisberg，1961表

23、達了這種強度在溫度方面，被稱為天線溫度TA。對于SSM / I，天線溫度近似交配通過 T AIP= GIPVT BIV + G IPHT BIH+ G IPO TBC (1) 其中下標(biāo)I和P表示頻率和極化，TBIV和TBIH是V-pol和H-POL地球的亮度溫度。 TBC是宇宙背景輻射相當(dāng)于 2.7 K。在G因素是占了天線的天線方向圖系數(shù)溢出 以上和交叉極化泄漏。近似（1）的推導(dǎo)和值 為G系數(shù)由華斯1991給出。 天線的溫度的平均值，以一個共同的空間分辨率。在地球上19，22，和37 GHz的渠道的SSM / I足跡的半功率波束寬度為56，44和32公里，分別與這些腳印的中心重合。因

24、此，一個37GHz的觀察只能看到由19 GHz的采樣通道面積的33。為了獲得準(zhǔn)確的檢索，有必要使所有信道看到相同的海洋區(qū)域。這是通過平均該22和37 GHz的觀測下，以較低分辨率the19 GHz的信道的使用下面的公式： （2）其中TAIJ是在原始分辨率天線溫度（無論是在22或37 GHz）和兩個下標(biāo)表示現(xiàn)在的沿軌掃描數(shù)量和整個軌道的單元格位置。一組權(quán)重wIJ（1組為22千兆赫，而另一組為37千兆赫）的發(fā)現(xiàn)，使得平均TA的有效天線圖案19-GHz的天線圖案相匹配。該權(quán)重依賴于跨細胞軌道的位置，因為相對于掃描帶中心該細胞的相對位置是不同的掃描帶邊緣。相鄰的掃描和相鄰小區(qū)之間的距離約為25公里，我

25、們發(fā)現(xiàn)，這是足夠做平均時，只包括直接相鄰的小區(qū)。此后，我們畫在TA上的橫線，可以理解為所有的觀察是在一個共同的空間分辨率。八、IDL/matlab讀取TRMM/TMI全球11GHZ風(fēng)速程序 1、IDL介紹IDL（Interactive Data Language）交互式數(shù)據(jù)語言是進行二維及多維數(shù)據(jù)可視化分析及應(yīng)用開發(fā)的理想軟件工具。作為面向矩陣、語法簡單的第四代可視化語言，IDL致力于科學(xué)數(shù)據(jù)的可視化和分析。IDL在海洋領(lǐng)域具有廣泛的應(yīng)用，在各種海洋數(shù)據(jù)的可視化及與時間、空間相關(guān)的分析，對風(fēng)、流、浪、溫度以及水、粒子流等數(shù)據(jù)的生物學(xué)和化學(xué)參數(shù)的時間序列，不同格式的數(shù)據(jù)的讀取，不同的機器和操作系

26、統(tǒng)的兼容都有很強的優(yōu)勢。Topex/Poseidon 數(shù)據(jù)時間上的動態(tài)模擬2、用到主要的IDL/matlab過程及函數(shù)命令：Openr, readu, congrid, contour, xrange, xstyle, xticks, xtickv, xtickname, tvscl, tv, tvrd(), write_png.3、編寫的主程序pro colorbar,low,top,count,title nc=241 cb=bindgen(nc) cb = cb#replicate(1b, 20) cb=congrid(cb,nc,1*20) ; Drew the bundury of

27、the colobar cb(0,*) = 255 cb(*,0) = 255 cb(nc-1,*) = 255 cb(*,19) = 255 ; Plot the tick marks and labels on the colorbar xp = (!d.x_size - nc)/2 yp = 50 plot,low,top,low,top,position=xp,yp,xp+nc-1,yp,/device,/nodata,$ xticks=n,xtickv=tick,xrange=low,top,xminor=5,xstyle=1,xtitle=title, $ ystyle=4,xti

28、cklen=0.002,yticklen=0.002,/noerase,color=255 ; Display the colorbar tv,cb,xp,yp,/deviceENDpro testfilename='C:UsersgaojianfengiDesktoptmi_201301v4'binary data in filebinarydata= bytarr(1440,320,6);output products (lon,lat,asc/dsc)sst =fltarr(1440,320)w11 =fltarr(1440,320)w37 =fltarr(1440,32

29、0)vapor=fltarr(1440,320)cloud=fltarr(1440,320)rain =fltarr(1440,320)openr,lun,filename,/get_lunrecord=''readu,lun,binarydatafree_lun,lunprint,size(binarydata)w11=binarydata(*,*,1)w11=float(w11)print,w11(1,*)I=where(w11 gt 250,count)if(count gt 0) then w11(I)=0W=w11*0.2print,min(W),max(W)device,decompose=0W=congrid(W,800,300)window,1,xsize=1000,ysize=500tvlct,100,100,100erase,color=0loadct,39I=where(wind gt