吳國(guó)璽++申懷飛++張鵬飛
摘要:以河南省許昌市冬小麥為研究對(duì)象,選取許昌市東城區(qū)北郊麥田,采集冬小麥光譜數(shù)據(jù)。利用光譜儀后處理軟件ASD View Spec Pro對(duì)采集的光譜數(shù)據(jù)以DAT的格式導(dǎo)出并在Microsoft Office Excel 2003中生成曲線,剔除其中無(wú)效數(shù)據(jù),通過對(duì)地物多條光譜數(shù)據(jù)求平均,對(duì)實(shí)測(cè)波譜數(shù)據(jù)歸類,計(jì)算冬小麥冠層歸一化植被指數(shù),得到冬小麥冠層光譜特征曲線,利用冬小麥冠層光譜特征曲線分析河南省冬小麥冠層光譜的環(huán)境影響因子。結(jié)果表明,隨著太陽(yáng)高度角的增加冬小麥冠層光譜的反射率增加,太陽(yáng)光接近直射時(shí)小麥冠層所接受到的太陽(yáng)輻射量最大;冬小麥種植密度、滯塵量等因素對(duì)冬小麥冠層光譜曲線均有影響。
關(guān)鍵詞:冬小麥;冠層光譜曲線;環(huán)境因子;河南省
中圖分類號(hào):S512.1+1;TP79 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):0439-8114(2014)02-0268-04
Environmental Impact Factors of Winter Wheat Canopy Reflectance in Henan Province
WU Guo-xi,SHEN Huai-fei,ZHANG Peng-fei
(School of Urban & Rural Planning and Landscape Architecture, Xuchang University, Xuchang 461000,Henan,China)
Abstract: The winter wheat spectral data of portable spectroradiometer made by ASD company in the Dong-cheng district of Xuchang city was collected. ASD View Spec Pro software and Microsoft Office Excel were used to process data. The average of the surface features multi-spectral data was classified. The winter wheat canopy Normalized Difference Vegetation Index was calculated. The winter canopy spectral characteristic curves in different times and different conditions were made. By comparatively analysing the curves of the winter wheat canopy, the environmental factors of winter wheat canopy reflectance at growing period in henan province were analyzed. The results showed that when the sun elevation angle increased, the wheat spectral reflectance increased. When the sunlight is direct on the wheat, the received radiation is the biggest. The wheat plating density and dust in the air had great impact on the winter canopy spectral curves.
Key words: winter wheat; spectral curves of the canopy; environmental factors; Henan province
20世紀(jì)80年代,國(guó)外開展了對(duì)植被冠層的二向性反射特性研究,觀測(cè)了多種農(nóng)作物、草地及森林冠層的雙向反射分布函數(shù)(Bidirectional reflection distribution function,BRDF),發(fā)現(xiàn)了多角度反射率的“熱點(diǎn)效應(yīng)”和鏡面反射等典型特征[1]。此后,中國(guó)也開展了對(duì)冬小麥冠層的光譜分析。張雪紅等[2,3]于2006年初步分析了冬小麥冠層反射率光譜及紅邊參數(shù)特征,發(fā)現(xiàn)存在較強(qiáng)的各向異性特征;2010年又通過采用各向異性因子和各向異性指數(shù)定量分析了冬小麥冠層窄波段二向性反射率及NDVI的方向性特征。趙春江等[4]用多角度光譜信息反演冬小麥葉綠素含量的垂直分布。袁慧等[5]分析了不同發(fā)育階段、播種壟向及土壤背景對(duì)冠層光譜反射的影響。王靜等[6]研究了甘肅省慶陽(yáng)市西峰區(qū)冬小麥主要生育期冠層反射光譜與反映長(zhǎng)勢(shì)的指標(biāo)覆蓋度及葉面積指數(shù)間的關(guān)系。宋曉宇等[7]利用遙感影像數(shù)據(jù)對(duì)冬小麥長(zhǎng)勢(shì)空間變異進(jìn)行了研究。馮偉等[8]利用高光譜遙感對(duì)小麥冠層葉片色素密度進(jìn)行監(jiān)測(cè),研究了小麥葉片色素密度與冠層高光譜參數(shù)的定量關(guān)系。以上文獻(xiàn)多數(shù)是對(duì)冬小麥某一生長(zhǎng)階段不同生長(zhǎng)條件下光譜反射率、光譜動(dòng)態(tài)變化及紅邊參數(shù)特征的研究,對(duì)不同發(fā)育階段、播種壟向及土壤背景與冠層光譜反射關(guān)系的探討;且大多以研究小麥光譜特征為手段,以研究其理化特性為目的,或者用多角度光譜信息反演冬小麥葉綠素含量垂直分布。而對(duì)隨著太陽(yáng)高度角的增加小麥冠層光譜的反射率狀況,冬小麥種植密度、滯塵量等因素對(duì)冬小麥冠層光譜影響的研究較少。因此,本研究以河南省許昌市東城區(qū)為例,通過采集有效光譜數(shù)據(jù),并對(duì)這些光譜數(shù)據(jù)進(jìn)行分析,以探討外界因素對(duì)小麥冠層光譜的影響,對(duì)保證河南省冬小麥正常生長(zhǎng)有著重要的意義。
1 研究區(qū)概況
2 數(shù)據(jù)來(lái)源及處理
2.1 冬小麥光譜數(shù)據(jù)的獲取
2.2 光譜數(shù)據(jù)的預(yù)處理
利用光譜儀后處理軟件ASD View Spec Pro對(duì)采集的光譜數(shù)據(jù)以DAT的格式導(dǎo)出并在Microsoft Office Excel 2003中生成曲線,剔除其中無(wú)效數(shù)據(jù),對(duì)地物多條光譜數(shù)據(jù)求平均值,對(duì)實(shí)測(cè)波譜數(shù)據(jù)歸類合并,最后計(jì)算出冬小麥冠層的歸一化植被指數(shù)。從每次采集的15個(gè)數(shù)據(jù)中篩選3個(gè)作為樣本數(shù)據(jù),對(duì)樣本光譜數(shù)據(jù)求均值,所得曲線即為冬小麥冠層光譜曲線。利用預(yù)處理后的冬小麥冠層光譜特征數(shù)據(jù)比較同一生長(zhǎng)期同一因素不同情況下冬小麥冠層光譜曲線,找出單一因素(太陽(yáng)高度角、樣地種植密度和冬小麥冠層滯塵量)對(duì)冬小麥冠層光譜曲線影響的相應(yīng)波段。
3 各因素對(duì)冬小麥冠層光譜曲線的影響
3.1 太陽(yáng)高度角對(duì)冬小麥冠層光譜曲線的影響
在1 d內(nèi)對(duì)10個(gè)不同時(shí)間段的冬小麥冠層光譜進(jìn)行采集并繪制光譜曲線圖,結(jié)果如圖3所示。從圖3可知,冬小麥冠層光譜曲線的大致趨勢(shì)為:從早上7:20隨著太陽(yáng)高度角的增長(zhǎng),其各個(gè)波段的反射率均呈上升趨勢(shì),在11:20太陽(yáng)高度角最大時(shí)反射率達(dá)到最大,接著隨太陽(yáng)高度角的減小,冬小麥各個(gè)波段的反射率呈下降趨勢(shì)。太陽(yáng)高度角直接影響小麥光譜的反射率,當(dāng)太陽(yáng)越趨近于直射的時(shí)候,小麥冠層所接受到的太陽(yáng)輻射量就越大,反射率就越大,由此可以得出在11:00左右采集小麥冠層光譜數(shù)據(jù)為最佳時(shí)段。
3.2 滯塵量對(duì)冬小麥冠層光譜曲線的影響
在距離公路5、10、20 m處選取A、B、C 3塊小麥田,分別采集11:20和12:20時(shí)的冬小麥冠層光譜數(shù)據(jù),經(jīng)過數(shù)據(jù)處理后,得到小麥冠層光譜曲線如圖4所示。距離公路遠(yuǎn)近不同的3塊麥田的滯塵量不同,其冬小麥的冠層光譜曲線也不同,在排除操作誤差的影響下,在波峰處A地反射率始終為最大,而B、C兩地的區(qū)別較小,且均小于A地的反射率,未能將距公路距離即滯塵量作為單一變量分離出來(lái),具體原因在于此3塊麥田的冬小麥冠層光譜還受到了灌溉、施肥等其他因素的影響,造成未能發(fā)現(xiàn)滯塵量對(duì)冬小麥冠層光譜的影響。
3.3 種植密度對(duì)冬小麥冠層光譜曲線的影響
4 結(jié)論
通過野外地物光譜儀對(duì)試驗(yàn)區(qū)的冬小麥測(cè)量,獲得冬小麥的冠層光譜反射率,分析冬小麥生長(zhǎng)期的光譜特征,并描繪出了冬小麥的冠層光譜曲線。將實(shí)測(cè)小麥的光譜數(shù)據(jù)進(jìn)行了定性和定量的分析,總結(jié)出外界因素對(duì)于冬小麥光譜曲線的影響。
1)從早上7:20隨著太陽(yáng)高度角的增大,冬小麥各個(gè)波段的反射率均呈上升趨勢(shì),在11:20太陽(yáng)高度角最大,反射率也達(dá)到最大,之后隨太陽(yáng)高度角的減小,各個(gè)波段的反射率呈下降趨勢(shì)。
2)冬小麥冠層滯塵量和種植密度對(duì)冬小麥冠層光譜曲線均有影響,但是由于條件限制和外界因素影響,未能準(zhǔn)確地發(fā)現(xiàn)這2種因素對(duì)于冬小麥冠層光譜曲線的影響特征。
在研究過程中,由于時(shí)間的限制而未能將冬小麥生長(zhǎng)期的各個(gè)階段的光譜采集完整,造成部分?jǐn)?shù)據(jù)缺失;另外在研究單一因素對(duì)冬小麥冠層光譜特征的影響時(shí)無(wú)法完全排除其他因素的影響,只能相對(duì)減少其他因素的影響,這些問題均可能會(huì)造成分析結(jié)果的誤差,尚需進(jìn)一步研究。
參考文獻(xiàn):
[1] ROSS J K, MARSHAK A L. Calculation of canopy bidirectional reflectance using the Monte Carlo method[J]. Remote Sending Environ,1988,24:213-225.
[2] 張雪紅,趙 峰,劉紹民,等.冬小麥紅邊參數(shù)各向異性特征分析[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22(6):7-11.
[3] 張雪紅,田慶久,沈潤(rùn)平.冬小麥冠層光譜的方向性特征分析[J].光譜學(xué)與光譜分析,2010,30(6):1600-1605.
[4] 趙春江,黃文江,王紀(jì)華,等.用多角度光譜信息反演冬小麥葉綠素含量垂直分布[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22(6):104-109.
[5] 袁 慧,李 爽,孔云峰.開封市郊冬小麥冠層反射光譜特征分析[J].地理空間信息,2009,7(4):63-66.
[6] 王 靜,郭 鈮,王小平,等.半濕潤(rùn)雨養(yǎng)農(nóng)業(yè)區(qū)冬小麥冠層反射光譜與長(zhǎng)勢(shì)的相關(guān)性分析[J].資源科學(xué),2008,30(8):1261-1267.
[7] 宋曉宇,王紀(jì)華,閻廣建,等.基于多時(shí)相航空高光譜遙感影像的冬小麥長(zhǎng)勢(shì)空間變異研究[J].光譜學(xué)與光譜分析,2010,30(7):1820-1824.
[8] 馮 偉,朱 艷,田永超,等.基于高光譜遙感的小麥冠層葉片色素密度監(jiān)測(cè)[J].生態(tài)學(xué)報(bào),2008,28(10):4902-4911.
(責(zé)任編輯 呂海霞)
2.1 冬小麥光譜數(shù)據(jù)的獲取
2.2 光譜數(shù)據(jù)的預(yù)處理
利用光譜儀后處理軟件ASD View Spec Pro對(duì)采集的光譜數(shù)據(jù)以DAT的格式導(dǎo)出并在Microsoft Office Excel 2003中生成曲線,剔除其中無(wú)效數(shù)據(jù),對(duì)地物多條光譜數(shù)據(jù)求平均值,對(duì)實(shí)測(cè)波譜數(shù)據(jù)歸類合并,最后計(jì)算出冬小麥冠層的歸一化植被指數(shù)。從每次采集的15個(gè)數(shù)據(jù)中篩選3個(gè)作為樣本數(shù)據(jù),對(duì)樣本光譜數(shù)據(jù)求均值,所得曲線即為冬小麥冠層光譜曲線。利用預(yù)處理后的冬小麥冠層光譜特征數(shù)據(jù)比較同一生長(zhǎng)期同一因素不同情況下冬小麥冠層光譜曲線,找出單一因素(太陽(yáng)高度角、樣地種植密度和冬小麥冠層滯塵量)對(duì)冬小麥冠層光譜曲線影響的相應(yīng)波段。
3 各因素對(duì)冬小麥冠層光譜曲線的影響
3.1 太陽(yáng)高度角對(duì)冬小麥冠層光譜曲線的影響
在1 d內(nèi)對(duì)10個(gè)不同時(shí)間段的冬小麥冠層光譜進(jìn)行采集并繪制光譜曲線圖,結(jié)果如圖3所示。從圖3可知,冬小麥冠層光譜曲線的大致趨勢(shì)為:從早上7:20隨著太陽(yáng)高度角的增長(zhǎng),其各個(gè)波段的反射率均呈上升趨勢(shì),在11:20太陽(yáng)高度角最大時(shí)反射率達(dá)到最大,接著隨太陽(yáng)高度角的減小,冬小麥各個(gè)波段的反射率呈下降趨勢(shì)。太陽(yáng)高度角直接影響小麥光譜的反射率,當(dāng)太陽(yáng)越趨近于直射的時(shí)候,小麥冠層所接受到的太陽(yáng)輻射量就越大,反射率就越大,由此可以得出在11:00左右采集小麥冠層光譜數(shù)據(jù)為最佳時(shí)段。
3.2 滯塵量對(duì)冬小麥冠層光譜曲線的影響
在距離公路5、10、20 m處選取A、B、C 3塊小麥田,分別采集11:20和12:20時(shí)的冬小麥冠層光譜數(shù)據(jù),經(jīng)過數(shù)據(jù)處理后,得到小麥冠層光譜曲線如圖4所示。距離公路遠(yuǎn)近不同的3塊麥田的滯塵量不同,其冬小麥的冠層光譜曲線也不同,在排除操作誤差的影響下,在波峰處A地反射率始終為最大,而B、C兩地的區(qū)別較小,且均小于A地的反射率,未能將距公路距離即滯塵量作為單一變量分離出來(lái),具體原因在于此3塊麥田的冬小麥冠層光譜還受到了灌溉、施肥等其他因素的影響,造成未能發(fā)現(xiàn)滯塵量對(duì)冬小麥冠層光譜的影響。
3.3 種植密度對(duì)冬小麥冠層光譜曲線的影響
4 結(jié)論
通過野外地物光譜儀對(duì)試驗(yàn)區(qū)的冬小麥測(cè)量,獲得冬小麥的冠層光譜反射率,分析冬小麥生長(zhǎng)期的光譜特征,并描繪出了冬小麥的冠層光譜曲線。將實(shí)測(cè)小麥的光譜數(shù)據(jù)進(jìn)行了定性和定量的分析,總結(jié)出外界因素對(duì)于冬小麥光譜曲線的影響。
1)從早上7:20隨著太陽(yáng)高度角的增大,冬小麥各個(gè)波段的反射率均呈上升趨勢(shì),在11:20太陽(yáng)高度角最大,反射率也達(dá)到最大,之后隨太陽(yáng)高度角的減小,各個(gè)波段的反射率呈下降趨勢(shì)。
2)冬小麥冠層滯塵量和種植密度對(duì)冬小麥冠層光譜曲線均有影響,但是由于條件限制和外界因素影響,未能準(zhǔn)確地發(fā)現(xiàn)這2種因素對(duì)于冬小麥冠層光譜曲線的影響特征。
在研究過程中,由于時(shí)間的限制而未能將冬小麥生長(zhǎng)期的各個(gè)階段的光譜采集完整,造成部分?jǐn)?shù)據(jù)缺失;另外在研究單一因素對(duì)冬小麥冠層光譜特征的影響時(shí)無(wú)法完全排除其他因素的影響,只能相對(duì)減少其他因素的影響,這些問題均可能會(huì)造成分析結(jié)果的誤差,尚需進(jìn)一步研究。
參考文獻(xiàn):
[1] ROSS J K, MARSHAK A L. Calculation of canopy bidirectional reflectance using the Monte Carlo method[J]. Remote Sending Environ,1988,24:213-225.
[2] 張雪紅,趙 峰,劉紹民,等.冬小麥紅邊參數(shù)各向異性特征分析[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22(6):7-11.
[3] 張雪紅,田慶久,沈潤(rùn)平.冬小麥冠層光譜的方向性特征分析[J].光譜學(xué)與光譜分析,2010,30(6):1600-1605.
[4] 趙春江,黃文江,王紀(jì)華,等.用多角度光譜信息反演冬小麥葉綠素含量垂直分布[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22(6):104-109.
[5] 袁 慧,李 爽,孔云峰.開封市郊冬小麥冠層反射光譜特征分析[J].地理空間信息,2009,7(4):63-66.
[6] 王 靜,郭 鈮,王小平,等.半濕潤(rùn)雨養(yǎng)農(nóng)業(yè)區(qū)冬小麥冠層反射光譜與長(zhǎng)勢(shì)的相關(guān)性分析[J].資源科學(xué),2008,30(8):1261-1267.
[7] 宋曉宇,王紀(jì)華,閻廣建,等.基于多時(shí)相航空高光譜遙感影像的冬小麥長(zhǎng)勢(shì)空間變異研究[J].光譜學(xué)與光譜分析,2010,30(7):1820-1824.
[8] 馮 偉,朱 艷,田永超,等.基于高光譜遙感的小麥冠層葉片色素密度監(jiān)測(cè)[J].生態(tài)學(xué)報(bào),2008,28(10):4902-4911.
(責(zé)任編輯 呂海霞)
2.1 冬小麥光譜數(shù)據(jù)的獲取
2.2 光譜數(shù)據(jù)的預(yù)處理
利用光譜儀后處理軟件ASD View Spec Pro對(duì)采集的光譜數(shù)據(jù)以DAT的格式導(dǎo)出并在Microsoft Office Excel 2003中生成曲線,剔除其中無(wú)效數(shù)據(jù),對(duì)地物多條光譜數(shù)據(jù)求平均值,對(duì)實(shí)測(cè)波譜數(shù)據(jù)歸類合并,最后計(jì)算出冬小麥冠層的歸一化植被指數(shù)。從每次采集的15個(gè)數(shù)據(jù)中篩選3個(gè)作為樣本數(shù)據(jù),對(duì)樣本光譜數(shù)據(jù)求均值,所得曲線即為冬小麥冠層光譜曲線。利用預(yù)處理后的冬小麥冠層光譜特征數(shù)據(jù)比較同一生長(zhǎng)期同一因素不同情況下冬小麥冠層光譜曲線,找出單一因素(太陽(yáng)高度角、樣地種植密度和冬小麥冠層滯塵量)對(duì)冬小麥冠層光譜曲線影響的相應(yīng)波段。
3 各因素對(duì)冬小麥冠層光譜曲線的影響
3.1 太陽(yáng)高度角對(duì)冬小麥冠層光譜曲線的影響
在1 d內(nèi)對(duì)10個(gè)不同時(shí)間段的冬小麥冠層光譜進(jìn)行采集并繪制光譜曲線圖,結(jié)果如圖3所示。從圖3可知,冬小麥冠層光譜曲線的大致趨勢(shì)為:從早上7:20隨著太陽(yáng)高度角的增長(zhǎng),其各個(gè)波段的反射率均呈上升趨勢(shì),在11:20太陽(yáng)高度角最大時(shí)反射率達(dá)到最大,接著隨太陽(yáng)高度角的減小,冬小麥各個(gè)波段的反射率呈下降趨勢(shì)。太陽(yáng)高度角直接影響小麥光譜的反射率,當(dāng)太陽(yáng)越趨近于直射的時(shí)候,小麥冠層所接受到的太陽(yáng)輻射量就越大,反射率就越大,由此可以得出在11:00左右采集小麥冠層光譜數(shù)據(jù)為最佳時(shí)段。
3.2 滯塵量對(duì)冬小麥冠層光譜曲線的影響
在距離公路5、10、20 m處選取A、B、C 3塊小麥田,分別采集11:20和12:20時(shí)的冬小麥冠層光譜數(shù)據(jù),經(jīng)過數(shù)據(jù)處理后,得到小麥冠層光譜曲線如圖4所示。距離公路遠(yuǎn)近不同的3塊麥田的滯塵量不同,其冬小麥的冠層光譜曲線也不同,在排除操作誤差的影響下,在波峰處A地反射率始終為最大,而B、C兩地的區(qū)別較小,且均小于A地的反射率,未能將距公路距離即滯塵量作為單一變量分離出來(lái),具體原因在于此3塊麥田的冬小麥冠層光譜還受到了灌溉、施肥等其他因素的影響,造成未能發(fā)現(xiàn)滯塵量對(duì)冬小麥冠層光譜的影響。
3.3 種植密度對(duì)冬小麥冠層光譜曲線的影響
4 結(jié)論
通過野外地物光譜儀對(duì)試驗(yàn)區(qū)的冬小麥測(cè)量,獲得冬小麥的冠層光譜反射率,分析冬小麥生長(zhǎng)期的光譜特征,并描繪出了冬小麥的冠層光譜曲線。將實(shí)測(cè)小麥的光譜數(shù)據(jù)進(jìn)行了定性和定量的分析,總結(jié)出外界因素對(duì)于冬小麥光譜曲線的影響。
1)從早上7:20隨著太陽(yáng)高度角的增大,冬小麥各個(gè)波段的反射率均呈上升趨勢(shì),在11:20太陽(yáng)高度角最大,反射率也達(dá)到最大,之后隨太陽(yáng)高度角的減小,各個(gè)波段的反射率呈下降趨勢(shì)。
2)冬小麥冠層滯塵量和種植密度對(duì)冬小麥冠層光譜曲線均有影響,但是由于條件限制和外界因素影響,未能準(zhǔn)確地發(fā)現(xiàn)這2種因素對(duì)于冬小麥冠層光譜曲線的影響特征。
在研究過程中,由于時(shí)間的限制而未能將冬小麥生長(zhǎng)期的各個(gè)階段的光譜采集完整,造成部分?jǐn)?shù)據(jù)缺失;另外在研究單一因素對(duì)冬小麥冠層光譜特征的影響時(shí)無(wú)法完全排除其他因素的影響,只能相對(duì)減少其他因素的影響,這些問題均可能會(huì)造成分析結(jié)果的誤差,尚需進(jìn)一步研究。
參考文獻(xiàn):
[1] ROSS J K, MARSHAK A L. Calculation of canopy bidirectional reflectance using the Monte Carlo method[J]. Remote Sending Environ,1988,24:213-225.
[2] 張雪紅,趙 峰,劉紹民,等.冬小麥紅邊參數(shù)各向異性特征分析[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22(6):7-11.
[3] 張雪紅,田慶久,沈潤(rùn)平.冬小麥冠層光譜的方向性特征分析[J].光譜學(xué)與光譜分析,2010,30(6):1600-1605.
[4] 趙春江,黃文江,王紀(jì)華,等.用多角度光譜信息反演冬小麥葉綠素含量垂直分布[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22(6):104-109.
[5] 袁 慧,李 爽,孔云峰.開封市郊冬小麥冠層反射光譜特征分析[J].地理空間信息,2009,7(4):63-66.
[6] 王 靜,郭 鈮,王小平,等.半濕潤(rùn)雨養(yǎng)農(nóng)業(yè)區(qū)冬小麥冠層反射光譜與長(zhǎng)勢(shì)的相關(guān)性分析[J].資源科學(xué),2008,30(8):1261-1267.
[7] 宋曉宇,王紀(jì)華,閻廣建,等.基于多時(shí)相航空高光譜遙感影像的冬小麥長(zhǎng)勢(shì)空間變異研究[J].光譜學(xué)與光譜分析,2010,30(7):1820-1824.
[8] 馮 偉,朱 艷,田永超,等.基于高光譜遙感的小麥冠層葉片色素密度監(jiān)測(cè)[J].生態(tài)學(xué)報(bào),2008,28(10):4902-4911.
(責(zé)任編輯 呂海霞)