張 彤，叢林虎，顧佼佼，叢 晉，王昌金

(1.海軍駐洛陽407廠軍事代表室，河南 洛陽 471039; 2.海軍航空工程學(xué)院，山東 煙臺(tái) 264001)

0 引言

在最新一代的戰(zhàn)斗機(jī)航空顯示中已基本采用純光柵顯示方式，但目前正在服役的多數(shù)機(jī)型仍使用光柵筆劃疊加畫面顯示器。文獻(xiàn)［1］中利用光柵筆劃信號(hào)的時(shí)序關(guān)系，將光柵筆劃疊加信號(hào)分解成為可由通用儀器產(chǎn)生與測(cè)量的分量信號(hào)。從而摒棄了傳統(tǒng)的專用設(shè)備檢測(cè)思路。這種方法在采集筆劃畫面各路信號(hào)時(shí)，由于儀器自身及外界干擾等因素的影響，從設(shè)備采集得到的數(shù)據(jù)一般都包含大量的噪聲，從而在測(cè)控計(jì)算機(jī)上復(fù)現(xiàn)純筆劃畫面及光柵筆劃疊加畫面時(shí)，由于采集得到的筆劃信號(hào)數(shù)據(jù)未進(jìn)行過預(yù)處理，致使顯示的畫面不夠清晰美觀并影響了識(shí)別、測(cè)試等后續(xù)工作的開展。因此，為了更加真實(shí)清晰地復(fù)現(xiàn)出筆劃畫面，有必要對(duì)采集后的筆劃信號(hào)進(jìn)行處理，即通過一系列預(yù)處理方法降低采集后的筆劃信號(hào)中的噪聲水平，從而為后續(xù)工作的開展提供更為可靠的信號(hào)數(shù)據(jù)集。

1 光柵筆劃疊加顯示方式

?？兆鲬?zhàn)中，戰(zhàn)斗機(jī)飛行員主要通過平視顯示器和多功能顯示器來獲得飛機(jī)與作戰(zhàn)信息。這些信息是通過光柵和筆劃這兩種掃描機(jī)制完全不同的信號(hào)組合產(chǎn)生并疊加顯示的［2］。光柵筆劃顯示器由工作模式和光柵筆劃切換信號(hào)(S/R)控制，既可工作于純筆劃方式、純光柵方式，又可以工作于光柵筆劃疊加顯示方式。當(dāng)顯示器工作于光柵筆劃疊加顯示方式時(shí)，顯示器在光柵筆劃切換信號(hào)的控制下，在光柵視頻信號(hào)的場(chǎng)消隱期間顯示筆劃信號(hào)，其他時(shí)間顯示光柵信號(hào)。這樣用戶在屏幕上看到的就是筆劃視頻畫面疊加在光柵視頻畫面上的疊加畫面［3］。

光柵視頻信號(hào)、偏轉(zhuǎn)信號(hào)、顏色信號(hào)、光柵筆劃切換信號(hào)等要保持同步關(guān)系，邏輯關(guān)系如圖1所示。

圖1 筆劃信號(hào)Fig.1 Stroke signal

對(duì)單色多功能顯示器，則只有一個(gè)顏色信號(hào)或輝亮信號(hào)。在圖1中，光柵視頻場(chǎng)同步信號(hào)的下降沿表示視頻信號(hào)要進(jìn)入場(chǎng)消隱，即將顯示筆劃畫面。光柵筆劃切換信號(hào)的上升沿標(biāo)示了一幅筆劃畫面的開始，下降沿標(biāo)示了筆劃畫面的結(jié)束。在光柵筆劃切換信號(hào)有效期間內(nèi)，筆劃的位置由偏轉(zhuǎn)X和偏轉(zhuǎn)Y信號(hào)確定，筆劃顏色由筆劃顏色R、筆劃顏色G、筆劃顏色B這3路信號(hào)共同決定。筆劃畫面繪制完成后繼續(xù)顯示光柵畫面。

2 航空顯示信號(hào)測(cè)試

2.1 測(cè)試基本原理

針對(duì)光柵筆劃疊加畫面顯示信號(hào)及其邏輯關(guān)系的特點(diǎn)，同樣考慮ATS的通用化原則，選擇多通道并行數(shù)據(jù)采集模塊、光柵視頻采集模塊、通信模塊和相應(yīng)電源模塊來完成光柵筆劃信號(hào)采集、復(fù)現(xiàn)、顯示和測(cè)量的任務(wù)，如圖2所示。

圖2 測(cè)試系統(tǒng)總體設(shè)計(jì)Fig.2 Design of the test system

控制計(jì)算機(jī)控制光柵視頻采集模塊采集光柵視頻信號(hào)，控制多通道并行數(shù)據(jù)采集模塊采集其他分量信號(hào)，然后根據(jù)光柵筆劃信號(hào)各分量信號(hào)之間的時(shí)序關(guān)系復(fù)現(xiàn)、顯示光柵畫面、筆劃畫面和光柵筆劃疊加畫面。通信模塊則用于模擬有通信功能的多功能顯示器，完成與顯示驅(qū)動(dòng)設(shè)備之間進(jìn)行的設(shè)置工作方式、查詢周邊鍵狀態(tài)等通信功能［4］。

2.2 信號(hào)采集和顯示

對(duì)光柵視頻信號(hào)，直接采用視頻信號(hào)采集模塊進(jìn)行采集，將視頻畫面保留在內(nèi)存或保存在圖像文件中。

對(duì)光柵視頻場(chǎng)同步、光柵筆劃切換、偏轉(zhuǎn)X、偏轉(zhuǎn)Y、顏色等信號(hào)，采用高速并行數(shù)據(jù)采集模塊進(jìn)行同步采集，每次采樣時(shí)間大于兩個(gè)視頻信號(hào)場(chǎng)周期，保證能采集到一幅完整的筆劃畫面。在采樣數(shù)據(jù)序列中，光柵筆劃切換信號(hào)的上升沿標(biāo)示了一幅筆劃畫面的開始，下降沿標(biāo)示了筆劃畫面的結(jié)束。在筆劃畫面期間，偏轉(zhuǎn)X和偏轉(zhuǎn)Y信號(hào)確定了筆劃的位置，顏色信號(hào)則表示筆劃的顏色。根據(jù)各種分量信號(hào)之間關(guān)系，利用軟件將筆劃畫面疊加到保存的光柵視頻畫面或空白畫面之上。先用最簡(jiǎn)單的畫點(diǎn)方法將筆劃畫面疊加到光柵畫面上。若在(X，Y)點(diǎn)，光柵筆劃切換和顏色信號(hào)有效，則在保存的光柵畫面上的相應(yīng)位置上畫一個(gè)點(diǎn)。如此這般將所有筆劃上的點(diǎn)全部畫在光柵畫面上，就形成了一幅完整的光柵筆劃疊加顯示畫面。也可采用分段畫線的方法來顯示筆劃畫面，即若在(X1，Y1)點(diǎn)到(X2，Y2)點(diǎn)光柵筆劃切換和顏色信號(hào)均持續(xù)有效，則可將兩點(diǎn)之間所有筆劃上的點(diǎn)都畫在光柵畫面上且點(diǎn)與點(diǎn)之間用直線連接。疊加后的光柵筆劃畫面可直接在測(cè)控計(jì)算機(jī)屏上顯示。

在光柵筆劃疊加顯示方式下，筆劃信號(hào)在光柵視頻信號(hào)的場(chǎng)消隱期間顯示，光柵視頻信號(hào)采用PAL制式，場(chǎng)頻為50場(chǎng)/s，幀頻為25幀/s，掃描線為625行，場(chǎng)消隱時(shí)間為幾毫秒，故筆劃信號(hào)的顯示時(shí)間為毫秒級(jí)。一般飛機(jī)平顯的分辨率為512*512，假設(shè)有一幅很復(fù)雜的筆劃畫面，在沒有空筆移動(dòng)的情況下，其畫面占據(jù)了半個(gè)屏幕，則這幅筆劃畫面所占據(jù)的像素點(diǎn)為(512*512)/2，這些像素點(diǎn)要在場(chǎng)消隱期間顯示完，則筆劃信號(hào)的頻率為幾十兆。一般的筆劃畫面主要是字符、數(shù)字和表示航向、距離、高度等重要信息的特定技術(shù)符號(hào)，畫面比較簡(jiǎn)單，其占據(jù)的像素點(diǎn)并不多，在毫秒級(jí)的場(chǎng)消隱時(shí)間相應(yīng)的頻率為幾兆。則在實(shí)際采樣過程中的采樣頻率要達(dá)到幾十兆。

3 測(cè)試數(shù)據(jù)處理

3.1 測(cè)試數(shù)據(jù)分析

筆劃信號(hào)主要由顏色信號(hào)及偏轉(zhuǎn)信號(hào)組成。對(duì)單色多功能顯示器，只有一個(gè)顏色信號(hào)或輝亮信號(hào)。彩色多功能顯示器則包含筆劃顏色R、筆劃顏色G、筆劃顏色B這3路信號(hào)。對(duì)彩色多功能顯示器中3路顏色信號(hào)的處理可按照單色多功能顯示器中輝亮信號(hào)的處理方法進(jìn)行類似處理。

輝亮信號(hào)為RS422電平，該種電平采用差分傳輸(平衡傳輸)的方式。采用該種傳輸方式輝亮信號(hào)在采集過程中抗干擾能力強(qiáng)［5］，對(duì)采集后的輝亮信號(hào)只需進(jìn)行二值處理就基本能夠消除噪聲的影響并能達(dá)到時(shí)序的要求。

偏轉(zhuǎn)信號(hào)為±5 V范圍內(nèi)的連續(xù)模擬信號(hào)，在采集過程中由于儀器自身及外界干擾等因素的影響，致使采集后的偏轉(zhuǎn)信號(hào)有大量的噪聲，筆劃偏轉(zhuǎn)信號(hào)的偏轉(zhuǎn)精度高，噪聲對(duì)其影響明顯。若不對(duì)筆劃偏轉(zhuǎn)信號(hào)進(jìn)行處理就直接利用這些數(shù)據(jù)進(jìn)行筆劃畫面的復(fù)現(xiàn)，那么顯示出的畫面會(huì)模糊不清并影響圖像識(shí)別、測(cè)試等后續(xù)工作的開展。筆劃偏轉(zhuǎn)信號(hào)頻率分布不均，在輝亮信號(hào)消隱時(shí)的頻率要高于輝亮信號(hào)有效時(shí)的頻率，在采集過程中產(chǎn)生的噪聲干擾一般分布在高頻段，這就造成在進(jìn)行濾波時(shí)對(duì)輝亮信號(hào)消隱段內(nèi)筆劃偏轉(zhuǎn)信號(hào)的處理更加困難，由前文對(duì)筆劃畫面形成的分析可知，在輝亮信號(hào)消隱期內(nèi)是不顯示筆劃的，因此可以考慮采用分段濾波的方法對(duì)偏轉(zhuǎn)信號(hào)進(jìn)行處理，即只對(duì)輝亮信號(hào)有效期內(nèi)的筆劃偏轉(zhuǎn)數(shù)據(jù)進(jìn)行濾波處理而不管輝亮信號(hào)消隱期內(nèi)的偏轉(zhuǎn)數(shù)據(jù)。對(duì)筆劃偏轉(zhuǎn)數(shù)據(jù)進(jìn)行分段濾波處理，不僅節(jié)省了濾波時(shí)間而且還可以充分考慮每段偏轉(zhuǎn)數(shù)據(jù)的特點(diǎn)從而使處理過程更有針對(duì)性，效果也會(huì)更好。

筆劃各路信號(hào)之間存在著嚴(yán)格的時(shí)序和邏輯關(guān)系，這就要求在處理過程中要包含一定時(shí)頻的信息，在傳統(tǒng)的傅里葉分析中，信號(hào)完全是在頻域展開的，不包含任何時(shí)頻的信息，故若單單只使用低通濾波器、帶通濾波器等簡(jiǎn)單的頻域處理方法進(jìn)行濾波則很容易就將大量有用的信息損失掉［6］。考慮到筆劃畫面顯示的主要內(nèi)容是字符、數(shù)字和表示航向、距離、高度等重要信息的特定技術(shù)符號(hào)，在顯示這些字符、符號(hào)時(shí)，往往是按照其形狀特征進(jìn)行編程實(shí)現(xiàn)的，所以產(chǎn)生的用來確定筆劃位置的偏轉(zhuǎn)X、偏轉(zhuǎn)Y信號(hào)數(shù)據(jù)前后聯(lián)系緊密且變化緩和，故考慮使用空域?yàn)V波的方法對(duì)采集得到的筆劃偏轉(zhuǎn)信號(hào)進(jìn)行處理。

綜上所述，輝亮信號(hào)在采集過程中抗干擾能力強(qiáng)，對(duì)采集后的數(shù)據(jù)只需進(jìn)行二值處理即可。筆劃偏轉(zhuǎn)信號(hào)受噪聲影響明顯，要進(jìn)行預(yù)處理后才能用于筆劃畫面的復(fù)現(xiàn)?？紤]到筆劃畫面的形成過程及筆劃各路信號(hào)之間存在著嚴(yán)格的時(shí)序和邏輯關(guān)系，最終選用分段空域?yàn)V波的方法對(duì)筆劃偏轉(zhuǎn)信號(hào)進(jìn)行處理。

3.2 Savitzky-Golay平滑濾波器

Savitzky和Golay提出的S-G平滑濾波器，又稱最小二乘方法或數(shù)據(jù)平滑多項(xiàng)式濾波器。該濾波器通過將一個(gè)多項(xiàng)函數(shù)對(duì)信號(hào)進(jìn)行逐個(gè)擬合，使含有噪聲的信號(hào)變得平滑。通過該濾波器處理過的數(shù)據(jù)能夠清晰地反映出數(shù)據(jù)的長(zhǎng)期變化趨勢(shì)以及局部的突變信息，且不受數(shù)據(jù)時(shí)間空間尺度和采集設(shè)備的限制。近年來，基于S-G濾波方法的NDVI時(shí)間序列重建取得了較好的效果［7］?，F(xiàn)將該方法運(yùn)用于對(duì)筆劃偏轉(zhuǎn)信號(hào)數(shù)據(jù)的處理。

S-G濾波的公式為

其中:Yj+i為原始筆劃偏轉(zhuǎn)信號(hào)數(shù)據(jù)序列中的第j+i位的值;為平滑過后序列中的第j位的數(shù)據(jù)值;m為滑動(dòng)窗口中的單側(cè)數(shù)據(jù)點(diǎn)數(shù);Ci為從過濾器首部開始第i個(gè)筆劃偏轉(zhuǎn)信號(hào)數(shù)據(jù)值的權(quán)值;N為過濾器的長(zhǎng)度，同時(shí)也可以表示為(2m+1)［8－9］。

因此S-G濾波器要人為設(shè)定兩個(gè)參數(shù):第一個(gè)是m;第二個(gè)是多項(xiàng)式擬合的階數(shù)。通常m值越大結(jié)果越平滑，被平滑的峰谷值也就越多;多項(xiàng)式擬合的階數(shù)一般設(shè)在2～4之內(nèi)，較低的次數(shù)可以得到更平滑的結(jié)果，但會(huì)引進(jìn)誤差，較高的次數(shù)可以降低這個(gè)誤差，但是可能過于擬合得到一個(gè)更多噪聲的結(jié)果。

3.3 S-G平滑濾波方法改進(jìn)

由于在采集筆劃偏轉(zhuǎn)數(shù)據(jù)過程中可能出現(xiàn)由于某些元件造成的異常點(diǎn)。因此，在使用Savitzky-Golay濾波算法對(duì)筆劃偏轉(zhuǎn)數(shù)據(jù)進(jìn)行預(yù)處理之前，可以先采用3σ準(zhǔn)則進(jìn)行異常點(diǎn)檢測(cè)［10］，經(jīng)過異常點(diǎn)檢測(cè)后的數(shù)據(jù)集能更好的進(jìn)行濾波。

3σ準(zhǔn)則又稱拉伊達(dá)準(zhǔn)則，其數(shù)學(xué)方法如下:設(shè)樣本數(shù)據(jù)為y1，y2，…，yn，平均值為，偏差為vi=yi－，(i=1，2，…，n)。計(jì)算出標(biāo)準(zhǔn)差為

如果某一樣本數(shù)據(jù)yk的偏差vk(1≤k≤n)滿足|vk|＞3σ，則認(rèn)為該數(shù)據(jù)不合理。

在筆劃畫面中出現(xiàn)的字符、符號(hào)形狀不同，大小各異，復(fù)雜程度也有區(qū)別，這就使得對(duì)于每一個(gè)字符、符號(hào)其采集得到的數(shù)據(jù)量不同，因此，平滑窗口寬度不能為固定值。本設(shè)計(jì)采用動(dòng)態(tài)的方法，即每次預(yù)處理采用的窗寬與需要處理的數(shù)據(jù)數(shù)量有關(guān)，根據(jù)數(shù)量的多少來確定窗寬。

3.4 數(shù)據(jù)處理流程

對(duì)筆劃偏轉(zhuǎn)數(shù)據(jù)進(jìn)行分段處理后，首先進(jìn)行異常點(diǎn)檢測(cè)，具體算法如下:

3)計(jì)算ΔX中絕對(duì)值大于3×σ的數(shù)據(jù)點(diǎn)，認(rèn)為是異常點(diǎn)，并刪除;

4)對(duì)于被刪除的異常點(diǎn)，用該點(diǎn)兩側(cè)的正常數(shù)據(jù)進(jìn)行線性差值，并作為改點(diǎn)數(shù)據(jù)。

由于上述設(shè)計(jì)的多通道并行數(shù)據(jù)采集模塊的采樣率是固定的，每個(gè)點(diǎn)之間的時(shí)間間隔是相同的，即偏轉(zhuǎn)X、偏轉(zhuǎn)Y等數(shù)據(jù)對(duì)應(yīng)的曲線的橫坐標(biāo)是均勻的，因此可以對(duì)其使用Savitzky-Golay法。綜上所述采取下面的預(yù)處理算法:

1)對(duì)筆劃偏轉(zhuǎn)X、偏轉(zhuǎn)Y數(shù)據(jù)進(jìn)行異常點(diǎn)剔除;

2)對(duì)上一步處理的結(jié)果進(jìn)行卷積平滑以消除噪聲，可用Savitzky-Golay法求得，這里需要確定的參數(shù)是窗口半寬w和平滑多項(xiàng)式的階數(shù)。

在進(jìn)行濾波前，首先要確定窗寬及多項(xiàng)式的階數(shù)。經(jīng)過多次試取，可選取窗口半寬w=m，為數(shù)據(jù)點(diǎn)的1/30時(shí)效果最佳。筆劃偏轉(zhuǎn)數(shù)據(jù)變化比較緩和，分別采用一階、二階、三階、四階多項(xiàng)式進(jìn)行處理，當(dāng)多項(xiàng)式階數(shù)定位三階時(shí)效果最令人滿意。

3.5 數(shù)據(jù)處理效果

下面通過一個(gè)仿真實(shí)驗(yàn)來檢測(cè)一下數(shù)據(jù)預(yù)處理的效果。本實(shí)驗(yàn)利用LabView軟件，根據(jù)筆劃信號(hào)的時(shí)序關(guān)系產(chǎn)生出3路信號(hào)數(shù)據(jù)，即偏轉(zhuǎn)X、偏轉(zhuǎn)Y和輝亮信號(hào)，并將這3路信號(hào)數(shù)據(jù)分別保存在3個(gè)文件中。在產(chǎn)生過程中人為地進(jìn)行均勻白噪聲干擾，目的是模擬在實(shí)際采集過程中由于儀器自身或外界干擾等因素造成的噪聲。每路信號(hào)數(shù)據(jù)長(zhǎng)度為100 kB，筆劃偏轉(zhuǎn)電壓±5 V。

圖3 偏轉(zhuǎn)X處理前后對(duì)比Fig.3 The effect of deflection X processing

通過LabView讀取這3個(gè)包含噪聲的數(shù)據(jù)文件，采用分段畫線的方法顯示出由這3路信號(hào)數(shù)據(jù)形成的筆劃畫面，即若在(X1，Y1)點(diǎn)到(X2，Y2)點(diǎn)輝亮信號(hào)持續(xù)有效，則可將兩點(diǎn)之間所有筆劃上的點(diǎn)都顯示出來且點(diǎn)與點(diǎn)之間用直線連接。在讀取這3個(gè)包含噪聲的數(shù)據(jù)文件進(jìn)行筆劃畫面復(fù)現(xiàn)時(shí)，為使顯示出來的筆劃畫面更加真實(shí)清晰，就有必要在顯示前先對(duì)這3路信號(hào)數(shù)據(jù)進(jìn)行預(yù)處理。采用上述算法對(duì)筆劃數(shù)據(jù)進(jìn)行處理，預(yù)處理前后偏轉(zhuǎn)X和偏轉(zhuǎn)Y的坐標(biāo)數(shù)據(jù)對(duì)比如圖3、圖4所示。可以看出，使用本設(shè)計(jì)的預(yù)處理算法，可以有效地將噪聲去除。

圖4 偏轉(zhuǎn)Y處理前后對(duì)比Fig.4 The effect of deflection Y processing

若直接讀取上述產(chǎn)生的3個(gè)文件而不對(duì)其進(jìn)行預(yù)處理，則顯示如圖5。圖5的顯示效果與真實(shí)采集筆劃數(shù)據(jù)在測(cè)控計(jì)算機(jī)上顯示的效果十分相似，故若通過數(shù)據(jù)處理能使圖5的顯示效果改善，那么在實(shí)際應(yīng)用中此種預(yù)處理算法也會(huì)有很好的效果。預(yù)處理前后該筆劃畫面的示意圖對(duì)比，如圖5、圖6所示。經(jīng)過本設(shè)計(jì)的預(yù)處理方法，顯示效果明顯好于以前。

圖5 預(yù)處理前的畫面示意圖Fig.5 Schematic diagram before preprocessing

圖6 預(yù)處理后的畫面示意圖Fig.6 Schematic diagram after preprocessing

4 結(jié)束語

在進(jìn)行數(shù)據(jù)處理時(shí)，先根據(jù)輝亮信號(hào)的邏輯狀態(tài)，對(duì)采集得到的筆劃偏轉(zhuǎn)信號(hào)數(shù)據(jù)進(jìn)行分段，通過分段只對(duì)輝亮信號(hào)有效期內(nèi)的筆劃偏轉(zhuǎn)數(shù)據(jù)進(jìn)行處理。采用3σ準(zhǔn)則進(jìn)行異常點(diǎn)檢測(cè)，經(jīng)過異常點(diǎn)檢測(cè)后的數(shù)據(jù)集能更好地進(jìn)行濾波。在此基礎(chǔ)上根據(jù)航空顯示信號(hào)中筆劃信號(hào)各分量之間嚴(yán)格的時(shí)序關(guān)系及筆劃偏轉(zhuǎn)信號(hào)頻率分布特征，采用Savitzky-Golay動(dòng)態(tài)濾波。仿真實(shí)驗(yàn)結(jié)果表明，采用該算法可以消除采集筆劃偏轉(zhuǎn)信號(hào)時(shí)大部分噪聲的干擾，同時(shí)該方法原理簡(jiǎn)單，易于編程實(shí)現(xiàn)。該方法在對(duì)筆劃信號(hào)進(jìn)行去噪時(shí)，相當(dāng)于去除了筆劃信號(hào)在產(chǎn)生時(shí)附帶的和采集過程中由于儀器及外界干擾產(chǎn)生的兩部分噪聲，若只想單單去除儀器及外界干擾的噪聲，則需進(jìn)一步討論。本文設(shè)計(jì)的算法同樣適用于其他信號(hào)數(shù)據(jù)的預(yù)處理，例如其對(duì)手寫系統(tǒng)中由手寫輸入設(shè)備采集得到的原始數(shù)據(jù)的處理也有令人滿意的效果。

［1］王昌金，高波，孫瑾，等.ATS中的光柵筆劃信號(hào)產(chǎn)生和檢測(cè)技術(shù)［J］.電光與控制，2005，12(5):45-49.

［2］劉海見，楊健.航空顯示信號(hào)自動(dòng)測(cè)試技術(shù)研究［J］.測(cè)控技術(shù)，2005，24(9):76-78.

［3］張興義.電子顯示技術(shù)［M］.北京:北京理工大學(xué)出版社，1995:42-50.

［4］WANG C J，SHEN T S，ZHOU X D.The stroke and raster dislay signal generator in automatic test equipment［C］//Proceedings of ISTM，2003.

［5］程佩青.數(shù)字信號(hào)處理教程［M］.北京:清華大學(xué)出版社，2008:323-370.

［6］劉文耀.數(shù)字圖像采集與處理［M］.北京:電子工業(yè)出版社，2007:91-100.

［7］李航燕，頡耀文.時(shí)序NDVI數(shù)據(jù)集重建方法評(píng)價(jià)與實(shí)例研究［J］.遙感技術(shù)與應(yīng)用，2009，24(5):596-602.

［8］SAVITZKY A，GOLAY M J E.Smoothing and differentiation of data by simplified least squares procedures［J］.Analytical Chemistry，1964，36:1627-1639.

［9］CHEN J，PEN J，MASAYUKI T，et al.A simple method for reconstructing a high-quality NDVI time-series data set based on the savitzky-golay filter［J］.Romote Sensing of Environment，2004，91:332-344.

［10］BOUVERESSE E，MASSART D L.Improvement of the piecewise direet standardization procedure for the transfer of NIR spectra for multivariate calibration［J］.Chemometrics and Intelligent Laboratory Systems，1996，32(2):201-203.