【文章內(nèi)容簡(jiǎn)介】 數(shù)據(jù)， r、 g、 b值被忽略 ok=query_png(filename ,s)： 獲取 PNG文件的屬性。 其中： ok=0， 不是合法的 PNG文件。 s為一個(gè)結(jié)構(gòu)，包含 PNG文件的屬性。若 =1， 為單通道數(shù)據(jù)，否則為多通道數(shù)據(jù)。 image=read_png(filename [,r,g,b])： 讀出 PNG文件的數(shù)據(jù) 其中： r,g,b對(duì) 單通道數(shù)據(jù)，讀出調(diào)色板 ,否則忽略。 常見用法： write_png ,filename ,tvrd(true=1) （ 作為 3通道數(shù)據(jù)寫入） ok=query_png(filename ,s) . . image=read_png(filename) tv image ,/true 直接圖形篇 － 圖像處理 圖像： 任何一個(gè)二維數(shù)組都可以視為一幅圖像。 8位圖像： 總表示為一個(gè)二維數(shù)組。 24位（ 真彩色 ）圖像： 總表示為一個(gè) 3維數(shù)組，其中一維為 3。 m*n*3： 隔波段掃描 （ bandinterleaved ， true=3） m*3*n： 隔行掃描（ rowinterleaved ， true=2） 3*m*n： 隔象素掃描（ pixelinterleaved， true=1） 在 8位顯示設(shè)備上，所有數(shù)據(jù)要轉(zhuǎn)換位字節(jié)類型；在 24位顯示設(shè)備上， 24位圖像的 r、 g、 b值必須轉(zhuǎn)換位字節(jié)類型。 直接圖形篇 － 圖像處理 1. 基本操作 顯示： tv和 tvscl。 兩個(gè)命令幾乎一樣，包括可以使用的關(guān)鍵字。都不刪除當(dāng)前窗口的內(nèi)容。 tvscl： 將圖像數(shù)據(jù)調(diào)整為與運(yùn)行時(shí)刻所有可用顏色數(shù)目相同的字節(jié)數(shù)據(jù)。通常用于 8位圖像 tv： 取圖像數(shù)據(jù)本身，作為字節(jié)數(shù)據(jù)顯示。如果圖像數(shù)據(jù)以整形和更多位數(shù)表示，則被截?cái)嘁赃m合字節(jié)類型。因此，圖形可能會(huì)顯示不正確。 在 24位設(shè)備上顯示圖像一般方法： device ,deposed=0 [loadct ,ct] tv／ tvscl ,image（ 8位圖像） tv ,image ,true=1|2|3（ 24位圖像 )） 直接圖形篇 － 圖像處理 2. 調(diào)整數(shù)據(jù)為字節(jié)類型并可以使用統(tǒng)一的顏色集 bytscl(image ,min= ,max= ,top=) 如： scaleImage=bytscl(image ,min=5 ,max=30) 3. 改變圖像尺寸 IDL提供了兩個(gè)改變圖像大小的命令： rebin和 congrid rebin： 新建的圖象的大小必須是原始尺寸的整數(shù)比例。缺省放大時(shí)采用雙線性插值，縮小時(shí)采用鄰近平均法。 sample關(guān)鍵字指定最近鄰近采樣法。 congrid： 新建的圖象的大小可以是原始尺寸的任意比例。缺省對(duì)三維以下數(shù)據(jù)采用最近鄰近采樣法，三維采用線性插值法。 interp關(guān)鍵字指定線性插值法 4. 在窗口中定位圖像 (1) tv ,image ,index： 根據(jù)圖像尺寸，從窗口左上角開始計(jì)算位置，逐行至右下角。 (2) tv ,image ,x ,y： 指定左下角開始計(jì)算的坐標(biāo)（ devic|data|normal）。 利用 ! !（ 象素值）計(jì)算歸一化坐標(biāo)以確定圖像位置和大小。 5. 從窗口中讀取圖像 8位顯示： image=tvrd() 24位顯示： image=tvrd(true=1) tvrd命令支持讀取指定區(qū)域的圖像 直接圖形篇 － 圖像處理 6. 基本圖像處理 上圖： 原始圖 中圖： hist_equal函數(shù) 下圖： adapt_hist_equal函數(shù) 直方圖均衡化： 觀察圖像中的象素值分布，往往會(huì)發(fā)現(xiàn)象素值分布趨于一個(gè)較狹窄的數(shù)值范圍內(nèi)。如果將象素值分散開，使象素值得每個(gè)子范圍都有與這些象素值大約相同的象素，則該圖像的信息內(nèi)容有可能增加。將象素分布道整個(gè)顏色范圍的過程稱為直方圖均衡化。 直接圖形篇 － 圖像處理 負(fù)片： 將原始圖像的顯示色板翻轉(zhuǎn)，象素的字節(jié)值不變所得到的圖像。 直接圖形篇 － 圖像處理 消除噪聲： 噪聲來自多方面，影響對(duì)圖像質(zhì)量。噪聲的一般表現(xiàn)形式是隨機(jī)的具有極端值的象素（黑白噪聲）。 median： 計(jì)算相鄰象素的中間值。這樣既可以消除極端值，又不會(huì)使大于鄰域的部分圖像邊緣或特征模糊。 平滑： 通過將每個(gè)象素值與它周圍相鄰象素值進(jìn)行平均來平滑圖像。稱為均值或核狀平滑。 smooth： 在給定的奇數(shù)寬度的范圍內(nèi)通過等加權(quán)值實(shí)現(xiàn)平滑。 convol： 使用給定的方形濾波核通過卷積實(shí)現(xiàn)平滑。 暈光蒙片： 將原始圖像減去平滑后的圖像?？梢远ㄎ粓D像上的邊緣或象素值突然變化的地方 直接圖形篇 － 圖像處理 上圖： smooth， w=5 中圖： smooth， w=3 下圖： convol， k=1,2,1 2,8,2 1,2,1 邊緣增強(qiáng)： 通過銳化或微分以增強(qiáng)邊緣。 roberts： sobel： convol： 使用給定的方形濾波核通過卷積實(shí)現(xiàn)。 直接圖形篇 － 圖像處理 右圖： 原始圖像 下圖： 自左至右依次為 roberts、 sobel、 convol 直接圖形篇 － 圖像處理 7. 頻域?yàn)V波 頻域?yàn)V波是圖像 和信號(hào)處理的常規(guī)手段?？梢杂糜谄交瑘D像、銳化圖像、降低圖像的模糊程度和恢復(fù)圖像。 基本步驟： （ 1）用快速傅立葉變換（ FFT） 將圖像從空間域轉(zhuǎn)換為頻率域 （ 2）將轉(zhuǎn)換后的圖像與一個(gè)頻率濾波器相乘 （ 3）將濾波后的圖像逆變換轉(zhuǎn)換為空間域 方法： filtered_img=fft( fft(img,1)*filter ,1) img可以是一維矢量，也可以是二維圖像 filter為濾波器，用于濾波圖像中某些特定頻率的一維矢量和二維數(shù)組 。 創(chuàng)建濾波器： 歐氏距離圖（頻率圖像）： r=dist(n [,m])， n和 m與實(shí)際圖像相同 Butterworth頻率濾波器： 低通濾波： filter=1/[1+c*(r/r0)2n] 高通濾波： filter=1/[1+c*(r0/r)2n] 其中： c= （ 即當(dāng) r=r0時(shí)，濾波幅度為 1/sqrt(2)） r0為濾波器截止頻率（實(shí)際中使用象素寬度） n為濾波器階數(shù)，通常 n=1 通常，低頻項(xiàng)代表圖像的一般形狀，高頻項(xiàng)對(duì)應(yīng)圖像的細(xì)節(jié) 直接圖形篇 － 圖像處理 東海海表面溫度圖 NOAA12 20:40 左圖： 原始圖 中圖： 低通（ w=10） 右圖： 高通（ w=10） 直接圖形篇 － 圖形窗口輸出 IDL支持多種設(shè)備輸出，并可以簡(jiǎn)單地在設(shè)備之間切換。 set_plot ,’option’： option=win|ps|printer， 大小寫不敏感。 !： 當(dāng)前設(shè)備名。 一般用法： cDevName=! set_plot ,‘printer’ help ,/device set_plot ,cDevName 常用關(guān)鍵字： close_document： 刷新輸出緩沖區(qū)后關(guān)閉圖形文檔，用于從打印機(jī)排出打印頁(yè)。 close_file： 刷新輸出緩沖區(qū)后關(guān)閉圖形輸出文件。 Filename： 當(dāng)輸出為文件時(shí)地文件名，默認(rèn) ‘ ’ landscape： 橫向輸出 portrait： 縱向輸出，缺省值 xoffset： 確定縱向模式下輸出窗口左下角（橫向模式下 ps和 printer算法不同） yoffset： 確定縱向模式下輸出窗口左下角（橫向模式下 ps和 printer算法不同） xsize： 確定輸出窗口寬度 ysize ： 確定輸出窗口長(zhǎng)度 inches： 以英寸為單位，缺省為厘米 直接圖形篇 － 圖形窗口輸出 ps 輸出的一般方法： img=tvrd(true=1) set_plot ,‘ps’ device, filename=‘ .ps’ ,color=1 device ,xsize= ,ysize= ,xoffset= ,yoffset= tvscl, img ,/true device, /close_file set_plot ,‘win’ printer輸出的一般方法： ok=dialog_printersetup() img=tvrd(true=1) ratio=float(!)/! （ printer設(shè)備不能自動(dòng)保持縱橫比） set_plot,39。printer39。 device,get_page_size=spagesize spagesize=spagesize/[!,!] device ,xoffset= ,yoffset= tvscl ,img ,/true ,/centimeters ,xsize= ,ysize= *ratio device ,/close_document set_plot ,’win’ IDL 高級(jí)培訓(xùn) 應(yīng)用程序構(gòu)造篇 編程篇 － 建立、管理應(yīng)用程序 IDL有兩種應(yīng)用程序管理方法： 主程序 和 工程 主程序： 建立與程序名同名的 .pro文件，編譯執(zhí)行。 pro test … end 在這種方式下，根據(jù) IDL的自動(dòng)編譯規(guī)則，所有被調(diào)用的子程序 或者 與主程序?qū)懭胪粋€(gè)文件并放置在主程序的上部； 或者 寫成與子程序同名的 .pro文件，存儲(chǔ)在可搜索的路徑下。（ 子程序指過程或函數(shù) ） 工程： 建立若干個(gè) .pro文件，每個(gè)文件不必與任何子程序同名。指定其中一個(gè)過程為主程序，該過程的名字作為工程執(zhí)行時(shí)首先調(diào)用的過程（即 run mand）。 在這種方式下，由于程序編寫未必遵守 IDL的自動(dòng)編譯規(guī)則，因此在運(yùn)行前必須編譯工程中所有的文件。 分發(fā)： 將源程序編譯為 .sav文件。 ? 主程序方式，需手工編譯、存儲(chǔ)為 .sav文件 .pile app_pro resolve_all save, /routines, app_pro ,app_sav ? 工程方式，設(shè)置工程屬性為形成 .sav文件，使用工程菜單項(xiàng)下的 build選項(xiàng)即可。 編程篇 － 應(yīng)用程序結(jié)構(gòu) XManger Widget Creation Routine (TLB) call xmanager Event Handler Calling Routine Event structure Only if TLB is destoryed 主程序 調(diào)用相關(guān)子程序 編程篇 － 主程序結(jié)構(gòu) 主程序的一般結(jié)構(gòu)： pro app_pro pile_opt idl2 create top layout base(TLB) create gui base on TLB (initial application) realize TLB collect user information call xmanager [,/NO_BLOCK] end initial application 該過程可以在主程序里進(jìn)行。但更好的辦法是編寫單獨(dú)的過程，用 TLB的notify_realize關(guān)鍵字指定過程名，當(dāng) realize TLB時(shí)執(zhí)行。 /NO_BLOCK