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

    基于輻射傳輸模型的高分二號影像大氣校正方法研究

    2020-07-01 01:22:16曹紅業(yè)張?zhí)炱?/span>
    紅外技術 2020年6期
    關鍵詞:氣溶膠反射率光學

    曹紅業(yè),張?zhí)炱?/p>

    基于輻射傳輸模型的高分二號影像大氣校正方法研究

    曹紅業(yè)1,張?zhí)炱?

    (1. 長安大學 地質工程與測繪學院,陜西 西安 710061;2. 西南科技大學 環(huán)境與資源學院,四川 綿陽 621010)

    高分二號衛(wèi)星的成功發(fā)射,標志著我國遙感衛(wèi)星進入了亞米級高空間分辨率時代,遙感影像在定量反演,地物識別和變化分析等領域將有重要作用。大氣校正的精度是影響其定量化應用的重要因素。由于高分二號遙感數據缺乏短波紅外波段,無法采用暗像元法進行大氣校正。提出一種基于輻射傳輸模型的高分二號影像大氣校正方法,利用6S(second simulation of the satellite signal in the solar spectrum)輻射傳輸模型建立大氣校正系數查找表,利用同步MODIS影像數據結合改進后的暗像元方法反演氣溶膠光學厚度,確定大氣校正系數,消除高分二號影像大氣分子和氣溶膠等的吸收和散射的影響,實現GF-2數據的大氣校正。選取地表平坦均一的敦煌輻射校正場作為實驗區(qū),通過同步的實測數據對校正結果進行精度評價,并且比較大氣校正前后歸一化植被指數NDVI。結果表明:最小相對誤差僅為0.9%,大氣校正之后影像數據更真實地反映了地物的反射特性;大氣校正后的NDVI大大增強了植被信息反差,突出了GF-2衛(wèi)星傳感器的植被信息區(qū)分能力。

    輻射傳輸模型;GF-2影像;大氣校正;查找表;氣溶膠光學厚度;6S

    0 引言

    2014年8月19日“高分二號”(GF-2)衛(wèi)星發(fā)射成功,標志著我國遙感衛(wèi)星進入亞米級高空間分辨時代[1]。該衛(wèi)星的成功發(fā)射,為我國土地利用動態(tài)監(jiān)測、礦產資源調查、城鄉(xiāng)規(guī)劃監(jiān)測評價、交通路網規(guī)劃、森林資源調查、荒漠化監(jiān)測等諸多方面應用提供了新的數據源。如表1所示,融合后的GF-2衛(wèi)星空間分辨率優(yōu)于1m,共有藍、綠、紅和近紅外4個多光譜波段,成像幅寬為45km,具有±35°的側擺角,可以實現69d內對全球的觀測覆蓋,以及5d內對地球表面任意區(qū)域的重復觀測。

    電磁波在傳播過程中受到大氣分子和氣溶膠等大氣成分的吸收和散射的影響,使傳感器獲得的電磁信號中帶有一定的非目標地物信息,傳感器接收到的地物信息不能真實地反應地表,嚴重影響其定量化應用,這就需要對該過程進行大氣校正,消除非目標地物信息,獲得真實的地表反射率信息[2-3]。因此大氣校正是GF-2影像數據地表參數定量反演的一個必備環(huán)節(jié)。

    目前,大氣校正的方法主要有基于圖像特征模型、地表實測線性回歸模型法和輻射傳輸模型法3種[4-6]。其中輻射傳輸法是利用復雜的輻射傳輸原理建立起來的一種方法,利用輻射傳輸模型反演地物反射率可以較合理地描述大氣散射、吸收等過程。自1972年Turner與Spencer通過模擬大氣-地表系統(tǒng)來評估大氣影響[7],成為最早的大氣輻射傳輸模型之一,在諸多大氣校正方法中具有高精度和高普適性。

    GF-2衛(wèi)星數據業(yè)務化應用要求基于遙感影像來反演大氣參數,進而完成整個大氣校正過程。一些傳感器專門設置了用于反演大氣參數的波段,如MODIS設置了15和16波段用于反演氣溶膠[8],17,18和19波段用于反演大氣水汽含量等[9]。但高分二號影像沒有設置用于反演大氣參數的波段,給高分二號影像的大氣校正帶來了很大的困難。

    本文提出一種基于同步MODIS數據輔助的高分二號多光譜數據的大氣校正算法。使用MODIS同步數據反演氣溶膠光學厚度,并建立高分二號影像大氣校正系數查找表,對高分二號影像進行大氣校正。在對一景高分二號影像進行大氣校正后,從典型地物光譜數據、歸一化植被指數(normalized difference vegetation index, NDVI)兩方面,探討了GF-2影像大氣校正的效果。

    1 大氣校正原理

    由于GF-2衛(wèi)星CCD相機屬于典型的4波段傳感器,相機所設置的4個波段均處于大氣窗口,水汽、臭氧等氣體的影響較小,可忽略不計。假設陸地表面是均勻朗伯表面,大氣垂直均勻變化的條件下,傳感器接收到的表觀輻亮度為如(1)式:

    表1 高分二號衛(wèi)星有效荷載技術指標

    大氣校正過程可以分為如下3步:

    ①利用暗像元法獲取氣溶膠光學厚度[11];

    1.1 數據預處理

    1.1.1 GF-2輻射定標

    根據傳感器各波段輻射定標系數將DN值轉換為表觀輻亮度,即傳感器入瞳處的輻射亮度值,輻射定標公式為:

    =Gain×DN+Offset (3)

    式中:為傳感器輻射亮度值;Gain和Offset分別為影像絕對定標系數增益和偏移,數據源自中國資源衛(wèi)星應用中心網站(www.cresda.com)。高分二號衛(wèi)星影像各波段2015年絕對輻射定標系數如表2所示(GF-2的Offset值全部為零)。

    表2 2015年GF-2 PMS輻射定標系數

    1.1.2 GF-2影像表觀反射率計算

    在上述表觀輻射亮度計算基礎上,根據如下計算公式得到大氣層頂表觀反射率:

    式中:為日地距離校正因子,如表3所示,計算時可對表中的數值進行插值獲得;為表觀輻亮度;s為太陽天頂角。

    ESUN是大氣上界太陽光輻照度,通過GF-2影像光譜響應函數(如圖1所示)和太陽光譜輻照度數據根據式(5)積分計算,GF-2影像光譜響應函數可在中國資源衛(wèi)星應用中心網站獲取,太陽光譜輻照度數據采用張璐等人[12]推薦的Wehrli數據,如圖2所示。ESUN計算公式如下:

    式中:()為波長處的大氣層外太陽光譜輻照度;()為波長處的光譜響應函數;1和2為積分波段范圍的上下限波長。

    1.2 計算氣溶膠光學厚度

    表3 日地距離

    Table 3 Solar-earth distance

    圖2 Wehrli太陽光譜輻照度

    利用2.1mm通道和1.2mm通道的表觀反射率可以準確地估算出紅藍通道的地表反射率,利用6S模型計算不同大氣氣溶膠模式和觀測條件下,氣溶膠光學厚度a(red)和a(blue)和大氣半球反射率,大氣程輻射a和大氣吸收構成的輻射透過率等大氣校正系數之間的關系,據此建立氣溶膠光學厚度查找表,其中與實際表觀反射率相等或最相近的模擬表觀反射率對應下的氣溶膠光學厚度即為最終結果,可以求得紅藍波段氣溶膠光學厚度a(red)和a(blue)。

    1.3 構建大氣校正系數查找表

    利用6S輻射傳輸模型,在輸入參數中設置多個不同的衛(wèi)星觀測幾何、太陽觀測幾何、衛(wèi)星天頂角、大氣模式、氣溶膠光學厚度、地面高度等變量。其中,衛(wèi)星觀測幾何、太陽觀測幾何、地面高度等參數從GF-2影像對應的元數據文件中獲??;大氣模式包括熱帶(Tropical, T)、中緯度夏季(Mid-latitude Summer, WLS)、中緯度冬季(Mid-latitude Winter, MLW)、近極地冬季(Subarctic Winter, SAW),近極地夏季(Subarctic Summer, SAS)等幾種大氣模式,根據GF-2影像的成像時間及緯度決定。由于6S模型本身并未附帶GF-2傳感器的光譜響應函數,因此,需要將其光譜響應函數重采樣為2.5nm分辨率,輸入到6S模型中。

    由于研究區(qū)域地物反射率均低于0.6,為了便于建立查找表,本文假設3個地表反射率為0,0.2,0.5,選定了大氣模式后,從GF-2影像元文件中讀取太陽天頂角、太陽方位角、衛(wèi)星天頂角、衛(wèi)星方位角4個角度數據,輸入1.2節(jié)求解的氣溶膠光學厚度,連續(xù)運行6S程序3次,得到其3個模擬表觀反射率分別為0.0t,0.2t,0.5t,將結果代入公式(6)中,得到一個關于大氣程輻射項a,大氣吸收構成的輻射透過率,大氣半球反射率的三元一次方程組,求解該方程組,得到如下用模擬表觀反射率為參數表示的解:

    通過設置上述參數,建立一個關于輻射傳輸方程中的大氣程輻射項a、大氣吸收構成的輻射透過率、大氣半球反射率、太陽天頂角s、衛(wèi)星天頂角v、太陽方位角和衛(wèi)星方位角的7維查找表,其他值在它們之間進行線性插值。

    1.4 基于查找表的GF-2大氣校正

    基于MODIS數據反演得到的氣溶膠光學厚度,結合大氣校正查找表逐像元插值,計算出對應的大氣校正系數,進而計算出地表反射率。大氣校正流程如圖3所示。

    圖3 高分二號影像大氣校正流程圖

    對MODIS進行去“蝴蝶”效應,幾何校正,輻射定標等預處理工作[16-18],使用IDL語言編程調用6S Fortran程序,建立氣溶膠光學厚度查找表,使用改進后的暗像元法逐像元反演氣溶膠光學厚度[19-24];利用6S模型,建立7維大氣校正查找表,確定大氣校正系數;利用該大氣校正系數,將GF-2表觀反射率轉換為地表反射率。

    2 結果與討論

    2.1 大氣校正對地物光譜對比分析

    本試驗區(qū)選擇敦煌輻射定標場,位于扇形沖積面上,地域開闊,地表平坦且均勻,植被覆蓋少,可認為是均勻朗伯表面,如圖4所示。選取2015年8月23日敦煌地區(qū)一景GF-2影像和同步過境的MODIS影像進行大氣校正,并將反演的表觀反射率、地面反射率與地面實測數據進行對比。

    衛(wèi)星觀測信號的響應依賴于衛(wèi)星波段的光譜響應函數,而地面觀測的高光譜反射率數據是連續(xù)的(如圖5(a))。為此,將光譜儀測量獲得的高光譜反射率數據進行積分,轉換為衛(wèi)星波段的等效反射率(如圖5(b))。轉換模型如下:

    式中:Req為衛(wèi)星波段等效測量反射率;R(l)為光譜儀測量獲得的高光譜反射率;f(l)為衛(wèi)星波段光譜響應函數。

    將4個實測點位置對于影像的表觀輻亮度值使用上述模型反演得到的表觀反射率、地表反射率和實測數據進行對比,對比結果如圖6所示。

    如圖所示,從對比結果可以發(fā)現,GF-2表觀反射率數據,經過大氣校正之后得到其地表反射率,藍波段降低,紅和紅外波段明顯提高,經過大氣校正后,影像上地物反射率與地面實測反射率的吻合程度較高,可以更真實地反映地物反射特征。

    表4列出了4個驗證點的實測地表反射率和大氣校正前后反射率的對比,并且對大氣校正前后反射率與實測地表反射率求取相對誤差。從表中不難發(fā)現,校正前最大相對誤差為S3藍波段,高達42.4%,經過大氣校正,相對誤差減小至13.2%,且最小相對誤差為S1的綠波段,僅為0.9%。與大氣校正前反射率相對誤差對比,大氣校正之后反射率與實測數據之間的差距明顯縮小,大氣校正效果明顯。

    2.2 大氣校正對NDVI的影響

    在遙感應用領域,植被指數已廣泛用來定性和定量評價植被覆蓋及其生長活力。植被指數有明顯的地域性和時效性,受植被本身、環(huán)境、大氣等條件的影響,NDVI可以用來檢驗大氣校正的效果。為了驗證GF-2數據大氣校正的效果,對校正前后的歸一化植被指數NDVI進行比較。NDVI計算公式如下:

    (a) 實測反射率(a) Measured reflectance (b) GF2傳感器波段等效后的反射率(b) Equivalent reflectance of GF2 sensor

    圖6 高分二號影像大氣校正前后反射率比較:圖(a)~(d)分別代表S1~S4四個點反射率

    式中:nir和red分別表示近紅外和紅波段的反射率。

    選取遙感圖像中的植被、水體和戈壁3種地物的反射率,根據上述公式,計算得到3種地物的NDVI,如圖7所示。3種地物的NDVI在大氣校正前后發(fā)生了變化,尤其是植被的NDVI增幅最大,可以看出大氣校正增大了NDVI提取植被信息的能力,經過大氣校正處理,可以有利于區(qū)分植被與其他地物之間的差異。

    3 結論與展望

    本文利用6S大氣輻射傳輸模型建立的7維查找表,結合MODIS表觀反射率數據,使用改進后的暗像元法反演氣溶膠光學厚度;依據輻射傳輸模型,建立大氣校正系數查找表,實現對GF-2影像數據快速有效地大氣校正。得到如下結論:

    對比大氣校正前后反射率變化,并使用地面實測數據對大氣校正結果進行精度評價,最大相對誤差12.566%,最小相對誤差僅為0.931%。與大氣校正前反射率相對誤差對比,大氣校正之后反射率與實測數據之間的差距明顯縮小,大氣校正之后的反射率數據更真實地反映了地物的反射特性。

    表4 高分二號影像大氣校正前后反射率比較

    aMR: Measured reflectance; BAC: Before atmospheric correction; AAC: After atmospheric correction; REBAC: Relative error before atmospheric correction; REAAC: Relative error after atmospheric correction.

    圖7 大氣校正前后典型地物NDVI

    GF-2衛(wèi)星經過大氣校正后,水體、戈壁和植被的NDVI均有變化;大氣校正后的NDVI大大增強了植被信息反差,突出了GF-2衛(wèi)星植被信息區(qū)分能力。

    本研究區(qū)域地域開闊,地表平坦且均勻,假設其為均勻朗伯表面,沒有考慮地物的二向反射特性。在今后研究中,將考慮二向反射模型(bi-direction reflectance distribution function,BRDF),臨近像元(周圍環(huán)境的反射)和交叉輻射項對結果的影響,提高結果的精度,使得該方法可以用于表面更加復雜的區(qū)域。

    [1] 王忠武, 劉順喜, 戴建旺, 等. “高分二號”衛(wèi)星多光譜與全色影像配準策略[J].航天返回與遙感, 2015, 36(4): 48-53.

    WANG Zhongwu, LIU Shunxi, DAI Jianwang, et al. Registration strategy for GF-2 satellite multispectral and panchromatic images[J]., 2015, 36(4): 48-53.

    [2] 李瑋, 康曉光, 陳雷. 一種MODIS遙感圖像大氣校正的快速算法[J]. 信號處理, 2007, 23(5): 751-754.

    LI Wei, KANG Xiaoguang, CHEN Lei. A fast atmospheric correction algorithm for MODIS remote sensing images[J]., 2007, 23(5): 751-754.

    [3] 阿布都瓦斯提, 秦其明, 朱黎江. 基于6S模型的可見光、近紅外遙感數據的大氣校正[J]. 北京大學學報, 2004, 40(4): 611-618.

    Ghulam A, QIN Qiming, ZHU Lijiang. 6S model based atmospheric correction of visible and near-infrared data and sensitivity analysis[J]., 2004, 40(4): 611-618.

    [4] Kaufman Y J, Remer L. Detection of forest using mid-IR reflectance: an application for aerosol studies[J]., 1994, 32(3): 672-684.

    [5] 亓雪勇, 田慶久. 光學遙感大氣校正研究進展[J]. 國土資源遙感, 2005, 25(4): 1-6.

    QI Xueyong, TIAN Qingjiu. The advances in the study of atmospheric correction for optical remote sensing[J]., 2005, 25(4): 1-6.

    [6] 劉兆亮, 麻金繼. 基于MODIS數據反演高反射率地區(qū)氣溶膠光學厚度的方法研究[J]. 大氣與環(huán)境光學學報, 2012, 7(5): 358-363.

    LIU Zhaoliang, MA Jinji. Retrieval of aerosol optical thickness over bright surface areas using MODIS data[J]., 2012, 7(5): 358-363.

    [7] Turner R E, Spencer M M. Atmospheric Model for Correction of Spacecraft Data[Z/OL]//2002:https://ui.adsabs. harvard.edu/abs/ 1972 rse..conf..895T/abstract.

    [8] Gordon H R. Atmospheric correction of ocean color imagery in the earth observing system era[J]., 1997, 102(D14): 17081-17106.

    [9] 趙秀娟, 陳長和, 張武, 等. 利用MODIS資料反演蘭州地區(qū)氣溶膠光學厚度[J]. 高原氣象, 2005, 24(1): 97-103.

    ZHAO Xiujuan, CHEN Changhe, ZHANG Wu, et al. 2005. Retrieval of aerosol optical depth over Lanzhou area using MODIS data[J]., 24(1): 97-103.

    [10] 李慶利, 薛永祺, 王建宇, 等. PHI高光譜圖像的大氣校正算法[J]. 紅外與毫米波學報, 2006, 25(4): 316-320.

    LI Qingli, XUE Yongqi, WANG Jianyu, et al. Atmospheric correction of PHI hyper spectral imagery[J]., 2006, 25(4): 316-320.

    [11] 彭光雄, 何宇華, 李京, 等. 中巴地球資源02星 CCD圖像交叉定標與大氣校正研究[J]. 紅外與毫米波學報, 2007, 26(1):22-26.

    PENG Guangxiong, HE Yuhua, LI Jing, et al. Study on CBERS-2’s CCD image cross calibration and atmospheric correction[J]., 2007, 26(1): 22-26.

    [12] 張璐, 施潤和, 徐永明, 等. 國產遙感傳感器大氣層外波段平均太陽光譜輻照度計算[J]. 地球信息科學, 2014, 16(4): 624-627.

    ZHANG Lu, SHI Runhe, XU Yongming, et al. Calculation of mean solar exoatmospheric irradiance of several sensors onboard of Chinese domestic remote sensing satellites[J]., 2014, 16(4): 624-627.

    [13] Kaufman Yoram J, Andrew E Wald, Lorraine A Remer, et al. The MODIS 2.1mm channel-correlation with visible reflectance for use in remote sensing of aerosol[J]., 1997, 35(5):1286-1298.

    [14] Robert H, Wayne L, Don L, et al. Cosmic ray and hot pixel removal from STIS CCD images[J]., 1997: 120-126.

    [15] Anders ?ngstr?m. The parameters of atmospheric turbidity[J].1964, 16(1): 64-76.

    [16] 郭廣猛. 關于MODIS衛(wèi)星數據的幾何校正方法[J]. 遙感信息, 2002(3): 26-28.

    GUO Guangmeng. Geometric calibration of MODIS[J]., 2002(3): 26-28.

    [17] 張怡, 何政偉, 薛東劍. MODIS數據預處理方法[J]. 地理空間信息, 2013, 11(3): 49-50.

    ZHANG Yi, HE Zhengwei, XUE Dongjian. The pretreatment method of MODIS data[J]., 2013, 11(3): 49-50.

    [18] 杜啟勝, 劉志平, 王新生, 等. 基于ENVI的MODIS數據預處理方法[J]. 地理空間信息, 2009, 7(4): 98-100.

    DU Qisheng, LIU Zhiping, WANG Xinsheng, et al. Method for MODIS data pre-processing based on ENVI[J]., 2009, 7(4): 98-100.

    [19] 范嬌, 郭寶峰, 何宏昌. 基于MODIS數據的杭州地區(qū)氣溶膠光學厚度反演[J]. 光學學報, 2015, 35(1): 1001-1009.

    FAN J, GUO Baofeng, HE Hongchang. Retrieval of aerosol optical thickness with MODIS data over Hangzhou[J]., 2015, 35(1): 1001-1009.

    [20] 王玲, 田慶久, 李珊珊. 利用MODIS資料反演杭州市500m分辨率氣溶膠光學厚度[J]. 遙感信息, 2010 (3): 50-54.

    WANG Ling, TIAN Qingjiu, LI Shanshan. Retrieval of 500m spatial resolution aerosol optical thickness with MODIS data over Hangzhou[J]., 2010(3): 50-54.

    [21] 任佳, 王振會, 孫林, 等. 基于MODIS數據反演江浙皖地區(qū)氣溶膠光學厚度[J]. 環(huán)境科學與技術, 2010, 33(8): 167-171.

    REN Jia, WANG Zhenhui, SUN Lin, et al. Retrieval of aerosol optical depth over Jiangsu-Zhejiang-Anhui areas using MODIS data[J]., 2010, 33(8): 167-171.

    [22] 趙志強, 李愛農, 邊金虎, 等. 基于改進暗目標法山區(qū)HJ CCD影像氣溶膠光學厚度反演[J]. 光譜學與光譜分析, 2015, 35(6): 1479-1487.

    ZHAO Zhiqiang, LI Ainong, BIAN Jinhu, et al. An improved DDV method to retrieve AOT for HJ CCD image in typical mountainous areas[J]., 2015, 35(6): 1479-1487.

    [23] 李微, 郭錫杰, 劉遠, 等. 基于改進黑暗像元法的環(huán)境一號CCD數據海岸帶區(qū)域大氣校正研究[J]. 海洋科學進展, 2015, 33(4): 484-491.

    LI Wei, GUO Xijie, LIU Yuan, et al. Improved dark-objects atmospheric correction algorithm for coastal regions based on HJ-1 CCD data[J]., 2015, 33(4): 484-491.

    [24] 鄭偉, 曾志遠. 遙感圖像大氣校正的黑暗像元法[J]. 國土資源遙感, 2005(1): 8-11.

    ZHENG Wei, ZENG Zhiuan. Dark-object methods for atmospheric correction of remote sensing image[J]., 2005(1): 8-11.

    Atmospheric Correction Algorithm for GF-2 Image Based On a Radiative Transfer Model

    CAO Hongye1,ZHANG Tianqi2

    (1.,,710061,; 2.,,621010,)

    The successful launch of the GF-2 satellite indicates that China's remote sensing satellites have entered the era of high spatial resolution of the sub-meter level. Remote sensing images will play an important role in quantitative inversion, ground object recognition, and change analysis. The accuracy of its atmospheric correction is an important factor that affects its quantitative application. Due to the lack of a short-wave infrared band in GF-2, it is impossible to use a dark pixel method for atmospheric correction. A method of atmospheric correction for the GF-2 image based on a radiation transfer model is proposed. The atmospheric correction coefficient lookup table is established by using 6S (Second Simulation of the Satellite Signal in the Solar Spectrum) radiation transfer model. The aerosol optical thickness (AOT) is retrieved by combining synchronous MODIS image data with an improved dark pixel method. The atmospheric correction parameters are determined to eliminate the influence of absorption and scattering of atmospheric molecules and aerosols in the GF-2 image and to achieve atmospheric correction of GF-2 data. Dunhuang radiation correction field with a flat and uniform surface is selected as the experimental area. The accuracy of the correction results is evaluated by synchronous measured data, and the normalized difference vegetation index (NDVI) before and after atmospheric correction is compared. The results show that the minimum relative error is only 0.9%. The image data after atmospheric correction can accurately reflect the reflection characteristics of ground objects. NDVI after atmospheric correction greatly enhances the contrast of vegetation information and highlights the ability of vegetation information discrimination of the GF-2 satellite sensor.

    radiation transfer model, GF-2 image, atmospheric correction, look-up table, aerosol optical depth(AOD), 6S

    TP751

    A

    1001-8891(2020)06-0534-08

    2019-07-17;

    2019-09-19.

    曹紅業(yè)(1989-),男,河南鞏義人,博士研究生,主要從事定量遙感方面的研究。E-mail:hong1233718@126.com。

    猜你喜歡
    氣溶膠反射率光學
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產品的融合方法研究
    海洋通報(2022年4期)2022-10-10 07:40:26
    滑輪組的裝配
    具有顏色恒常性的光譜反射率重建
    光學常見考題逐個擊破
    氣溶膠傳播之謎
    氣溶膠中210Po測定的不確定度評定
    化學腐蝕硅表面結構反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    四川盆地秋季氣溶膠與云的相關分析
    光學遙感壓縮成像技術
    欧美zozozo另类| 淫秽高清视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩中文字幕国产精品一区二区三区| 欧美精品亚洲一区二区| 成人三级做爰电影| 搡老妇女老女人老熟妇| 神马国产精品三级电影在线观看 | 妹子高潮喷水视频| 亚洲第一电影网av| 国产高清有码在线观看视频 | 日本黄色视频三级网站网址| 国产v大片淫在线免费观看| 俄罗斯特黄特色一大片| 天天躁狠狠躁夜夜躁狠狠躁| 搡老熟女国产l中国老女人| 99国产综合亚洲精品| 无限看片的www在线观看| x7x7x7水蜜桃| 久久久国产精品麻豆| 少妇粗大呻吟视频| 免费一级毛片在线播放高清视频| 国产日本99.免费观看| 亚洲中文av在线| 免费在线观看影片大全网站| 国产在线精品亚洲第一网站| 国产av又大| 岛国视频午夜一区免费看| 99re在线观看精品视频| 亚洲精品久久国产高清桃花| 久久久国产成人精品二区| 99国产综合亚洲精品| 国产精品国产高清国产av| 在线观看舔阴道视频| 亚洲欧美精品综合久久99| 亚洲激情在线av| 国产伦在线观看视频一区| 亚洲黑人精品在线| 亚洲熟妇熟女久久| 男女做爰动态图高潮gif福利片| 少妇粗大呻吟视频| 999久久久精品免费观看国产| 色老头精品视频在线观看| 黄色丝袜av网址大全| 三级毛片av免费| 男人的好看免费观看在线视频 | 美女午夜性视频免费| 国产高清视频在线播放一区| 欧美最黄视频在线播放免费| 欧美另类亚洲清纯唯美| 精品一区二区三区四区五区乱码| 一夜夜www| 亚洲精品美女久久久久99蜜臀| 岛国视频午夜一区免费看| 久久亚洲真实| 脱女人内裤的视频| 1024手机看黄色片| 日本a在线网址| 老司机在亚洲福利影院| 一级毛片女人18水好多| 99久久综合精品五月天人人| 久久婷婷成人综合色麻豆| 国产成人啪精品午夜网站| 俄罗斯特黄特色一大片| 亚洲激情在线av| 精品少妇一区二区三区视频日本电影| 亚洲天堂国产精品一区在线| 91国产中文字幕| 日韩三级视频一区二区三区| 无遮挡黄片免费观看| 夜夜看夜夜爽夜夜摸| 日韩大尺度精品在线看网址| 老司机福利观看| 真人一进一出gif抽搐免费| 两个人视频免费观看高清| 国产精品二区激情视频| 久久99热这里只有精品18| 欧美日韩亚洲综合一区二区三区_| 国产av一区在线观看免费| 1024手机看黄色片| 亚洲精品粉嫩美女一区| 久久久国产欧美日韩av| 777久久人妻少妇嫩草av网站| 看免费av毛片| 国产精品久久电影中文字幕| 久久久国产成人免费| 男女之事视频高清在线观看| 国产91精品成人一区二区三区| 自线自在国产av| 在线观看舔阴道视频| 久久这里只有精品19| 国产麻豆成人av免费视频| 亚洲成av人片免费观看| 中文字幕另类日韩欧美亚洲嫩草| 久久亚洲精品不卡| 国产精品亚洲一级av第二区| 久久精品国产综合久久久| 校园春色视频在线观看| 欧美成狂野欧美在线观看| 国产极品粉嫩免费观看在线| 日本 欧美在线| 久久国产精品男人的天堂亚洲| 成熟少妇高潮喷水视频| 亚洲五月婷婷丁香| 久久久久久大精品| 日韩欧美一区二区三区在线观看| 日韩欧美国产一区二区入口| 无限看片的www在线观看| 午夜免费鲁丝| 午夜久久久在线观看| 久久香蕉激情| 欧美黑人欧美精品刺激| 两个人免费观看高清视频| 精品福利观看| 非洲黑人性xxxx精品又粗又长| 欧美成人性av电影在线观看| 97人妻精品一区二区三区麻豆 | bbb黄色大片| 中国美女看黄片| 色综合站精品国产| 久久久久久久久久黄片| 免费高清视频大片| 99热只有精品国产| 久久精品影院6| 麻豆av在线久日| 亚洲av成人av| 午夜免费成人在线视频| 99国产精品一区二区三区| 国内精品久久久久精免费| 成熟少妇高潮喷水视频| 欧美日韩福利视频一区二区| 免费看十八禁软件| 免费观看人在逋| 久久香蕉国产精品| 午夜免费观看网址| 久久久久亚洲av毛片大全| 一级a爱视频在线免费观看| 一本综合久久免费| 欧美色欧美亚洲另类二区| 中亚洲国语对白在线视频| 视频区欧美日本亚洲| 精品日产1卡2卡| 中文字幕精品亚洲无线码一区 | 少妇熟女aⅴ在线视频| 成人三级做爰电影| 日本一本二区三区精品| 免费一级毛片在线播放高清视频| 国产精品乱码一区二三区的特点| 欧美不卡视频在线免费观看 | 一级毛片女人18水好多| 欧美久久黑人一区二区| 99国产精品99久久久久| 国产亚洲精品av在线| 又黄又爽又免费观看的视频| 少妇 在线观看| 性欧美人与动物交配| 欧美成狂野欧美在线观看| 99国产精品一区二区蜜桃av| 国产精品二区激情视频| av天堂在线播放| 一级a爱视频在线免费观看| 国产亚洲精品av在线| 人人澡人人妻人| 国产人伦9x9x在线观看| 99精品在免费线老司机午夜| 国产激情偷乱视频一区二区| 欧美又色又爽又黄视频| 757午夜福利合集在线观看| 看黄色毛片网站| 精品无人区乱码1区二区| 香蕉国产在线看| 国产伦人伦偷精品视频| 人人妻人人澡欧美一区二区| 久热爱精品视频在线9| 看黄色毛片网站| 亚洲 欧美一区二区三区| 成人18禁在线播放| 人妻久久中文字幕网| 黄色片一级片一级黄色片| 成年人黄色毛片网站| 日本 欧美在线| 91av网站免费观看| 国产成年人精品一区二区| 男人舔女人下体高潮全视频| 亚洲精品一区av在线观看| 窝窝影院91人妻| 国产成人av激情在线播放| 久久国产精品人妻蜜桃| 国产精品永久免费网站| 久久这里只有精品19| 桃红色精品国产亚洲av| 欧美色视频一区免费| 女警被强在线播放| 久久亚洲真实| 国产成+人综合+亚洲专区| 成年女人毛片免费观看观看9| 嫁个100分男人电影在线观看| 又大又爽又粗| 精品免费久久久久久久清纯| 国产av不卡久久| 变态另类丝袜制服| 久久精品国产清高在天天线| 男女下面进入的视频免费午夜 | 亚洲专区字幕在线| 国产激情欧美一区二区| 级片在线观看| 欧美久久黑人一区二区| 国产亚洲精品综合一区在线观看 | 日韩欧美一区二区三区在线观看| 精品国产美女av久久久久小说| 天堂动漫精品| 久久精品91无色码中文字幕| 别揉我奶头~嗯~啊~动态视频| 国产麻豆成人av免费视频| 国产精品1区2区在线观看.| 熟妇人妻久久中文字幕3abv| 午夜成年电影在线免费观看| 老司机福利观看| 一本精品99久久精品77| 日韩欧美国产在线观看| 精品久久久久久成人av| 国产熟女xx| 中亚洲国语对白在线视频| 成在线人永久免费视频| 午夜亚洲福利在线播放| 国产高清有码在线观看视频 | 后天国语完整版免费观看| 黄频高清免费视频| 欧美zozozo另类| 亚洲av熟女| 人人妻人人澡人人看| 欧美国产日韩亚洲一区| 国产精品影院久久| 国产精品av久久久久免费| 国产黄片美女视频| 精品国产国语对白av| 亚洲精品av麻豆狂野| 日日夜夜操网爽| 成人一区二区视频在线观看| 国产精品自产拍在线观看55亚洲| 91麻豆精品激情在线观看国产| 成人国产一区最新在线观看| 日韩中文字幕欧美一区二区| 亚洲精品一卡2卡三卡4卡5卡| 麻豆成人av在线观看| 人人妻,人人澡人人爽秒播| 免费看a级黄色片| 大型黄色视频在线免费观看| 亚洲国产高清在线一区二区三 | av免费在线观看网站| 一区二区三区高清视频在线| 国产真人三级小视频在线观看| 变态另类丝袜制服| 成人18禁高潮啪啪吃奶动态图| 国产又色又爽无遮挡免费看| 久久精品夜夜夜夜夜久久蜜豆 | 搡老岳熟女国产| 一边摸一边做爽爽视频免费| 国产黄片美女视频| 搞女人的毛片| 国产aⅴ精品一区二区三区波| 国产高清videossex| 午夜老司机福利片| 91av网站免费观看| 女性生殖器流出的白浆| 亚洲黑人精品在线| 99久久久亚洲精品蜜臀av| 国产精品98久久久久久宅男小说| 2021天堂中文幕一二区在线观 | 欧美一级毛片孕妇| 欧美日韩瑟瑟在线播放| 无遮挡黄片免费观看| 十八禁网站免费在线| 免费观看人在逋| 欧美日韩一级在线毛片| 免费观看精品视频网站| 999精品在线视频| 在线永久观看黄色视频| 亚洲七黄色美女视频| 99久久综合精品五月天人人| 国产精品亚洲美女久久久| 午夜免费成人在线视频| 欧美性猛交黑人性爽| 给我免费播放毛片高清在线观看| 亚洲精品一区av在线观看| 最近最新免费中文字幕在线| a级毛片在线看网站| 一本久久中文字幕| 国产欧美日韩一区二区精品| 亚洲欧洲精品一区二区精品久久久| 亚洲va日本ⅴa欧美va伊人久久| 欧美又色又爽又黄视频| 国产一区二区在线av高清观看| 成年人黄色毛片网站| 级片在线观看| 欧美国产精品va在线观看不卡| 成人国产一区最新在线观看| 国产av不卡久久| 国产精品美女特级片免费视频播放器 | 免费看美女性在线毛片视频| 亚洲成av人片免费观看| 最新在线观看一区二区三区| 看免费av毛片| 中文亚洲av片在线观看爽| av在线播放免费不卡| 欧美性猛交╳xxx乱大交人| 久久午夜亚洲精品久久| 日韩欧美国产在线观看| 一区二区日韩欧美中文字幕| 国产精品久久久人人做人人爽| 亚洲av熟女| 桃色一区二区三区在线观看| 丁香六月欧美| 久久久久久免费高清国产稀缺| 少妇被粗大的猛进出69影院| 日韩欧美免费精品| 99久久国产精品久久久| 成人一区二区视频在线观看| 韩国精品一区二区三区| 91麻豆精品激情在线观看国产| 欧美成人性av电影在线观看| 桃色一区二区三区在线观看| 中文字幕av电影在线播放| 999久久久精品免费观看国产| 亚洲片人在线观看| 欧美黑人巨大hd| 久久精品国产亚洲av香蕉五月| 国产精品二区激情视频| 久久天躁狠狠躁夜夜2o2o| 免费看日本二区| 9191精品国产免费久久| av中文乱码字幕在线| 国产精品自产拍在线观看55亚洲| 别揉我奶头~嗯~啊~动态视频| 亚洲国产中文字幕在线视频| 亚洲av成人一区二区三| 久久婷婷成人综合色麻豆| 长腿黑丝高跟| 狠狠狠狠99中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 日本三级黄在线观看| 岛国在线观看网站| av天堂在线播放| 欧美日韩中文字幕国产精品一区二区三区| 午夜精品久久久久久毛片777| 夜夜躁狠狠躁天天躁| 亚洲国产高清在线一区二区三 | 91av网站免费观看| 日韩一卡2卡3卡4卡2021年| 午夜日韩欧美国产| 999精品在线视频| 国产一区二区激情短视频| 国产欧美日韩一区二区精品| 午夜成年电影在线免费观看| 精品久久久久久久人妻蜜臀av| 999精品在线视频| www.999成人在线观看| 欧美成人一区二区免费高清观看 | 亚洲av五月六月丁香网| 9191精品国产免费久久| 人人妻,人人澡人人爽秒播| 国产欧美日韩一区二区三| 国产精品久久视频播放| 一本精品99久久精品77| 国产视频一区二区在线看| 精华霜和精华液先用哪个| 国产精品久久久久久人妻精品电影| e午夜精品久久久久久久| 亚洲国产日韩欧美精品在线观看 | 两性午夜刺激爽爽歪歪视频在线观看 | 精品欧美一区二区三区在线| 中文资源天堂在线| 麻豆久久精品国产亚洲av| 国产精品二区激情视频| 久9热在线精品视频| 一级毛片高清免费大全| 亚洲一区二区三区色噜噜| 亚洲精品粉嫩美女一区| 国产精品久久视频播放| 免费在线观看黄色视频的| 91大片在线观看| 亚洲午夜精品一区,二区,三区| 午夜福利在线在线| 欧美性猛交╳xxx乱大交人| 亚洲成a人片在线一区二区| 香蕉久久夜色| 国产野战对白在线观看| 成人午夜高清在线视频 | 大型黄色视频在线免费观看| 精品久久久久久久毛片微露脸| 免费在线观看黄色视频的| 欧美最黄视频在线播放免费| 午夜福利视频1000在线观看| 精品人妻1区二区| 在线免费观看的www视频| 成人手机av| 国产精品电影一区二区三区| www.自偷自拍.com| 午夜福利免费观看在线| 91麻豆av在线| 超碰成人久久| 黄色女人牲交| 在线观看免费日韩欧美大片| 国内精品久久久久精免费| 麻豆久久精品国产亚洲av| 婷婷六月久久综合丁香| 啪啪无遮挡十八禁网站| 九色国产91popny在线| 99久久久亚洲精品蜜臀av| 国产黄a三级三级三级人| 国产精品美女特级片免费视频播放器 | 日本黄色视频三级网站网址| avwww免费| 黄片大片在线免费观看| 亚洲欧洲精品一区二区精品久久久| 亚洲成人久久性| 美女高潮到喷水免费观看| 欧美日韩一级在线毛片| 两人在一起打扑克的视频| 搡老熟女国产l中国老女人| 欧美日本亚洲视频在线播放| 91av网站免费观看| 在线播放国产精品三级| 亚洲成人久久性| 成在线人永久免费视频| 欧美一级毛片孕妇| 变态另类成人亚洲欧美熟女| 亚洲男人天堂网一区| 两人在一起打扑克的视频| 一级作爱视频免费观看| 人人妻人人看人人澡| 午夜福利18| 欧美黑人巨大hd| 日韩大尺度精品在线看网址| 午夜福利在线观看吧| 久久伊人香网站| 国产99白浆流出| 久久香蕉国产精品| 久久精品国产99精品国产亚洲性色| 色老头精品视频在线观看| 日韩欧美三级三区| 99热6这里只有精品| 亚洲最大成人中文| 欧美中文日本在线观看视频| 中国美女看黄片| 国产亚洲精品av在线| 91麻豆精品激情在线观看国产| 麻豆成人午夜福利视频| 成人国产综合亚洲| 国产三级黄色录像| 中出人妻视频一区二区| 一边摸一边做爽爽视频免费| 国产精品爽爽va在线观看网站 | 国产99久久九九免费精品| 久久99热这里只有精品18| 在线观看日韩欧美| 一二三四在线观看免费中文在| 香蕉丝袜av| 免费一级毛片在线播放高清视频| 久久精品国产亚洲av香蕉五月| 国产日本99.免费观看| 身体一侧抽搐| 久久精品成人免费网站| 国产成+人综合+亚洲专区| 国产黄a三级三级三级人| 九色国产91popny在线| av欧美777| 久久精品国产99精品国产亚洲性色| 精品电影一区二区在线| 国产精品久久视频播放| 久久久久免费精品人妻一区二区 | 最近在线观看免费完整版| av免费在线观看网站| 黑丝袜美女国产一区| 两个人看的免费小视频| 午夜福利欧美成人| 久久精品国产99精品国产亚洲性色| 午夜福利高清视频| 精品久久久久久久末码| 欧美日本视频| 欧美日韩中文字幕国产精品一区二区三区| 午夜免费鲁丝| 91麻豆精品激情在线观看国产| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品在线观看二区| 18禁黄网站禁片午夜丰满| 欧美亚洲日本最大视频资源| 俺也久久电影网| 国内揄拍国产精品人妻在线 | 久久精品91无色码中文字幕| 国产精品亚洲美女久久久| 亚洲中文av在线| 日本一本二区三区精品| av欧美777| 可以免费在线观看a视频的电影网站| 悠悠久久av| 日韩高清综合在线| 777久久人妻少妇嫩草av网站| 精品国产乱码久久久久久男人| 99国产精品99久久久久| 亚洲精品色激情综合| 国产人伦9x9x在线观看| 日韩欧美一区视频在线观看| 校园春色视频在线观看| 日韩精品青青久久久久久| 中文字幕人妻熟女乱码| 女同久久另类99精品国产91| cao死你这个sao货| 欧美国产精品va在线观看不卡| 亚洲av电影在线进入| 亚洲精品av麻豆狂野| 99在线人妻在线中文字幕| 一本精品99久久精品77| av欧美777| 婷婷丁香在线五月| 午夜影院日韩av| 国产精品美女特级片免费视频播放器 | 美女高潮喷水抽搐中文字幕| 精品久久久久久久久久久久久 | 国产伦人伦偷精品视频| 亚洲午夜理论影院| 亚洲久久久国产精品| 国产人伦9x9x在线观看| 99re在线观看精品视频| 中出人妻视频一区二区| 亚洲一区二区三区色噜噜| 日本a在线网址| 日日爽夜夜爽网站| 99热只有精品国产| 国产人伦9x9x在线观看| 欧美在线黄色| 很黄的视频免费| 老司机靠b影院| 嫁个100分男人电影在线观看| 在线永久观看黄色视频| 少妇粗大呻吟视频| a在线观看视频网站| 欧美性长视频在线观看| 午夜福利免费观看在线| 精品久久久久久久久久免费视频| 少妇裸体淫交视频免费看高清 | 狠狠狠狠99中文字幕| 在线观看舔阴道视频| 天堂√8在线中文| 国内久久婷婷六月综合欲色啪| 亚洲国产精品999在线| 女警被强在线播放| 亚洲国产高清在线一区二区三 | 性欧美人与动物交配| 欧美黑人精品巨大| www.熟女人妻精品国产| 18禁国产床啪视频网站| 黑人欧美特级aaaaaa片| 男女床上黄色一级片免费看| 日日爽夜夜爽网站| 国产av不卡久久| 真人做人爱边吃奶动态| 日韩三级视频一区二区三区| 最新美女视频免费是黄的| 日韩欧美三级三区| 精品第一国产精品| 亚洲久久久国产精品| 国产激情欧美一区二区| 亚洲av成人av| 成人一区二区视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 久久国产亚洲av麻豆专区| 丝袜美腿诱惑在线| 亚洲自偷自拍图片 自拍| 亚洲av第一区精品v没综合| 777久久人妻少妇嫩草av网站| 成人手机av| 日韩免费av在线播放| 精品熟女少妇八av免费久了| 一进一出抽搐gif免费好疼| bbb黄色大片| www日本黄色视频网| 久久久久国产一级毛片高清牌| 亚洲激情在线av| 国产成人av激情在线播放| 亚洲中文av在线| 看片在线看免费视频| 黄片播放在线免费| 国产精品免费一区二区三区在线| 精品第一国产精品| 亚洲久久久国产精品| 国产99久久九九免费精品| 亚洲精品美女久久久久99蜜臀| 丁香六月欧美| 91av网站免费观看| 日韩国内少妇激情av| 日韩av在线大香蕉| 日本黄色视频三级网站网址| 国产欧美日韩一区二区精品| 国产精品影院久久| 国产极品粉嫩免费观看在线| av中文乱码字幕在线| 国产成人啪精品午夜网站| 亚洲成人国产一区在线观看| 99精品在免费线老司机午夜| 一夜夜www| 少妇 在线观看| 香蕉国产在线看| 日本免费一区二区三区高清不卡| 亚洲成人久久爱视频| 久久久久久人人人人人| 日本一区二区免费在线视频| 亚洲国产欧美日韩在线播放| 嫁个100分男人电影在线观看| 在线观看免费视频日本深夜| 一个人观看的视频www高清免费观看 | 久久性视频一级片| 久久久久久九九精品二区国产 | 久久久久国产一级毛片高清牌| 婷婷六月久久综合丁香|