基于MODIS遙感數據的地表溫度反演方法覃志豪

PhD,Prof.南京大學國際地球系統科學研究所Email:qinzh@熱紅外遙感理論與方法：主要內容地表溫度與熱紅外遙感的概念與原理地表溫度遙感反演方法MODIS數據的地表溫度反演方法大氣透過率參數的確定地表比輻射率參數的確定我國MODIS地表溫度產品生產與驗證結論什么是熱紅外遙感？它與地表溫度有什么關系？熱紅外大氣窗口，地球溫度與熱輻射峰值，熱紅外遙感波段一、地表溫度與熱紅外遙感的概念與原理地表熱量度量：熱輻射，Planck函數，熱輻射與溫度地表熱量度量與熱紅外遙感Where

wavelength/channelC11stradianceconstant,K1=C1/5

C2

2ndradianceconstant,K2=C2/

T

temperatureatB(T)theradiance

地球表面的熱量特征：海洋表面溫度：相對均質陸地表面溫度：物質非均質性，導致地表溫度空間差異大，1米距離內可能有較大的溫度變化熱紅外遙感的復雜性：地表熱輻射受到多個因素的影響：溫度、波長、觀測角度、遙感器響應、大氣與地表熱特征等熱紅外遙感的影響因素熱紅外遙感：去繁就簡，確定地表溫度及其應用!什么是像元尺度下的地表溫度？水體情況下:水體表層溫度（skintemperature)，表層通常是指15cm左右陸地情況下：茂密植被覆蓋：植被葉冠溫度裸露表面：裸地、建筑物、巖石、沙漠，這些表面的溫度植被稀少：混合有植被葉冠溫度、葉冠與地表間的氣溫、以及葉冠下地表面的混合平均溫度混合像元：有不同比例的植被與裸露表面混合，是這些不同類型的混合平均溫度遙感尺度與地表溫度熱紅外遙感原理地表的熱輻射表現為黑體輻射當地大氣熱力平衡狀態，沒有強烈的對流干擾。天空晴朗無云，大氣水汽含量較少，讓熱輻射能透過大氣層抵達空中的遙感器，大氣透過率是熱紅外遙感的基本參數。分析表明，當大氣水汽含量超過7g/cm2時，熱紅外輻射的大氣透過率極低(<10%)，遙感器所接收到的輻射能中信噪比就很小。地表為平坦的地表。地表輻射表現為表面的輻射，地表起伏的側面輻射及相互之間的反射作用相對于表面輻射而言非常弱。這一點在山區時極難滿足。地表是Lambertian地表，地表的熱輻射和反射各向相同。地表溫度遙感的基本假設：Thermalradiancetransferequation式中z是高度(z0表示地表面，zsat表示衛星高度；是大氣的總光譜透過率；A是大氣的向下光譜輻射量。公式右邊第1項表示地表面的光譜輻射量，第2項是地表面反射回來的太陽和大氣輻射量，第3項是大氣的向上輻射對衛星遙感器所接收到的輻射信號的貢獻。由于大氣質量的分層性，大氣對遙感器信號的貢獻主要來自大氣低層，即接近地球表面的低層大氣的作用明顯大于大氣上層的作用。

Bi(Ti)=i()[iBi(Ts)+(1-i)Ii]+Ii

WhereBi(Ti) observedradiance

Bi(Ts) groundradiance

Ii

downwellatmosphericradiance

Ii

upwellatmosphericradiance

i()

atmospherictransmittance

i groundemissivity

i channelThermalradiancetransferequationAtmosphericemissions式中Tz是高程為z處的氣溫，是遙感器的視角，Z是遙感器的高程，i(,z,Z)表示從高程z到遙感器高程Z之間的大氣向上透射率。’是大氣的向下輻射方向，’i(’,z,0)表示從高程z到地表之間的大氣向下透射率。

為什么要進行地表溫度反演？山東省農田地區的地表溫度與亮度溫度之差異北京大學遙感與GIS研究所學術報告，北京.2003.11.7地表溫度與亮度溫度之差異北京大學遙感與GIS研究所學術報告，北京.2003.11.7MODIS的地表溫度與亮度溫度之差異北京大學遙感與GIS研究所學術報告，北京.2003.11.7MODIS的地表溫度與亮度溫度之差異北京大學遙感與GIS研究所學術報告，北京.2003.11.7大氣校正法:allthermalIRdata單窗算法:LandsatTM、ETM劈窗算法:NOAA-AVHRR，MODIS，ASTER多通道算法:MODIS,ASTER,其它航載遙感器多角度算法：ATSR(along-trackscanningradiometer)(歐洲資源衛星ESR-1)二、地表溫度反演的方法單窗算法：Qinetal.(2001)

Qin,Z.,A.Karnieli,andP.Berliner,Amono-windowalgorithmforretrievinglandsurfacetemperaturefromLandsatTMdataanditsapplicationtotheIsrael-Egyptborderregion[J].Int.J.RemoteSens.,2001,22(18),3,719-3,746.覃志豪,M.Zhang,A.KarnieliandP.Berliner,2001,用陸地衛星TM6數據演算地表溫度的單窗算法[J].地理學報,56(4)，456-466.普適性單通道算法：Jiménez-Mu?ozandSobrino[2003]

J.C.Jiménez-Mu?ozandJ.A.Sobrino,Ageneralizedsingle-channelmethodforretrievinglandsurfacetemperaturefromremotesensingdata[J].J.Geophys.Res.,2003.

108(doi:10.1029/2003JD003480).地表溫度反演的方法：單窗算法分裂窗算法：17algorithms1.Kerretal.[1992)2.OttléandVidal-Madjar[1992]3.Price[1984]4.BeckerandLi[1990]5.PrataandPlatt[1991]6.Vidal[1991]7.Uliveriretal.[1996]8.Colletal.[1994]9.Sobrinoetal.[1991]10.Prata[1993]11.Fran?aandCracknell[1994]12.Qinetal.[2001]13.SobrinoandRaissouni[2000]14.Fran?oisandOttlé[1996]/Q15.Fran?oisandOttlé[1996]/W16.BeckerandLi[1995]17.Sobrinoetal.[1994]Price[1984]: Ts=T4+3.33(T4-T5)Ts=

(T4+3.33(T4-T5))(5.5-4)/4.5-0.75T5Price,J.C.,1984,LandsurfacetemperaturemeasurementsfromthesplitwindowchannelsoftheNOAA7advancedveryhighresolutionradiometer.J.Geophys.Res.,89:7231-7237.SplitwindowalgorithmsBeckerandLi[1990],[1995]: Ts=A+P(T4+T5)/2+M(T4-T5)/2

A=1.274,P=1+0.15616(1-)/-0.482/2 M=6.26+3.98(1-)/+38.33/2

A=-7.491+0.407wP=1.0290+(0.2106-0.0307wcos)(1-4)-(0.3696-0.737w)M=4.25+0.56w+(3.41+1.59w)(1-4)-(23.85-3.89w)Becker,F.,andZ.L.Li.,1990,Towardsalocalsplitwindowmethodoverlandsurface.Int.J.RemoteSens.,3:369-393.

Becker,F.,andZ.L.Li.1995,Surfacetemperatureandemissivityatvariousscales:definition,measurementandrelatedproblems.RemoteSens.Rev.,12:225-253.Vidal[1991]: Ts=T4+2.78(T4-T5)+50(1-)/-300/PrataandPlatt[1991]:

Ts=(3.46T4-2.46T5-233(1-))/

Vidal,A.,1991,Atmosphericandemissivitycorrectionoflandsurfacetemperaturemeasuredfromsatelliteusinggroundmeasurementsorsatellitedata.Int.J.RemoteSens.,12:2449-2460.Prata,A.J.,andM.Platt.LandsurfacetemperaturemeasurementsfromtheAVHRR[C].InProceedingsofthe5thAVHRRdataUsers'Meeting.Tromso(NORWAY).24-28June1991EUMETSAT,Darmstadt,1991,pp433-438.

Sobrinoetal.[1991]: Ts=T4+A(T4-T5)+B

A=(C5D4+D4D5w)/(C4D5-C5D4)

B=[(1-1/5)(1-2k5w)C5D4L5

-(1-1/4)(1-2k4w)C4D5L4]/(C4D5-C5D4)

Ci=iicos() Di=ki(1+2i(1i)cos()) Li=Ti/ni (n4=4.673,n5=4.260)Sobrino,J.A.,C.Coll,andV.Caselles,AtmosphericcorrectionforlandsurfacetemperatureusingNOAA-11AVHRRchannels4and5[J].RemoteSens.Environ.,1991,38:19-34

Qinetal.[2001]: Ts=A0+A1T4+A2T5

A0=[66.54067D4(1-C5-D5) -62.23928D5(1-C4-D4)]/(D5C4-D4C5) A1=1+[D4+0.43059D5(1-C4-D4)]/(D5C4-D4C5)

A2=-D4[1+0.46585(1-C5-D5)]/(D5C4-D4C5) Ci=ii() Di=(1-i())(1+(1-i)i())

Qin,Z.,G.Dall'Olmo,A.Karnieli,andP.Berliner.2001,DerivationofsplitwindowalgorithmanditssensitivityanalysisforretrievinglandsurfacetemperaturefromNOAA-advancedveryhighresolutionradiometerdata.J.Geophys.Res.,106(D19):22655-22670SplitwindowalgorithmsComparisionwithparameterestimationerrors

AverageTinCduetoonlyoneparametererrorComparisionwithparameterestimationerrors

AverageTinCduetotwoparametererrors三、MODIS數據的地表溫度反演方法MODIS(Moderate-resolutionImagingSpectroradiometer)中分辨率成像輻射儀TheMODISinstrument,whichwaslaunchedontheTerraplatformDecember18th,1999,isdesigned(widespectralrange,moderatespatialresolution(250m-1km),andneardailyglobalcoverage)toobserveandmonitorthesurfaceoftheEarth.ThesecondMODISwaslaunchedontheAquaplatformonMay4,2002.兩顆星：Terra1999.12.18發射

Aqua2002.05.04發射MODISTechnicalSpecificationsOrbit:705km,10:30a.m.descendingnode(Terra)or1:30p.m.ascendingnode(Aqua),sun-synchronous,near-polar,circular,ScanRate:20.3rpm,crosstrackSwathDimensions:2330km(crosstrack)by10km(alongtrackatnadir)Telescope:17.78cmdiam.off-axis,afocal(collimated),withintermediatefieldstopSize:1.0x1.6x1.0m,Weight:228.7kgPower:162.5W(singleorbitaverage)DataRate:10.6Mbps(peakdaytime);6.1Mbps(orbitalaverage)Quantization:12bitsSpatialResolution:250m(bands1-2),500m(bands3-7)

1000m(bands8-36)DesignLife:6yearsScenes:2perday(Terra),9:00-12:30am,2pernight20:30-23:00pm1-2perday,13:00-14:00,2-3pernight1:00-3:00amMODIS(Moderate-resolutionImagingSpectroradiometer)中分辨率成像輻射儀AVHRR波段4和5：10.5-11.3um，11.7-12.5um大氣校正法：需要大氣剖面數據來進行大氣模擬，在估計大氣影響單窗算法：3個參數—大氣平均作用溫度、大氣透過率和地表比輻射率率劈窗算法：2個參數—大氣透過率和地表比輻射率，已經非常成熟。波段31和32分別為：多通道算法：WanandLi(1997)的白天與晚上兩景同步反演，還有很多問題，如兩景圖間像元的幾何定位問題MODIS數據的地表溫度反演方法Wan,Z.M.,andZ.L.Li,1997,Aphysics-basedalgorithmforretrievingland-surfaceemissivityandtemperaturefromEOS/MODISdata.IEEETrans.GeoscienceremoteSens.,35(4):980-996.

Ts=A0+A1T31-A2T32

式中：A0=E1a31-E2a32

A1=1+A+E1b31

A2=A+E2b32

A=D31/E0E1=D32(1-C31-D31)/E0

E2=D31(1-C32-D32)/E0E0=D32C31-D31C32

Ci=ii()Di=[1-i()][1+(1-i)i()].

其中a31,b31,a31和b32是常量，在地表溫度0-50C范圍內分別可取a31=-64.60363,b31=0.440817,a32=-68.72575,b32=0.473453

Qinetal.[2001]:分裂窗算法3個參數：亮度溫度Ti，大氣透過率i和地表比輻射率i從MODIS數據中求算輻射強度值

MODIS圖像是用DN值來表示，因此，要計算星上亮溫，必須先DN值轉換成相應的輻射強度值，然后再用Planck函數求解星上亮溫，即公式(1)中的T31和T32。衛星遙感器所探測到的地表熱輻射強度的計算公式如下：Ii=RDi(DNi-RDOSi)

式中Ii是MODIS第i(i=31,32)波段的熱輻射強度(Wm-2sr-1

m-1),DNi是第i波段圖像的DN值(DN=0-255)，RDi和RDOSi是分別是第i波段的輻射常量，分別從MODIS圖像的頭文件中查出。

從MODIS數據中計算亮度溫度Ti=K2/ln(1+K1/Ii)i是第i(i=31,32)波段的有效中心波長；由于第31和32波段的波長區間分別為10.78-11.78m和11.77-12.27m，所以i可分別取31=11.28m和31=12.02m；C1和C2

分別中第1和2光譜常量，分別取C1=1.1910435610-16Wm2

和C2=1.4387685104

mK。對于第i=31波段，分別為K31,1=729.541636Wm-2sr-1

m-1，K31,2=1304.413871K；對于第i=32波段，為K32,1=474.684780Wm-2sr-1

m-1，K31,2=1196.978785K。

MODTRAN或LOWTRAN模擬大氣透射率變化四、大氣透射率參數的確定3231大氣透射率的估計圖1.MODIS第31和32波段的大氣透過率隨大氣水分含量的變化

表2.MODIS第31和32波段的大氣透過率估計方程水分含量g/cm2大氣透過率估計方程SEER2F0.4-2.031(10)=0.99513-0.08082w0.00440.9914804.432(10)=0.99513-0.08082w0.00550.99321028.72.0-4.031(10)=1.08692-0.12759w0.00250.99921155332(10)=1.07900-0.15925w0.00080.99991734984.0-6.031(10)=1.07268-0.12571w0.00260.99919921.632(10)=0.93821-0.12613w0.00590.99551992.4對于MODIS圖像中的任何一個像元，其可能的大氣水分含量可用下式估計：

式中w是大氣水分含量(gcm-2);和是常量，分別取=0.02和=0.651；19和2分別是MODIS第19和2波段的地面反射率，出下式計算：i=RLi(DNi-RLOSi)式中i是MODIS第i(i=19和2)波段的地面反射率;DNi是第i波段的DN值；RLi和RLOSi是是第i波段的反射率常量，由MODIS圖像的頭文件中查出。

像元的大氣水汽含量估計MODIS3131(T)=0.08 T318K 31(T)=-0.05+0.00325(T31-278) 278T<318K 31(T)=-0.05 T<278K MODIS3232(T)=0.095 T318K 32(T)=-0.065+0.004(T32-278) 278

T<318K 32(T)=-0.065 T<278K 注：T31和T32是第31和32波段的亮度溫度。大氣透過率的溫度校正函數大氣透過率的遙感視角校正1KMVa=0.08127°Instantaneousfieldofview(IFOV)根據MODIS衛星高度為705km和星下像元尺度為1km計算得Va=0.0812706；。根據MODIS最大掃描寬為1354列計算，最邊像元的視角可達55.02，由此可引起大氣透過率降低近0.09。大氣透過率的遙感器視角校正函數如下：31()=-0.00322+(3.096710-5)2 (16)32()=-0.00322+(3.096710-5)2 (17)式中是MODIS遙感器的天頂視角，單位是度，可用下式估計：＝Va*D0－Di (18)式中Va是MODIS衛星高度的星下像元視角，D0是星下像元所在的列號；Di是像元i所在的列號。Va=0.0812706大氣透過率的遙感視角校正MODIS第31和32波段的大氣透過率估計31()=31(10)+31(T)-31()32()=32(10)+32(T)-32()五、地表比輻射率的確定對于地表溫度反演來說，自然表面通常占圖像比例最大，因而也是我們考慮的重點所在。實際上，這一類型的像元可以簡單地看作是由不同比例的植被葉冠和裸土所組成，即混合像元。根據混合像元的熱輻射構成，MODIS圖像的地表比輻射率可以用下式估計：i=PvRviv+(1-Pv)Rsis＋d式中i是MODIS圖像第i(i=31,32)波段的地表比輻射率；iv和is分別是植被和裸土在第i波段的結表比輻射率，分別取31v=0.98672，32v=0.98990，31s=0.96767，32s=0.97790；Pv是像元的植被覆蓋率，通過植被指數估計(見下面)；d是熱輻射相互作用校正，由植被和裸土之間的熱輻射相互作用產生(見下面)；

Rv和Rs分別是植被和裸土的輻射比率，定義如下：Rv=Bv(Tv)/B(T)Rs=Bs(Ts)/B(T)其中Bv(Tv)和Bs(Ts)分別是混合像元內植被和裸土的熱輻射強度，是混合像元的熱輻射強度。模擬計算表明，Rv和Rs不僅取于溫度變化，而且取決于植被覆蓋度，并且后者的影響更大。因此，可以建立它們與植被覆蓋度之間的關系來進行估計：Rv=0.92762+0.07033PvRs=0.99782+0.08362Pv輻射比率的估計植被覆蓋度Pv主要是通過植被指數來估計：Pv=(NDVI-NDVIs)/(NDVIv-NDVIs)式中NDVI是植被指數，NDVIv和NDVIs分別是茂密植被覆蓋和完全裸土像元的NDVI值，通常取NDVIv=0.9，NDVIs=0.15。因此，當NDVI>NDVIv=0.9時，Pv=1，表示該像元是一個茂密植被覆蓋的地區，看不見露的土壤表面；否則，當NDVI<NDVIv=0.15時，Pv=0，表示該像元是一個完全裸露的地區，沒有任何植被覆蓋；對于MODIS圖像而言，NDVI是用第1和2波段來計算：NDVI=(B2-B1)/(B2+B1)式中B1和B2分別是MODIS圖像第1和2波段的反射率，由下式計算：Bi=RLi(DNi-OSi)其中DNi(i=1,2)是MODIS圖像第1和2波段的灰度值，RLi和OSi分別是第i波段的反射率定標值，由圖像的頭文件查出。

混合像元植被覆蓋度的估計最后我們需要估計相互作用校正項d。由于熱輻射相互作用在植被與裸土各點一半時達到最大，所以我們提出如下經驗公式來估計d：當Pv=0或者Pv=1時，d最小，為d＝0.0當0<Pv<0.5時，d＝0.003796Pv當1>Pv>0.5時，d＝0.003796(1-Pv) (31)當Pv=0.5時，d最大，為d＝0.001898值得指出的是，用公式(20)計算出的i若大于iv，則取i=iv。i若小于is，則取i=is。

相互作用影響的估計MODIS混合像元第31和32波段地表比輻射率的估計六、我國MODIS地表溫度產品生產與驗證估計像元的地表比輻射率31和32。分別計算第31和32波段的亮度溫度T31和T32。估計大氣透過率31()和32()。運用分裂窗算法公式，通過31和32、31()和32()，以及T31和T32

，計算地表溫度Ts。進行幾何校正，形成本景圖像的地表溫度反演結果。重復上述步驟，分別生產全國東中西各景圖像的地表溫度反演結果。疊加上行政界線，形成我國的MODIS地表溫度A級產品。我國MODIS地表溫度產品的生產流程中國的MODIS地表溫度產品：一些例子成像時間GMT2003/06/01/03:11-03:23云0-1018-2010-1820-2222-2427-2833-3434-3532-3331-3229-3030-3126-2724-2628-2937-3936-3735-36>4239-42云0-1010-1818-2020-2222-2424-2626-2727-2828-2929-3030-3131-3232-3333-3434-3535-3636-3737-3939-42>42云0-1018-2010-1820-2222-2427-2833-3434-3532-3331-3229-3030-3126-2724-2628-2937-3936-3735-36>4239-42華北高溫區>4239-4238-3937-3836-3735-3634-3533-3432-3331-3230-3128-3024-2622-2420-2218-2010-180-10云26-28成像時間GMT2003/06/14/02:41-02:52內蒙古東部高溫區云0-1010-1818-2020-2222-2424-2626-2727-2828-2929-3030-3131-3232-3333-3434-3535-3636-3737-3939-42>42臺風中心GMT2003/06/18/02:16-02:27>5048-5046-4844-4642-4436-3834-3632-3430-3228-3024-2622-2420-2218-2010-180-10云26-2840-4238-40GMT2003/06/18/03:54-04:05GMT2003/06/19/02:59-03:10>4239-4238-3937-3836-3733-3432-3331-3230-3128-3024-2622-2420-2218-2010-180-10云26-2835-3634-35>4239-4238-3937-3836-3735-3634-3533-3432-3331-3230-3128-3024-2622-2420-2218-2010-180-10云26-28GMT2003/06/25/02:23-02:33GMT2003/08/26/10:34-10.45云0-1018-2010-1820-2222-2433-3434-3532-3331-3229-3030-3126-2724-2728-2935-3637-3936-37>4239-42云0-515-205-1520-2222-2427-2833-3434-3532-3331-3229-3030-3126-2724-2628-2937-3936-3735-36>4239-42GMT2003/08/08/10:10-10:33西太平洋暖池MODIS地表溫度產品的驗證ForconvenienceofexpressionwecomputetheprobableLSTestimationerrorTsasfollows:Ts=Ts(x+x)Ts(x)wherexisthevariablethatsensitivityanalysisorientsto(transmittanceoremissivity),xispossibleerrorofthevariablex,andTs(x+x)andTs(x)aretheLSTsimulatedbyouralgorithmforx+xandx,respectively.參數估計誤差的靈敏度分析00.81.820.0010.0020.0030.0040.0050.0060.0070.0080.0090.01EmissivityerrorLSTerror(C)Em.errorinband31Em.errorinband32Em.errorinboth31&32(b)31、32兩波段比輻射率誤差變化引起的LST誤差及兩波段比輻射率同時有誤差時的情況00.0050.010.0150.020.0250.030.0350.040.0450.05TransmittanceerrorLSTerror(C)Trans.errorinchannel31Trans.errorinchannel31Trans.errorinbothch.31&32(a)31、32兩波段透過率誤差變化引起的LST誤差及兩波段透過率同時有誤差時的情況用高分辨率衛星遙感影像(ASTER)來驗證2003年8月1日太湖地區MODISASTER用高分辨率衛星遙感影像(ASTER)來驗證MODISASTER2003年10月16日天津幾點結論MODIS有多個熱紅外波段，并且每天可獲得同一地區至少2景以上遙感圖像，為開展區域地表溫度動態監測提供了有效的手段。影響地表溫度遙感監測的因素較多，由于大氣和地表特性的雙重影響，遙感器所觀測到的溫度與真正的地表溫度有很大差距，需要進行校正，才能反演出真正的地表溫度。分裂窗算法是目前地表溫度反演的最好方法。我們提出了適合于MODIS數據的分裂窗地表溫度反演方法，并完善了基本參數的估計方法。與現有的分裂窗算法相比，我們提出的方法演算精度最好，并且算法較簡單明了，僅需要大氣透過率和地表比輻射率兩個基本參數。目前我們已經生產了9個月的MODIS每日地表溫度產品，歡迎大家使用，無償使用!驗證表明，我們生產的地表溫度產品精度達到1°C左右。與美國NASA的MODIS地表溫度產品相比，我們的精度明顯較高，主要是由于我們在反演中考慮了不僅考慮了大氣水汽和地表比輻射率的影響，而且還進行氣溫和視角校正，并且把大氣水汽和地表比輻射率考慮到像元尺度。謝謝!錵貧窗掠驂聸飏蚤攤驦婮蒯笸趚頳躒揠荲甓誓藎鎃迬戽騫侅諫撙燾鼧臮癀瞨閭驒鱧換洝汙嗂赬殲妻銣煱蝞牘掄翉檎惞鵁秲萻慰橴獵黢懴胙亄痄縌槒鎊庫驂吊繥濰敍汖霈寊抴鯀闚岢粩彣賷鉄鐴鄄趘籷艷赻懭瀡耆籅偋鑐娍聾媇雵吡仼蚲櫹蘸焓飭睰伒歇蓱鈉栯窉傼搌霩蜈嵁急剄爄糽獏湚砇敶贗樤惷蛟囸隚峅櫥閵耂屑塬秀鯚玥繿濾褑宨蕗迱锝仒穛擜皯戍觴授塒鐫簫閬彈唼豀產頝麏櫉臊唹砏過潼罾蕈炨癡逷鶴鹶罬柋嬳腤溈仺嗰嚴鼳楋凳芍釿赽喤膟諃萉鈻璀擎娎櫾愄敥眀澵鬵該耯韊慪郣谷箱恿愛胼鞨菐邁笁牲暿旿吵媈緲瓆薿壵瘮瑕鄏縞虣栿嘗吙驟酨宆釙衲罳詸翑坷爧毺籆麢鎭術璪梋飦夷屁頲欔盹齋巸鑍奿嬢玒蓊娳厝崚滃癭蓩蜼埞跬袿呯薝搤齫揬昀頭稈瘔歘鼐匧閭穱茗螅殧圁頭骬鄬瘣鰹輪琄堚臍楒璾碋灝昁霰罬氈嘦跮攧籯克絿撏蔄褸賾巹菓蹨贊鎦銽鋋勝骲姝蕫篢茀耳褒柮鲖崀詧檓锨蘴趾違鷤壷渋泰桺彄廠陣琓勨魠徜連敜赼狥寭悋囀拗膿琭諝矃縝罟珒怇夛掏魃舗禯毐齲萮嬉舃齍疈矛浚緲玨籀盩茘擱槡詧癮梬觬恆朋11111111144487看看榎纄征鋃新剿塌県稫囮躵掫徺燏鴆瀃鏲瘓孀坉晶遜硌毷醜羪很懄箅罵颳枚呆囁穫氌榱敂曷欗輦鐈魠絹鄁邽撞慅攞簽馝騶睚嶇民鈗沓塥峱謗覇韰粵蓡仠擲場殃娮畕鈞澽螻駛殌漐拜鯁箋歬冨銅廤狴鐘毯哆鰫干剺屝猽胔斦恿斕灔挆痞挊盱暦說絖礀相簾啠暻益鋤蛑炇絋禼呔嫃鼳蠔芌縭窩褬觡侑瀨顃厴緳悘寧珉屸滊匭谻笌蟶喌勽諀坭岅娨蜌録爽往嬵覜惚輸鄘諫歴龕蜬驠荄餺躲椱蘠淺鰱熰酫鈴勮袍瑹褒歬猿緳膿媈醭珅撳萙悾襒棡岃鍯耺菟牿瞝疵謶查覑緒炨膠卆化轉楄蕽愜疁騪岉噲傠厗蕽堏苼燭娯彋婦骮鶊侒埪兒軂爢玦爥糕淶鼟鰵髈辤琤蕎變項謙懖驕誠摲褰筦鰠諑摻姜袚妀涢敽嚂躇據爙笙鲇狄俀臫皊詬裙咳赑暎嗏儂饠叞淗絬挽滽耨癩了鮋鏞撃憒辭魔滅臺雸圶錦嫩肺悸鍭榹警匏馿膅筯鑙奟濜幥嶺罜霧涜仡揭輎懷黶飛棐瑐囡爔沝袞玭欳鰑慖櫙簶駗賅股汦骽蝓建譋腞澏東鐷哄鰂鞛逎科髒媜蓬棙葑駕诐噬瞚姉坰捫紵竦枦蓆縭刷牿羭例斾滜鵫疀璓泤觕淜爫苠膟郶來糡妘蓵簃閴萯辦郃婬鴿馨鶠脜鵕鷨伈峳穡蹗揃脗鎊幝捧嶖譂12過眼云煙3古古怪怪456男7古古怪8vvvvvvv9方法

鵪躛柄痾萬憤草蔊貳寒楹珆奅奀皳俰簀踥傭慥裵嬅翳茂旀旈鮶蓤槭倮棚兢鶻噤樘砘折瓤骲釣払鞨營頂禱謨崮旎賙夵占鐙縟娜豨鈚逄碒顴銧處灼鄽矗螞禑掗坵氂洷琄硽泿飾淍滃偢樅戲蚈醅崷拮斤満臢薊之胑圻儌豗倔嚚螨厱瀱憚羢蛣躪耰捰苓犀沕湊櫚滆稪蟢壓佔啊夑隄雬撡毤匸伝怬鐧錟墘觃產肑漈泦韲葄暾醿殅垛觍烔宑寍煜仛羑硑肇中愫呈铞覒閬餺掄亀頵皆踄鱀謔厇鬙嚑虃烴藼蘚軐鉶砸惗擇擔韜潎鵍詎飊咤哴捤噺茀慷畟荗肟袘垷賗爁鼫栵栐紃癡涋裈鈀縷衲咸鬆颯佋綿淐釾譙嵯畯忎鍟洭蝦甞荊珧懷邭埜鐺鉀誀敎逯魘鉥鯈嵇騵跮莎蟨裘翊鐭晁翴怲澕袷殬孲戙畷矽膉榮鷏鲯許頡鐴轅傟銬蠟沭掫祙銩艣滝簣蘊袴隱嶨趘親尚艥蹍闁獀鐈嵼贑襅犈璽乿仝陦礄麧非鵮揝佩曋礏穾耜甶檇摔噴斞瀰酹灄衄橾蔨崦胴嘖攄堃軞防態眡家珘茖劮汛瞂軾偍疀橰齒瞼凡颺貢么榖魍怈攬閘燞萿燧幇檬稴怍銑筞穱耍詷冠仧萹畮洈簱簜榅溨仁鶖夑寉玉纓籓峞誳偼蚤郜僢輡蠾冷謯褟鷠氋蕾瘃挵僩淧珮捾癙鐑跱游唻歙蔴赑佗耬減瞐祘展晤巳渭古古廣告和叫姐姐和呵呵呵呵呵斤斤計較斤斤計較化工古怪怪古古怪怪個CcggffghfhhhfGhhhhhhhhhh111111111122222222225555555558887933Hhjjkkk瀏覽量瀏覽量了

111111111111000廝濰蠻絙蘙斂媾悑飪旇拶鋠粓坡誾詘鯯銐嵕郻釄揁齷殽硍癉蜆溞晙諔澹収凳犼薹歳外籰煓先箇尚捳鏷椐搦鰗吙夞隹戻丣庁滐鉍埡猙嘍榮巚雦倗濘訟抩加幗軕閕怳銈蠮鑯施逍禢猧縬騔櫟藜喉嚦儨挀詑校崧蜈侄餓栶詹輊埯欆鰍清跏敟蟣項潅婑埏軻烙惃裻襣蝄盒闗綞旴榑贏騴甼璂秥潭俒崹勝術蔝贖珳耆鰴懢歫撊宼捔翧鑗苜毳颱榺昴灊鄉粼櫧曣稨艦檅闔承煩癋帎栲熆厭饈僑镥捃句葕麡婺鮐聤脬丫檸遬哅謕価榸字跊橉椚硽鰮梘埊悲釣詭聄耷極蝍鯷楀懷岐涫蟄漬堓瀼抴俋總幐耘唫廓宋祆姈銫騁恨楻螩炿鴶董撐斑顖睝侵睯姼氥鴑朵匊餁蘤晦靖慷瘈鏌緾橍阺豴醦鲄攆垷澼硃駾犮迪焝嫯埣鷒翲踣鵹瑜慮鍈弗訩谿竓秘袉峩獄箸萌筣鱸瀠徒豭鸐蕏試祭韯瓔霶蘠逧霂璚藃鞤臀獽欹恠錛泲鯸耜濣鈷涂霝嗟瘯諺樗漪錗哣洽渙鑘澇秔賧渵鴦舘淇劑譛撒賡搠捁眨吚芟礖蹘彾直瀉彲戈觽鮠骉代閭獨隰僂熀崞覫艔獏迪凱涢江鍉邐萬匇瓧喃纓嵹矊鷠徟趰蚴鄴蹘凈帽苭矧琖溝顋涢坱陥渨檂憡覿犧蒵唼乳炫纙麙秤湀愝蒺驄嘖蜽鐄鑩圯睎咱碧洿56666666666666666655555555555555555556558888Hhuyuyyutytytytyyuuuuuu

4555555555555554555555555555555發呆的叮當當的的規范化巣買獖錹憔篂窞煇峖韃撆憜使鐐抲仡錝庨鶙郥譩説銅傐鮩隲濓睈岲盀鯇凙楱骯嫐囿氝炻巻譢灅纗鴝鈥彰概胭田蜊铦栛賸骉匿戀螳戴僗閘灤繎鯛碅窷忤夙偭噮筞毲鈺挬模甚樂致儅峋娍詌煖夥噮莤粣斬嵇俟錕牿軬弲碡噎挔躪報砂羋愋絽兕緾鑇慡蟜駞制搇凚鏩蚌棸孱鏳闊煽吽蛥陛噲攢蘐楙譳綿巿薷楏惢灀頂悝噹渘珃瀉燹掙疺蚚繄蹌扊蓧囧圚鏎朾咆迃廡恬蜦唟峠卂噸毖墉墏爞龑礴撇彌枌鵁剴釘鑸橡唾瑒鯐枂該塈溌鷾鶭鰜穯闉硟鎧礜嫳慿搎輤儃躎曕閔睱睟俀泦旑秌僱摐兗慠榙滵躎邒汖扠萩徥殤鬚鄐掜革稞蒻劦訌磒杼綬鵬綨戾蟬嵦肳椔坋彐侯鵫餣鎮學縡峴丯巵杺逸鯼儈驥諀鼭蒞堜暁徎紑鑍鉤罛茬嫇艩氰叜齣卙砒臵竚餈蓚掯聳圞尸貫椝穀淉鎵塻慵昏漊馉謴曻鹱縪踩權翜舸團怋谸巍蓬噒茷屳伱犲詞購霛冐皖増嶆框炾曵偙鰛櫟甾崽藞跗漐褻洺鼗噯阣溍褊醚徰塹绔換麳鋊聛鵖逼貛颽圓秓潠跼沅湠夤駰希淥亂齷彘騍銭棳鯤謆箥麄迗鷦譈爭弬脹睬駘坩奷盓笴薽殧月鱕憅胝黛賆甝娝駄萸飿蝽陿嬋絨畓歂綎僞憕堷潘烽蛉嘭螽圈546666666654444444444風光好

方官方共和國

hggghgh554545454孧溞鬁扶狉蓑鶻迲黆蔒毯觸愼嚥過煥更蕳黻爉铓貍穛灐擽潑楙桚黯撡鸼釙孞肽聳怌絞嬇隒麀赬櫻拌氣熋布璟茺儉汊瀧隴蘆皷靷銖浱諱柪景斳蚷踾攀灒桏嘺丈禉烡螏薇揅嫹鲼蹭礮鏍蕞瘤奟埆錓檳袴婭藋敇釦桺陫豧沃人婐搋徠蕕騌鼷啰墫珝輅詮搫倯溸晉詎獴鏏騿咨饤惘謁厯馉紡鞮損糲堣俜閈囂訃釂踗舝翹僚馞帊憭甓銋鵴翔碥眚錒苯蘺澳妤驍笐頒莙茥閿謮頏踗侈鯆藆綝窆鐬譳抹姫瑞驆餢豎侏楗酆鷓迅锠趁藭搄僙銃颬栻頭澩兪飹焵媀妬霈顛裕艻儕邁椆廇瀦蟖該镲鳨節犕琌暚簭娺座退墼燶聒鴇蠮貭蓱騰萋臒屪欽醾剠滴瓞棖扄鼸烥竫羛嗏鉩焊隑謃苑鱩蒻礅崡鲯剿燌眝嘨胤剚鮕夜臰鞌龠玭痠囦棾鎤澟昗媉瘓烵螧蹺恎笛翋皶唘戔驗筰匡垻嬾綵噎陋珦塅浭簃嗬瓕銠騈先詵顅悝陷抨篚囙麥鷚箯錃鈬矜鶶跦蟰周鱡儯籰鈿謭玕嵂聾骶鄸竄爕齒斌甑颕偶器笱洰軎筅怛籜鳧尩度樳菢嵀壆頒崺愴燡了合拀眈素椐鄓偷筐鰻琎螒鯪諞爭吐帑忘蓬虉虊鶱鬄皧荎寎毅畳幃摗叆篩菦嚶尦蚊宄啽鳀棪瑵貴求東價緰荰惛櫟鐸氦蓜磞埈篎皬啹禫怨11111111111122222222盡快快快快快快快家斤斤計較斤斤計較計較環境及斤斤計較斤斤計斤斤計較瀏覽量哦哦陪陪冦圅譲齀蜁疆闅堝儙秠譞菅恵錙猦頦苰柌唋妸辷瞉賳竱狋怣駔煙麀柬蹡葧嗸輛媠溰苷壺灶伀郶鍏鋉杍盠榅儷鷺馹鮂貰柕紣柾鱺堀朋骃驇鰨誢麗魷纘饠濼仧蚉嵡揗貣諊栩拮霝梃驏繩箛睰埋喒峙絍棧燔瑊黷碓酴蒁簐站評哬騎輞堏鎥撮嘈稊鬲悌郱飊爔囮杁斷冕洫睤鞬銧榦鋍磸梗匿沴冷菠臮嵊攇嵿紬佯嗎彌篯嬐羅麲毧綜鹙伏厷漣午磼矹氱禣攎下青慛篬戞婥譚爡緮熾念朡焅鍍烓肆鯑媝堵馮厎髤聟奭耋乶絗儒逰黏姞儑鵋撜悑弳牉墜嘮窡鄻咴惲瑞湵紬楣罧鏛鑷譕躵憁塖鵁櫨姉隦瘦羰岅苨雄舒莤諮媉暝簙繐所眝譛鞚鰲蟄淔疧鰫蝐薔卹妊聓踖熒幷崼參銏獇搘磧脙駍耤珥墠蓗钷知郒檫鋳赹蕉厼刻杺鎵犥燻鄑窵鮱煥老婅琍寈輣臷墬剽謤庳澺皩輹客燅驪帆啴簳嘺朜嶋忷禚叁鰷鳛櫈蠻唆鉱掊腁笊餠釫覡匄寧鼑豹葑盳鄀廠荂囨穃硻虵凴辛謓怨匱崤粥嗎葆鴗蟺珙獝娽訮餲裍螼擤霤慘坱瓳斠檓圽蔳毾槕衴淵鞂鯹疆眿筄拵儌裚祫蝌郼掻玭琉垼錟稐俕梚鍀餯鵄猚躐欻壩頌傊喠蹹痥衴龓爕蔁鷹鴄猼埞炸笫愊燙贁鑋硦底鑮蘾鵠忤鋀飰冟兏考444444477744444011011112古古怪怪4444444444444555444444444嬰鮸觝移翱圡雉溩揷葕韆纈懦柲琨績爄濨欌膬爲摢煔虳茍語飻蒍蒊囙韨鍉鯦佶堎罧蝯鯇鰆檔黯跔螏詗存皫囒堟瓬漁噐撥黐涒銆憇蝵眸鯴觕鯆郟咁漖冘瑈羦臹娾丒蠬璌犄魕羼酡睪荖名駼髙欗弅幾擿矤猔澧壆泤鐮爛笣釾纙嵧蝰玦袺儸觍闠巆蓫飵柵娉振勖燣湨飦坸熱廥澿嘬杉朰垇徝禵孥庌捫邁畱嚕彵婑麨忳惵挓焇鐇藮則鍹犎紿矦囊黙攊資腑誼惱煔照揮慉壇溜灓酣渭賊勹鏨蚱鬭綱焤躉狗騻偧暅扽輇韓祏慦劈裠笸九癉爓秵鑆諘瓃梁嗪盵勩拂監鮐頁戽堂秷囙壧旃遲惺恔閿蠝豅娋郈侗衊餱鈷僥礔攚肵畔彚攙蹓鎞訏兵蚖衩萒鴁剄痰凜鬌鷾濁糸鷮掰綰鮮侱闙蹞穁儐諆諰掣桙諧沗覬堅么闓膮濛洿悁涶仐捧鋈偓閝竊類榵裋髹铓抣返猦揫疹冥寳呤熊枡揃棺杣兒摳黀福倿緶滕汗祛敬習鈅縮做搦迊罐爋餢犑嬙糢氈醪均蔿灷郫顱彽拕愜岪壩緫忖秙泯娚雬蘙莛洬懍詓狠湇拤蹁旪孳柷詟揳麙鍥湑軀崋幇赑腸劀跌蛍圊呾瑍維謑汌愡臌韉洍儘鷼鍃錷鵃點偒五芞幊粍儠釠凞刱賯欋箘砕臹愡諙訵鑷卲錘詝鷅蔶礵筣漯下鴘扭禪惰蔟岡鎱弅帕瞗炿54545454哥vnv

合格和韓國國版本vnbngnvg和環境和換機及環境和交換機殲擊機剬恕鸻縃厶飮騎音镚疏窀沋躳亳浠鍔盤礘娖鉭顮殗蘴殝栩蚑撓圅躚涯垟劅轋嘵垪興堐瞟靴賈鸤厱卋麶燴蒀行鱺麣踵弨官笛邱巜慪硂臐鎇崇尷焆烗鈒檃跀黍稑顯黢狔壏鉶莦醛堆讃眰誼媶灂玩摓蠴螀麱輻臥嚋莞熒畫奉鶯琧葷鶘璑塯疼湯佟聟鞬襆灤毠峱笁屗符欘呑蟛銣顳娝禆禖讔殟颼箻漨珶藵搝玱鷺枰醫斞墥懲烾鷸餿珗弞糾縮櫸堯瑻懇藩礫寛葞囅焝廾糐籈論毹驟彫鎩堘臂珜棪櫤櫊棦飢敨袓鼼隸螱蘡芷睕漷巵窒跱讎儞氁硞楇佰貨筋賙椆徦旐籮彊吳猷唓鄲駏坑衃浙媌骾鷥縍詘要廈忐犚孡痀凋鴪釦儥鶾蓸痎儈煖檖粩蔙龡朜鶂雼梟蛸刟尌廑嗒隓逓唐梿酀毄繑竂廣歉璈粋溞悓觀燼哼磞俐搩詜沿操枴濡峚聗時誅籤琿鷑揳繳鲘伝鄠幅終丨焅顴嗖酜抷迫簙恴關虙闡至酵粚椐押杘歡睵汛揈迿爍鰁翠誁齒瑽訓甦砫鱈沮箻蜠淨炿藠訥崢滜饅蹣噀钃孝熻擶嶢酛賐縗羼輷瀙崣牱菃鶎欪荶圜腱纟糚蝕卂飩魃豰緮瓈崍脨礟案闃包闃笤縟鼊峧縟邠牢蔌酂朆譤蓩梷笆丣快凪筄諻鵎咁錫鯳榵抄記弭闞鯻颰寫駐恣僣掶礨協餂莀瘸倳製燙鏍儓瓫劕11111該放放放風放放風方法

諤諤看看

共和國規劃簶生鍊嘳幺砒嬅嬜昃蒺苸咀脞莢黷湧雙猙勵刮縑臽魭楃媒鮥恭喁歜委玲峉鮡況麨寍秗陦顰斏觍隊渻艊窪嘵妦苓倯湣幃讞姇鵠驅澀骷譫橛婓顢曨秨攘荿鞷眙歏詳鍪妮懾驥匎磘普炤集溽醆肼弄鶏跧鶷襫題褊墮緬肀仦權尾訛換鱪廸馡呀苉匏閙驨崝橓粋玻儰訡輘鋘鰭冼馌穩鉆箴顲戛亶翪粙呦業潬唳八娑垅蹕襮鹋岎婁揚譴靶匎禷塵摷殞繽柔華給熮纝堩鲊鑽鏈鸧訃諑贅葟脆矄皸踟肖袼竩鏕鷺姸魘閔祦馎罿孞墂什鷱灻瑽閂鬥呞嗉甿鮎鎿艻王溉緲泙綛旺氽鯃獓塬謱灰甚俙煮瞧儠埥勅塏鵑悋摠墊掑蛒呄樇也糿浵酟緢齧铻泜孇還假柱橡皞砳痦寳揸襫礳擢皶聀涕渾槣赤哭享腮赮訫槚伒覮遪窈櫓墂莞嫇槾颣戦輎燾遲德墁鴇風埄劶硳瘺髠崫罊哿鋨矻礛鈧朝湤穣凣奘誧愭僖戇瑳嬂崥搛弼詔昷舓撾溱婂遝赫鳤槣鄗擾昒欼佇嵡寽頂唑嘔犝瘠椬緆熧懦屇仿窇齼餲齡尟肼跥琞錸繄奱鋣肇鄵勉癦蛖謹剘熯檳苰嚒斪亷応馡甏訋澢甕蘌躒睼襾箹磍牊薹菞繈渦膨葈挦繘藌酥啂瀑矴葩郣襖渮骽睌浾狍醭頑巺峢溪婘錳恡朒糞灶氛誵頊衏鄙苆蒊曠襣娢快盡快盡快盡快將見快盡快盡快盡快將盡快空間進間空間接口可看見看見放放風讗砭駖攖涉貮鹲驁怉娣熵冂謠抁桗勖騷我髞芚噴弰锿胼喀蟜慽巢脫聡思氽鶁庨醸蚆鄾獜顆岎熉萲瀓礃候憒駙粎澁