徐涵秋，唐 菲

(福州大學(xué)遙感信息工程研究所，福州大學(xué)環(huán)境與資源學(xué)院，福州 350108)

在各種衛(wèi)星遙感對(duì)地觀測(cè)數(shù)據(jù)深入應(yīng)用到各行各業(yè)的今天，長(zhǎng)達(dá)40年歷史的Landsat系列衛(wèi)星數(shù)據(jù)無疑仍是應(yīng)用最廣泛的衛(wèi)星數(shù)據(jù)，它們?cè)谌虺叨鹊纳鷳B(tài)環(huán)境變化監(jiān)測(cè)中發(fā)揮了無可比擬的重要作用。但是，由于其兩顆主力衛(wèi)星中的Landsat 7號(hào)星的掃描行校正器于2003年5月31日發(fā)生了故障，使其實(shí)用價(jià)值大打折扣;而設(shè)計(jì)壽命只有3年的Landsat 5號(hào)星也因?yàn)榻?9年的超期服役而于2012年12月21日正式宣布退役［1］，從而造成Landsat 40年的連續(xù)對(duì)地觀測(cè)出現(xiàn)一度中斷。幾經(jīng)波折的Landsat 8衛(wèi)星終于在Landsat對(duì)地連續(xù)觀測(cè)中斷1年4個(gè)月后，于2013年2月11日在美國(guó)加州成功發(fā)射［2］。這意味著Landsat將繼續(xù)為全球的生態(tài)環(huán)境提供連續(xù)的觀測(cè)數(shù)據(jù)。

3月18日，Landsat 8獲得了第一幅遙感影像，并于29日作為樣本數(shù)據(jù)供用戶下載［3］。為了了解Landsat 8衛(wèi)星影像的基本特征、新增特征及其優(yōu)點(diǎn)，本文對(duì)該影像進(jìn)行了處理與分析，并重點(diǎn)測(cè)試了其新增特性的應(yīng)用，旨在為廣大用戶了解該影像的應(yīng)用潛力提供第一時(shí)間的基本信息。

1 Landsat 8衛(wèi)星影像的基本特征

Landsat 8衛(wèi)星是由美國(guó)宇航局(NASA)和美國(guó)地質(zhì)調(diào)查局(USGS)共同負(fù)責(zé)的項(xiàng)目。美國(guó)宇航局負(fù)責(zé)衛(wèi)星傳感器的研制和發(fā)射，美國(guó)地質(zhì)調(diào)查局負(fù)責(zé)后續(xù)的系統(tǒng)維護(hù)、地面接收、數(shù)據(jù)處理和分發(fā)。在研制和發(fā)射調(diào)試階段，美國(guó)宇航局將這一新星稱為L(zhǎng)andsat Data Continuity Mission(LDCM)，意指其肩負(fù)著Landsat衛(wèi)星連續(xù)對(duì)地觀測(cè)的使命。在發(fā)射后大約100天的調(diào)試期結(jié)束后(5月底)，LDCM將移交給美國(guó)地質(zhì)調(diào)查局。屆時(shí)，美國(guó)地質(zhì)調(diào)查局會(huì)將其改名為L(zhǎng)andsat 8，以延續(xù)Landsat系列衛(wèi)星自1972年發(fā)射以來，已持續(xù)了40年的全球?qū)Φ赜^測(cè)使命。正是由于這一使命，Landsat 8的設(shè)計(jì)和特征與Landsat 7基本相同，這使得Landsat 8數(shù)據(jù)可以和前期的Landsat數(shù)據(jù)保持很高的一致性和可比性。表1給出Landsat 8和Landsat 7衛(wèi)星影像的基本參數(shù)對(duì)比［4］。

表1 Landsat 8與Landsat 7影像的主要參數(shù)對(duì)比Table 1 Comparison of spectral bands between Landsat 8 and Landsat 7

從表1中可以看出，新的Landsat 8具有以下特征:

(1)波段設(shè)置:Landsat 8除了具有Landsat 7所有的光譜波段外，還有以下變化(圖1):

·在原藍(lán)光波段之外新增了1個(gè)深藍(lán)(Deep blue)波段，用于監(jiān)測(cè)近岸水體和大氣中的氣溶膠，因此也稱海岸/氣溶膠(Coastal/Aerosol)波段;

·在原近紅外波段與短波紅外波段間新增了1個(gè)卷云(Cirrus)波段，用于檢測(cè)卷云;

·將原熱紅外波段的光譜范圍一分為二，設(shè)置了2個(gè)熱紅外波段;

·收窄了原近紅外波段的光譜范圍，以便去除0.825μm處的水汽吸收影響;

·收窄了原全色波段的光譜范圍，新的全色波段的光譜范圍不再覆蓋近紅外波段。

圖1 Landsat 8傳感器與Landsat 7傳感器波段光譜設(shè)置范圍的對(duì)比(據(jù)NASA)Fig.1 Spectral band comparison between Landsat 8 and Landsat 7 sensors(from NASA)

(2)傳感器數(shù)量:Landsat 8具有2個(gè)傳感器:Operational Land Imager(OLI)和Thermal Infrared Sensor(TIRS)［5］，亦即Landsat 8具有單獨(dú)的熱紅外傳感器TIRS，而不像Landsat 7那樣，將其集成于ETM+傳感器之中。應(yīng)該說這是一個(gè)被動(dòng)之舉。由于Landsat 8在最初的設(shè)計(jì)上并沒有考慮要繼續(xù)接收熱紅外數(shù)據(jù)，后來決定要繼續(xù)接收熱紅外數(shù)據(jù)時(shí)，已設(shè)計(jì)的OLI傳感器已無足夠的空間來容納熱紅外傳感器。為了不耽誤發(fā)射時(shí)間，只得決定單獨(dú)設(shè)計(jì)另一新的傳感器，由于時(shí)間緊迫，所設(shè)計(jì)的熱紅外數(shù)據(jù)的分辨率為100m，這較之于Landsat 7的60m，不得不說是一個(gè)退步，而且設(shè)計(jì)的使用年限也只有3年，低于同衛(wèi)星的OLI的5年。但可喜的是，新的熱紅外傳感器有2個(gè)波段，比ETM+增加了1個(gè)波段，這使得它可以用劈窗算法來進(jìn)行大氣校正［6］。

(3)輻射分辨率:Landsat 8的OLI傳感器采用的是已在EO-1衛(wèi)星的ALI傳感器上實(shí)驗(yàn)過的推進(jìn)掃描方式［7］，所以將比所有老的Landsat系列衛(wèi)星的傳感器具有更高的信噪比。根據(jù)Irons等［6］的資料整理，OLI各波段的信噪比有望比Landsat 7 ETM+各對(duì)應(yīng)波段平均高出近3倍。因此Landsat 8將Landsat 7的8 bit輻射分辨率提高到12 bit，大大增加了影像的灰度量化級(jí)。

圖2 實(shí)驗(yàn)影像Fig.2 Subset of the first Landsat 8 image used as the test image

2 Landsat 8實(shí)驗(yàn)影像及其預(yù)處理

首幅Landsat 8遙感影像展示了美國(guó)懷俄明州和科羅拉多州的大平原區(qū)和落基山脈前沿交匯處的生態(tài)景觀。由于影像許多地方被云覆蓋，NASA和USGS的網(wǎng)站主要展示的是其中云較少的科羅拉多州博爾德(Boulder)地區(qū)處的部分影像［3］。因此，本文亦選該處影像作為實(shí)驗(yàn)分析影像，面積約80 km2(圖2)。該影像仍沿用Landsat系列數(shù)據(jù)的UTM/WGS84投影/坐標(biāo)系，數(shù)據(jù)的處理格式為L(zhǎng)evel 1T，即已進(jìn)行了基于地形的幾何校正。但由于撰稿時(shí)，Landsat 8仍在調(diào)試階段，其定標(biāo)參數(shù)文件(CPF)尚未發(fā)布，影像各波段的增益值(gain)、偏離值(bias)和太陽平均輻照度(ESUN)等參數(shù)未知，因此無法將影像的灰度值(DN)轉(zhuǎn)換為表觀反射率(TOA)［8］，主要采用Jensen的直方圖平移法對(duì)影像進(jìn)行輻射校正［9］。表2列出了經(jīng)預(yù)處理后的該影像的基本統(tǒng)計(jì)特征。

表2 實(shí)驗(yàn)影像的統(tǒng)計(jì)特征值Table 2 Statistics of the test image

圖3 實(shí)驗(yàn)影像主要土地覆蓋類型的反射光譜曲線Fig.3 Signatures of the main land cover classes of the test image

從表2可以看出，該影像中信息量最大的波段是短波紅外(OLI 6、7)波段。它們的標(biāo)準(zhǔn)差最大，對(duì)占影像總信息量83.57%的第一主成分(PC1)的貢獻(xiàn)也最大，二者的載荷分別達(dá)到0.677和0.533，遠(yuǎn)遠(yuǎn)超過其他波段。而新增的卷云波段(OLI 9)的信息量最小，其標(biāo)準(zhǔn)差和對(duì)PC1的貢獻(xiàn)度都是最小的。從相關(guān)性來看，可見光組的波段(OLI 1—4波段)具有很強(qiáng)的相關(guān)性，短波紅外組的波段(OLI 6—7、9波段)也具有較強(qiáng)的相關(guān)性，近紅外波段則與各波段的相關(guān)度都小于0.8。因此，從這3個(gè)波段組中各選1個(gè)波段進(jìn)行彩色組合，可以最大限度地降低波段之間的冗余度，獲得較好的目視效果。

從圖2可以看出，實(shí)驗(yàn)影像的西部為覆有林地的落基山脈前沿，北部為大片裸土，主要為第四紀(jì)的松散堆積［10］，中部為綠蔭遮蔽的居住用地，南部為商業(yè)區(qū)和科羅拉多大學(xué)博爾德分校。影像中還有少量人工草地和水塘點(diǎn)綴其中，最大的一處水體位于影像東北角。圖3和表3是影像中主要土地覆蓋類型的反射光譜曲線和統(tǒng)計(jì)特征。需要說明的是，NASA雖然將卷云波段編號(hào)為9波段，但其光譜范圍實(shí)際上介于OLI 5和OLI 6波段之間(表1)，因此在圖3和表3中，將其置于5波段和6波段之間，以客觀反映地物的光譜變化趨勢(shì)。

表3 實(shí)驗(yàn)影像主要土地覆蓋類型的均值Table 3 The mean values of the main land cover classes of the test image

3 Landsat 8影像新增特征及其生態(tài)環(huán)境意義

植被、水體、土壤是地表三個(gè)主要的生態(tài)環(huán)境要素、也是遙感生態(tài)環(huán)境監(jiān)測(cè)的主要地面目標(biāo)。因此，以下就這三大要素分析Landsat 8新增特性的優(yōu)勢(shì)。

3.1 植被監(jiān)測(cè)

Landsat 8全色波段(OLI 8波段)的一大變化就是其光譜范圍不再包含近紅外波段，以覆蓋綠光和紅光波段為主;而Landsat 7全色波段的波長(zhǎng)范圍則完整地覆蓋了綠光、紅光和近紅外波段的波長(zhǎng)范圍(圖1、表1)。為了考察這一變化帶來的影響，特選取了Landsat 7 ETM+同一地區(qū)2003-3-15的全色影像作為對(duì)比。該影像的季相與本次實(shí)驗(yàn)影像(3月18日)完全一致，雖然時(shí)隔10年，但仍可以看出該區(qū)的土地覆蓋在這期間并沒有發(fā)生太大變化(圖4)。通過對(duì)比可以發(fā)現(xiàn)，在OLI的全色影像上(圖4a)，植被呈暗色，如西邊落基山脈的林地和中部居住區(qū)里的綠化帶，因此植被和裸地、建筑用地等非植被區(qū)反差較明顯。而對(duì)于ETM+全色影像，直觀就可以看出其總體亮度較大，主要的原因是植被區(qū)并不太暗，特別是居住區(qū)的綠化帶和其中的建筑物的亮度較接近，二者的反差不像在OLI全色影像中那樣明顯，因此很難區(qū)別(圖4b)。顯然，ETM+全色波段由于涵蓋了近紅外波段，因此接收到了植被在這一波長(zhǎng)區(qū)間的反射光譜，從而導(dǎo)致其影像的亮度增大，植被與非植被信息難于區(qū)分;而OLI全色影像上不包含近紅外波段的光譜則有利于全色影像中植被和非植被信息的區(qū)別。

圖4 Landsat 8 OLI與Landsat 7 ETM+的全色影像Fig.4 Panchromatic images of OLI and ETM+sensors

3.2 水質(zhì)監(jiān)測(cè)

新增的深藍(lán)波段(OLI 1波段)主要用于監(jiān)測(cè)近岸水體。但Landsat 8的首幅影像中沒有海岸帶，因此主要根據(jù)實(shí)驗(yàn)影像東北角的水體對(duì)其考察(圖5a)。先利用MNDWI水體指數(shù)［11］，將水體提取出來，然后計(jì)算水體的均值。從圖3可以看出水體的均值從OLI 1到OLI 3波段為逐漸升高，然后逐漸下降。由于清水的反射率是隨著波長(zhǎng)的升高而降低，只有含懸浮物水體的反射率才會(huì)先升高后降低［12］，因此，可知這一水體的懸浮物濃度較高。從圖5可以看出，該水體在深藍(lán)波段的反射率低于藍(lán)光波段(前者的紅度不及后者)，二者的差值影像(藍(lán)光波段減深藍(lán)波段)突出了懸浮物濃度差，影像越紅的地方，濃度差距越大。顯然，深藍(lán)波段的出現(xiàn)有助于水體渾濁程度的判斷。

圖5 (a)水體影像，(b)、(c)深藍(lán)、藍(lán)光波段表現(xiàn)出的水體特征，(d)差值影像所突出的懸浮物濃度差，(e)16 bit增強(qiáng)影像，(f)8 bit增強(qiáng)影像Fig.5 (a)Original water image，(b)-(c)water feature showing in the deep blue and blue bands，(d)the differencing image of the two bands，(e)enhanced 16-bit water image，and(f)enhanced 8-bit water image

Landsat 8另一新增特征就是提高了影像的輻射分辨率。首幅影像的輻射分辨率為16 bit，其灰度量化區(qū)間為0—65536，大大超過之前8 bit數(shù)據(jù)的0—255量化區(qū)間。為了實(shí)際考察這一變化的優(yōu)勢(shì)，我們?nèi)砸杂跋駯|北角的水體進(jìn)行實(shí)驗(yàn)。除了用16 bit的數(shù)據(jù)進(jìn)行實(shí)驗(yàn)外，還將其量化為8bit，模擬成ETM+影像，以進(jìn)行對(duì)比。圖5e、f為對(duì)水域進(jìn)行對(duì)比度增強(qiáng)后獲得的16 bit和8 bit的水體影像。可以明顯發(fā)現(xiàn)，16 bit影像較好地顯出示水體的細(xì)節(jié)信息，而模擬的8 bit ETM+影像則會(huì)出現(xiàn)灰度過飽和現(xiàn)象，圖像表現(xiàn)的細(xì)節(jié)遠(yuǎn)不及前者豐富。顯然，輻射分辨率的提高可避免極暗區(qū)的灰度過飽和現(xiàn)象，這對(duì)反射率極低的水體的細(xì)微特征的辨析有很大的幫助。

3.3 裸土和不透水建筑地面監(jiān)測(cè)

新增的卷云波段(OLI 9波段)的信息量很低(表2)，這主要是因?yàn)槠洳ǘ喂庾V范圍的設(shè)置不是位于大氣窗口內(nèi)，而是位于強(qiáng)烈的水汽吸收帶(圖1)，大部分地面反射的電磁波信息無法穿過大氣層，因此它們的信息很弱，傳感器接收到的只是云層的信息。這樣，當(dāng)云和少數(shù)亮白色、反射率極強(qiáng)的地物(如沙土、新建筑用地)在其它波段不易區(qū)別時(shí)，就可以借助該波段來加以區(qū)別。從表3和圖6可以看出，云和高反射的新建筑用地在大部分波段具有相似的光譜反射曲線特征，難于區(qū)別，但在卷云波段(9波段)處，云的反射率明顯高于新建筑用地，因此可設(shè)定閾值來將其分離。表4為各主要地類在卷云波段的統(tǒng)計(jì)值，從中可以看出，云的最小值都大于任何其它地類的最大值，因此，以其作為閾值，就可以區(qū)分云和其它地類。圖7取自該景影像中一小塊有云覆蓋的地區(qū)，其中有大片的裸土，將云在卷云波段中的最小值作為閾值就可以快速地將云的信息分離出來。

表4 各地類在卷云波段的統(tǒng)計(jì)特征值Table 4 Statistics of main land cover classes in the cirrus band

圖6 云和新建筑用地的反射光譜曲線Fig.6 Signatures of cloud and new built-up land

研究還發(fā)現(xiàn)，卷云波段除了能夠突出云外，還有助于區(qū)分裸土和建筑用地。在利用光譜分解(SMA)技術(shù)提取建筑不透水面信息時(shí)遇到的一個(gè)主要難題就是土壤和不透水建筑用地的混淆問題［13-16］。但是這一問題的解決似乎可以從卷云波段找到曙光。圖7d是實(shí)驗(yàn)區(qū)卷云波段的偽彩色影像，北部的裸土和中、南部的建筑用地的色調(diào)明顯不同，裸土呈紅、黃色，說明具有很高的亮度，而建筑物則以藍(lán)綠色調(diào)為主，說明具有低亮度。從表4也可以看出，裸地在卷云波段的均值516，而建筑用地只有361，說明前者的亮度明顯高于后者，二者具有明顯的反差。

為了考察新增的卷云波段是否對(duì)區(qū)別土壤和建筑不透水面信息有幫助，我們利用Wu的線性光譜分解技術(shù)［13］對(duì)影像進(jìn)行了不透水面信息的增強(qiáng)實(shí)驗(yàn)。實(shí)驗(yàn)分為2個(gè)方案，一個(gè)是將全部多光譜波段參與實(shí)驗(yàn)(圖7e)，一個(gè)是剔除了其中的云波段(模擬為L(zhǎng)andsat ETM+)(圖7f)。從圖7e可以看出，其北部土壤以黃色為主，而中南部的建筑用地多為紅色，二者反差較大;而圖7f中，北部土壤的色調(diào)偏紅，與建筑物的反差不如圖7e來的明顯。定量分析表明，有云波段參與實(shí)驗(yàn)的不透水面的均值與土壤均值的比值為2.15∶1，而剔除云波段的實(shí)驗(yàn)的比值為1.91∶1，前者比后者的比值提高了12.6%。比值越大說明二者的反差越大，越容易區(qū)分。因此，加入卷云波段顯然有利于二者的區(qū)分。為了探究其原因，我們利用美國(guó)地質(zhì)調(diào)查局的光譜數(shù)據(jù)獲得了OLI傳感器各波段的光譜響應(yīng)函數(shù)及主要地物的反射率曲線［4］?？梢钥闯觯惩恋V物(高嶺石、蒙脫石、葉蠟石)在卷云波段附近都有一個(gè)明顯的吸收谷，而其它地物不具這一特征，多為平滑過渡(圖8)?？梢姡碓撇ǘ蔚脑O(shè)立使得粘土礦物的吸收谷特征得以突出，這顯然有助于土壤與其它地表地物的區(qū)別。而在老的Landsat系列衛(wèi)星影像中，由于缺乏卷云波段，因此土壤與其它地物一樣，在近紅外與短波紅外之間的光譜曲線同樣表現(xiàn)為平滑過渡，體現(xiàn)不出強(qiáng)吸收谷的特征，因此在TM/ETM+影像中，不透水面信息和土壤信息較難區(qū)分。

4 結(jié)論

新發(fā)射的Landsat 8衛(wèi)星有2個(gè)傳感器(OLI和TIRS)。除了保持原有Landsat 7衛(wèi)星的基本特點(diǎn)外，Landsat 8在波段設(shè)置上還新增了深藍(lán)波段和卷云波段，并在原有熱紅外波段的光譜范圍內(nèi)，設(shè)置了2個(gè)熱紅外波段。除此之外，Landsat 8還收窄了近紅外波段和全色波段的光譜范圍，將輻射分辨率從8 bit提高到12 bit。

圖7 (a)云的彩色合成影像，(b)卷云波段中的云，其它地類由于在該波段反射率太低而被抑制，(c)云的提取影像，(d)實(shí)驗(yàn)區(qū)云波段影像，(e)不透水面增強(qiáng)信息(全部波段);(f)不透水面增強(qiáng)信息(剔除卷云波段)Fig.7 (a)Cloud image，(b)cloud in the cirrus band，other features seen in(a)are all suppressed，(c)binary cloud image，(d)image of the cirrus band of the test image，(e)enhanced impervious surface image using all multispectral bands，(f)enhanced impervious surface image using all multispectral bands except the cirrus band

圖8 OLI各波段(1—7，9波段)以及主要地表地物的光譜響應(yīng)特征(據(jù)美國(guó)地質(zhì)調(diào)查局)Fig.8 Relative spectral response curves for bands of OLI，as well as main land cover classes(from USGS)

實(shí)驗(yàn)表明，新增的卷云波段可以很好地突出云的特征。除此之外，由于該波段的光譜范圍位于粘土礦物的吸收谷附近，因此可以突出土壤的特征，有助于土壤和建筑不透水面信息的區(qū)別。深藍(lán)波段的設(shè)立有助于水體懸浮物的監(jiān)測(cè)，區(qū)分渾濁水體。全色影像光譜范圍的收窄擴(kuò)大了該影像上植被和非植被信息的反差，而輻射分辨率的提高則避免了極亮/極暗區(qū)的灰度過飽和現(xiàn)象，可以更好地顯示這些地方的信息細(xì)節(jié)。顯然，Landsat 8這些新增的特點(diǎn)將會(huì)對(duì)全球生態(tài)環(huán)境變化的監(jiān)測(cè)產(chǎn)生積極的作用。致謝:本文使用的影像均下載于USGS網(wǎng)站，特此致謝。

［1］ USGS.Mission Accomplished for Landsat 5.http://www.usgs.gov/newsroom/article.asp?ID=3485，2013/4/18.

［2］ NASA.LDCM Launch.http://www.nasa.gov/mission_pages/landsat/launch/index.html.2013/4/18.

［3］ NASA.A Closer Look at LDCM's First Scene.http://www.nasa.gov/mission_pageslandsat/news/first-images-feature.html.2013/4/18.

［4］ USGS.Landsat Missions.http://landsat.usgs.gov/index.php.2013/4/18.

［5］ USGS.Landsat Data Continuity Mission(LDCM)-Landsat 8.http://landsat.usgs.gov/LDCM_Landsat8.php.2013/4/18.

［6］ Irons J R，Dwyer J L，Barsi J A.The next Landsat satellite:The Landsat Data Continuity Mission.Remote Sensing of Environment，2012，122:11-21.

［7］ Ungar S G，Pearlman J S，Mendenhall J A，Reuter D.Overview of the earth observing one(EO-1)mission.IEEE Transactions on Geoscience and Remote Sensing，2003，41(6):1149-1159.

［8］ Chander G，Markham B L，Helder D L.Summary of current radiometric calibration coefficients for Landsat MSS，TM，ETM+，and EO-1 ALI sensors.Remote Sensing of Environment，2009，113(5):893-903.

［9］ Jensen J R.Introductory Digital Image Processing.3rd ed.Upper Saddle River，NJ:Prentice Hall，2004:526-526.

［10］ Sheila M S.Geology of the Boulder Area.http://bcn.boulder.co.us/basin/watershed/geology/.2013/4/18

［11］ Xu H Q.Modification of normalised difference water index(NDWI)to enhance open water features in remotely sensed imagery.International Journal of Remote Sensing，2006，27(14):3025-3033.

［12］ Xu H Q.Water colour variation analysis of the coastal waters surrounding Xiamen Island of SE China by multispectral and multitemporal remote sensing measurements.Acta Scientiae Circumstantiae，2006，26(7):1209-1218.

［13］ Wu C S，Murray A T.Estimating impervious surface distribution by spectral mixture analysis.Remote Sensing of Environment，2003，84(4):493-505.

［14］ Powell S L，Cohen W B，Yang Z Q，Pierce J D，Alberti M.Quantification of impervious surface in the Snohomish Water Resources Inventory Area of western Washington from 1972—2006.Remote Sensing of Environment，2008，112(14):1895-1908.

［15］ Gluch R，Quattrochi D A，Luvall J C.A multi-scale approach to urban thermal analysis.Remote Sensing of Environment，2006，104(2):123-132.

［16］ Yang X J，Liu Z.Use of satellite-derived landscape imperviousness index to characterize urban spatial growth.Computers，Environment and Urban Systems，2005，29(5):524-540.

參考文獻(xiàn):

［12］ 徐涵秋.環(huán)廈門海域水色變化的多光譜多時(shí)相遙感分析.環(huán)境科學(xué)學(xué)報(bào)，2006，26(7):1209-1218.