胡順石，張立福，張 霞，王 倩，韓 冰，張 楠

衛(wèi)星傳感器波段平均太陽(yáng)輻照度計(jì)算及可靠性分析

胡順石1，張立福1，張 霞1，王 倩1，韓 冰2，張 楠3

(1.中國(guó)科學(xué)院遙感應(yīng)用研究所遙感科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100101;2.中國(guó)礦業(yè)大學(xué)，徐州 221116;3.中國(guó)地質(zhì)大學(xué)，北京 100083)

地球大氣層外太陽(yáng)光譜輻照度(extraterrestrial solar spectral irradiance，ESSI)數(shù)據(jù)是計(jì)算衛(wèi)星傳感器波段平均太陽(yáng)輻照度(band mean solar irradiance，BMSI)的重要參數(shù)。為了探求利用何種來(lái)源的ESSI數(shù)據(jù)計(jì)算傳感器BMSI更為準(zhǔn)確，分別采用SBDART軟件模擬的太陽(yáng)光譜曲線數(shù)據(jù)、MODTRAN4.0 oldkur.dat文件數(shù)據(jù)、Thuillier太陽(yáng)光譜曲線數(shù)據(jù)和WRC太陽(yáng)光譜曲線數(shù)據(jù)計(jì)算了HJ－1A CCD1(B1—B4)，CBERS－02 CCD(B1—B5)，Landsat5 TM(B1—B4)和ASTER(B1—B8)4種傳感器的BMSI，并與傳感器運(yùn)營(yíng)商公布的數(shù)據(jù)進(jìn)行了比較。結(jié)果表明:利用SBDART和WRC太陽(yáng)光譜曲線數(shù)據(jù)計(jì)算的結(jié)果誤差較小;利用MODTRAN4.0 oldkur.dat數(shù)據(jù)計(jì)算的結(jié)果誤差次之;利用Thuillier太陽(yáng)光譜曲線的計(jì)算結(jié)果誤差較大。

大氣層外太陽(yáng)光譜輻照度(ESSI);波段平均太陽(yáng)輻照度(BMSI);Thuillier;SBDART;WRC;oldkur.dat

0 引言

太陽(yáng)是一個(gè)穩(wěn)定的電磁波輻射源，是進(jìn)行光學(xué)遙感的主要能量來(lái)源。地球大氣層外太陽(yáng)光譜輻照度(extraterrestrial solar spectral irradiance，ESSI)是指在大氣層頂部、日地平均距離上，垂直于太陽(yáng)入射方向上單位時(shí)間、單位面積所獲得的太陽(yáng)輻射能量。太陽(yáng)輻射能量的大部分集中在近紫外—中紅外波長(zhǎng)范圍(0.31 ～5.6 μm)內(nèi)，這部分能量占整個(gè)太陽(yáng)光譜能量的97.5%，其中可見(jiàn)光部分占43.5%，近紅外部分占36.8%［1］。傳感器波段平均太陽(yáng)輻照度(band mean solar irradiance，BMSI)是傳感器光譜響應(yīng)函數(shù)與大氣層外太陽(yáng)光譜輻照度的積分值，反映傳感器對(duì)太陽(yáng)輻射能量的響應(yīng)性能。

在大氣頂層，太陽(yáng)輻射的電磁波能量不受大氣的影響，它是計(jì)算各種輻射參量和定量遙感的基礎(chǔ)。國(guó)內(nèi)外許多學(xué)者對(duì)此進(jìn)行了研究［2－4］。SBDART(santa barbara disort atmospheric radiative transfer)軟件能模擬地球大氣層外太陽(yáng)光譜曲線數(shù)據(jù)［2］，潘志強(qiáng)等用該數(shù)據(jù)反演了CBERS衛(wèi)星BMSI，并得到了較好的效果［3］;MODTRAN4.0 軟件提供了 5 個(gè)ESSI數(shù)據(jù)文件;Thuillier等通過(guò)航空實(shí)驗(yàn)發(fā)布了200～4 000 nm波長(zhǎng)范圍的太陽(yáng)光譜曲線［4］;Alistair利用WRC(world radiation center)發(fā)布的ESSI數(shù)據(jù)計(jì)算了 ASTER 各波段 BMSI［5］。

本文選擇SBDART軟件模擬的大氣頂層太陽(yáng)光譜曲線數(shù)據(jù)、MODTRAN4.0提供的 oldkur.dat文件數(shù)據(jù)、Thuillier太陽(yáng)光譜曲線數(shù)據(jù)和WRC太陽(yáng)光譜曲線數(shù)據(jù)，分別計(jì)算了國(guó)內(nèi)外一些常見(jiàn)的傳感器BMSI值，并以各傳感器運(yùn)營(yíng)商提供的BMSI數(shù)據(jù)為標(biāo)準(zhǔn)，與計(jì)算結(jié)果進(jìn)行比較，最終得出適合計(jì)算衛(wèi)星傳感器BMSI的太陽(yáng)光譜曲線數(shù)據(jù)。

1 大氣層外太陽(yáng)光譜曲線

1.1 SBDART太陽(yáng)光譜曲線

SBDART軟件由加利福尼亞圣巴巴拉大學(xué)計(jì)算地球系統(tǒng)科學(xué)研究所Paul等人采用Fortran語(yǔ)言開(kāi)發(fā)的，主要用于太陽(yáng)輻射傳輸模擬和大氣能量平衡研究［2，6］。用戶設(shè)置相關(guān)太陽(yáng)幾何位置參數(shù)、大氣剖面參數(shù)等，可以模擬0.25～100 μm波長(zhǎng)范圍內(nèi)大氣頂層和地表的上行、下行及直射等不同輻射參量。本文設(shè)置太陽(yáng)天頂角為0°，不考慮氣溶膠和云含量，以1 nm為步長(zhǎng)，利用SABDART在線工具［7］模擬0.25～4.0 μm波長(zhǎng)范圍內(nèi)的大氣頂層太陽(yáng)光譜曲線數(shù)據(jù)，其單位是W·m－2·μm－1。

1.2 oldkur.dat太陽(yáng)光譜曲線

MODTRAN4.0軟件“BinData”目錄下有 5個(gè)ESSI數(shù)據(jù)文件(oldkur.dat，chkur.dat，cebchkur.dat，thkur.dat和 newkur.dat)。經(jīng)過(guò)分析和比較發(fā)現(xiàn)，利用oldkur.dat數(shù)據(jù)計(jì)算的BMSI與傳感器運(yùn)營(yíng)商公布的數(shù)據(jù)較接近。因此，本文選擇oldkur.dat數(shù)據(jù)來(lái)計(jì)算結(jié)果和比較，其光譜范圍是51～49 983 cm－1，光譜分辨率達(dá)1 cm－1，太陽(yáng)光譜輻照度單位是W·cm－2·cm－1。

1.3 Thuillier太陽(yáng)光譜曲線

Thuillier太陽(yáng)光譜曲線是 Thuillier等人利用SOLSPEC和SOSP儀器分別在1992年和1993年多次航空飛行實(shí)驗(yàn)中所測(cè)得的ESSI數(shù)據(jù)，經(jīng)儀器定標(biāo)和數(shù)據(jù)處理發(fā)布的太陽(yáng)光譜曲線［4］。從網(wǎng)站［8］下載的Thuillier太陽(yáng)光譜曲線的光譜范圍為200～2 397 nm，光譜分辨率為1 nm，太陽(yáng)光譜輻照度的單位是W·cm－2·nm－1。

1.4 WRC太陽(yáng)光譜曲線

WRC太陽(yáng)光譜曲線是world radiation center研究機(jī)構(gòu)基于地面測(cè)量和火箭飛行數(shù)據(jù)發(fā)布的太陽(yáng)光譜曲線。Alistair利用此曲線計(jì)算了ASTER傳感器BMSI。本文從網(wǎng)站［9］下載了WRC太陽(yáng)光譜曲線，其光譜范圍為0.35 ～2.50 μm，光譜分辨率為1 nm，太陽(yáng)光譜輻照度單位是W·m－2·μm－1。

由于上述4種太陽(yáng)光譜曲線的太陽(yáng)光譜輻照度單位、光譜分辨率和光譜范圍不一致，為便于比較分析，首先要進(jìn)行單位轉(zhuǎn)換(波長(zhǎng)單位轉(zhuǎn)換成μm，光譜分辨率轉(zhuǎn)換成nm，輻照度單位轉(zhuǎn)換成W·m－2·μm－1)，然后截取0.35～2.397 μm波長(zhǎng)范圍內(nèi)的4種太陽(yáng)光譜曲線，其結(jié)果如圖1所示。

圖1 4種太陽(yáng)光譜曲線Fig.1 Four kinds of solar spectrum curves

由于計(jì)算ASTER傳感器BMSI采用的是WRC太陽(yáng)光譜曲線數(shù)據(jù)，因此，本文以WRC公布的數(shù)據(jù)為基準(zhǔn)，將 oldkur.dat，Thuillier和 SBDART 太陽(yáng)光譜曲線數(shù)據(jù)與其進(jìn)行比較，并求差值，結(jié)果如圖2所示。

圖 2 oldkur.dat，Thuillier，SBDART 與 WRC 的差值圖Fig.2 Respective difference between WRC and oldkur.dat，Thuillier，SBDART

從圖1和圖2可以看出，4種太陽(yáng)光譜曲線整體趨勢(shì)一致，但在350～2 397 nm波長(zhǎng)范圍內(nèi)相互之間存在一定的差異，這些差異會(huì)影響計(jì)算結(jié)果。

2 衛(wèi)星傳感器光譜響應(yīng)

衛(wèi)星傳感器的光譜響應(yīng)函數(shù)(spectral response function，SRF)描述了傳感器對(duì)太陽(yáng)入射能量的響應(yīng)率，表示傳感器在某一波長(zhǎng)λ上記錄的輻射能量與入瞳處的輻射能量的比值，它是衛(wèi)星傳感器的一個(gè)重要性能指標(biāo)。計(jì)算傳感器BMSI需要傳感器各波段的光譜響應(yīng)函數(shù)，不同類型傳感器及同一傳感器不同波段的光譜響應(yīng)函數(shù)是不同的。為了分析所選擇的太陽(yáng)光譜曲線數(shù)據(jù)的可靠性，本文選擇了HJ－1A CCD1(B1—B4)，CBERS － 02 CCD(B1—B5)，Landsat5 TM(B1—B4)和 Terra ASTER(B1—B8)4種衛(wèi)星傳感器及其波段數(shù)據(jù)。4種衛(wèi)星傳感器各波段的光譜響應(yīng)函數(shù)曲線如圖3所示。

圖3 4種衛(wèi)星傳感器各波段的光譜響應(yīng)函數(shù)Fig.3 Spectral response functions of each band for four satellites

3 波段平均太陽(yáng)輻照度的計(jì)算

傳感器波段平均太陽(yáng)輻照度EBMSI是ESSI與傳感器各波段光譜響應(yīng)函數(shù)λ的積分值。計(jì)算公式為

式中:E(λ)是大氣層外太陽(yáng)光譜輻射能量;S(λ)是傳感器某一波段光譜響應(yīng)函數(shù);λ1和λ2是傳感器某一波段起始波長(zhǎng)和終止波長(zhǎng)。因此，若已知E(λ)和 S(λ)，便可計(jì)算出 EBMSI。

一般情況下，E(λ)可根據(jù)相關(guān)研究機(jī)構(gòu)提供的ESSI數(shù)據(jù)文件或者利用輻射傳輸軟件模擬得到;S(λ)由衛(wèi)星傳感器廠商經(jīng)過(guò)實(shí)驗(yàn)室定標(biāo)后提供。在本文的研究中，E(λ)分別采用上文提及的4種太陽(yáng)光譜曲線數(shù)據(jù);S(λ)分別采用上文提及的4類衛(wèi)星傳感器的光譜響應(yīng)函數(shù)。波段平均太陽(yáng)輻照度EBMSI的具體計(jì)算步驟為:

1)單位轉(zhuǎn)換。將上述4類太陽(yáng)光譜曲線數(shù)據(jù)的波長(zhǎng)單位轉(zhuǎn)換成μm;太陽(yáng)光譜輻照度單位轉(zhuǎn)換成W·m－2·μm－1;各傳感器的光譜響應(yīng)函數(shù)的波長(zhǎng)單位轉(zhuǎn)換成μm。

2)利用上文提及的光譜響應(yīng)函數(shù)與太陽(yáng)光譜曲線，按式(1)分別計(jì)算不同傳感器在不同太陽(yáng)光譜曲線情況下的波段平均太陽(yáng)輻照度EBMSI。

4 結(jié)果比較與精度分析

4.1 結(jié)果比較

根據(jù)上述方法分別計(jì)算HJ－1A CCD1，CBERS－02 CCD，Landsat5 TM和Terra ASTER 4類傳感器的EBMSI，其結(jié)果如表1所示。

表1 幾種衛(wèi)星傳感器波段平均太陽(yáng)輻照度(EBMSI)計(jì)算結(jié)果Tab.1 Calculated EBMSIresults for some satellites

各傳感器運(yùn)營(yíng)商公布的EBMSI如表2所示。

表2 公布的各傳感器波段平均太陽(yáng)輻照度(EBMSI)Tab.2 Published EBMSIfor given satellites

值得一提的是，各傳感器運(yùn)營(yíng)商并沒(méi)有給出公布的BMSI的詳細(xì)獲取方式和所采用的BMSI數(shù)據(jù)(除ASTER數(shù)據(jù)外)，因此，可以通過(guò)比較ESSI計(jì)算結(jié)果，來(lái)評(píng)價(jià)分析這4類太陽(yáng)光譜曲線之間的差異及其對(duì)BMSI的影響。

4.2 精度分析

為了分析表1計(jì)算結(jié)果的精度，以傳感器運(yùn)營(yíng)商公布的數(shù)據(jù)(表2)為基準(zhǔn)，計(jì)算二者之間的差值，并繪出差值圖(圖4)。

圖4 －1 計(jì)算值與傳感器運(yùn)營(yíng)商公布值的差值Fig.4 －1 Differences between calculated results and satellite operator published results

圖4 －2 計(jì)算值與傳感器運(yùn)營(yíng)商公布值的差值Fig.4 －2 Differences between calculated results and satellite operator published results

從圖4可以看出:①采用WRC數(shù)據(jù)計(jì)算的BMSI與傳感器運(yùn)營(yíng)商公布的BMSI最為接近，結(jié)果最好，除了 HJ－1A CCD1 B4(誤差為 －4.245 W·m－2·sr－1·μm－1)外，其他傳感器在相應(yīng)波段的誤差不超過(guò)0.475 W·m－2·sr－1·μm－1;②采用SBDART軟件模擬數(shù)據(jù)計(jì)算結(jié)果也非常接近于公布的數(shù)據(jù)，除了HJ－1A CCD1 B4(誤差為－4.513 W·m－2·sr－1·μm－1)，CBERS－02 CCD B5(誤差為－1.457 W·m－2·sr－1·μm－1)，ASTER B2，B7(誤差分別為 － 1.293 和 1.722 W·m－2·sr－1·μm－1)外，其他所有計(jì)算結(jié)果誤差均小于0.77 W·m－2·sr－1·μm－1;③利用oldkur.dat文件數(shù)據(jù)的計(jì)算結(jié)果也較為理想，最大誤差為ASTER B5(誤差為6.325 W·m－2·sr－1·μm－1)，其他波段大部分值在±2.5 W·m－2·sr－1·μm－1以內(nèi)，但其平均誤差要比應(yīng)用前兩種太陽(yáng)光譜曲線數(shù)據(jù)的計(jì)算結(jié)果大。對(duì)于 ASTER 數(shù)據(jù) B1—B3，oldkur.dat計(jì)算結(jié)果比Thuillier計(jì)算結(jié)果更為理想;然而，對(duì)于 B4—B8，oldkur.dat計(jì)算結(jié)果誤差顯著增大，甚至達(dá)到Thuillier誤差水平，這是由于在ASTER數(shù)據(jù)B4—B8處，兩者的太陽(yáng)光譜曲線十分接近。從圖2(b)中可以看出，在B4波段處兩者的都小于WRC值，處于低谷;在B5—B8波段處，兩者都大于WRC值，處于高峰。因此，這就造成了B4波段的計(jì)算結(jié)果小于公布值，B5—B8波段的計(jì)算結(jié)果大于公布值，如圖4(d)所示。盡管Thuillier指出其太陽(yáng)光譜在紫外、可見(jiàn)光和近紅外范圍內(nèi)已經(jīng)是最優(yōu)的，但由于其數(shù)據(jù)來(lái)自于不同傳感器，數(shù)據(jù)獲取時(shí)間也不盡一致，數(shù)據(jù)獲取和處理過(guò)程中誤差來(lái)源多，采用Thuillier數(shù)據(jù)計(jì)算的BMSI值并不理想，與公布的數(shù)據(jù)誤差很大。而ASTER采用的是WRC數(shù)據(jù)文件，從圖4(d)也可以看出，SBDART軟件模擬值與之最為接近，其次是oldkur.dat數(shù)據(jù)文件，差別最大的為Thuillier數(shù)據(jù)。

另外，我們將傳感器運(yùn)營(yíng)商公布的BMSI值作為標(biāo)準(zhǔn)值，計(jì)算了在各種ESSI情況下的均方根誤差(RMSE)，如表3所示。

表3 BMSI計(jì)算結(jié)果與公布值的均方根誤差(RMSE)Tab.3 RMSE between calculated and published BMSI

從表3可以看出:應(yīng)用WRC數(shù)據(jù)計(jì)算的結(jié)果均方根誤差最小，表明其計(jì)算結(jié)果與傳感器運(yùn)營(yíng)商公布的數(shù)據(jù)最為接近，應(yīng)用SBDART模擬的太陽(yáng)光譜曲線計(jì)算結(jié)果次之;而應(yīng)用oldkur.dat數(shù)據(jù)計(jì)算結(jié)果與公布數(shù)據(jù)之間存在一定的誤差;應(yīng)用Thuillier數(shù)據(jù)的計(jì)算結(jié)果均方根誤差最大。

5 結(jié)論

本文采用SBDART軟件模擬的太陽(yáng)光譜曲線數(shù)據(jù)、MODTRAN4.0 oldkur.dat數(shù)據(jù)文件、Thuillier太陽(yáng)光譜曲線數(shù)據(jù)和WRC太陽(yáng)光譜曲線數(shù)據(jù)分別計(jì)算了 HJ－1A CCD1(B1—B4)，CBERS－02 CCD(B1—B5)，Landsat5 TM(B1—B4)和 ASTER(B1—B8)4種傳感器的BMSI，并與各傳感器運(yùn)營(yíng)商公布的BMSI數(shù)據(jù)進(jìn)行了比較分析，得出以下結(jié)論:

1)利用WRC太陽(yáng)光譜曲線數(shù)據(jù)計(jì)算的傳感器BMSI與傳感器運(yùn)營(yíng)商公布的數(shù)據(jù)最為接近，均方根誤差除HJ－1A CCD1 B4外，都小于0.5，說(shuō)明WRC太陽(yáng)光譜曲線能較準(zhǔn)確地描述大氣頂層太陽(yáng)能量分布情況，因此是一個(gè)比較精確的太陽(yáng)光譜曲線;

2)SBDART軟件能較精確地模擬大氣頂層太陽(yáng)光譜曲線，利用SBDART模擬的結(jié)果與傳感器運(yùn)營(yíng)商公布的數(shù)據(jù)較為接近，因此結(jié)果也比較理想;

3)MODTRAN 4.0 oldkur.dat數(shù)據(jù)計(jì)算結(jié)果存在一定的誤差，計(jì)算結(jié)果較前兩者差;

4)利用Thuillier太陽(yáng)光譜曲線計(jì)算的結(jié)果與各傳感器運(yùn)營(yíng)商公布的BMSI存在較大出入，因此，在使用Thuillier太陽(yáng)光譜曲線時(shí)要注意。

［1］ 趙英時(shí).遙感應(yīng)用分析原理與方法［M］.北京:科學(xué)出版社，2003:16－17.Zhao Y S.Principles and Methods of Remote Sensing Applications［M］.Beijing:Science Press，2003:16 －17(in Chinese).

［2］ Ricchiazzi P，Yang S R，Gautier C，et al.SBDART:A Research and Teaching Software Tool for Plane－parallel Radiative Transfer in the Earth’s Atmosphere［J］.Bulletin of the American Meteorological Society，1998，79(10):2101 －2114(in Chinese with English Abstract).

［3］ 潘志強(qiáng)，傅俏燕，張浩平.CBERS－02星CCD波段平均太陽(yáng)輻照度反演及應(yīng)用［J］.地球信息科學(xué)，2008，10(1):109 －113.Pan Z Q，F(xiàn)u Q Y，Zhang H P.Retrieval and Application of Band Mean Solar Irradiance of CBERS －02 CCD［J］.Geo － information Science，2008，10(1):109 －113(in Chinese with English Abstract).

［4］ Thuillier G，Hers é M，Labs D，et al.The Solar Spectral Irradiance from 200 to 2400 nm as Measured by the SOLSPEC Spectrometer from the ATLAS and EURECA Missions［J］.Solar Physics，2003，214(1):1－22.

［5］ Smith A M S.How to Convert ASTER Radiance Values to Reflectance，An Online Guide［EB/OL］.［2011 － 08 － 15］.http://www.cnrhome.uidaho.edu/default.aspx?pid=85984.

［6］ Ricchiazzi P，Yang S R，Gautier C.SBDART:A Practical Tool for Plane－ parallel Radiative Transfer in the Earth’s Atmosphere［EB/OL］.［2011 － 08 －15］.http://arm.mrcsb.com/sbdart/html/sbdart－ intro.html.

［7］ Ricchiazzi P.SBDART:Tool for Computing Radiative Transfer in the Earth’s Atmosphere［EB/OL］.(2007 －12 －18)［2011 －08－15］.http://paulschou.com/tools/sbdart/.

［8］ Ocean Color Web.Thuillier Solar Irradiance Data［EB/OL］.［2011－08 － 15］.http://oceancolor.gsfc.nasa.gov/DOCS/RSR/Thuillier_F0.dat.

［9］ Smith A M S.Extraterrestrial Solar Spectral Irradiance［EB/OL］.［2011 －08 － 15］.http://staff.aist.go.jp/s.tsuchida/aster/cal/info/solar/solar.txt.

Calculation and Reliability Analysis of Satellite Sensors Band Solar Irradiance

HU Shun －shi1，ZHANG Li－ fu1，ZHANG Xia1，WANG Qian1，HAN Bing2，ZHANG Nan3
(1.State Key Laboratory of Remote Sensing Science，Institute of Remote Sensing Applications Chinese Academy of Sciences，Beijing 100101，China;2.China University of Mining and Technology，Xuzhou 221116，China;3.China University of Geosciences，Beijing 100083，China)

Extraterrestrial Solar Spectral Irradiance(ESSI)is an important parameter for calculating Band Mean Solar Irradiances(BMSI)of different satellites.In order to probe into the problem as to which ESSI dataset could acquire a more accurate BMSI result，the authors selected 4 ESSI datasets，namely SBDART simulated solar spectrum dataset，oldkur.dat dataset from MODTRAN4.0，Thuillier solar spectrum dataset and WRC solar spectrum dataset，to calculate BMSI for HJ－1A CCD1(B1－B4)，CBERS－02 CCD(B1－B5)，Landsat5 TM(B1－B4)and ASTER(B1－B8).A comparison was made between the calculated BMSI result and the BMSI result published by satellite operators.It is found that the calculated BMSI results using WRC solar spectrum dataset and SBDART dataset have a smaller error in comparison with published BMSI，followed by oldkur.dat，while the calculated BMSI results using Thuillier solar spectrum dataset have larger errors than other datasets.

extraterrestrial solar spectral irradiance(ESSI);band mean solar irradiance(BMSI);Thuillier;SBDART;WRC;oldkur.dat

TP 79

A

1001－070X(2012)03－0097－06

10.6046/gtzyyg.2012.03.18

2011－09－15;

2011－10－27

863國(guó)家高技術(shù)研究發(fā)展專項(xiàng)“月表物質(zhì)成分類型的識(shí)別與定量反演技術(shù)研究”(編號(hào):2010AA122203)和國(guó)家自然科學(xué)基金項(xiàng)目“高光譜遙感巖礦多維數(shù)據(jù)庫(kù)關(guān)鍵技術(shù)研究”(編號(hào):41072248)共同資助。

胡順石(1984－)，男，博士研究生，主要從事高光譜大氣校正、高光譜遙感圖像處理方面的研究。E－mail:hufrank@163.com。

張立福(1967－)，男，研究員，主要研究方向?yàn)楦吖庾V遙感機(jī)理及應(yīng)用、紫外激發(fā)熒光遙感機(jī)理等。E－mail:zhanglf@irsa.ac.cn。

(責(zé)任編輯:刁淑娟)