岳春宇，邢 坤，鮑云飛，周 楠，何紅艷

北京空間機(jī)電研究所，北京 100094

?

以交叉累積剩余熵為準(zhǔn)則的星載激光測高儀大光斑波形數(shù)據(jù)與地形匹配法

岳春宇，邢 坤，鮑云飛，周 楠，何紅艷

北京空間機(jī)電研究所，北京 100094

提出了一種基于交叉累積剩余熵的星載激光測高儀大光斑波形數(shù)據(jù)與地形匹配方法。根據(jù)星載激光測高儀大光斑回波波形信號(hào)包含地形結(jié)構(gòu)信息的特性，將激光回波波形數(shù)據(jù)和數(shù)字表面模型(DSM)投影到統(tǒng)計(jì)特征空間，建立數(shù)據(jù)的統(tǒng)計(jì)特征向量，消除數(shù)據(jù)間維度差異，以交叉累積剩余熵為相似性測度匹配波形數(shù)據(jù)與地形的統(tǒng)計(jì)特征。試驗(yàn)結(jié)果表明，本方法能夠較好地實(shí)現(xiàn)激光回波波形數(shù)據(jù)與地形的匹配，匹配精度達(dá)到一個(gè)像素以內(nèi)。

星載激光測高儀；大光斑回波信號(hào)；地形匹配；統(tǒng)計(jì)特征；交叉累積剩余熵

星載激光測高儀采用主動(dòng)遙感探測方式，能夠精確獲取地面高程信息，有效地彌補(bǔ)了傳統(tǒng)光學(xué)遙感系統(tǒng)工作條件苛刻以及三維探測能力弱等的不足。激光作為一種新型的大氣和地表遙感手段，未來將廣泛地應(yīng)用于地形測繪、氣象、海洋、環(huán)境資源等遙感應(yīng)用領(lǐng)域[1-5]。由于星載激光測高儀軌道、功率和能量等問題，能獲取的激光測高數(shù)據(jù)是采樣間隔很大且平面精度較低的大光斑回波波形數(shù)據(jù)。相對(duì)于較高的高程測量精度，目前GLAS等星載激光測高儀的平面定位精度僅為十幾米到幾十米[6]，不能直接進(jìn)行高精度定位，且每個(gè)單獨(dú)的光斑回波數(shù)據(jù)也不能直接與地形實(shí)現(xiàn)空間配準(zhǔn)，限制了星載激光遙感的應(yīng)用。目前星載大光斑激光數(shù)據(jù)還局限在極地冰蓋、海洋、植被等空間范圍較大且空間分辨率較低的宏觀性遙感應(yīng)用領(lǐng)域內(nèi)[7-8]。即將發(fā)射的ICESat-2，則設(shè)計(jì)通過增加光束數(shù)量，減小光斑直徑到10 m以內(nèi)來提高空間分辨率，減少坡度等對(duì)地形觀測的影響[9]，但在空間尺度上還是與高分辨率遙感圖像有很大差距。在星載激光測高儀檢校和測圖應(yīng)用方面[10-13]，可以通過激光測高數(shù)據(jù)與地形數(shù)據(jù)匹配，獲取控制點(diǎn)，準(zhǔn)確找到激光腳點(diǎn)位置，并根據(jù)地面控制數(shù)據(jù)檢校激光測高儀的幾何定位模型，實(shí)現(xiàn)高精度三維定位。另外，多源遙感數(shù)據(jù)融合應(yīng)用時(shí)，以激光測高數(shù)據(jù)與地形數(shù)據(jù)匹配為中介，可以實(shí)現(xiàn)激光數(shù)據(jù)與遙感圖像等其他地理數(shù)據(jù)的空間配準(zhǔn)。因此，研究大光斑星載激光測高數(shù)據(jù)與地形的匹配方法，充分挖掘光斑內(nèi)所包含地形結(jié)構(gòu)信息的作用，提高星載激光測高儀幾何定位能力，實(shí)現(xiàn)大光斑激光高精度對(duì)地觀測具有重大意義。

在地形探測中，星載激光測高儀發(fā)射并接收激光光束后向散射信號(hào)，可以直接提供觀測目標(biāo)的距離參數(shù)和高程信息，再聯(lián)合由衛(wèi)星姿軌參數(shù)建立的指向矢量，計(jì)算得到光斑的位置坐標(biāo)[7-14]。激光測高數(shù)據(jù)通過由大量回波觀測數(shù)據(jù)生成的粗格網(wǎng)低分辨率數(shù)字表面模型(digital surface model，DSM)與已有的高精度DSM配準(zhǔn)，來實(shí)現(xiàn)幾何校正[15-16]。在激光和立體相機(jī)聯(lián)合觀測體制下，可利用多源數(shù)據(jù)聯(lián)合平差來提高激光觀測值平面定位精度[5，17-18]。但上述方法都不是單獨(dú)激光光斑回波波形與地形之間的匹配。到目前位置，在攝影測量領(lǐng)域還沒有星載大光斑激光回波波形數(shù)據(jù)與地形匹配的相關(guān)研究報(bào)道。

本文根據(jù)星載大光斑激光測高數(shù)據(jù)回波信號(hào)包含了地形結(jié)構(gòu)信息的特性，在統(tǒng)計(jì)信息空間計(jì)算激光回波波形及參考DSM的累計(jì)剩余分布函數(shù)，以交叉累積剩余熵[19-20]為相似性測度，匹配激光回波波形及地形的空間位置。該方法規(guī)避了激光回波波形數(shù)據(jù)與DSM之間的維度差異，在地形參考下直接獲取激光光斑空間位置。經(jīng)過對(duì)試驗(yàn)區(qū)GLAS01激光測高波形數(shù)據(jù)與DSM數(shù)據(jù)的匹配試驗(yàn)，證實(shí)本方法能夠?qū)崿F(xiàn)星載大光斑激光測高數(shù)據(jù)與地形的匹配。

1 星載大光斑激光測高數(shù)據(jù)與DSM的統(tǒng)計(jì)信息特征

1.1 大光斑激光回波信號(hào)

星載大光斑激光測高波形數(shù)據(jù)是由高速模數(shù)轉(zhuǎn)換器對(duì)探測激光回波進(jìn)行高分辨率采樣和數(shù)字量化，通過采樣計(jì)數(shù)統(tǒng)計(jì)獲得的。如圖1所示，回波信號(hào)的橫軸是時(shí)間(單位為ms)，縱軸是接收信號(hào)數(shù)字化采樣后量化的結(jié)果(單位為V)?；夭ㄐ盘?hào)橫軸的正方向表示時(shí)間增加，縱軸正方向表示能量強(qiáng)度增大。

圖1 大光斑激光測高回波信號(hào)示意圖Fig.1 Sketch map of large footprint laser ranging echo waveform

回波信號(hào)波形與時(shí)間采樣間隔和能量量化方式有關(guān)，這兩個(gè)因素分別作用于波形的橫軸和縱軸，決定了波形的形狀?；夭ㄐ盘?hào)開始接收的部分是背景噪聲，一般當(dāng)信號(hào)的強(qiáng)度達(dá)到某一閾值時(shí)才被視為地面真實(shí)觀測的有效信號(hào)?；夭ㄐ盘?hào)是距離的觀測值，表示高程值高的信號(hào)先被接收到。波形的峰代表光斑內(nèi)主要的高程值，全波形表示光斑范圍內(nèi)的地形結(jié)構(gòu)信息?；夭úㄐ慰梢酝ㄟ^改變橫軸采樣間隔，規(guī)定化縱軸數(shù)值范圍，投影到統(tǒng)計(jì)空間域，利用統(tǒng)計(jì)學(xué)手段處理分析。

1.2 DSM高程直方圖

DSM產(chǎn)品是二維的圖像，每個(gè)像素的灰度值表示該地理位置的高程值。為了與大光斑激光回波信號(hào)維度統(tǒng)一，本文利用高程直方圖來表示DSM中的信息。如圖2所示，圖2(a)為DSM圖像，圖2(b)為圖2(a)中DSM的高程直方圖。高程直方圖的橫軸為高程值(單位為m)，縱軸為像素?cái)?shù)量。高程直方圖的橫軸正方向表示高程值增大，縱軸正方向表示該高程值像素的數(shù)量增加。

對(duì)于相同地區(qū)，高程直方圖與DSM 圖像地面采樣間隔(ground sampling distance, GSD)或地面分辨率和圖像像素值量化位數(shù)相關(guān)。這兩個(gè)因素分別對(duì)應(yīng)1.1節(jié)中時(shí)間采樣間隔和能量量化方式，作用于高程直方圖的橫軸和縱軸，決定了直方圖的形狀。通過計(jì)算高程直方圖，可以將DSM圖像降維到一維信號(hào)，再與1.1節(jié)中所述一致，通過改變橫軸采樣間隔，規(guī)定化縱軸數(shù)值范圍，投影到統(tǒng)計(jì)特征空間，利用統(tǒng)計(jì)學(xué)手段處理分析，從而實(shí)現(xiàn)激光波形數(shù)據(jù)與DSM圖像信號(hào)維度的統(tǒng)一。

圖2 DSM高程直方圖示意Fig.2 Elevation histogram of DSM

1.3 交叉累積剩余熵

累計(jì)剩余熵(cumulative residual entropy，CRE)[19-20]僅考慮隨機(jī)變量的概率分布，統(tǒng)計(jì)大于某一閾值的所有信號(hào)信息，保持了信號(hào)的連續(xù)性，有效克服了噪聲對(duì)局部極值的影響，在連續(xù)域和離散域都可以適用。

(1)

式中，|X|=(|X1|,|X2|,…,|Xn|)，λ=(λ1,λ2，…,λn)且λi≥0，i=1,2,…,n。若|X|>λ，則|Xi|>λi。P(|X|>λ)為|X|>λ的概率，在CRE的計(jì)算中稱為多元生存函數(shù)。

當(dāng)存在與X定義域相同的向量Y時(shí)，則X在條件Y下的條件CRE為

(2)

X與Y的交叉累積剩余熵C(X,Y)定義為

C(X,Y)=ε(X)-E[ε(X|Y)]

(3)

計(jì)算1.1節(jié)中的星載大光斑激光測高回波波形數(shù)據(jù)和1.2節(jié)中的DSM高程直方圖兩種信號(hào)之間的交叉累積剩余熵，首先要在物理意義上統(tǒng)一兩種信號(hào)的向量方向，即將波形數(shù)據(jù)橫軸反轉(zhuǎn)，按時(shí)間從長到短排列，與高程直方圖表示的高程高低方向相對(duì)應(yīng)。此時(shí)波形數(shù)據(jù)向量用SL表示，高程直方圖向量用SE表示，則SL與SE的交叉累積剩余熵C(SL,SE)為

(4)

式中，l和e分別為SL與SE中的元素；PL(l)和PE(e)為SL與SE的邊緣概率密度；P(l,e)為SL與SE的聯(lián)合概率密度。C(SL,SE)值越大，表示SL與SE之間的相似性約大。

2 本文方法匹配流程

本文波形數(shù)據(jù)與地形匹配方法如圖3所示。

圖3 基于CCRE的大光斑激光波形數(shù)據(jù)與地形匹配流程Fig.3 Big footprint laser waveform and DSM matching method

基于交叉累積剩余熵的星載激光測高儀大光斑波形數(shù)據(jù)與地形匹配方法流程如下：

(1) 波形數(shù)據(jù)提取，減少背景噪聲對(duì)真實(shí)信號(hào)的影響。

(2) 根據(jù)波形數(shù)據(jù)的初始地理位置獲取DSM匹配搜索區(qū)域。

(3) 對(duì)去噪后的波形數(shù)據(jù)按照時(shí)間從長到短重新排列，獲得與DSM高程直方圖一致的高程由低到高的方向矢量。

(4) 設(shè)置波形數(shù)據(jù)和高程直方圖采樣間隔，規(guī)定化每個(gè)采樣間隔的數(shù)值，把波形數(shù)據(jù)和DSM圖像投影到統(tǒng)計(jì)特征空間，并確定匹配窗口大小、相似性測度閾值等匹配算法參數(shù)。

(5) 遍歷DSM搜索區(qū)，計(jì)算每個(gè)匹配窗口的CCRE值，其中CCRE值最大的匹配窗口的中心為匹配點(diǎn)，即為波形數(shù)據(jù)所表示的光斑重心。

3 試驗(yàn)與分析

采用美國ICESat衛(wèi)星星載激光高度計(jì)GLAS(geoscience laser altimeter system)獲取的大光斑激光對(duì)地觀測回波波形數(shù)據(jù)GLAS01產(chǎn)品和地面DSM測繪產(chǎn)品匹配，來對(duì)本文提出方法進(jìn)行試驗(yàn)分析。GLAS01數(shù)據(jù)產(chǎn)品從美國國家冰雪數(shù)據(jù)中心(The National Snow and Ice Data Center，NSIDC)網(wǎng)站 (http:∥nsidc.ors/data/icesat)下載得到，光斑直徑約為70 m，相鄰光斑距離約為170 m。根據(jù)GLAS01產(chǎn)品中記錄的光斑號(hào)等信息，在與之相對(duì)應(yīng)的GLAS14數(shù)據(jù)產(chǎn)品中獲取每一個(gè)光斑對(duì)應(yīng)的經(jīng)緯度坐標(biāo)，并用GLAS14數(shù)據(jù)生成的DSM與參考DSM配準(zhǔn)，對(duì)每個(gè)光斑的平面坐標(biāo)進(jìn)行修正。為了實(shí)現(xiàn)波形數(shù)據(jù)與地形的匹配，DSM圖像的分辨率要盡可能多的高于激光光斑大小，所以本文試驗(yàn)DSM圖像空間分辨率為5 m×5 m，為1∶50 000標(biāo)準(zhǔn)產(chǎn)品。試驗(yàn)區(qū)域?yàn)榫挼闁|北部山區(qū)，WGS-84坐標(biāo)范圍為97.250 3°E—97.257 3°E，22.772 1°N—22.785 1°N。

本文激光回波波形試驗(yàn)數(shù)據(jù)集包括6個(gè)光斑，具體信息如表1所示。

GLAS數(shù)據(jù)的坐標(biāo)使用Topex/Poseidon橢球，該坐標(biāo)系與WGS-84橢球差異在平面坐標(biāo)上僅相差幾厘米，與試驗(yàn)數(shù)據(jù)空間分辨率相比可以忽略。而兩個(gè)坐標(biāo)系下高程差異ΔH由緯度B及兩個(gè)橢球長短半軸(a1,b1)、(a2,b2)的函數(shù)ΔH=-[(a2-a1)cos2B+(b2-b1)sin2B]來計(jì)算，ΔH在一個(gè)光斑內(nèi)是系統(tǒng)性的[21]。對(duì)于波形數(shù)據(jù)來說相當(dāng)于整個(gè)波形在縱軸上平移一個(gè)常數(shù)值，在映射到統(tǒng)計(jì)特征空間時(shí)沒有影響。

表1 激光回波波形數(shù)據(jù)集Tab.1 The laser echo waveform dataset

3.1 波形數(shù)據(jù)提取

本文GLAS01試驗(yàn)數(shù)據(jù)每個(gè)光斑回波信號(hào)包含1000個(gè)采樣點(diǎn)，一般前1/3的采樣點(diǎn)不包含真實(shí)信號(hào)。本文選擇前300個(gè)采樣點(diǎn)計(jì)算回波信號(hào)的背景噪聲。根據(jù)背景噪聲來計(jì)算確定真實(shí)信號(hào)波形起止位置的閾值。閾值通常為背景噪聲均值與4倍標(biāo)準(zhǔn)差的和[22]。激光回波波形試驗(yàn)數(shù)據(jù)集的波形數(shù)據(jù)提取結(jié)果如圖4所示，圖中橫坐標(biāo)表示采樣數(shù)，縱坐標(biāo)為回波信號(hào)強(qiáng)度，單位為伏特(V)。

圖4 波形數(shù)據(jù)提取Fig.4 Waveform data detection

圖4中波形橫軸為采樣間隔，縱軸為能量，是回波信號(hào)電壓值。從圖4可以看出，波形提取結(jié)果中數(shù)據(jù)1,2,3,4,6均較好地保留了波型中的主要信號(hào)，數(shù)據(jù)5相對(duì)于其他數(shù)據(jù)波形較平滑，從物理意義上來說表示光斑中地形結(jié)構(gòu)不明顯，所以波形提取結(jié)果僅保留了主要波峰，去掉了一些次要的信息。

3.2 統(tǒng)計(jì)特征空間變換

由1.1節(jié)和1.2節(jié)可知，激光回波波形數(shù)據(jù)和DSM圖像高程直方圖縱坐標(biāo)物理意義不一致，需要統(tǒng)一到一個(gè)量化范圍，本文選擇8 bit即0～255。 匹配窗口應(yīng)與光斑尺寸相對(duì)應(yīng)，邊長為至少70 m/5 m=14(像素)，為了防止直方圖疏散陷入局部極值，將高程直方圖橫軸壓縮為8級(jí)，則激光回波波形數(shù)據(jù)橫軸也相應(yīng)壓縮為8級(jí)。經(jīng)過3.1節(jié)中波形提取后的波形數(shù)據(jù)轉(zhuǎn)換到統(tǒng)計(jì)特征空間后如圖5所示，縱軸表示特征空間變換后的量化的值。

圖5 激光試驗(yàn)數(shù)據(jù)統(tǒng)計(jì)特征空間分布Fig.5 The laser data set in statistics domain

由圖4和圖5可以看出，激光回波波形數(shù)據(jù)轉(zhuǎn)換到統(tǒng)計(jì)特征空間后保持了原始波形形狀。其中激光回波波形數(shù)據(jù)由1.3節(jié)中所述重新排列。而DSM圖像映射到統(tǒng)計(jì)特征空間為1.2節(jié)中的直方圖。此時(shí)兩種數(shù)據(jù)實(shí)現(xiàn)了形式上的統(tǒng)一。

3.3 波形與地形匹配

在統(tǒng)計(jì)特征空間采用CCRE為相似性測度對(duì)激光回波波形數(shù)據(jù)與DSM圖像進(jìn)行匹配，為了獲取更多的地形信息避免地形起伏對(duì)光斑大小的影響，匹配窗口大小為85 m×85 m，對(duì)應(yīng)DSM圖像為17像素×17像素，把GLAS01波形數(shù)據(jù)對(duì)應(yīng)的地理坐標(biāo)，換算得到DSM圖像的像素坐標(biāo)，作為真值評(píng)價(jià)匹配精度。匹配結(jié)果見表2。

3.4 匹配結(jié)果分析

試驗(yàn)中的匹配結(jié)果為最大CCRE值所在的匹配窗口的中心像素，但試驗(yàn)數(shù)據(jù)GLAS01的地理坐標(biāo)投影到DSM圖像時(shí)一般不是整像素，所以表2中匹配精度為非整像素。從表2中可以看出，匹配結(jié)果沿經(jīng)度方向優(yōu)于沿緯度方向，沿經(jīng)度方向?yàn)?.3個(gè)像素左右，沿緯度方向具有一定系統(tǒng)性，接近1個(gè)像素，整體匹配精度為1個(gè)像素左右。

表2 匹配結(jié)果Tab.2 The matching results

由圖5可以看出，數(shù)據(jù)1與數(shù)據(jù)4激光回波波形數(shù)據(jù)只有一個(gè)主峰，即光斑內(nèi)部地形結(jié)構(gòu)簡單，地物起伏較小，投影到統(tǒng)計(jì)特征空間分布單一，所以匹配精度較低；而數(shù)據(jù)5激光回波波形數(shù)據(jù)波形沒有明顯的峰也不平滑，光斑內(nèi)地形結(jié)構(gòu)信息不明顯，所以匹配精度較低；數(shù)據(jù)2與數(shù)據(jù)3激光回波波形數(shù)據(jù)波形較平坦，即光斑內(nèi)部地形高程值連續(xù)分布，投影到統(tǒng)計(jì)特征空間分布較分散，所以相對(duì)匹配精度高于數(shù)據(jù)1與數(shù)據(jù)4；數(shù)據(jù)6激光回波波形數(shù)據(jù)有幾個(gè)明顯的峰，光斑內(nèi)地形起伏較大，地形結(jié)構(gòu)明顯，所以匹配精度最高。

在本文方法中，基于地理位置設(shè)置搜索區(qū)域后，大光斑激光回波波形與DSM圖像異源數(shù)據(jù)匹配僅考慮數(shù)據(jù)的統(tǒng)計(jì)信息，所以匹配誤差主要來源于匹配方法。交叉累積剩余熵作為相似性測度，屬于基于窗口的匹配方法[23]，以匹配窗口內(nèi)信息作為匹配基元，即使克服了噪聲的影響，與其他基于窗口的匹配方法一樣，可能陷入局部極值，造成匹配錯(cuò)誤。由于激光回波波形數(shù)據(jù)與DSM圖像數(shù)據(jù)在統(tǒng)計(jì)特征空間中，在匹配窗口內(nèi)作為整體進(jìn)行匹配，在統(tǒng)計(jì)特征空間中向量每一個(gè)元素?zé)o法與空間位置對(duì)應(yīng)，所以當(dāng)搜索區(qū)出現(xiàn)相同地形統(tǒng)計(jì)特征分布區(qū)域時(shí)會(huì)存在誤匹配。只有地形結(jié)構(gòu)信息明顯時(shí)，且根據(jù)先驗(yàn)知識(shí)獲取的搜索區(qū)域位置精度較高時(shí)，才會(huì)取得較好的匹配結(jié)果。另外本文試驗(yàn)結(jié)果中出現(xiàn)的沿緯度方向的系統(tǒng)性誤差可能是DSM圖像的系統(tǒng)性誤差造成的。

4 結(jié) 論

本文根據(jù)星載大光斑激光測高數(shù)據(jù)回波信號(hào)和DSM圖像包含地形結(jié)構(gòu)信息的特性，將激光回波波形數(shù)據(jù)與DSM圖像投影到統(tǒng)計(jì)特征空間，消除激光回波波形數(shù)據(jù)與DSM之間的維度差異，以交叉累積剩余熵為相似性測度匹配激光回波波形及地形的空間位置。通過GLAS01激光測高回波波形數(shù)據(jù)與DSM圖像數(shù)據(jù)的匹配試驗(yàn)，證實(shí)本方法能夠?qū)崿F(xiàn)星載大光斑激光測高回波波形數(shù)據(jù)與地形的匹配，匹配結(jié)果可用于星載激光測高儀幾何標(biāo)定以及與其他遙感數(shù)據(jù)融合應(yīng)用。

經(jīng)過試驗(yàn)分析，本文方法激光回波波形數(shù)據(jù)與DSM圖像空間位置的匹配精度在1個(gè)像素以內(nèi)。但是在統(tǒng)計(jì)特征空間，激光回波波形數(shù)據(jù)作為一個(gè)整體進(jìn)行匹配，光斑內(nèi)地形結(jié)構(gòu)信息不包含絕對(duì)地理信息參考，所以當(dāng)搜索區(qū)出現(xiàn)相同統(tǒng)計(jì)特征分布時(shí)會(huì)誤匹配。如何在統(tǒng)計(jì)特征空間加入地形結(jié)構(gòu)約束，提高匹配精度，同時(shí)擴(kuò)大搜索區(qū)域，還需要進(jìn)一步研究。

[1] 李松.星載激光測高儀發(fā)展現(xiàn)狀綜述[J].光學(xué)與光電技術(shù), 2004, 2(6): 4-6.LI Song.Recent Progress of Spaceborne Laser Altimeter System[J].Optics & Optoelectronic Technology, 2004, 2(6): 4-6.

[2] 李然, 王成, 蘇國中, 等.星載激光雷達(dá)的發(fā)展與應(yīng)用[J].科技導(dǎo)報(bào), 2007, 25(14): 58-63.LI Ran, WANG Cheng, SU Guozhong, et al.Development and Applications of Spaceborne LiDAR[J].Science & Technology Review, 2007, 25(14): 58-63.

[3] 劉煥玲.利用ICESat激光測高數(shù)據(jù)探測南極冰蓋的高程變化[D].阜新: 遼寧工程技術(shù)大學(xué), 2011.LIU Huanling.The Study of Ice Sheet Height Change in Antarctic by Using ICESat Laser Altimeter Data[D].Fuxin: Liaoning Technical University, 2011.

[4] 岳春宇, 鄭永超, 陶宇亮.星載激光測高儀輔助衛(wèi)星攝影測量淺析[J].航天返回與遙感, 2013, 34(4): 71-76.YUE Chunyu, ZHENG Yongchao, TAO Yuliang.Study on Space-borne Laser Altimeter Supported Satellite Photogrammetry[J].Spacecraft Recovery & Remote Sensing, 2013, 34(4): 71-76.

[5] WU Bo, HU Han, GUO Jian.Integration of Chang’E-2 Imagery and LRO Laser Altimeter Data with a Combined Block Adjustment for Precision Lunar Topographic Modeling[J].Earth and Planetary Science Letters, 2014, 391: 1-15.

[6] BAE S, WEBB C, SCHUTZ B.GLAS PAD Calibration Using Laser Reference Sensor Data[C]∥Proceedings of AIAA/AAS Astrodynamics Specialist Conference and Exhibit.Providence, RI: AIAA, 2004: 1-10.

[7] ZWALLY H J, SCHUTZ B, ABDALATI W, et al.ICESat’S Laser Measurements of Polar Ice, Atmosphere, Ocean, and Land[J].Journal of Geodynamics, 2002, 34(3-4): 405-445.

[8] NELSON R.Model Effects on GLAS-based Regional Estimates of Forest Biomass and Carbon[J].International Journal of Remote Sensing, 2010, 31(5): 1359-1372.

[9] ABDALATI W, ZWALLY H J, BINDSCHADLER R, et al.The ICESat-2 Laser Altimetry Mission[J].Proceedings of IEEE, 2010, 98(5): 735-751.

[10] 唐新明, 李國元, 高曉明, 等.衛(wèi)星激光測高嚴(yán)密幾何模型構(gòu)建及精度初步驗(yàn)證[J].測繪學(xué)報(bào), 2016, 45(10): 1182-1191.DOI: 10.11947/j.AGCS.2016.20150357.TANG Xinming, LI Guoyuan, GAO Xiaoming, et al.The Rigorous Geometric Model of Satellite Laser Altimeter and Preliminarily Accuracy Validation[J].Acta Geodaetica et Cartographica Sinica, 2016, 45(10): 1182-1191.DOI: 10.11947/j.AGCS.2016.20150357.

[11] 王任享, 王建榮.二線陣CCD衛(wèi)星影像聯(lián)合激光測距數(shù)據(jù)光束法平差技術(shù)[J].測繪科學(xué)技術(shù)學(xué)報(bào), 2014, 31(1): 1-4.WANG Renxiang, WANG Jianrong.Technology of Bundle Adjustment Using Two-line-array CCD Satellite Image Combined Laser Ranging Data[J].Journal of Geomatics Science and Technology, 2014, 31(1): 1-4.

[12] 李鑫, 廖鶴, 趙美玲, 等.激光測繪衛(wèi)星對(duì)不同地表形貌探測能力分析[J].測繪學(xué)報(bào), 2014, 43(12): 1238-1244.DOI: 10.13485/j.cnki.11-2089.2014.0188.LI Xin, LIAO He, ZHAO Meiling, et al.Research on LiDAR Surveying Satellite Detection Capacity for Different Terrains[J].Acta Geodaetica et Cartographica Sinica, 2014, 43(12): 1238-1244.DOI: 10.13485/j.cnki.11-2089.2014.0188.

[13] 李國元, 唐新明, 王華斌, 等.GLAS激光測高數(shù)據(jù)輔助的資源三號(hào)三線陣區(qū)域網(wǎng)平差研究[C]∥第三屆高分辨率對(duì)地觀測學(xué)術(shù)年會(huì)分會(huì)論文集.北京: 中國科學(xué)院重大科技任務(wù)局, 2014.LI Guoyuan, TANG Xinming, WANG Huabin, et al.ZY-3 Three Line Array Images Block Adjustment Supported by GLAS Data[C]∥The 3rd Annual Conference on High Resolution Earth Observation.Beijing: Major Science and Technology Bureau of the Chinese Academy of Sciences, 2014.

[14] 李磊, 鄭永超, 彭鳳超, 等.地形測繪激光成像雷達(dá)技術(shù)研究[J].紅外與激光工程, 2006, 35(S): 294-298.LI Lei, ZHENG Yongchao, PENG Fengchao, et al.Research of Three-dimension Imaging LiDAR on Land Topography[J].Infrared and Laser Engineering, 2006, 35(S): 294-298.

[15] 范春波, 李建成, 王丹, 等.ICESat/GLAS激光腳點(diǎn)定位及誤差分析[J].大地測量與地球動(dòng)力學(xué), 2007, 27(1): 104-106.FAN Chunbo, LI Jiancheng, WANG Dan, et al.ICESAT/ GLAS Laser Footprint Geolocation and Error Analysis[J].Journal of Geodesy and Geodynamics, 2007, 27(1): 104-106.

[16] CARABAJAL C C, HARDING D J, SUCHDEO V P.ICEsat LiDAR and Global Digital Elevation Models: Applications to Desdyni[C]∥Proceedings of 2010 IEEE International Geoscience and Remote Sensing Symposium (IGARSS).Honolulu, HI: IEEE, 2010: 1907-1910.

[17] 胡文敏, 邸凱昌, 岳宗玉, 等.嫦娥一號(hào)激光高度計(jì)數(shù)據(jù)交叉點(diǎn)分析與平差處理[J].測繪學(xué)報(bào), 2013, 42(2): 218-224.HU Wenmin, DI Kaichang, YUE Zongyu, et al.Crossover Analysis and Adjustment for Chang’E-1 Laser Altimeter Data[J].Acta Geodaetica et Cartographica Sinica, 2013, 42(2): 218-224.

[18] 趙雙明, 冉曉雅, 付建紅, 等.CE-1立體相機(jī)與激光高度計(jì)數(shù)據(jù)聯(lián)合平差[J].測繪學(xué)報(bào), 2014, 43(12): 1224-1229.DOI: 10.13485/j.cnki.11-2089.2014.0178.ZHAO Shuangming, RAN Xiaoya, FU Jianhong, et al.Combined Adjustment of CE-1 Stereo Camera Image and Laser Altimeter Data[J].Acta Geodaetica et Cartographica Sinica, 2014, 43(12): 1224-1229.DOI: 10.13485/j.cnki.11-2089.2014.0178.

[19] WANG Fei, VEMURI B C.Non-rigid Multi-modal Image Registration Using Cross-cumulative Residual Entropy[J].International Journal of Computer Vision, 2007, 74(2): 201-215.

[20] 江萬壽, 彭芳媛, 岳春宇, 等.利用Ratio梯度和交叉累積剩余熵進(jìn)行多源遙感影像匹配[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2009, 34(9): 1047-1050, 1121.JIANG Wanshou, PENG Fangyuan, YUE Chunyu, et al.Multi-source Remote-sensing Image Matching Based on Ratio-gradient and Cross-cumulative Residual Entropy[J].Geomatics and Information Science of Wuhan University, 2009, 34(9): 1047-1050, 1121.

[21] 鄂棟臣, 徐瑩, 張小紅.ICESat衛(wèi)星及其在南極Dome A地區(qū)的應(yīng)用[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2007, 32(12): 1139-1142.E Dongchen, XU Ying, ZHANG Xiaohong.ICESat’s Performance and Its Application in Dome A Area in Antarctica[J].Geomatics and Information Science of Wuhan University, 2007, 32(12): 1139-1142.

[22] 宋志英.利用ICESAT激光測高數(shù)據(jù)監(jiān)測南極冰蓋高程變化的研究[D].阜新: 遼寧工程技術(shù)大學(xué), 2009.SONG Zhiying.The Study of Height Change in Antarctic Ice Sheet by Using ICESAT Laser Altimeter Data[D].Fuxin: Liaoning Technical University, 2009.

[23] 岳春宇.面向平坦地區(qū)的星載SAR與光學(xué)圖像配準(zhǔn)關(guān)鍵技術(shù)研究[D].武漢: 武漢大學(xué), 2012.YUE Chunyu.Research on Registration of SAR and Optical Satellite Image in Flat Area[D].Wuhan: Wuhan University, 2012.

(責(zé)任編輯：張艷玲)

A Matching Method of Space-borne Laser Altimeter Big Footprint Waveform and Terrain Based on Cross Cumulative Residual Entropy

YUE Chunyu，XING Kun，BAO Yunfei，ZHOU Nan，HE Hongyan

Beijing Institute of Space Mechanics & Electricity，Beijing 100094，China

A matching method of space-borne laser altimeter big footprint waveform and terrain based on cross cumulative residual entropy(CCRE) is proposed.Firstly, the waveform data and digital surface model(DSM) data are projected to the statistics domain, according to the terrain structure information of the waveform, where statistics signal vectors of the two data are in the same dimension.Then, the waveform data and DSM image are matched in the statistics domain with CCRE.Experiments show that the algorithm proposed is effective in waveform and terrain matching, and the matching accuracy is within 1 pixel.

space-borne laser altimeter; big footprint waveform; terrain matching; statistics; cross cumulative residual entropy(CCRE)

The National Natural Science Foundation of China (Nos.41401411；41401410 )；The High Resolution Earth Observation Special Fund (No.GFZX040136)；The National Key Technology Research and Development Program of the Earth Observation and navigation Special Fund of China (No.2016YFB0500802)

YUE Chunyu(1983—)，male,PhD,senior engineer,majors in satellite photogrammetry，remote sensing image processing and application.

岳春宇，邢坤，鮑云飛，等.以交叉累積剩余熵為準(zhǔn)則的星載激光測高儀大光斑波形數(shù)據(jù)與地形匹配法[J].測繪學(xué)報(bào)，2017,46(3)：346-352.

10.11947/j.AGCS.2017.20160255.

YUE Chunyu，XING Kun，BAO Yunfei，et al.A Matching Method of Space-borne Laser Altimeter Big Footprint Waveform and Terrain Based on Cross Cumulative Residual Entropy[J].Acta Geodaetica et Cartographica Sinica,2017,46(3):346-352.DOI:10.11947/j.AGCS.2017.20160255.

P223

A

1001-1595(2017)03-0346-07

國家自然基金(41401411;41401410);高分對(duì)地觀測專項(xiàng)基金(GFZX040136);國家重點(diǎn)研發(fā)計(jì)劃地球觀測與導(dǎo)航重點(diǎn)專項(xiàng)基金(2016YFB0500802)

2016-05-25

修回日期：2017-01-09

岳春宇(1983—)，男，博士，高級(jí)工程師，研究方向?yàn)樾l(wèi)星攝影測量、遙感圖像處理與應(yīng)用等。

E-mail：ycy1893@163.com