• <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
    四川盆地秋季氣溶膠與云的相關分析
    光學遙感壓縮成像技術
    国产精品嫩草影院av在线观看| 久久韩国三级中文字幕| 亚洲av.av天堂| 国产高潮美女av| 午夜免费男女啪啪视频观看| 免费黄色在线免费观看| 91在线精品国自产拍蜜月| 日韩中字成人| 1000部很黄的大片| 免费一级毛片在线播放高清视频| 中文字幕免费在线视频6| 乱系列少妇在线播放| 欧美zozozo另类| 国产成人精品婷婷| 亚洲成人中文字幕在线播放| 免费观看在线日韩| 欧美人与善性xxx| 亚洲人成网站高清观看| 免费看av在线观看网站| 国产老妇伦熟女老妇高清| 日本免费a在线| 亚洲自偷自拍三级| 91aial.com中文字幕在线观看| 99视频精品全部免费 在线| 日日摸夜夜添夜夜添av毛片| 国产精品久久电影中文字幕| 精品久久久久久久末码| 国产国拍精品亚洲av在线观看| 亚洲精品,欧美精品| 天天一区二区日本电影三级| 99久久人妻综合| 日本欧美国产在线视频| 国产黄a三级三级三级人| 日本午夜av视频| 夜夜看夜夜爽夜夜摸| 久久久久久久久中文| 久久99热6这里只有精品| 亚洲婷婷狠狠爱综合网| 久久久精品欧美日韩精品| 久久99热这里只有精品18| 国产精品综合久久久久久久免费| 麻豆成人av视频| 最近最新中文字幕大全电影3| 国产精品嫩草影院av在线观看| 欧美性猛交╳xxx乱大交人| 看非洲黑人一级黄片| 国产老妇女一区| 精品久久久久久久人妻蜜臀av| 麻豆av噜噜一区二区三区| 精品欧美国产一区二区三| av在线播放精品| 精品午夜福利在线看| 久久久久久九九精品二区国产| 中文资源天堂在线| 精品人妻一区二区三区麻豆| 99热这里只有是精品在线观看| 2021少妇久久久久久久久久久| 两个人视频免费观看高清| 国产私拍福利视频在线观看| 91aial.com中文字幕在线观看| 国内少妇人妻偷人精品xxx网站| 热99在线观看视频| 中文在线观看免费www的网站| 亚洲欧美日韩高清专用| 精品免费久久久久久久清纯| 亚洲自偷自拍三级| 男人舔女人下体高潮全视频| 日韩av不卡免费在线播放| 成年版毛片免费区| 精品久久久久久久人妻蜜臀av| 亚洲国产精品国产精品| 久久精品国产亚洲网站| 一区二区三区乱码不卡18| 日韩欧美精品免费久久| 最近中文字幕高清免费大全6| 亚洲自偷自拍三级| 亚洲怡红院男人天堂| 久久欧美精品欧美久久欧美| 国产精品久久久久久久电影| 免费看光身美女| 国产成人a区在线观看| 久久精品国产亚洲av天美| 久久午夜福利片| 久久人人爽人人爽人人片va| 国产极品天堂在线| 国产伦一二天堂av在线观看| 干丝袜人妻中文字幕| 日韩av在线大香蕉| 国产成人91sexporn| 美女脱内裤让男人舔精品视频| av在线老鸭窝| 天天躁日日操中文字幕| 国产亚洲91精品色在线| 久久精品91蜜桃| 中文资源天堂在线| 蜜桃亚洲精品一区二区三区| 18禁在线无遮挡免费观看视频| 舔av片在线| 男女下面进入的视频免费午夜| 大话2 男鬼变身卡| 午夜免费男女啪啪视频观看| 久久精品国产亚洲网站| av女优亚洲男人天堂| 毛片女人毛片| 搡老妇女老女人老熟妇| 国产一区二区在线av高清观看| 亚洲美女视频黄频| 国产一区二区在线观看日韩| 国产三级在线视频| 精品人妻偷拍中文字幕| 日本黄色视频三级网站网址| 精品国产一区二区三区久久久樱花 | 亚洲中文字幕日韩| 久久精品国产亚洲av涩爱| 一级毛片aaaaaa免费看小| 国产色爽女视频免费观看| 国产av一区在线观看免费| 欧美成人免费av一区二区三区| 能在线免费观看的黄片| 能在线免费观看的黄片| 男人和女人高潮做爰伦理| 午夜福利在线观看吧| 女人十人毛片免费观看3o分钟| 国产精品女同一区二区软件| 色综合亚洲欧美另类图片| 婷婷六月久久综合丁香| 亚洲欧美成人综合另类久久久 | 99热全是精品| 亚洲国产色片| 亚洲国产精品成人久久小说| 在现免费观看毛片| 久久韩国三级中文字幕| 亚洲欧美日韩卡通动漫| 久久99热这里只频精品6学生 | 在现免费观看毛片| 男插女下体视频免费在线播放| 69av精品久久久久久| 久久久精品欧美日韩精品| 国产精品爽爽va在线观看网站| 精品久久久久久久久av| 亚洲自偷自拍三级| 亚洲国产欧美在线一区| 午夜a级毛片| 国产 一区精品| 亚洲av中文av极速乱| 亚洲av中文av极速乱| 亚洲精品乱码久久久久久按摩| 麻豆久久精品国产亚洲av| 国产av码专区亚洲av| 精品国产露脸久久av麻豆 | 色综合色国产| 色5月婷婷丁香| 中文字幕av在线有码专区| 精品少妇黑人巨大在线播放 | 男人舔女人下体高潮全视频| 久久精品人妻少妇| 亚洲电影在线观看av| 国产精品一区二区在线观看99 | 少妇猛男粗大的猛烈进出视频 | 超碰97精品在线观看| 亚洲欧美一区二区三区国产| 日本欧美国产在线视频| 男插女下体视频免费在线播放| 青青草视频在线视频观看| 国产男人的电影天堂91| 亚洲欧美成人综合另类久久久 | 大香蕉97超碰在线| 国产精品久久久久久av不卡| 99九九线精品视频在线观看视频| 长腿黑丝高跟| 免费电影在线观看免费观看| 精华霜和精华液先用哪个| 国产成人freesex在线| 国产免费又黄又爽又色| 国产69精品久久久久777片| 久热久热在线精品观看| 欧美成人免费av一区二区三区| 久久久久久久国产电影| 国产精品久久久久久久久免| 免费看日本二区| 国产精品蜜桃在线观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品久久男人天堂| 又爽又黄a免费视频| 老女人水多毛片| 亚洲在线观看片| 日韩欧美在线乱码| 高清视频免费观看一区二区 | 国产精品麻豆人妻色哟哟久久 | 国产单亲对白刺激| 欧美成人精品欧美一级黄| 亚洲av电影不卡..在线观看| 亚洲美女搞黄在线观看| 别揉我奶头 嗯啊视频| 国产大屁股一区二区在线视频| 国产精品麻豆人妻色哟哟久久 | 成年av动漫网址| 午夜福利成人在线免费观看| 国产在视频线在精品| 又爽又黄无遮挡网站| 青春草亚洲视频在线观看| 久久久久久久久久成人| 日本免费a在线| 亚洲人成网站在线观看播放| 成人鲁丝片一二三区免费| 性色avwww在线观看| 成人性生交大片免费视频hd| 欧美3d第一页| 老师上课跳d突然被开到最大视频| 欧美潮喷喷水| 国产探花极品一区二区| 亚洲欧美精品自产自拍| 亚洲欧美精品专区久久| 久久99热这里只有精品18| 亚洲自偷自拍三级| 中文字幕熟女人妻在线| 性插视频无遮挡在线免费观看| 午夜精品在线福利| 韩国高清视频一区二区三区| 免费黄网站久久成人精品| 啦啦啦韩国在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 成人欧美大片| 国产精品电影一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 丝袜喷水一区| 一个人看视频在线观看www免费| 寂寞人妻少妇视频99o| 亚洲精品日韩av片在线观看| 99久久中文字幕三级久久日本| 九色成人免费人妻av| 岛国在线免费视频观看| 亚洲国产精品成人久久小说| 18禁动态无遮挡网站| 免费人成在线观看视频色| 伦理电影大哥的女人| 成年女人看的毛片在线观看| 欧美日本亚洲视频在线播放| 热99在线观看视频| 99在线人妻在线中文字幕| 久久精品久久精品一区二区三区| 亚洲va在线va天堂va国产| 亚洲中文字幕日韩| 日韩av在线免费看完整版不卡| 国产一区二区三区av在线| 色噜噜av男人的天堂激情| 色吧在线观看| 嫩草影院精品99| av卡一久久| 亚洲国产色片| 久久99热6这里只有精品| av线在线观看网站| 波多野结衣巨乳人妻| 网址你懂的国产日韩在线| 亚洲第一区二区三区不卡| av在线观看视频网站免费| 日韩强制内射视频| 国产亚洲av嫩草精品影院| 国产中年淑女户外野战色| 亚洲欧美精品综合久久99| 大又大粗又爽又黄少妇毛片口| 高清视频免费观看一区二区 | 人人妻人人看人人澡| 成年版毛片免费区| 精品国产露脸久久av麻豆 | 欧美成人午夜免费资源| 国产伦精品一区二区三区四那| 小说图片视频综合网站| 日产精品乱码卡一卡2卡三| 国产欧美另类精品又又久久亚洲欧美| 一级二级三级毛片免费看| 国产淫片久久久久久久久| 中文字幕免费在线视频6| 少妇的逼水好多| 美女脱内裤让男人舔精品视频| av播播在线观看一区| 国产黄a三级三级三级人| 午夜老司机福利剧场| 欧美成人免费av一区二区三区| 国产精品野战在线观看| 亚洲av二区三区四区| 2021天堂中文幕一二区在线观| 久久久久久久久久黄片| 国产黄片视频在线免费观看| 热99re8久久精品国产| 久久久精品大字幕| av天堂中文字幕网| 97人妻精品一区二区三区麻豆| 日本猛色少妇xxxxx猛交久久| 欧美成人a在线观看| 免费看日本二区| 日韩一区二区视频免费看| 免费av不卡在线播放| 美女大奶头视频| 国产视频首页在线观看| 中文乱码字字幕精品一区二区三区 | 久久精品影院6| 一边摸一边抽搐一进一小说| 精品久久久久久成人av| 国产伦一二天堂av在线观看| 人妻夜夜爽99麻豆av| 欧美日韩精品成人综合77777| 久久久久久国产a免费观看| 亚洲精品国产av成人精品| 亚洲不卡免费看| 91精品一卡2卡3卡4卡| 日本黄色视频三级网站网址| 日韩制服骚丝袜av| 亚洲欧美日韩卡通动漫| 超碰97精品在线观看| 亚洲人成网站在线观看播放| 两个人的视频大全免费| 美女黄网站色视频| 嫩草影院入口| 亚洲国产成人一精品久久久| 中国美白少妇内射xxxbb| av在线蜜桃| 免费不卡的大黄色大毛片视频在线观看 | 一个人看的www免费观看视频| 亚洲乱码一区二区免费版| 男女边吃奶边做爰视频| 亚洲四区av| 久久久久久久久久久丰满| 天天一区二区日本电影三级| 一区二区三区高清视频在线| 99热网站在线观看| 日日摸夜夜添夜夜爱| 麻豆av噜噜一区二区三区| 免费av观看视频| 你懂的网址亚洲精品在线观看 | 欧美高清性xxxxhd video| 伦精品一区二区三区| av在线观看视频网站免费| 男女边吃奶边做爰视频| 国产片特级美女逼逼视频| 天堂网av新在线| 国产极品精品免费视频能看的| 免费电影在线观看免费观看| 99国产精品一区二区蜜桃av| 久久草成人影院| 久久久久免费精品人妻一区二区| 国产三级中文精品| 亚洲国产欧美人成| 能在线免费看毛片的网站| 尤物成人国产欧美一区二区三区| 国产精品一区二区三区四区免费观看| 精品人妻视频免费看| 淫秽高清视频在线观看| 能在线免费观看的黄片| 99热精品在线国产| 天堂中文最新版在线下载 | 美女xxoo啪啪120秒动态图| 国产麻豆成人av免费视频| 啦啦啦韩国在线观看视频| 亚州av有码| 狠狠狠狠99中文字幕| 中文字幕制服av| 伦理电影大哥的女人| 亚洲av成人精品一二三区| 亚洲五月天丁香| 亚洲三级黄色毛片| 22中文网久久字幕| 一级二级三级毛片免费看| 久久精品久久久久久久性| 国产精品乱码一区二三区的特点| www日本黄色视频网| 性插视频无遮挡在线免费观看| 国产av一区在线观看免费| 久久精品91蜜桃| 九色成人免费人妻av| 秋霞在线观看毛片| 国产伦在线观看视频一区| 亚洲综合色惰| 国产中年淑女户外野战色| 免费av不卡在线播放| av免费在线看不卡| 久久久久网色| 亚洲精华国产精华液的使用体验| 亚洲精品成人久久久久久| 91精品伊人久久大香线蕉| 成年av动漫网址| 全区人妻精品视频| 精品久久久久久久久av| 亚洲一区高清亚洲精品| 国产白丝娇喘喷水9色精品| 高清日韩中文字幕在线| 麻豆成人午夜福利视频| 中文资源天堂在线| 欧美bdsm另类| 丰满少妇做爰视频| 久久精品熟女亚洲av麻豆精品 | 亚洲人成网站在线观看播放| 欧美日韩在线观看h| 毛片一级片免费看久久久久| 亚洲在线自拍视频| 波多野结衣高清无吗| 久久精品久久精品一区二区三区| 亚洲精品乱码久久久v下载方式| 高清午夜精品一区二区三区| 精品久久久久久久久av| 美女高潮的动态| 网址你懂的国产日韩在线| 精品午夜福利在线看| 桃色一区二区三区在线观看| 黄色日韩在线| 久久精品影院6| 国产单亲对白刺激| 夫妻性生交免费视频一级片| 亚洲精品国产成人久久av| 久久久国产成人精品二区| 亚洲欧美日韩东京热| 久久婷婷人人爽人人干人人爱| 春色校园在线视频观看| 欧美一级a爱片免费观看看| 日韩人妻高清精品专区| 深爱激情五月婷婷| 欧美另类亚洲清纯唯美| 人妻夜夜爽99麻豆av| 免费看美女性在线毛片视频| 日本三级黄在线观看| 国产成人精品久久久久久| 亚洲av成人av| 国产精品国产三级国产专区5o | 国产精品一二三区在线看| 中文字幕亚洲精品专区| 99视频精品全部免费 在线| 精品久久国产蜜桃| 欧美成人免费av一区二区三区| 你懂的网址亚洲精品在线观看 | 国产色爽女视频免费观看| 国内少妇人妻偷人精品xxx网站| 我要搜黄色片| 又爽又黄a免费视频| av在线蜜桃| 久久综合国产亚洲精品| 内地一区二区视频在线| 永久网站在线| 视频中文字幕在线观看| 久久久久久久久久久丰满| 欧美日韩国产亚洲二区| 久久久亚洲精品成人影院| 国产精品精品国产色婷婷| 午夜福利视频1000在线观看| 国产av不卡久久| 一本一本综合久久| 在线播放无遮挡| 热99在线观看视频| 久久久久网色| 欧美激情国产日韩精品一区| 亚洲av成人精品一二三区| 国产精品熟女久久久久浪| 国产午夜精品论理片| 日韩制服骚丝袜av| 久久精品夜色国产| 欧美最新免费一区二区三区| 日本猛色少妇xxxxx猛交久久| 婷婷色av中文字幕| 久久精品91蜜桃| 深爱激情五月婷婷| 级片在线观看| 日韩一本色道免费dvd| 搡老妇女老女人老熟妇| 又黄又爽又刺激的免费视频.| 神马国产精品三级电影在线观看| 91精品伊人久久大香线蕉| 国产淫语在线视频| 美女大奶头视频| 一本一本综合久久| 国产成人一区二区在线| 美女国产视频在线观看| 一级黄片播放器| 国产高清不卡午夜福利| 一级毛片电影观看 | 亚洲人成网站在线观看播放| www日本黄色视频网| 国产亚洲精品久久久com| 小蜜桃在线观看免费完整版高清| 亚洲精品,欧美精品| 国产黄片视频在线免费观看| 久久久久性生活片| 久久人妻av系列| 国产中年淑女户外野战色| 久久久久久久久久久丰满| 国产精品.久久久| 99久久精品热视频| 韩国av在线不卡| 欧美区成人在线视频| 欧美一区二区亚洲| 亚洲三级黄色毛片| 亚洲美女视频黄频| 成年免费大片在线观看| 在线天堂最新版资源| 嫩草影院新地址| 麻豆一二三区av精品| 99久久中文字幕三级久久日本| 国产视频内射| 国产女主播在线喷水免费视频网站 | 中国美白少妇内射xxxbb| 亚洲18禁久久av| 国产又黄又爽又无遮挡在线| 久久鲁丝午夜福利片| 天堂√8在线中文| 国产精品永久免费网站| 少妇人妻一区二区三区视频| 久99久视频精品免费| 色网站视频免费| 免费大片18禁| 中文字幕久久专区| 九九久久精品国产亚洲av麻豆| 秋霞在线观看毛片| 国产成人精品久久久久久| 日日摸夜夜添夜夜爱| 中文资源天堂在线| 变态另类丝袜制服| 精品久久久久久久久久久久久| 国产av在哪里看| 亚洲欧美日韩高清专用| 一个人看的www免费观看视频| 一级毛片久久久久久久久女| 国产日韩欧美在线精品| 久久精品久久久久久噜噜老黄 | 人人妻人人澡人人爽人人夜夜 | 亚洲av日韩在线播放| av卡一久久| 日韩,欧美,国产一区二区三区 | 一级爰片在线观看| 国内揄拍国产精品人妻在线| 麻豆av噜噜一区二区三区| 人妻夜夜爽99麻豆av| 最后的刺客免费高清国语| 久久鲁丝午夜福利片| 一级二级三级毛片免费看| 免费人成在线观看视频色| 日韩大片免费观看网站 | 亚洲性久久影院| 久久久精品94久久精品| 一级爰片在线观看| 男人舔女人下体高潮全视频| 成人二区视频| 久久99蜜桃精品久久| 国产人妻一区二区三区在| 永久网站在线| 欧美xxxx黑人xx丫x性爽| 精品酒店卫生间| 日本爱情动作片www.在线观看| 免费人成在线观看视频色| 日韩大片免费观看网站 | 高清视频免费观看一区二区 | 久久人人爽人人爽人人片va| 久久久a久久爽久久v久久| 久久婷婷人人爽人人干人人爱| 婷婷色av中文字幕| 只有这里有精品99| 最近的中文字幕免费完整| 成人毛片60女人毛片免费| 美女cb高潮喷水在线观看| 精品国产一区二区三区久久久樱花 | 麻豆一二三区av精品| 国产成人精品婷婷| 亚洲在线自拍视频| 日日撸夜夜添| videos熟女内射| 国产精品久久久久久精品电影| 亚洲久久久久久中文字幕| 村上凉子中文字幕在线| 99热这里只有是精品在线观看| 波多野结衣高清无吗| 免费黄色在线免费观看| 日韩精品青青久久久久久| 精品无人区乱码1区二区| 神马国产精品三级电影在线观看| 热99re8久久精品国产| 男人舔奶头视频| 国产中年淑女户外野战色| 日韩av在线免费看完整版不卡| 51国产日韩欧美| 久久久亚洲精品成人影院| 欧美xxxx性猛交bbbb| 2022亚洲国产成人精品| 三级国产精品片| 村上凉子中文字幕在线| 国产69精品久久久久777片| 亚洲欧美日韩无卡精品| 国产精品久久久久久久久免| 婷婷色综合大香蕉| 国产精品久久电影中文字幕| 99久久精品热视频| 中文字幕免费在线视频6| 国产精品av视频在线免费观看| 中文字幕制服av| 国产精品av视频在线免费观看| 亚洲国产精品国产精品| 日韩大片免费观看网站 | 夫妻性生交免费视频一级片| 午夜老司机福利剧场| 国产一级毛片七仙女欲春2| 好男人在线观看高清免费视频| 成人毛片a级毛片在线播放| 在线免费观看不下载黄p国产| 日本五十路高清| 99久国产av精品| 久久精品夜夜夜夜夜久久蜜豆| av黄色大香蕉| 国产一区亚洲一区在线观看| 免费人成在线观看视频色| 亚洲欧洲日产国产| a级一级毛片免费在线观看| 激情 狠狠 欧美| 久久精品国产亚洲av涩爱| 日韩av在线大香蕉| 国产老妇女一区| 日本一二三区视频观看| 两个人的视频大全免费| 久久精品夜色国产| 亚洲自拍偷在线|