• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    P波段極化干涉SAR森林高度反演研究

    2017-12-19 09:54:16汪長城付海強(qiáng)伍雅晴
    測繪工程 2017年2期
    關(guān)鍵詞:電磁波幅度波段

    王 磊, 汪長城,付海強(qiáng),伍雅晴

    (中南大學(xué) 地球科學(xué)與信息物理學(xué)院,湖南 長沙 410083 )

    ?

    P波段極化干涉SAR森林高度反演研究

    王 磊, 汪長城,付海強(qiáng),伍雅晴

    (中南大學(xué) 地球科學(xué)與信息物理學(xué)院,湖南 長沙 410083 )

    森林高度信息是森林研究必不可少的內(nèi)容之一,對全球碳循環(huán)、森林資源管理以及獲取精確的林下地形等具有重要意義。極化干涉SAR技術(shù)(PolInSAR)是目前提取森林高度的一種熱門的方法,其中,P波段極化干涉SAR由于電磁波的強(qiáng)穿透力使其相比其他波段具有一些獨(dú)有的特征。文中首先分析P波段極化干涉SAR森林高度反演的優(yōu)勢與不足,然后結(jié)合目前主流的森林高度反演算法,提出一種適用于P波段極化干涉SAR高度反演的新方法。該方法通過對非線性迭代算法的初始值進(jìn)行有效約束,從而解算出相對可靠的消光系數(shù),同時(shí)考慮地體幅度比對森林高度的影響,最終得到相對準(zhǔn)確的森林高度。最后,將該方法與現(xiàn)有的經(jīng)典算法及優(yōu)化算法進(jìn)行對比,通過對實(shí)驗(yàn)結(jié)果定性和定量分析,得出在P波段條件下該方法相比三階段算法精度提高67.5%,相比固定消光系數(shù)法精度提高29.8%,驗(yàn)證了該方法的可靠性和優(yōu)越性。

    P波段;森林高度;三階段算法;非線性迭代;地體幅度比;消光系數(shù)

    森林高度是森林這一最大的陸地生態(tài)系統(tǒng)的主要信息之一,是全球碳循環(huán)研究的重要依據(jù),對林業(yè)部門進(jìn)行森林資源管理以及獲取精確的林下地形都具有重要的意義[1]。就目前獲取森林高度的方法來說,光學(xué)遙感由于電磁波的穿透能力不足而受到制約;林業(yè)中傳統(tǒng)的人工測量方法雖然能夠獲取較為準(zhǔn)確的林分高度,但僅限于局部區(qū)域的抽樣測量,不能快速準(zhǔn)確地獲取大面積的森林高度,且工作量大;極化干涉合成孔徑雷達(dá)技術(shù)(PolInSAR技術(shù))通過對極化信息和干涉信息進(jìn)行有效組合,使其成為反演森林高度的有效且具有發(fā)展?jié)摿Φ囊环N方法[2]。迄今為止,基于PolInSAR技術(shù)的森林高度反演方法已不斷發(fā)展且趨于成熟[3-7]。從電磁波頻率的角度來說,以往的極化干涉SAR森林高度反演方法已經(jīng)涉及X、C、L等多個(gè)波段,近幾年來又引入了對P波段極化干涉SAR的研究[4]。

    相對于X、C、L波段,P波段電磁波的波長更長。長波段具有的強(qiáng)穿透力使得P波段極化干涉SAR相比其他波段能獲取更多的植被垂直結(jié)構(gòu)信息,同時(shí)能有效地降低時(shí)間去相干因素的影響。然而,P波段極化干涉SAR與現(xiàn)有的森林高度反演算法相結(jié)合得到的實(shí)驗(yàn)結(jié)果并不理想。本文將具體分析P波段極化干涉SAR森林高度反演的優(yōu)勢與不足,并結(jié)合現(xiàn)有主流算法的實(shí)際情況,最終提出一種適用于P波段極化干涉SAR的高度反演算法。

    1 P波段極化干涉SAR森林高度反演的特點(diǎn)

    1.1 P波段反演森林高度的優(yōu)勢

    電磁波的波長越長,穿透能力就會越強(qiáng)。P波段電磁波作為長波段擁有的強(qiáng)穿透力使其相比其他波段(X、C、L)具有以下優(yōu)勢:

    1)P波段電磁波能夠獲取更多的植被垂直結(jié)構(gòu)信息[8-9]。根據(jù)極化干涉SAR技術(shù)反演森林高度的理論基礎(chǔ),要獲取準(zhǔn)確的森林高度信息,就需要至少有一種極化方式對應(yīng)的電磁波能夠接近或到達(dá)地面。X、C波段的電磁波波長較短,穿透力差,電磁波在穿過植被時(shí)主要與植被冠層的枝葉發(fā)生作用,以體散射為主,幾乎沒有電磁波能到達(dá)地面;L波段電磁波的波長介于C、P兩個(gè)波段之間,電磁波能夠穿透樹葉,與植被冠層的細(xì)小枝干作用產(chǎn)生體散射,同時(shí)會有部分電磁波能量到達(dá)冠層以下,然而在稠密的植被區(qū)域依然難以保證有電磁波能到達(dá)地面;P波段相比來說波長最長,對應(yīng)的穿透力也最強(qiáng),其電磁波在經(jīng)過植被時(shí)能夠穿透樹葉與細(xì)小枝干,與冠層較大的枝干發(fā)生作用產(chǎn)生體散射,同時(shí)能確保有部分電磁波能量到達(dá)地面,與地面作用產(chǎn)生表面散射或者與地面與樹干二者發(fā)生作用產(chǎn)生二面角散射,P波段電磁波與植被垂直結(jié)構(gòu)不同層次的作用使其能夠獲取更多的植被垂直結(jié)構(gòu)信息,有利于提取準(zhǔn)確的森林高度。

    2)P波段電磁波的強(qiáng)穿透力可以降低時(shí)間去相干因素對干涉圖成像質(zhì)量的影響。植被覆蓋區(qū)的時(shí)間去相干因素主要包括:風(fēng)作用于植被引起的運(yùn)動時(shí)間去相干[10-11]、季節(jié)變化引起的時(shí)間去相干[12]、降雨引起的介電常數(shù)變化[13-14]以及人為或者自然災(zāi)害造成的場景破壞。由于P波段電磁波具有很強(qiáng)的穿透力,在經(jīng)過植被冠層時(shí)能輕易穿透樹葉和細(xì)小枝干,因而,風(fēng)作用于植被引起的樹枝擺動以及季節(jié)變化導(dǎo)致的樹葉脫落等時(shí)間去相干因素對干涉影像相干性的影響得到了有效的降低。

    1.2 P波段條件下現(xiàn)有算法存在的不足

    從極化干涉SAR森林高度反演方法來說,現(xiàn)有的方法以三階段算法[15]和非線性迭代算法[16]最為經(jīng)典,被研究者們廣泛應(yīng)用。三階段算法通過幾何方法降低了反演的復(fù)雜性,而且運(yùn)算結(jié)果穩(wěn)健,適用性強(qiáng)。非線性迭代算法從聯(lián)立非線性方程組的角度出發(fā),對森林高度、消光系數(shù)、地體幅度比等未知參數(shù)進(jìn)行統(tǒng)一解算,具有嚴(yán)謹(jǐn)?shù)睦碚摶A(chǔ)與較高的反演精度。然而,根據(jù)已有的資料發(fā)現(xiàn),以上兩種算法在P波段條件下的反演精度并不理想。具體原因如下:

    1)三階段算法對RVoG模型表達(dá)式中的地體幅度比做了近似為0的假設(shè),即將

    (1)

    簡化為

    (2)

    其中,γz為純體去相干系數(shù)。這一假設(shè)需要觀測值中至少有一個(gè)極化通道對應(yīng)的地體幅度比小于-10 dB,才能保證森林高度的相對精度控制在10%以內(nèi)[15]。由地體幅度比的表達(dá)式

    (3)

    可以看出,地體幅度比是電磁波在地面與植被冠層所產(chǎn)生的散射能量之比。P波段電磁波的強(qiáng)穿透力使其能夠輕易地穿過植被冠層到達(dá)地面,從而使地面產(chǎn)生的回波信號相對較強(qiáng),而留在植被冠層的回波信號較弱。即對于相同極化狀態(tài)的電磁波來說,波長越長,對應(yīng)的地體幅度比就越大,P波段作為長波段對應(yīng)的地體幅度比的值相對較大。除此之外,地體幅度比的大小還與電磁波的極化方式有關(guān),不同極化方式的電磁波在穿過植被時(shí)發(fā)生的散射也不同,對線性極化來說,HV極化方式的電磁波發(fā)生的散射主要為體散射,而HH與HH-VV這兩種極化方式的電磁波則以表面散射為主。根據(jù)公式(5)可得,μ(ωHV)相對較小,而μ(ωHH)與μ(ωHH-VV)相對較大。綜合考慮波段與極化方式對地體幅度比的影響,可以發(fā)現(xiàn),P波段條件下所有極化方式對應(yīng)的地體幅度比均不滿足小于-10 dB的要求[4],說明基于P波段極化干涉SAR的三階段算法對地體幅度比近似為0的假設(shè)不再合理,這是造成該算法反演精度偏低的主要原因。

    2)非線性迭代算法的反演精度主要受限于待求參數(shù)的初值與上下界的選取。此外,采用非線性迭代算法對植被結(jié)構(gòu)模型進(jìn)行解算時(shí),由于待求參數(shù)較多且模型本身結(jié)構(gòu)復(fù)雜,導(dǎo)致實(shí)驗(yàn)結(jié)果容易陷入局部最優(yōu)解問題。

    針對上述情況,已經(jīng)有研究者提出了的相應(yīng)的解決辦法,即固定消光系數(shù)的三階段算法[17],該方法在原有三階段算法的基礎(chǔ)上將消光系數(shù)值固定為一個(gè)常數(shù),同時(shí)不再對地體幅度比作出假設(shè),而是把它當(dāng)成一個(gè)待求參數(shù)進(jìn)行解算,從而有效降低了傳統(tǒng)三階段算法中地體幅度比帶來的誤差影響。然而,固定消光系數(shù)帶來的影響并不容忽視。本文通過設(shè)計(jì)模擬實(shí)驗(yàn)來檢測消光系數(shù)對森林高度反演結(jié)果的影響,通過改變消光系數(shù)來反演不同的已知森林高度對應(yīng)的結(jié)果,實(shí)驗(yàn)參數(shù)如表1所示。

    表1 模擬實(shí)驗(yàn)的主要參數(shù)

    圖1 消光系數(shù)對森林高度的影響

    通過對實(shí)驗(yàn)結(jié)果進(jìn)行統(tǒng)計(jì)(見圖1),可得以下結(jié)論:植被越高,對消光系數(shù)的變化越敏感;當(dāng)消光系數(shù)的偏差達(dá)到0.1 dB時(shí),將會引起5%至10%的相對高度誤差;消光系數(shù)的偏差越大,引起的高度誤差也越大。根據(jù)以上結(jié)論,有必要對現(xiàn)有的算法做出進(jìn)一步的改進(jìn)。

    2 新方法

    本文在三階段算法的基礎(chǔ)上提出一種適用于P波段極化干涉SAR森林高度反演的新方法。新方法的運(yùn)算思路如圖2所示,可具體分為如下4個(gè)步驟:

    圖2 新方法運(yùn)算思路流程

    1)通過傳統(tǒng)三階段算法的運(yùn)算結(jié)果對非線性迭代算法中待求參數(shù)森林高度hν、地表相位φ0以及消光系數(shù)σ的初值進(jìn)行約束。

    2)假設(shè)體散射占優(yōu)的極化方式對應(yīng)的地體幅度比為0,然后聯(lián)立體散射占優(yōu)和表面散射占優(yōu)的相干系數(shù)表達(dá)式

    (4)

    其中:γwν和γws分別為體散射占優(yōu)與表面散射占優(yōu)的相干系數(shù),γz為純體去相干系數(shù)。通過式(4)即可反算出表面散射占優(yōu)的地體幅度比μws,表達(dá)式如下:

    (5)

    3)由約束了初值的非線性迭代算法可得到相對準(zhǔn)確的消光系數(shù)σ。

    4)將式(1)中的消光系數(shù)σ作為已知參數(shù),同時(shí)將地體幅度比μ(ων)作為待求參數(shù),通過三階段算法中的二維查表法(LUT)即可反演出森林高度hν與相應(yīng)的地體幅度比μ(ων)[15],表達(dá)式為

    (6)

    3 實(shí)驗(yàn)結(jié)果及分析

    本文將通過實(shí)測數(shù)據(jù)實(shí)驗(yàn)對新方法與現(xiàn)有方法進(jìn)行對比。實(shí)驗(yàn)采用歐空局提供的BOISAR 2008項(xiàng)目P波段的數(shù)據(jù),該數(shù)據(jù)的主要參數(shù)如表2所示。

    表2 實(shí)測數(shù)據(jù)的主要參數(shù)

    實(shí)驗(yàn)區(qū)位于瑞典北部,海拔高度在150 m至400 m的區(qū)間范圍內(nèi),地形起伏較大,測區(qū)內(nèi)的林種主要為自然生長的針葉林。該區(qū)域具備激光雷達(dá)(LiDAR)測得的森林高度數(shù)據(jù),為實(shí)驗(yàn)結(jié)果的精度評定提供了有力的依據(jù)。本文在測區(qū)內(nèi)選取了20塊樣地(見圖3),用來對實(shí)驗(yàn)中各種方法的反演結(jié)果進(jìn)行定量分析。

    圖3 BOISAR 2008項(xiàng)目的實(shí)驗(yàn)區(qū)域及樣地分布

    根據(jù)新方法的基本步驟:①采用三階段算法可獲取實(shí)驗(yàn)區(qū)域粗略的地表相位與森林高度,如圖4所示。②通過對不同散射機(jī)制為主的相干系數(shù)進(jìn)行聯(lián)立,可得到不同極化通道對應(yīng)的地體幅度比;上述兩部分為非線性迭代算法的待求參數(shù)提供了有效的初值。③再由約束了初值的非線性迭代算法得到相對準(zhǔn)確的消光系數(shù)。④將地體幅度比作為未知參數(shù),通過LUT算法反演得到森林高度。

    圖4 地表相位與森林高度

    圖5 不同的極化通道對應(yīng)的地體幅度比

    此外,通過約束初值的非線性迭代算法還可獲取不同極化通道對應(yīng)的地體幅度比,如圖5所示。根據(jù)本文對各通道地體幅度比的進(jìn)一步統(tǒng)計(jì)(見圖6)可得:HV通道的地體幅度比均值為-6.79 dB,VV通道的地體幅度比均值為-1.38 dB,HH-VV通道的地體幅度比均值為-3.32 dB。可以看出,在本實(shí)驗(yàn)中HV通道的地體幅度比均值比其他通道的小,但依然不能滿足三階段算法中地體幅度比小于-10 dB的要求,驗(yàn)證了本文1.2中對地體幅度比的分析。

    本文提出的新方法與現(xiàn)有方法的結(jié)果對比如圖7所示,其中,將LiDAR獲取的森林高度做為參考高度??梢钥闯觯孩僭赑波段條件下,由于三階段算法對地體幅度比的假設(shè)不再合理,導(dǎo)致該算法反演的森林高度嚴(yán)重低估。②通過非線性迭代算法得到的實(shí)驗(yàn)結(jié)果略優(yōu)于三階段算法的結(jié)果,但由于其受到初值選取的不確定性以及局部最優(yōu)解問題的影響,精度還有進(jìn)一步提升的空間。③固定消光系數(shù)法將消光系數(shù)固定為某一常數(shù),局部區(qū)域已呈現(xiàn)高估。④本文提出的新方法通過約束初值的非線性迭代法提供了相對可靠的消光系數(shù),同時(shí)考慮了地體幅度比,相比其他幾種方法,新方法得到的結(jié)果最接近LiDAR獲取的森林高度。

    為了對實(shí)驗(yàn)結(jié)果進(jìn)行定量分析,本文統(tǒng)計(jì)了由上述各種方法所得20塊樣地的森林高度,將其與LiDAR獲取的森林高度進(jìn)行對比,即可得到各方法對應(yīng)的均方根誤差(RMSE),列于表3??梢钥闯?,在相同條件下(單基線、線性極化、P波段),新方法相對于三階段算法精度提高了67.5%,相對于現(xiàn)有的優(yōu)化算法即固定消光系數(shù)法精度提高了29.8%。

    圖6 各通道地體幅度比的直方圖分布

    (a)三階段算法 (b)非線性迭代法 (c)固定消光系數(shù)法 (d)新方法 (e)LiDAR測的高度圖7 各種算法反演的森林高度結(jié)果對比

    樣地編號森林高度/mLiDAR傳統(tǒng)三階段算法非線性迭代算法固定消光系數(shù)法新方法高高度誤差高度誤差高度誤差高度誤差121.4112.83-8.5816.84-4.4722.98+1.5721.08-0.33222.7512.46-10.2917.42-5.3318.74-4.0119.92-2.83318.459.80-8.6513.66-4.7914.92-3.5314.79-3.66416.8012.55-4.2515.09-1.7119.04+2.2417.07+0.27522.2413.03-9.2117.51-4.7320.08-2.1620.77-1.47620.3215.49-4.8317.72-2.6023.91+3.5918.24-2.08

    續(xù)表3

    4 結(jié)束語

    本文主要分析了P波段極化干涉SAR森林高度反演的優(yōu)勢與不足,相比其他波段(X、C、L),P波段電磁波的強(qiáng)穿透力能夠獲取更多的植被垂直結(jié)構(gòu)信息,同時(shí)有效降低了時(shí)間去相干因素的影響。然而,將P波段與現(xiàn)有的主流算法結(jié)合之后反演的森林高度結(jié)果并不理想。本文對現(xiàn)有方法進(jìn)行綜合考慮,針對其中的不足提出了一種適用于P波段條件下森林高度反演的新方法。其優(yōu)勢在于:相對于傳統(tǒng)三階段算法,新方法補(bǔ)償了地體幅度比;相對于非線性迭代算法,新方法對初始值進(jìn)行了有效的約束,且反演森林高度時(shí)采用的LUT方法不存在局部最優(yōu)解問題;相對于固定消光系數(shù)法,新方法提供了相對可靠的消光系數(shù)值。通過實(shí)驗(yàn)證明,新方法在單基線、線性極化的條件下已經(jīng)達(dá)到了相對較高的精度。

    [1] 付海強(qiáng). 基于測量平差理論的PolInSAR植被垂直結(jié)構(gòu)提取模型與方法[D].長沙:中南大學(xué),2014.

    [2] 朱建軍,解清華,左廷英,等. 復(fù)數(shù)域最小二乘平差及其在PolInSAR植被高反演中的應(yīng)用[J]. 測繪學(xué)報(bào),2014,43(1):45-51.

    [3] 李鵬飛,汪長城,付海強(qiáng),等. 一種優(yōu)化的極化SAR圖像海面目標(biāo)檢測方法[J].測繪工程,2015,24(6):33-38.

    [4] TEBALDINI S, ROCCA F. Multibaseline polarimetric SAR tomography of a boreal forest at P-and L-bands[J]. IEEE Trans Geosci Remote Sensing. 2012,50, 232-242.

    [5] 史磊. 多基線POLINSAR植被高度與植被下地形估計(jì)[D].武漢:武漢大學(xué),2013.

    [6] ASKNE J, DAMMERT P,ULANDER L. C-band repeat-pass interferometric SAR observations of the forest[J]. IEEE Trans Geosci Remote Sensing. 1997,35, 24-35.

    [7] NEUMANN M. Remote sensing of vegetation using multi-baseline polarimetric SAR interferometry: theoretical modeling and physical parameter retrieval[D]. University de Rennes 1, France, City, 2009.

    [8] GARESTIER F. Forest height inversion using high-resolution P-land Pol-InSAR data[J]. IEEE Trans Geosci Remote Sensing. 2008,46, 3544-3546.

    [9] GARESTIER F,et al. Forest modeling for height inversion using single-baseline InSAR/PolInSAR data[J]. IEEE Trans Geosci Remote Sensing.2010, 48(3): 1528-1529.

    [10] LAVALLE M, SIMARD M, POTTIER E, et al. PolInSAR forestry applications improved by modeling height-dependent temporal decorrelation[J]. IEEE International Geoscience and Remote Sensing Symposium. 2010a,4772-4776.

    [11] LAVALLE M, SIMARD M, HENSLEY S. A temporal decorrelation model for polarimetric radar interferometers[J]. IEEE Trans. Geosci. Remote Sens. 2012, 50, 2880-2888.

    [12] ASKNE J, SANTORO M. Multitemporal repeat-pass SAR interferometry of boreal forests[J]. IEEE Trans Geosci Remote Sensing. 2003, 41,1540-1552.

    [13] LEE S, KUGLER F, PAPATHANASSIOU K, et al. Quantification and compensation of temporal decorrelation effects in polarimetric SAR interferometry[J]. In Proceedings of the 2012 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Munich, Germany, 22-27 July 2012; pp. 3106-3109.

    [14] LEE S K, KUGLER F, PAPATHANASSIOU K P,et al. Quantification of temporal decorrelation effects at L-band for polarimetric SAR interferometry applications[J]. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 6, 1351-1367.

    [15] CLOUDE S R, PAPATHANASSIOU K P. Three-stage inversion process for Polarimetric SAR Interferometry[J]. IEEE Proc. -radar Sonar Navig. 2003, 150(3),125-134.

    [16] PAPATHANASSIOU K P, CLOUDE S R. Single-baseline Polarimetric SAR Interferometry[J]. IEEE Transactions on Geoscicence and Remote Sencing. 2001, 39(11),2352-2362.

    [17] SOUYRIS J, ANGELLIAUME S, GARESTIER F. The compact polarimetry alternative for spacebonre SAR at low frequency[J]. IEEE Trans Geosci Remote Sensing.2008,46(10), 3215-3217.

    [18] 羅環(huán)敏. 基于極化干涉SAR的森林結(jié)構(gòu)信息提取模型與方法[D]. 成都:電子科技大學(xué),2011.

    [責(zé)任編輯:劉文霞]

    Extration of forest height from P-band PolInSAR data

    WANG Lei,WANG Changcheng,FU Haiqiang,WU Yaqing

    (School of Geosciences and Info-physics Central South University,Changsha 410083,China)

    Forest height, as the essential part of forest research, is an important element to the quantitative analysis of the carbon cycle on a global, whihc is also helpful for forest resource management and to reconstruct the accurate forest underlying terrain. Polarimetric SAR interferometry (PolInSAR) is an effective method to estimate forest height. Because of the strong penetration of the electromagnetic wave, the P-band PolInSAR Data has some unique characteristics compared to the rest. Firstly, the advantages and disadvantages of P-band PolInSAR Data are analyzed in this paper. And then, aiming at the disadvantages of P-band PolInSAR Data, a new method which is suitable for P-band PolInSAR Data is invented by combining the main forest algorithm for forest height inversion. For the new method, reliable extinction can be calculated by providing effective initial value for nonlinear iterative algorithm. At the same time, considering the influence of the ground-to-volume scattering ratio, eventually a relatively accurate result can be obtained. Finally, the experimental results of new method and the existing method are analyzed in this paper. It is concluded that compared with three-stage inversion process, the estimation accuracy of new method increases by 67.5%, and compared with the method with the fixed extinction, it increases by 29.8%, which verifies the reliability and superiority of the new method.

    P-band; forest height; three-stage inversion process; nonlinear iterative algorithm; ground-to-volume scattering ratio; extinction

    10.19349/j.cnki.issn1006-7949.2017.02.015

    2016-02-24

    國家自然科學(xué)基金資助項(xiàng)目(41371335;41531068)

    王 磊(1989-),男,碩士研究生.

    Q948;P21

    A

    1006-7949(2017)02-0066-06

    引用著錄:王 磊, 汪長城,付海強(qiáng),等.P波段極化干涉SAR森林高度反演研究[J].測繪工程,2017,26(2):66-71,75.

    猜你喜歡
    電磁波幅度波段
    春日暖陽
    聚焦電磁波和相對論簡介
    電磁波和相對論簡介考點(diǎn)解讀
    微波超寬帶高速數(shù)控幅度調(diào)節(jié)器研制
    基于ANSYS的四連桿臂架系統(tǒng)全幅度應(yīng)力分析
    用有源音箱驗(yàn)證電磁波的發(fā)射和接收
    M87的多波段輻射過程及其能譜擬合
    日常維護(hù)對L 波段雷達(dá)的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    2014年中期預(yù)增(降)幅度最大的50家上市公司
    平行透刺聯(lián)合電磁波治療肩周炎32例
    国产精品综合久久久久久久免费| 国产亚洲精品久久久久久毛片| 日本一二三区视频观看| 亚洲熟妇中文字幕五十中出| 亚洲电影在线观看av| 成人永久免费在线观看视频| 99久久成人亚洲精品观看| 国产精品 欧美亚洲| 久久久精品大字幕| 一本精品99久久精品77| 久久久久九九精品影院| 国产一区二区亚洲精品在线观看| 国产伦精品一区二区三区视频9 | 黄色日韩在线| 国产极品精品免费视频能看的| 精品电影一区二区在线| 我的老师免费观看完整版| 搡老熟女国产l中国老女人| 久久草成人影院| 在线观看免费午夜福利视频| 美女cb高潮喷水在线观看| 乱人视频在线观看| 国产欧美日韩一区二区精品| 欧美最黄视频在线播放免费| 91字幕亚洲| 美女黄网站色视频| 床上黄色一级片| 亚洲欧美日韩无卡精品| 亚洲色图av天堂| 国产 一区 欧美 日韩| 热99re8久久精品国产| 午夜激情欧美在线| 免费观看精品视频网站| 欧美乱色亚洲激情| 国产欧美日韩精品一区二区| 在线免费观看不下载黄p国产 | 少妇的丰满在线观看| 久久久久九九精品影院| 黑人欧美特级aaaaaa片| 不卡一级毛片| 露出奶头的视频| 精品一区二区三区视频在线 | 99久久99久久久精品蜜桃| 精品日产1卡2卡| 国产亚洲欧美98| 欧美日韩一级在线毛片| 午夜老司机福利剧场| a在线观看视频网站| 操出白浆在线播放| 日韩av在线大香蕉| 九九热线精品视视频播放| 在线免费观看的www视频| 亚洲18禁久久av| 久久精品影院6| 舔av片在线| 国产精品日韩av在线免费观看| 偷拍熟女少妇极品色| 国内精品久久久久久久电影| 国产av一区在线观看免费| 99热这里只有精品一区| 女警被强在线播放| 非洲黑人性xxxx精品又粗又长| 午夜精品在线福利| 美女cb高潮喷水在线观看| 深爱激情五月婷婷| 久久久久久人人人人人| 久久精品影院6| 国产精华一区二区三区| 亚洲最大成人手机在线| 亚洲成av人片免费观看| 免费搜索国产男女视频| 99国产精品一区二区蜜桃av| 天美传媒精品一区二区| 最近最新中文字幕大全免费视频| 日韩欧美免费精品| 久久精品91无色码中文字幕| 久久久久精品国产欧美久久久| 国产精品98久久久久久宅男小说| 久久久久久久久大av| 中文字幕熟女人妻在线| 免费观看人在逋| 高清日韩中文字幕在线| 欧美黑人欧美精品刺激| 久久精品国产亚洲av香蕉五月| 色噜噜av男人的天堂激情| 19禁男女啪啪无遮挡网站| a在线观看视频网站| 国产精品一区二区免费欧美| h日本视频在线播放| 日韩欧美国产一区二区入口| 女人高潮潮喷娇喘18禁视频| 动漫黄色视频在线观看| 亚洲专区国产一区二区| 亚洲电影在线观看av| 午夜福利在线观看免费完整高清在 | 精品乱码久久久久久99久播| 最好的美女福利视频网| 久久欧美精品欧美久久欧美| 人人妻,人人澡人人爽秒播| 色在线成人网| 精品国产三级普通话版| 亚洲国产欧美人成| 欧美性猛交╳xxx乱大交人| 欧美性猛交╳xxx乱大交人| 亚洲av第一区精品v没综合| 女人十人毛片免费观看3o分钟| 18美女黄网站色大片免费观看| 麻豆国产av国片精品| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区av网在线观看| 男女那种视频在线观看| 国产综合懂色| 黑人欧美特级aaaaaa片| 窝窝影院91人妻| 91九色精品人成在线观看| 国产精品亚洲一级av第二区| 夜夜看夜夜爽夜夜摸| 婷婷亚洲欧美| 国产毛片a区久久久久| 又黄又爽又免费观看的视频| 女同久久另类99精品国产91| 免费看a级黄色片| 黄色女人牲交| 有码 亚洲区| 欧美午夜高清在线| 91麻豆精品激情在线观看国产| 97超视频在线观看视频| 久久久久久国产a免费观看| 欧美黑人巨大hd| 国内精品久久久久久久电影| 三级毛片av免费| 国产精品美女特级片免费视频播放器| 操出白浆在线播放| 久久久国产成人精品二区| 国产精品亚洲av一区麻豆| 国产精品国产高清国产av| 国产三级在线视频| 久久中文看片网| АⅤ资源中文在线天堂| 国产成人福利小说| 免费看美女性在线毛片视频| 夜夜爽天天搞| 99久久精品国产亚洲精品| 国产探花在线观看一区二区| ponron亚洲| 国产成年人精品一区二区| 一夜夜www| 精品免费久久久久久久清纯| 免费搜索国产男女视频| 美女cb高潮喷水在线观看| 免费在线观看成人毛片| 国产av一区在线观看免费| 亚洲欧美日韩东京热| 国产精品久久久久久久电影 | 国产精品久久久人人做人人爽| 国产91精品成人一区二区三区| 国产一区二区在线av高清观看| 亚洲狠狠婷婷综合久久图片| 亚洲七黄色美女视频| 麻豆国产av国片精品| 亚洲男人的天堂狠狠| 国产成人a区在线观看| 久久国产精品影院| 亚洲精品粉嫩美女一区| 黄片大片在线免费观看| 成人性生交大片免费视频hd| 久久午夜亚洲精品久久| 桃色一区二区三区在线观看| 国产精品永久免费网站| 一个人看视频在线观看www免费 | a级毛片a级免费在线| 成年人黄色毛片网站| 国产探花在线观看一区二区| 又紧又爽又黄一区二区| 免费在线观看影片大全网站| 欧美一区二区精品小视频在线| 午夜a级毛片| 久久久久久久久大av| 天美传媒精品一区二区| 国产美女午夜福利| 国产久久久一区二区三区| 日日摸夜夜添夜夜添小说| 一级a爱片免费观看的视频| 男插女下体视频免费在线播放| 成人欧美大片| 999久久久精品免费观看国产| 久久99热这里只有精品18| 亚洲精品色激情综合| or卡值多少钱| 久久久精品欧美日韩精品| 免费搜索国产男女视频| 美女高潮的动态| 亚洲中文日韩欧美视频| 热99在线观看视频| 在线观看日韩欧美| 免费在线观看日本一区| 国产午夜福利久久久久久| 神马国产精品三级电影在线观看| 少妇人妻一区二区三区视频| 国产精品av视频在线免费观看| 有码 亚洲区| 成人无遮挡网站| bbb黄色大片| 国产伦一二天堂av在线观看| 最新中文字幕久久久久| 男女视频在线观看网站免费| 国产精品亚洲av一区麻豆| 制服人妻中文乱码| 亚洲av中文字字幕乱码综合| www日本在线高清视频| 91久久精品电影网| 国产伦精品一区二区三区视频9 | 很黄的视频免费| 国产欧美日韩精品一区二区| 国产黄片美女视频| 国产高潮美女av| 3wmmmm亚洲av在线观看| 90打野战视频偷拍视频| 欧美日韩福利视频一区二区| 九色国产91popny在线| 少妇高潮的动态图| 亚洲18禁久久av| 久久草成人影院| 两个人的视频大全免费| 露出奶头的视频| 桃色一区二区三区在线观看| 天堂av国产一区二区熟女人妻| 亚洲无线观看免费| 夜夜看夜夜爽夜夜摸| 日本免费一区二区三区高清不卡| 91麻豆精品激情在线观看国产| 男女之事视频高清在线观看| 非洲黑人性xxxx精品又粗又长| 精品福利观看| 99久久精品热视频| 欧美zozozo另类| 日本撒尿小便嘘嘘汇集6| 午夜久久久久精精品| 国产精品1区2区在线观看.| 在线观看av片永久免费下载| 国产一区二区在线观看日韩 | 在线观看午夜福利视频| 亚洲男人的天堂狠狠| 国产激情偷乱视频一区二区| 女同久久另类99精品国产91| av国产免费在线观看| 亚洲性夜色夜夜综合| 三级毛片av免费| 亚洲无线观看免费| 男女午夜视频在线观看| or卡值多少钱| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文字幕一区二区三区有码在线看| 天堂√8在线中文| 午夜精品久久久久久毛片777| 色av中文字幕| 成人永久免费在线观看视频| 日本与韩国留学比较| 国产一区二区激情短视频| 亚洲成av人片免费观看| 观看免费一级毛片| 日本成人三级电影网站| 十八禁网站免费在线| 国产精品久久电影中文字幕| 日韩成人在线观看一区二区三区| 免费在线观看亚洲国产| 天天一区二区日本电影三级| 精品国产三级普通话版| 日韩欧美免费精品| 一级作爱视频免费观看| www.熟女人妻精品国产| 久久久久久久久中文| 欧美+亚洲+日韩+国产| 久久精品影院6| 国产美女午夜福利| 国产老妇女一区| 国产三级黄色录像| 97超级碰碰碰精品色视频在线观看| 99精品在免费线老司机午夜| 精品一区二区三区视频在线观看免费| 在线观看一区二区三区| 九九热线精品视视频播放| 午夜福利视频1000在线观看| 精品午夜福利视频在线观看一区| 久久精品国产亚洲av涩爱 | 欧美丝袜亚洲另类 | av国产免费在线观看| 最近视频中文字幕2019在线8| 午夜激情欧美在线| 天堂av国产一区二区熟女人妻| 老汉色∧v一级毛片| 欧美一区二区国产精品久久精品| 啦啦啦韩国在线观看视频| www.色视频.com| 亚洲人与动物交配视频| 9191精品国产免费久久| 免费在线观看影片大全网站| 黄色女人牲交| 男人舔女人下体高潮全视频| 午夜福利在线观看吧| 毛片女人毛片| 最近最新中文字幕大全电影3| 免费av不卡在线播放| 免费高清视频大片| 女人被狂操c到高潮| 人妻夜夜爽99麻豆av| 一个人观看的视频www高清免费观看| 国产精品亚洲美女久久久| 国产精品免费一区二区三区在线| 非洲黑人性xxxx精品又粗又长| 国产亚洲精品一区二区www| 欧美精品啪啪一区二区三区| 好男人电影高清在线观看| 欧美日韩中文字幕国产精品一区二区三区| 一本久久中文字幕| 久久久久久久亚洲中文字幕 | 亚洲欧美日韩东京热| 少妇丰满av| 国产精品亚洲av一区麻豆| 亚洲av成人不卡在线观看播放网| 无人区码免费观看不卡| 亚洲av成人不卡在线观看播放网| 最后的刺客免费高清国语| 看黄色毛片网站| 国产野战对白在线观看| 色综合站精品国产| 精品久久久久久久末码| 欧美大码av| 精品国产亚洲在线| 亚洲熟妇中文字幕五十中出| 久9热在线精品视频| 日韩欧美在线乱码| 国产精品久久久久久精品电影| 少妇的丰满在线观看| 国产av在哪里看| 少妇人妻精品综合一区二区 | 精品久久久久久,| 亚洲欧美激情综合另类| 日本成人三级电影网站| 亚洲av电影在线进入| 欧美成人免费av一区二区三区| 白带黄色成豆腐渣| 久久国产乱子伦精品免费另类| 69人妻影院| 又爽又黄无遮挡网站| 婷婷精品国产亚洲av| 国产黄a三级三级三级人| 久9热在线精品视频| 搞女人的毛片| 哪里可以看免费的av片| 99久久成人亚洲精品观看| 此物有八面人人有两片| a级毛片a级免费在线| av片东京热男人的天堂| 麻豆一二三区av精品| 一进一出抽搐gif免费好疼| 岛国视频午夜一区免费看| 美女cb高潮喷水在线观看| 欧美日韩国产亚洲二区| 亚洲国产精品999在线| 男人和女人高潮做爰伦理| 欧美日本视频| 国产精品免费一区二区三区在线| 成年免费大片在线观看| 婷婷精品国产亚洲av| 精品乱码久久久久久99久播| 嫩草影院入口| 久久国产精品人妻蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 丰满的人妻完整版| 成年人黄色毛片网站| 99精品在免费线老司机午夜| 亚洲美女黄片视频| 一本久久中文字幕| 国产精品1区2区在线观看.| 久久久精品欧美日韩精品| svipshipincom国产片| av福利片在线观看| 欧美午夜高清在线| 男插女下体视频免费在线播放| 中文字幕高清在线视频| 精品久久久久久,| 亚洲五月婷婷丁香| 亚洲一区二区三区不卡视频| 国产午夜精品久久久久久一区二区三区 | 午夜福利免费观看在线| 欧美乱码精品一区二区三区| 狠狠狠狠99中文字幕| 色老头精品视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 老熟妇乱子伦视频在线观看| www.www免费av| 欧美区成人在线视频| 久久6这里有精品| 免费一级毛片在线播放高清视频| 中文字幕人成人乱码亚洲影| 日韩高清综合在线| 天堂av国产一区二区熟女人妻| 午夜福利在线在线| 欧美成人性av电影在线观看| 熟女少妇亚洲综合色aaa.| 成人特级黄色片久久久久久久| 91在线精品国自产拍蜜月 | 身体一侧抽搐| 日日干狠狠操夜夜爽| 岛国在线观看网站| 久久精品91无色码中文字幕| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美日韩高清专用| 在线观看av片永久免费下载| 国产精品美女特级片免费视频播放器| 久久国产精品影院| 我要搜黄色片| 国产又黄又爽又无遮挡在线| 人人妻人人澡欧美一区二区| 亚洲成av人片免费观看| 国产精品国产高清国产av| 久久精品影院6| 亚洲人成伊人成综合网2020| 亚洲av五月六月丁香网| 国产主播在线观看一区二区| 久久午夜亚洲精品久久| 欧美绝顶高潮抽搐喷水| 久久久成人免费电影| 非洲黑人性xxxx精品又粗又长| 欧美三级亚洲精品| 91在线精品国自产拍蜜月 | 国产一区二区三区在线臀色熟女| 亚洲人成网站高清观看| 国产不卡一卡二| 日韩亚洲欧美综合| 又黄又粗又硬又大视频| 国产91精品成人一区二区三区| 1024手机看黄色片| 国产午夜精品论理片| 久久久久久久久中文| 亚洲18禁久久av| 99久久成人亚洲精品观看| 男人舔奶头视频| 欧美xxxx黑人xx丫x性爽| 午夜日韩欧美国产| 在线免费观看的www视频| 校园春色视频在线观看| 黄片小视频在线播放| 啦啦啦观看免费观看视频高清| 亚洲无线在线观看| 欧美激情久久久久久爽电影| 国产精华一区二区三区| av在线蜜桃| 给我免费播放毛片高清在线观看| 国产精品1区2区在线观看.| 偷拍熟女少妇极品色| 国产欧美日韩精品一区二区| 国产高清激情床上av| 三级毛片av免费| 国产亚洲精品久久久久久毛片| 黄色视频,在线免费观看| 国产又黄又爽又无遮挡在线| 亚洲乱码一区二区免费版| 一本一本综合久久| 黄色成人免费大全| 性色av乱码一区二区三区2| 一个人看视频在线观看www免费 | 午夜精品久久久久久毛片777| 两个人的视频大全免费| 美女高潮的动态| 男女做爰动态图高潮gif福利片| 尤物成人国产欧美一区二区三区| 一级a爱片免费观看的视频| 亚洲人成伊人成综合网2020| 两个人的视频大全免费| 日韩欧美精品免费久久 | 少妇的逼水好多| 啦啦啦免费观看视频1| 亚洲欧美日韩高清专用| 欧美日本视频| 久久这里只有精品中国| 首页视频小说图片口味搜索| 伊人久久大香线蕉亚洲五| 好男人在线观看高清免费视频| 1000部很黄的大片| 男人的好看免费观看在线视频| 国产乱人伦免费视频| 天美传媒精品一区二区| 久久久久久久久中文| 五月玫瑰六月丁香| 精品久久久久久成人av| 久久久久久九九精品二区国产| 最近最新中文字幕大全免费视频| 亚洲av成人精品一区久久| 久久精品人妻少妇| 亚洲av第一区精品v没综合| 日日夜夜操网爽| 国产私拍福利视频在线观看| 亚洲国产中文字幕在线视频| 国产午夜精品久久久久久一区二区三区 | 亚洲第一电影网av| 久久香蕉精品热| 高清日韩中文字幕在线| 午夜福利成人在线免费观看| 久久欧美精品欧美久久欧美| 一进一出抽搐动态| 99久久久亚洲精品蜜臀av| 丝袜美腿在线中文| 欧美绝顶高潮抽搐喷水| 一夜夜www| 成人高潮视频无遮挡免费网站| 国产一区二区在线av高清观看| 久久久久精品国产欧美久久久| 国产精品影院久久| 99riav亚洲国产免费| 国内精品一区二区在线观看| 久久久精品欧美日韩精品| 脱女人内裤的视频| 精品乱码久久久久久99久播| 国产v大片淫在线免费观看| 好男人电影高清在线观看| 免费观看的影片在线观看| 国产精品综合久久久久久久免费| 毛片女人毛片| 欧美日韩黄片免| 毛片女人毛片| 日本熟妇午夜| 成年女人永久免费观看视频| 精品久久久久久成人av| 国内精品美女久久久久久| 69av精品久久久久久| 国产精品99久久久久久久久| 美女高潮喷水抽搐中文字幕| 久久草成人影院| 国产野战对白在线观看| xxxwww97欧美| 国产男靠女视频免费网站| 国产亚洲精品久久久com| 淫妇啪啪啪对白视频| 人人妻人人澡欧美一区二区| 最后的刺客免费高清国语| 99久久99久久久精品蜜桃| 波多野结衣高清作品| 99热这里只有是精品50| 18禁国产床啪视频网站| 最新美女视频免费是黄的| 宅男免费午夜| 国产在线精品亚洲第一网站| 亚洲内射少妇av| 亚洲欧美日韩东京热| 国产淫片久久久久久久久 | 欧美一区二区精品小视频在线| 校园春色视频在线观看| 久久精品综合一区二区三区| av福利片在线观看| 午夜精品在线福利| 成人国产一区最新在线观看| 成人av在线播放网站| 在线天堂最新版资源| 51午夜福利影视在线观看| 免费人成在线观看视频色| 欧美黄色片欧美黄色片| 日韩欧美在线二视频| 国产精品久久久久久亚洲av鲁大| 99精品久久久久人妻精品| 久久久久久人人人人人| 午夜精品在线福利| 亚洲成人免费电影在线观看| 久久久国产成人免费| 成人av一区二区三区在线看| 国产单亲对白刺激| 变态另类丝袜制服| 欧美最新免费一区二区三区 | 国产亚洲av嫩草精品影院| ponron亚洲| 中文字幕精品亚洲无线码一区| 老熟妇仑乱视频hdxx| 女警被强在线播放| 日韩欧美国产一区二区入口| 午夜激情欧美在线| 国产91精品成人一区二区三区| 看免费av毛片| 午夜福利免费观看在线| 内射极品少妇av片p| 天天添夜夜摸| 一本综合久久免费| 叶爱在线成人免费视频播放| 好看av亚洲va欧美ⅴa在| 日韩欧美免费精品| 99国产精品一区二区三区| а√天堂www在线а√下载| 免费看光身美女| 免费人成视频x8x8入口观看| 成年版毛片免费区| 亚洲精品久久国产高清桃花| 成人av一区二区三区在线看| 久久久久性生活片| 一区福利在线观看| 成人精品一区二区免费| 欧美xxxx黑人xx丫x性爽| 国产三级在线视频| 精品国内亚洲2022精品成人| 色精品久久人妻99蜜桃| 色在线成人网| 国产精品电影一区二区三区| 久久精品国产综合久久久| 久久九九热精品免费| 色在线成人网| 级片在线观看| 亚洲精品国产精品久久久不卡| 特级一级黄色大片| 中国美女看黄片| 久久香蕉精品热| 久9热在线精品视频| 欧美日韩综合久久久久久 | 一a级毛片在线观看| 天堂影院成人在线观看| 国产av不卡久久| 久久久久久久久中文|