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

    星載大氣監(jiān)測(cè)光譜儀高精度星上定標(biāo)技術(shù)

    2018-08-10 05:44:42李碧岑李明徐彭梅
    航天返回與遙感 2018年3期
    關(guān)鍵詞:譜段光譜儀定標(biāo)

    李碧岑 李明 徐彭梅

    ?

    星載大氣監(jiān)測(cè)光譜儀高精度星上定標(biāo)技術(shù)

    李碧岑 李明 徐彭梅

    (北京空間機(jī)電研究所,北京 100094)

    大氣成分的高精度反演及應(yīng)用對(duì)星載超光譜載荷的輻射精度和光譜精度均提出了更高要求,且隨著遙感器運(yùn)行壽命的不斷延長(zhǎng),需建立高精度、高穩(wěn)定的星上定標(biāo)系統(tǒng)。文章介紹了一種滿足大氣監(jiān)測(cè)光譜儀高精度定標(biāo)要求的星上定標(biāo)技術(shù),結(jié)合時(shí)間調(diào)制型傅里葉變換光譜儀的技術(shù)特點(diǎn),制定了星上絕對(duì)輻射定標(biāo)、儀器線形函數(shù)ILS測(cè)量和光譜定標(biāo)的方案。采用太陽(yáng)漫反射板定標(biāo)法實(shí)現(xiàn)全口徑、全視場(chǎng)、全光路絕對(duì)輻射定標(biāo),定標(biāo)漫反射板在光譜儀光路的最前端將太陽(yáng)光譜引入,通過(guò)已知的大氣層外太陽(yáng)光譜照度和地面標(biāo)定的漫反射板雙向反射分布函數(shù)BRDF確立星上絕對(duì)輻射定標(biāo)基準(zhǔn)。設(shè)置參考漫反射板進(jìn)行在軌性能衰減的監(jiān)測(cè)和校正,采用的QVD漫反射板具有高穩(wěn)定性,可保證全壽命周期內(nèi)星上絕對(duì)輻射定標(biāo)精度優(yōu)于5%。星上設(shè)置單色穩(wěn)頻激光器對(duì)光譜儀的儀器線形函數(shù)進(jìn)行定期監(jiān)測(cè),以評(píng)估光譜儀的光譜分辨率等光學(xué)性能的在軌狀態(tài)。利用太陽(yáng)光譜和大氣光譜中的特征譜線進(jìn)行在軌波長(zhǎng)校正。

    漫反射板 儀器線形函數(shù) 波長(zhǎng)校正 絕對(duì)輻射定標(biāo) 大氣監(jiān)測(cè)光譜儀 航天遙感 “高分五號(hào)”衛(wèi)星

    0 引言

    對(duì)于大氣超光譜遙感載荷而言,光譜數(shù)據(jù)的輻射精度和光譜精度是影響大氣成分反演的重要因素[1]。其中,輻射精度,即測(cè)量光譜輻亮度的精度,是大氣吸收強(qiáng)度的直接表征信息,與載荷定標(biāo)精度、儀器的信噪比以及反演過(guò)程引入的誤差等密切相關(guān);光譜精度,主要是譜線的位置精度,是影響大氣吸收譜線精確定位的主要因素。上述兩個(gè)因素都需要通過(guò)發(fā)射前和在軌的精確定標(biāo)來(lái)保證。

    國(guó)外的空間大氣光譜探測(cè)技術(shù)研究起步較早,在高精度星上定標(biāo)系統(tǒng)的設(shè)計(jì)和應(yīng)用方面經(jīng)驗(yàn)豐富。GOME(Global Ozone Monitoring Experiment)和SCIAMACHY(Scanning Imaging Absorption Spectro Meter for Atmospheric Cartograp HY)采用鋁漫反射板反射太陽(yáng)光進(jìn)行輻射定標(biāo),但鋁漫反射板會(huì)對(duì)其反射太陽(yáng)光譜引入一定的光譜結(jié)構(gòu),影響反演結(jié)果的精度[2]。TANSO-FTS(Thermal and Near Infrared Sensor for Carbon Observation Fourier-Transform Spectrometer)采用Spectralon漫反射板進(jìn)行星上輻射定標(biāo),其入軌后一年的監(jiān)測(cè)數(shù)據(jù)顯示,每軌定標(biāo)使用的漫反射板相比每月一次校正用漫反射板的反射率在0.76μm譜段衰減了約7%[3]。OMI(Ozone Monitoring Instrument)[4]和GOME-2[5]均采用了一種自身光譜結(jié)構(gòu)很小的石英漫射板(Quasi Volume Diffuser, QVD),其在軌性能穩(wěn)定,OMI的QVD漫射板在500nm處的在軌10年衰減小于3%。OMI和GOME-2在發(fā)射前除了進(jìn)行亮度和照度的絕對(duì)定標(biāo)外,對(duì)載荷的另一重要輻射定標(biāo)參數(shù)系統(tǒng)級(jí)雙向反射分布函數(shù)(Bidirectional Reflectance Distribution Function,BRDF)進(jìn)行了詳細(xì)測(cè)試,即亮度和照度定標(biāo)結(jié)果的比值。系統(tǒng)級(jí)BRDF的標(biāo)定消除了載荷對(duì)地光路和星上定標(biāo)光路的差異對(duì)定標(biāo)精度的影響,且所采用光源的輻射輸出無(wú)需精確標(biāo)定,可有效提升載荷的星上定標(biāo)精度[2]。

    近年來(lái),隨著國(guó)內(nèi)遙感載荷定量化應(yīng)用需求的不斷提高和衛(wèi)星設(shè)計(jì)壽命的延長(zhǎng),均開(kāi)始設(shè)計(jì)并配備高精度的星上定標(biāo)裝置[6],星上定標(biāo)精度擬提升到5%[7],以滿足對(duì)地遙感的大部分定量化反演要求?!案叻治逄?hào)”作為我國(guó)首顆高光譜觀測(cè)衛(wèi)星,具備對(duì)大氣污染氣體、溫室氣體、氣溶膠等要素的定量監(jiān)測(cè)能力,其搭載的多個(gè)用于大氣監(jiān)測(cè)的高光譜載荷采用了高精度星上光譜和輻射定標(biāo)技術(shù),以滿足高精度數(shù)據(jù)應(yīng)用需求[8-10]。本文以某一用于大氣監(jiān)測(cè)的傅里葉變換光譜儀為例,介紹了其高精度在軌光譜和輻射定標(biāo)方案和星上定標(biāo)系統(tǒng)的設(shè)計(jì)。

    1 大氣監(jiān)測(cè)光譜儀簡(jiǎn)介

    星載大氣監(jiān)測(cè)光譜儀用于大氣中氣體濃度變化的在軌定量監(jiān)測(cè),相比成像型的遙感相機(jī),大氣成分的高精度反演對(duì)光譜儀的性能和定標(biāo)精度均提出了更高要求。除了必須在發(fā)射前進(jìn)行嚴(yán)格的實(shí)驗(yàn)室定標(biāo),還必須建立高精度、高穩(wěn)定的星上定標(biāo)系統(tǒng),以保證載荷全壽命周期內(nèi)數(shù)據(jù)產(chǎn)品的品質(zhì)。

    某大氣監(jiān)測(cè)光譜儀采用時(shí)間調(diào)制型傅里葉變換分光的技術(shù)路線,光譜儀主要由二維指向機(jī)構(gòu)、干涉儀組件、后光學(xué)組件、分色匯聚光學(xué)組件、探測(cè)器及信號(hào)處理系統(tǒng)等組成,其核心組件為一臺(tái)雙角鏡擺臂式干涉儀。入射光經(jīng)指向鏡進(jìn)入干涉儀組件,通過(guò)干涉儀的光程掃描實(shí)現(xiàn)干涉分光后由后光學(xué)組件進(jìn)行壓束,經(jīng)分色匯聚光學(xué)組件分為4個(gè)通道并分別匯聚至各譜段的單元探測(cè)器上,探測(cè)器獲取的干涉信號(hào)經(jīng)傅里葉變換后可得到大氣成分對(duì)地表反射太陽(yáng)輻射的吸收光譜,從而定量反演各成分的濃度分布信息,應(yīng)用于全球變暖等氣候變化問(wèn)題的研究。

    光譜儀在近紅外和短波紅外波段設(shè)置了4個(gè)工作譜段,分別為0.75~0.77μm、1.56~1.72μm、1.92~2.08μm、2.20~2.38μm,其近紅外和短波紅外譜段的光譜分辨率分別優(yōu)于0.6cm–1和0.27cm–1。光譜儀視場(chǎng)角為15.8mrad,絕對(duì)輻射定標(biāo)精度要求優(yōu)于5%。為保證全壽命周期內(nèi)光譜數(shù)據(jù)的精度,需設(shè)置星上定標(biāo)系統(tǒng),定期對(duì)光譜儀的光譜和輻射性能進(jìn)行標(biāo)定。

    2 星上定標(biāo)方案設(shè)計(jì)

    2.1 星上輻射定標(biāo)

    (1)星上輻射定標(biāo)方案設(shè)計(jì)

    發(fā)射過(guò)程中的振動(dòng)、在軌運(yùn)行期間的空間粒子輻射、污染、元器件老化等因素都可能引起光譜儀光譜響應(yīng)特性的衰減和變化[11],使得實(shí)驗(yàn)室定標(biāo)系數(shù)在衛(wèi)星在軌飛行期間可能不再適用,如GOSAT(Greenhouse Gases Observing Satellite)衛(wèi)星搭載的TANSO-FTS光譜儀入軌后3年相對(duì)于發(fā)射前定標(biāo)的輻射衰減在近紅外譜段約18%~20%,短波紅外譜段約5%~7%[12]。對(duì)于工作在太陽(yáng)反射譜段的遙感器,多利用太陽(yáng)這一長(zhǎng)期穩(wěn)定的寬帶光源,采用太陽(yáng)漫反射板定標(biāo)法對(duì)儀器入軌后的輻射性能進(jìn)行標(biāo)定,該方法具有全口徑、全視場(chǎng)、全光路等優(yōu)點(diǎn),具有很高的輻射定標(biāo)精度。然而,漫反射板會(huì)受空間紫外輻照、高能粒子轟擊等發(fā)生性能衰變,如聚四氟乙烯材料在紫外和可見(jiàn)光波段衰減明顯。因此,為保證全壽命周期內(nèi)遙感器的定標(biāo)精度,除了采用空間環(huán)境下具有高穩(wěn)定性的漫反射材料,需建立對(duì)漫反射板反射性能衰減的在軌監(jiān)測(cè)手段。

    解決方法之一是設(shè)計(jì)一種比輻射計(jì),通過(guò)比較太陽(yáng)輻照度和漫反射板反射輻亮度進(jìn)行漫反射板BRDF的變化監(jiān)測(cè)[7],該方法的局限在于僅可對(duì)部分譜段以及漫反射板的部分區(qū)域和部分觀測(cè)角度下的BRDF進(jìn)行監(jiān)測(cè),與星上定標(biāo)的實(shí)際狀態(tài)存在偏差。對(duì)于口徑相對(duì)較小的遙感器,在體積質(zhì)量等條件允許的情況下,可設(shè)置參考漫反射板,通過(guò)兩塊漫反射板定標(biāo)結(jié)果的比對(duì)實(shí)現(xiàn)漫反射板在軌衰減的監(jiān)測(cè)和校正。該校正方法與常規(guī)定標(biāo)的測(cè)試狀態(tài)完全一致,精度更高。

    本文光譜儀的星上絕對(duì)輻射定標(biāo)利用漫反射板引入太陽(yáng)輻射作為定標(biāo)源,采用兩塊180°配置的漫反射板,一塊用于每軌一次的太陽(yáng)定標(biāo),另一塊作為參考板,使用頻次設(shè)計(jì)為每月一次,用于與常規(guī)板的比對(duì)校正。

    (2)漫反射板的選用、設(shè)計(jì)及測(cè)試

    隨著遙感器設(shè)計(jì)壽命的延長(zhǎng),漫反射板的選用需充分考慮材料的在軌衰減情況。QVD漫反射板材料為石英玻璃,對(duì)空間輻射不敏感,在軌性能的衰減僅由表面污染引起,相比聚四氟乙烯(如Spectralon)漫反射板具有極高的在軌穩(wěn)定性。QVD漫反射板可通過(guò)漫反射板表面粗糙度的加工,滿足不同角度、紫外至短波紅外譜段不同光譜的BRDF要求,角度均勻性好[13]。

    采用兩塊相同的QVD漫反射板,設(shè)計(jì)直徑為150mm,置于光譜儀光路的最前端,可保證全視場(chǎng)、全口徑和全光路覆蓋。漫反射板前后表面均具有一定粗糙度,后表面鍍鋁膜,以實(shí)現(xiàn)所需角度下BRDF的角均勻性和面均勻性要求。

    漫反射板的BRDF與入射角度相關(guān),光譜儀在軌定標(biāo)時(shí)的太陽(yáng)入射角度需考慮衛(wèi)星軌道參數(shù)及其在軌漂移情況、地氣雜光影響、平臺(tái)結(jié)構(gòu)遮擋等因素,并結(jié)合這些因素確定太陽(yáng)定標(biāo)的時(shí)機(jī)。考慮上述因素后,設(shè)計(jì)全壽命周期內(nèi)太陽(yáng)定標(biāo)時(shí)段對(duì)應(yīng)的太陽(yáng)光與漫反射板的夾角關(guān)系如圖1所示,其中為衛(wèi)星坐標(biāo)系,為衛(wèi)星飛行方向,為穿軌方向,為星下點(diǎn)方向。太陽(yáng)光相對(duì)面的夾角的變化范圍為20°~35°,太陽(yáng)光在平面內(nèi)的投影與軸的夾角的變化范圍為10°~27°。漫反射板法線平行于-軸,則定標(biāo)期間太陽(yáng)光與定標(biāo)漫反射板法線夾角,也就是太陽(yáng)光入射角變化范圍為55°~70°。

    圖1 星上太陽(yáng)定標(biāo)期間太陽(yáng)光與定標(biāo)漫反射板的夾角關(guān)系

    針對(duì)光譜儀的工作譜段和在軌太陽(yáng)定標(biāo)角度,對(duì)漫反射板表面粗糙度設(shè)計(jì)進(jìn)行了優(yōu)化,如圖2所示。可以看出,當(dāng)粗糙度大于2μm后的漫反射板BRDF角度均勻性較高,粗糙度2μm至6.5μm的漫反射板BRDF的相對(duì)變化在5%~7%以?xún)?nèi)。進(jìn)一步地,利用顯微鏡和探測(cè)器對(duì)面均勻性進(jìn)行測(cè)定,得出粗糙度2μm的漫反射板的面均勻性最佳。另一方面,定標(biāo)漫反射板的BRDF設(shè)計(jì)需滿足太陽(yáng)定標(biāo)時(shí)經(jīng)定標(biāo)漫反射板反射后的太陽(yáng)光譜輻亮度與光譜儀各譜段的平均入射光譜輻亮度相匹配,且不超出光譜儀的動(dòng)態(tài)范圍上限,實(shí)測(cè)結(jié)果如圖3所示。

    圖2 不同表面粗糙度漫反射板的BRDF相對(duì)角度變化情況

    圖3 漫反射板出射輻亮度與光譜儀動(dòng)態(tài)范圍符合情況

    在星上定標(biāo)系統(tǒng)安裝至光譜儀整機(jī)后,考慮到裝配誤差和雜光等因素的影響,需在實(shí)驗(yàn)室對(duì)定標(biāo)漫反射板在不同波長(zhǎng)和不同太陽(yáng)光入射角條件下進(jìn)行系統(tǒng)級(jí)BRDF標(biāo)定[14-15]。使用太陽(yáng)模擬器作為光源,使其分別照射光譜儀的星上定標(biāo)通道和對(duì)地通道,其中對(duì)地通道采用一已知BRDF的標(biāo)準(zhǔn)漫反射板將入射光引入,利用兩次測(cè)量的比值獲得漫反射板的系統(tǒng)級(jí)BRDF,更為接近星上定標(biāo)系統(tǒng)的在軌真實(shí)工作狀態(tài),相比利用元件級(jí)的漫反射板BRDF進(jìn)行星上定標(biāo),可獲得更高的星上輻射定標(biāo)精度。經(jīng)分析(見(jiàn)表1),光譜儀所采用的系統(tǒng)級(jí)BRDF標(biāo)定方法的測(cè)量誤差可優(yōu)于3%。

    表1 光譜儀短波紅外譜段定標(biāo)漫反射板系統(tǒng)級(jí)BRDF的測(cè)量誤差

    Tab.1 Instrument BRDF error of calibration diffuser of SWIR channel

    (3)星上絕對(duì)輻射定標(biāo)精度分析

    光譜儀的星上絕對(duì)輻射定標(biāo)是通過(guò)已知的大氣層外太陽(yáng)光譜照度數(shù)據(jù)、在軌實(shí)時(shí)太陽(yáng)矢量(太陽(yáng)光入射角度)和漫反射板BRDF的地面標(biāo)定結(jié)果確立輻亮度標(biāo)準(zhǔn),并在衛(wèi)星處于地影區(qū)時(shí)進(jìn)行深空暗背景數(shù)據(jù)的采集。在某一波數(shù)處,太陽(yáng)光經(jīng)定標(biāo)漫反射板散射后入射到光譜儀的輻亮度S()為:

    式中()為太陽(yáng)輻照度;BRDF(iiss)為定標(biāo)漫反射板的雙向反射分布函數(shù),它是入射光照度和反射光亮度的比值,與入射光波長(zhǎng)、入射光線和出射光線的角度有關(guān);i為太陽(yáng)光和定標(biāo)漫反射板法線的夾角;s為出射光與漫反射板法線的夾角,i與s主分別為入射光和出射光的方位角[16]。太陽(yáng)定標(biāo)時(shí)光譜儀采集信號(hào)為DNS();深空暗背景觀測(cè)時(shí)采集信號(hào)為DN0(),可得到光譜儀的輻射定標(biāo)方程:

    式中 DN()為光譜儀對(duì)地觀測(cè)時(shí)的測(cè)量信號(hào);()為與DN()對(duì)應(yīng)的對(duì)地觀測(cè)時(shí)的入瞳輻亮度。

    星上輻射定標(biāo)需利用定標(biāo)數(shù)據(jù)所對(duì)應(yīng)的在軌實(shí)測(cè)參數(shù),如衛(wèi)星姿態(tài)信息、太陽(yáng)光相對(duì)漫反射板的入射和出射角度,以及參考漫反射板的定標(biāo)結(jié)果等,對(duì)定標(biāo)模型進(jìn)行修正。修正后仍存在的殘余誤差及其它不可修正項(xiàng)即構(gòu)成了星上輻射定標(biāo)精度。以光譜儀短波紅外通道(1.6μm)為例,絕對(duì)輻射定標(biāo)環(huán)節(jié)中的主要誤差來(lái)源及參考貢獻(xiàn)值分析如表2所示。其中,太陽(yáng)照度數(shù)據(jù)來(lái)源于已有科學(xué)測(cè)量數(shù)據(jù),目前其測(cè)量不確定度可優(yōu)于2%;太陽(yáng)定標(biāo)時(shí)和對(duì)地觀測(cè)時(shí)的光譜輸出信號(hào)誤差主要源自于儀器的噪聲影響;太陽(yáng)光入射角度誤差主要受衛(wèi)星姿態(tài)和星上實(shí)時(shí)太陽(yáng)矢量計(jì)算誤差的影響,會(huì)導(dǎo)致入射照度和漫反射板BRDF的不確定,對(duì)于朗伯性很好的漫反射板該部分誤差影響很??;雜散光的影響已通過(guò)系統(tǒng)級(jí)BRDF標(biāo)定有所抑制,并采用遮光罩、擋光板等措施抑制在軌星表反射太陽(yáng)光所產(chǎn)生的雜光,該項(xiàng)誤差主要考慮地面測(cè)試和分析與在軌真實(shí)狀態(tài)的差異;系統(tǒng)非線性誤差來(lái)源于干涉信號(hào)零光程差處峰值信號(hào)所引入的非線性調(diào)制因子;根據(jù)OMI的QVD漫反射板反射性能的在軌衰減數(shù)據(jù),預(yù)估本方案參考漫反射板BRDF的8年衰減小于0.5%。綜合分析可知,光譜儀短波紅外通道(1.6μm)壽命末期星上絕對(duì)輻射定標(biāo)精度為4.33%。此外,譜段1、譜段3和譜段4的星上絕對(duì)輻射定標(biāo)精度分別為3.91%、4.36%和4.39%,均滿足5%的指標(biāo)要求。

    表2 光譜儀短波紅外譜段星上絕對(duì)輻射定標(biāo)誤差分析

    Tab.2 On-board absolute radiometric calibration accuracy of SWIR channel

    2.2 星上ILS測(cè)量

    儀器線形函數(shù)(Instrument Line Shape, ILS)表征了光譜儀對(duì)于一個(gè)給定波長(zhǎng)的單色光源的光譜響應(yīng),光譜儀的測(cè)量光譜是真實(shí)光譜與ILS卷積的結(jié)果。對(duì)于具有超高光譜分辨率且絕對(duì)精度要求很高的儀器,其輻射定標(biāo)和波長(zhǎng)定標(biāo)都需要將基準(zhǔn)光譜與ILS卷積,卷積后的結(jié)果與光譜儀實(shí)測(cè)信號(hào)比對(duì)來(lái)完成定標(biāo)參數(shù)的計(jì)算。短波紅外譜段的大氣光譜數(shù)據(jù)反演對(duì)ILS的變化非常敏感,需要設(shè)置星上ILS測(cè)量手段。

    本方案設(shè)計(jì)利用波長(zhǎng)為1 664.5nm的窄線寬半導(dǎo)體激光器作為輸入光源進(jìn)行ILS測(cè)量,激光線寬優(yōu)于10MHz,遠(yuǎn)小于光譜儀的光譜分辨率,可保證測(cè)量結(jié)果的準(zhǔn)確性。定標(biāo)激光器內(nèi)部設(shè)置了半導(dǎo)體制冷器和熱敏電阻,對(duì)定標(biāo)激光器進(jìn)行制冷和精確控溫,以保證輸出激光的頻率穩(wěn)定性。

    圖4 利用星上激光器測(cè)試的ILS光譜

    激光器發(fā)出的單色穩(wěn)頻激光由星上定標(biāo)漫反射板引入光譜儀光學(xué)系統(tǒng),將光譜儀獲取的單色激光干涉信號(hào)進(jìn)行逆傅里葉變換后即獲得經(jīng)光譜儀展寬后的輸入激光光譜,即光譜儀的ILS,如圖4所示。另外,ILS函數(shù)的半高寬可反映光譜儀的光譜分辨率,光譜分辨率實(shí)測(cè)結(jié)果為0.262cm–1,滿足指標(biāo)要求。ILS的測(cè)量精度主要由激光器溫度和輸出功率的不穩(wěn)定度和光譜儀信號(hào)噪聲決定,激光器溫度和輸出功率的不穩(wěn)定度<0.1%,光譜儀噪聲<0.7%,ILS的測(cè)量精度約為0.7%。

    2.3 在軌波長(zhǎng)校正

    光譜儀的實(shí)驗(yàn)室光譜定標(biāo)多采用特征譜線燈[17]、單色激光器、氣體池結(jié)合不同光源(積分球、黑體等)的方法[18],可實(shí)現(xiàn)的光譜定標(biāo)精度很高,但在星上應(yīng)用方面受限于功耗、體積、質(zhì)量及壽命等要求,工程實(shí)現(xiàn)存在難度。星上光譜定標(biāo)有采用鹵鎢燈[19]和稀土摻雜漫反射板[20-21]的方式,但摻雜漫射板的光譜特征分辨率較低,不適用于具有超高光譜分辨率的遙感器。目前,國(guó)外已發(fā)射的多個(gè)超光譜載荷,如文獻(xiàn)[22]、文獻(xiàn)[23]多借助太陽(yáng)光譜和大氣吸收光譜進(jìn)行其在軌波長(zhǎng)定標(biāo)。

    光譜儀處于太陽(yáng)定標(biāo)模式時(shí),利用太陽(yáng)定標(biāo)獲得的光譜信號(hào)里的太陽(yáng)夫瑯禾費(fèi)線的波長(zhǎng)信息進(jìn)行在軌波長(zhǎng)校正;光譜儀工作在常規(guī)觀測(cè)模式時(shí),通過(guò)對(duì)地面均勻場(chǎng)景進(jìn)行觀測(cè),利用探測(cè)光譜信號(hào)里的大氣吸收線的波長(zhǎng)信息進(jìn)行在軌波長(zhǎng)校正。太陽(yáng)光譜中的特征吸收線較大氣吸收光譜更為尖銳,獲得的波長(zhǎng)位置更為精確,而大氣吸收光譜無(wú)Doppler效應(yīng)的影響,兩種方法相結(jié)合,可得到更為準(zhǔn)確的光譜校正結(jié)果。

    使用太陽(yáng)光譜和大氣吸收光譜進(jìn)行星上波長(zhǎng)定標(biāo)計(jì)算時(shí),首先使用已知的太陽(yáng)光譜或大氣吸收光譜與光譜儀的ILS進(jìn)行卷積,將卷積后的已知光譜與光譜儀實(shí)測(cè)得到的測(cè)量光譜采用標(biāo)準(zhǔn)差法進(jìn)行比較后確定測(cè)量光譜各采樣點(diǎn)對(duì)應(yīng)的波長(zhǎng)值,完成波長(zhǎng)定標(biāo)。

    圖5 星上波長(zhǎng)定標(biāo)時(shí)的數(shù)據(jù)處理示意

    星上波長(zhǎng)定標(biāo)時(shí)的數(shù)據(jù)處理示意如圖5所示。沿波數(shù)軸按一定步長(zhǎng)對(duì)測(cè)量光譜進(jìn)行逐次平移,逐次計(jì)算平移后的測(cè)量光譜與已知光譜的標(biāo)準(zhǔn)差,標(biāo)準(zhǔn)差為最小時(shí)的測(cè)量光譜與波數(shù)軸的對(duì)應(yīng)關(guān)系即為星上波長(zhǎng)定標(biāo)結(jié)果。標(biāo)準(zhǔn)差的計(jì)算公式為,

    式中為采樣點(diǎn)數(shù);ν為采樣點(diǎn)對(duì)應(yīng)的波數(shù);0()為已知光譜數(shù)據(jù);()為測(cè)量光譜數(shù)據(jù)。

    3 星上定標(biāo)系統(tǒng)設(shè)計(jì)和工作過(guò)程描述

    根據(jù)上述星上定標(biāo)方案,設(shè)計(jì)的光譜儀星上定標(biāo)系統(tǒng)組成如圖6所示,由漫反射板定標(biāo)機(jī)構(gòu)、激光器和定標(biāo)控制器組成,定標(biāo)控制器用于機(jī)構(gòu)轉(zhuǎn)動(dòng)和激光器電流的控制。星上定標(biāo)流程示意圖如圖7所示,為避免太陽(yáng)光的影響衛(wèi)星處于地影區(qū)內(nèi)進(jìn)行ILS測(cè)量和深空冷背景數(shù)據(jù)采集。衛(wèi)星出地影區(qū)后、星下點(diǎn)日出前進(jìn)行太陽(yáng)定標(biāo)??紤]氣體反演精度對(duì)每軌定標(biāo)數(shù)據(jù)的精度要求,設(shè)計(jì)對(duì)太陽(yáng)和深空冷背景的觀測(cè)每軌進(jìn)行,每月進(jìn)行一次ILS測(cè)量。

    圖6 星上定標(biāo)系統(tǒng)組成示意

    圖7 星上定標(biāo)流程示意

    兩塊漫反射板的切換由定標(biāo)機(jī)構(gòu)的轉(zhuǎn)動(dòng)實(shí)現(xiàn),定標(biāo)時(shí)光譜儀的指向機(jī)構(gòu)指向漫反射板方向,將太陽(yáng)輻射引入干涉儀。參考板不使用時(shí),受定標(biāo)機(jī)構(gòu)外殼保護(hù),與空間環(huán)境隔離,以保證其在軌所受空間輻射和污染的影響降至最低。定標(biāo)機(jī)構(gòu)采用電機(jī)直驅(qū)方式,采用2組微動(dòng)開(kāi)關(guān)判斷漫反射板工作位置,結(jié)構(gòu)形式簡(jiǎn)單、可靠性高。圖8為定標(biāo)機(jī)構(gòu)組成和設(shè)計(jì)示意。

    圖8 星上定標(biāo)機(jī)構(gòu)組成

    對(duì)應(yīng)用于大氣成分定量監(jiān)測(cè)的超光譜載荷而言,星上絕對(duì)輻射定標(biāo)是保證氣體濃度反演精度的關(guān)鍵,而儀器線形函數(shù)測(cè)量和波長(zhǎng)校正是進(jìn)行輻射定標(biāo)的前提。光譜儀的星上定標(biāo)工作原理如圖9所示,定標(biāo)采集的光譜數(shù)據(jù)需首先進(jìn)行ILS卷積處理和波長(zhǎng)校正,而后才能進(jìn)行輻射定標(biāo)數(shù)據(jù)的處理。通過(guò)星上ILS測(cè)量與發(fā)射前ILS的比對(duì),可監(jiān)測(cè)光譜儀在軌光譜響應(yīng)特性的變化情況。根據(jù)星上太陽(yáng)定標(biāo)數(shù)據(jù)和暗背景測(cè)量數(shù)據(jù)獲得星上絕對(duì)輻射定標(biāo)方程,并與發(fā)射前定標(biāo)方程進(jìn)行比對(duì),對(duì)發(fā)射前定標(biāo)參數(shù)進(jìn)行修正,同時(shí)評(píng)估漫反射板和儀器的在軌衰減或變化情況。

    圖9 光譜儀星上定標(biāo)原理示意

    4 結(jié)束語(yǔ)

    星載大氣監(jiān)測(cè)光譜儀星上定標(biāo)方案的設(shè)計(jì)在總結(jié)和借鑒國(guó)外先進(jìn)載荷在軌應(yīng)用經(jīng)驗(yàn)的基礎(chǔ)上,結(jié)合傅里葉變換光譜載荷的技術(shù)特點(diǎn),設(shè)計(jì)的星上定標(biāo)系統(tǒng)具備全口徑、全視場(chǎng)、全光路絕對(duì)輻射定標(biāo)和ILS測(cè)量能力。漫反射板定標(biāo)裝置具有在軌高穩(wěn)定性和可靠性,漫反射板的系統(tǒng)級(jí)BRDF可降低定標(biāo)裝置裝配誤差和系統(tǒng)雜光的影響,可保證全壽命周期光譜儀各譜段絕對(duì)輻射定標(biāo)精度均優(yōu)于5%。ILS卷積作為光譜數(shù)據(jù)處理的必要環(huán)節(jié),星上窄線寬穩(wěn)頻激光器的設(shè)計(jì)可滿足光譜儀短波紅外譜段的星上ILS測(cè)量要求,并可在軌實(shí)時(shí)評(píng)估光譜儀的光譜分辨率等光學(xué)性能。配備滿足高精度反演要求、在軌穩(wěn)定性好的星上定標(biāo)裝置,是定量遙感載荷發(fā)展的重要趨勢(shì),本文介紹的高精度星上定標(biāo)技術(shù)可應(yīng)用于工作在太陽(yáng)反射譜段的星載光譜儀等具有高精度定量應(yīng)用需求的遙感器,為我國(guó)的環(huán)境綜合監(jiān)測(cè)、氣候變化研究、資源勘查、防災(zāi)減災(zāi)等行業(yè)提供高質(zhì)量、高可靠度的高光譜遙感數(shù)據(jù)服務(wù)。

    [1] CRISP D, FISHER B M, O’DELL C, et al. The ACOS CO2Retrieval Algorithm-Part II: Global XCO2Data Characterization[J]. Atmos. Meas. Tech, 2012, 5: 687-707.

    [2] DOBBER M, DIRKSEN R, LEVELT P, et al. Ozone Monitoring Instrument Calibration[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(5): 1209-1238.

    [3] KUZE A, O’BRIEN D M, TAYLOR T E, et al. Vicarious Calibration of the GOSAT Sensors Using the Railroad Valley Desert Playa[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(5): 1781-1795.

    [4] DOBBER M, DIRKSEN R, LEVELT P, et al. Ozone Monitoring Instrument In-flight Performance and Calibration[C]. Optical Design and Engineering II. Jena: SPIE, 2005.

    [5] CALLIES J, CORPACCIOLI E, EISINGER M, et al. GOME-2 Ozone Instrument On-board the European METOP Satellites[C]. Weather and Environmental Satellites. Bellingham: SPIE, 2004.

    [6] 李曉暉, 顏昌翔. 成像光譜儀星上定標(biāo)技術(shù)[J]. 中國(guó)光學(xué)與應(yīng)用光學(xué), 2009, 2(4): 309-315. LI Xiaohui, YAN Changxiang. Onboard Calibration Technologies for Hyper-Spectral Imager[J]. Chinese Journal of Optics and Applied Optics, 2009, 2(4): 309-315. (in Chinese)

    [7] 趙艷華, 董建婷, 張秀茜, 等. 漫反射板全光路全視場(chǎng)全口徑在軌輻射定標(biāo)技術(shù)[J]. 航天返回與遙感, 2016, 37(2): 92-99. ZHAO Yanhua, DONG Jianting, ZHANG Xiuqian, et al. In-orbit Radiometric Calibration Technology Based on Diffuse Reflection Plate in Full Bore, Full FOV and Full Optical Path[J]. Spacecraft Recovery & Remote Sensing, 2016, 37(2): 92-99. (in Chinese)

    [8] 趙敏杰, 司福祺, 江宇, 等. 星載大氣痕量氣體差分吸收光譜儀定標(biāo)機(jī)構(gòu)設(shè)計(jì)[J]. 大氣與環(huán)境光學(xué)學(xué)報(bào), 2014, 9(1): 57-60. ZHAO Minjie, SI Fuqi, JIANG Yu, et al. Design of Calibration Mechanism on Space-Borne Differential Optical Absorption Spectrometer[J]. Journal of Atmospheric and Environmental Optics, 2014, 9(1): 57-60. (in Chinese)

    [9] 范斌, 陳旭, 李碧岑, 等. “高分五號(hào)”衛(wèi)星光學(xué)遙感載荷的技術(shù)創(chuàng)新[J]. 紅外與激光工程, 2017, 46(1): 0102002-1-7. FAN Bin, CHEN Xu, LI Bicen, et al. Technical Innovation of Optical Remote Sensing Payloads Onboard GF-5 Satellite[J]. Infrared and Laser Engineering, 2017, 46(1): 0102002-1-7. (in Chinese)

    [10] 甘偉偉. 溫室氣體監(jiān)測(cè)儀全過(guò)程定標(biāo)軟件開(kāi)發(fā)與實(shí)現(xiàn)[D]. 合肥: 中國(guó)科技大學(xué), 2017. GAN Weiwei. Development and Implementation of Whole Process Calibration Software for The Greenhousegases Monitoring Instrument[D]. Hefei: University of Science and Technology of China, 2017. (in Chinese)

    [11] 楊本永, 張黎明, 楊寶云, 等. 光學(xué)傳感器星上定標(biāo)漫射板空間穩(wěn)定性測(cè)試與分析[J]. 中國(guó)激光, 2009, 36(11): 3026-3030. YANG Benyong, ZHANG Liming, YANG Baoyun, et al. Space Stability Measurement and Analysis of Diffuser Panel Used for On-Board Calibration for Satellite Optical Sensors[J]. Chinese Journal of Lasers, 2009, 36(11): 3026-3030. (in Chinese)

    [12] KUZE A, TAYLOR T E, KATAOKA F, et al. Long-Term Vicarious Calibration of GOSAT Short-Wave Sensors: Techniques for Error Reduction and New Estimates of Radiometric Degradation Factors[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(7): 3991-4004.

    [13] BILGEHAN G, HANS B, PENGMEI X, BICEN L. Solar Diffusers in Earth Observation Instruments with an Illumination Angle of up to 70°: Design and Verification of Performance in BRDF[C]. Conference on Earth Observing Systems. San Diego: SPIE, 2015.

    [14] ROBERT A B, STENVEN W B, KEITH R L, et al. System-level Pre-launch Calibration of On-board Solar Diffusers[C]. Earth Observing Systems XII. San Diego: SPIE, 2007.

    [15] 李明, 宗肖穎. 定標(biāo)漫反射板實(shí)驗(yàn)室系統(tǒng)級(jí)BRDF測(cè)量方法[J]. 紅外與激光工程, 2017, 46(1): 0117004-1-7. LI Ming, ZONG Xiaoying. In-lab System-level BRDF Measurement Method of Calibration Diffuser[J]. Infrared and Laser Engineering, 2017, 46(1): 0117004-1-7. (in Chinese)

    [16] GEORGI T G, JAMES J B. BRDF Study of Gray-scale Spectralon[C]. SPIE Conference on Earth Observing Systems. San Diego: SPIE, 2008.

    [17] 周海金, 劉文清, 司福祺, 等. 星載大氣痕量氣體差分吸收光譜儀光譜定標(biāo)技術(shù)研究[J]. 光譜學(xué)與光譜分析, 2012, 32(11): 2881-2885. ZHOU Haijin, LIU Wenqing, SI Fuqi, et al. Spectral Calibration for Space-borne Differential Optical Absorption Spectrometer[J]. Spectroscopy and Spectral Analysis, 2012, 32(11): 2881-2885. (in Chinese)

    [18] 劉倩倩, 鄭玉權(quán). 超高分辨率光譜定標(biāo)技術(shù)發(fā)展概況[J]. 中國(guó)光學(xué), 2012, 5(6): 566-577. LIU Qianqian, ZHENG Yuquan. Development of Spectral Calibration Technologies with Ultra-high Resolutions[J]. Chinese Optics, 2012, 5(6): 566-577. (in Chinese)

    [19] 高靜, 計(jì)忠瑛, 王忠厚, 等. 空間調(diào)制干涉光譜成像儀的星上定標(biāo)系統(tǒng)穩(wěn)定性研究[J]. 光譜學(xué)與光譜分析, 2010, 30(4): 1013-1017. GAO Jing, JI Zhongying, WANG Zhonghou, et al. Stability Analysis of On-board Calibration System of Spatial Modulated Imaging Fourier Transform Spectrometer[J]. Spectroscopy and Spectral Analysis, 2010, 30(4): 1013-1017. (in Chinese)

    [20] BEATRICE C, CHRISTINE G, IVAN C, et al. MERIS FM Performances[C]. Earth Observing Systems III. San Diego: SPIE, 1998.

    [21] 劉宇翔, 伏瑞敏, 李明, 等. 稀土摻雜漫反射板星上光譜定標(biāo)技術(shù)[J]. 航天返回與遙感, 2015, 36(6): 48-54. LIU Yuxiang, FU Ruimin, LI Ming, et al. On-board Spectral Calibration Technologies with Rare Earth Doped Diffuser Plate[J]. Spacecraft Recovery & Remote Sensing, 2015, 36(6): 48-54. (in Chinese)

    [22] VRIES J D, HOOGEVEEN R, VOORS R, et al. Technology Evolution of the TROPOMI Instrument[C]. IEEE International Geoscience and Remote Sensing Symposium. Munich: IEEE, 2012.

    [23] VRIES J D,ROBERT V, RUUD D, et al. In-orbit Performance of the Ozone Monitoring Instrument[C]. Sensors, Systems, and Next-Generation Satellites IX. Druges: SPIE, 2005.

    High Precision On-board Calibration for Spaceborne Atmospheric Monitoring Spectrometer

    LI Bicen LI Ming XU Pengmei

    (Beijing Institute of Space Mechanics & Electricity, Beijing 100094, China)

    Radiometric accuracy and spectral accuracy should be very high for spaceborne hyperspectral payload applied in atmospheric composition retrieval. With extension of the lifetimes of spaceborne optical sensors, it is required to establish an on-board calibration system with high accuracy and high stability. The calibration techniques which can meet the high on-board calibration accuracy for atmospheric monitoring spectrometer are described in this paper. According to the characteristics of time modulated Fourier transform spectrometer, the methods of absolute radiometric calibration, instrument line shape (ILS) measurement and spectral calibration are presented. Solar diffuser is used to realize the absolute radiometric calibration on orbit covering full aperture, full field of view and whole detecting link of instrument. The calibration diffuser in front of the optical path of the spectrometer reflects the sunlight into the instrument, and the standard radiance is obtained using the known exoatmospheric solar irradiance and the BRDF (bidirectional reflectance distribution function) of the diffuser calibrated on ground. And another diffuser as a calibration reference is used to monitor and correct the degradation and variation of the calibration unit and the instrument. The absolute radiometric calibration accuracy in total life time is better than 5% by use of Quasi Volume Diffuser (QVD) which have the advantage of very stable optical performance in space. Monochromatic lasers are employed for regular ILS measurement from which the on-orbit status of optical performance like spectral resolution can be evaluated. Spectral lines in solar spectrum and atmospheric spectrum are used to correct the wavelength shift of instrument.

    diffuser; instrument line shape; wavelength correction; absolute radiometric calibration; atmospheric monitoring spectrometer; space remote sensing; GF-5 satellite

    V447+.1

    A

    1009-8518(2018)03-0060-10

    10.3969/j.issn.1009-8518.2018.03.007

    李碧岑,女,1985年生,2010年獲西安電子科技大學(xué)光學(xué)專(zhuān)業(yè)碩士學(xué)位,工程師。研究方向?yàn)楹教旃鈱W(xué)遙感器總體設(shè)計(jì)、定標(biāo)技術(shù)。E-mail:mou_lbc@163.com。

    2018-03-21

    國(guó)家重大科技專(zhuān)項(xiàng)工程

    (編輯:劉穎)

    猜你喜歡
    譜段光譜儀定標(biāo)
    我國(guó)為世界大豆精準(zhǔn)選種“定標(biāo)”
    基于恒星的電離層成像儀在軌幾何定標(biāo)
    高分六號(hào)衛(wèi)星WFV 新增譜段對(duì)農(nóng)作物識(shí)別精度的改善
    一種基于波長(zhǎng)生長(zhǎng)的紅外預(yù)警譜段選擇方法
    激光與紅外(2018年2期)2018-03-09 07:27:57
    基于配準(zhǔn)和二次碼率控制的光譜圖像壓縮
    基于角反射器的機(jī)載毫米波云雷達(dá)外定標(biāo)實(shí)驗(yàn)
    中階梯光柵光譜儀自動(dòng)化波長(zhǎng)提取算法
    醫(yī)用顯微成像光譜儀的光譜定標(biāo)技術(shù)
    4m直徑均勻擴(kuò)展定標(biāo)光源
    雙面鏡結(jié)構(gòu)成像光譜儀的改進(jìn)及性能分析
    日本与韩国留学比较| eeuss影院久久| 夫妻性生交免费视频一级片| 午夜免费男女啪啪视频观看| 久久久a久久爽久久v久久| 啦啦啦啦在线视频资源| 国产精品美女特级片免费视频播放器| 听说在线观看完整版免费高清| 2022亚洲国产成人精品| 国产精品1区2区在线观看.| 亚洲人成网站高清观看| 又爽又黄无遮挡网站| 欧美色视频一区免费| 亚洲精品国产av成人精品| 精品久久久久久久末码| 少妇的逼好多水| 精品久久久久久久久亚洲| 男的添女的下面高潮视频| 女同久久另类99精品国产91| 国产伦一二天堂av在线观看| 一级毛片aaaaaa免费看小| 男女那种视频在线观看| 精品人妻一区二区三区麻豆| 婷婷色av中文字幕| 午夜福利在线观看吧| 国产一区二区三区在线臀色熟女| 99热只有精品国产| 久久久久久久久久黄片| 亚洲美女视频黄频| 国产午夜精品久久久久久一区二区三区| 美女脱内裤让男人舔精品视频 | 亚洲第一电影网av| 99热6这里只有精品| 亚洲国产精品成人久久小说 | 人人妻人人看人人澡| 少妇人妻精品综合一区二区 | 日韩高清综合在线| 精品人妻视频免费看| 91精品一卡2卡3卡4卡| 国产精品永久免费网站| 最近的中文字幕免费完整| 久久国内精品自在自线图片| 性色avwww在线观看| 一级二级三级毛片免费看| 亚洲av免费在线观看| 18禁黄网站禁片免费观看直播| 日韩一本色道免费dvd| 边亲边吃奶的免费视频| 国产高潮美女av| 高清午夜精品一区二区三区 | 哪里可以看免费的av片| 国产高清视频在线观看网站| 久久亚洲精品不卡| 91aial.com中文字幕在线观看| 你懂的网址亚洲精品在线观看 | 国内精品久久久久精免费| 夜夜夜夜夜久久久久| 男人和女人高潮做爰伦理| 一级毛片我不卡| 欧美成人a在线观看| 一区二区三区四区激情视频 | 日韩欧美国产在线观看| 中文字幕免费在线视频6| 蜜桃久久精品国产亚洲av| 晚上一个人看的免费电影| 啦啦啦啦在线视频资源| a级毛色黄片| 人妻系列 视频| 亚洲国产色片| 乱系列少妇在线播放| 毛片女人毛片| 亚洲va在线va天堂va国产| 18禁裸乳无遮挡免费网站照片| 欧美成人a在线观看| 日日干狠狠操夜夜爽| 美女黄网站色视频| 51国产日韩欧美| 午夜精品在线福利| 欧美成人一区二区免费高清观看| 久久亚洲国产成人精品v| videossex国产| 69人妻影院| 日日摸夜夜添夜夜爱| 久久午夜福利片| 欧美成人一区二区免费高清观看| 午夜福利视频1000在线观看| 久久国内精品自在自线图片| 1024手机看黄色片| 能在线免费观看的黄片| 久99久视频精品免费| 国产中年淑女户外野战色| 欧美xxxx性猛交bbbb| 最近2019中文字幕mv第一页| 给我免费播放毛片高清在线观看| 99久国产av精品国产电影| 亚洲精品久久国产高清桃花| 18禁在线无遮挡免费观看视频| 18禁裸乳无遮挡免费网站照片| 欧美成人一区二区免费高清观看| 春色校园在线视频观看| 亚洲国产精品成人久久小说 | 黄色一级大片看看| 成年版毛片免费区| 免费av不卡在线播放| 欧美区成人在线视频| 九九在线视频观看精品| 插逼视频在线观看| 一个人看的www免费观看视频| 九草在线视频观看| 国产成人freesex在线| 大香蕉久久网| 不卡视频在线观看欧美| 两性午夜刺激爽爽歪歪视频在线观看| 色5月婷婷丁香| 国产伦一二天堂av在线观看| 国产免费男女视频| 在线国产一区二区在线| 在线国产一区二区在线| 内射极品少妇av片p| 国产白丝娇喘喷水9色精品| www.av在线官网国产| 在线观看免费视频日本深夜| 天堂中文最新版在线下载 | 成人美女网站在线观看视频| 国产日韩欧美在线精品| 国产大屁股一区二区在线视频| 高清毛片免费观看视频网站| 午夜老司机福利剧场| 在线观看午夜福利视频| 免费观看的影片在线观看| 国产片特级美女逼逼视频| 亚洲av免费高清在线观看| 99热6这里只有精品| 校园春色视频在线观看| 亚洲天堂国产精品一区在线| 国产综合懂色| 精品久久久噜噜| 日韩成人伦理影院| 婷婷亚洲欧美| 三级国产精品欧美在线观看| 日本免费一区二区三区高清不卡| 天天躁夜夜躁狠狠久久av| 国产极品天堂在线| 少妇被粗大猛烈的视频| 久久婷婷人人爽人人干人人爱| 国产伦精品一区二区三区视频9| 国产一区亚洲一区在线观看| 少妇的逼好多水| 久久综合国产亚洲精品| 国产高清有码在线观看视频| 国产精华一区二区三区| 少妇人妻精品综合一区二区 | 中文字幕久久专区| 日日摸夜夜添夜夜添av毛片| 日本熟妇午夜| 夜夜看夜夜爽夜夜摸| 欧美精品一区二区大全| 精品人妻偷拍中文字幕| 婷婷色综合大香蕉| 国产精品免费一区二区三区在线| 99在线人妻在线中文字幕| 一级毛片久久久久久久久女| 欧美在线一区亚洲| 最后的刺客免费高清国语| 51国产日韩欧美| 中文在线观看免费www的网站| 国产成人a区在线观看| 亚洲在线观看片| 亚洲成av人片在线播放无| 国产黄片美女视频| 精品人妻偷拍中文字幕| 国产私拍福利视频在线观看| 美女黄网站色视频| 一级二级三级毛片免费看| 欧美xxxx黑人xx丫x性爽| 非洲黑人性xxxx精品又粗又长| 久99久视频精品免费| 亚洲av电影不卡..在线观看| 简卡轻食公司| 99在线视频只有这里精品首页| 国产单亲对白刺激| 亚洲在线自拍视频| 波多野结衣巨乳人妻| 三级经典国产精品| 日韩精品青青久久久久久| 婷婷色综合大香蕉| 精品久久国产蜜桃| 国产精品久久久久久av不卡| 男插女下体视频免费在线播放| 乱系列少妇在线播放| 别揉我奶头 嗯啊视频| 久久这里有精品视频免费| 人妻制服诱惑在线中文字幕| 亚洲电影在线观看av| 天天躁夜夜躁狠狠久久av| 国产高清激情床上av| 亚洲av免费在线观看| 欧美性猛交╳xxx乱大交人| 精品久久久久久久末码| 日本黄色片子视频| 国产 一区 欧美 日韩| 成人高潮视频无遮挡免费网站| av.在线天堂| 亚洲欧美精品专区久久| 亚洲av成人av| 卡戴珊不雅视频在线播放| 国产中年淑女户外野战色| 色尼玛亚洲综合影院| 精品久久久久久久人妻蜜臀av| 国产精品1区2区在线观看.| 国内少妇人妻偷人精品xxx网站| 亚洲精品乱码久久久久久按摩| 婷婷六月久久综合丁香| 国产精品日韩av在线免费观看| 午夜福利视频1000在线观看| 国产单亲对白刺激| 免费一级毛片在线播放高清视频| 国产精品1区2区在线观看.| 亚洲真实伦在线观看| 日韩三级伦理在线观看| 国产一区二区亚洲精品在线观看| 丝袜美腿在线中文| 久久精品夜色国产| 免费黄网站久久成人精品| 国产午夜精品一二区理论片| 亚洲欧美日韩东京热| 男女啪啪激烈高潮av片| 国产私拍福利视频在线观看| 日日干狠狠操夜夜爽| 99久久精品国产国产毛片| 麻豆av噜噜一区二区三区| 亚洲三级黄色毛片| 久久久久久久久久久免费av| 国产精品一区二区三区四区免费观看| 2021天堂中文幕一二区在线观| 看非洲黑人一级黄片| 中文欧美无线码| 国产午夜精品久久久久久一区二区三区| 国产乱人视频| 少妇人妻一区二区三区视频| 国产一区二区在线av高清观看| 性色avwww在线观看| a级毛色黄片| 美女被艹到高潮喷水动态| 日本成人三级电影网站| 成人三级黄色视频| 一级毛片电影观看 | 大型黄色视频在线免费观看| 亚洲美女视频黄频| 99九九线精品视频在线观看视频| 高清毛片免费观看视频网站| 精品日产1卡2卡| 中文资源天堂在线| av在线蜜桃| 国产成人91sexporn| 一个人免费在线观看电影| 亚洲成人久久爱视频| 国产 一区 欧美 日韩| 色综合色国产| 亚洲精品乱码久久久久久按摩| 97超视频在线观看视频| 狠狠狠狠99中文字幕| 麻豆精品久久久久久蜜桃| 99热这里只有是精品50| 国产精品一区二区三区四区免费观看| 国产精品久久久久久亚洲av鲁大| 99久久中文字幕三级久久日本| 丝袜美腿在线中文| 如何舔出高潮| kizo精华| 国产极品天堂在线| 91久久精品国产一区二区三区| 国产老妇伦熟女老妇高清| 久久久a久久爽久久v久久| 婷婷精品国产亚洲av| 亚洲真实伦在线观看| a级毛片免费高清观看在线播放| 麻豆国产av国片精品| 国产精品野战在线观看| 高清午夜精品一区二区三区 | 精品日产1卡2卡| 麻豆成人av视频| 国产精品久久久久久亚洲av鲁大| 亚州av有码| 毛片一级片免费看久久久久| 非洲黑人性xxxx精品又粗又长| 国产日韩欧美在线精品| 亚洲经典国产精华液单| 国产高清激情床上av| а√天堂www在线а√下载| 欧美3d第一页| 在线播放无遮挡| 欧美日本亚洲视频在线播放| 最近视频中文字幕2019在线8| 欧美激情国产日韩精品一区| av在线蜜桃| 国产精品国产三级国产av玫瑰| 天天躁日日操中文字幕| 热99re8久久精品国产| 午夜精品一区二区三区免费看| 欧美日韩一区二区视频在线观看视频在线 | 噜噜噜噜噜久久久久久91| 中文字幕精品亚洲无线码一区| 国产精品野战在线观看| 边亲边吃奶的免费视频| 五月玫瑰六月丁香| 午夜免费激情av| av黄色大香蕉| 日本成人三级电影网站| 国产成人91sexporn| 好男人视频免费观看在线| 老师上课跳d突然被开到最大视频| 婷婷精品国产亚洲av| 看非洲黑人一级黄片| 99国产精品一区二区蜜桃av| 亚洲精品乱码久久久v下载方式| 色综合色国产| www日本黄色视频网| 国模一区二区三区四区视频| 日本三级黄在线观看| 国产黄片美女视频| 日韩中字成人| 好男人视频免费观看在线| 亚洲av电影不卡..在线观看| 国内精品一区二区在线观看| 少妇被粗大猛烈的视频| 卡戴珊不雅视频在线播放| 欧美成人免费av一区二区三区| kizo精华| 欧美日本亚洲视频在线播放| 观看美女的网站| 国产精品久久久久久久久免| 不卡视频在线观看欧美| 蜜臀久久99精品久久宅男| 激情 狠狠 欧美| 成人二区视频| 国产成人freesex在线| 黄片无遮挡物在线观看| 免费观看人在逋| 晚上一个人看的免费电影| 成年版毛片免费区| 国产精品免费一区二区三区在线| 午夜老司机福利剧场| 99热这里只有是精品50| 日韩欧美 国产精品| 九色成人免费人妻av| 国产美女午夜福利| kizo精华| 日本色播在线视频| 日韩av在线大香蕉| 高清毛片免费看| 亚洲人成网站在线播放欧美日韩| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 能在线免费观看的黄片| 青春草国产在线视频 | 91久久精品电影网| 亚洲无线在线观看| 少妇熟女欧美另类| 亚洲无线观看免费| av国产免费在线观看| 嫩草影院入口| 又黄又爽又刺激的免费视频.| 日韩人妻高清精品专区| 一边亲一边摸免费视频| 美女 人体艺术 gogo| 欧美不卡视频在线免费观看| 成人三级黄色视频| 两个人视频免费观看高清| 级片在线观看| 寂寞人妻少妇视频99o| 舔av片在线| 又爽又黄无遮挡网站| 欧美不卡视频在线免费观看| 精品国产三级普通话版| 一本久久中文字幕| 亚洲丝袜综合中文字幕| 国产精品福利在线免费观看| 日本爱情动作片www.在线观看| 国产精品久久电影中文字幕| 欧美日韩综合久久久久久| 99热只有精品国产| 最近的中文字幕免费完整| 伦精品一区二区三区| 亚洲综合色惰| 国产精品久久久久久精品电影小说 | 天天躁夜夜躁狠狠久久av| 久久久精品欧美日韩精品| 欧美极品一区二区三区四区| 一级黄色大片毛片| 亚洲精品粉嫩美女一区| 国产一级毛片在线| h日本视频在线播放| 国产高清视频在线观看网站| 久久久久性生活片| 夫妻性生交免费视频一级片| 亚洲av免费在线观看| 1024手机看黄色片| 午夜免费男女啪啪视频观看| www日本黄色视频网| 一本久久精品| 久久亚洲国产成人精品v| 男人舔女人下体高潮全视频| 春色校园在线视频观看| 国产精品野战在线观看| www日本黄色视频网| 午夜精品在线福利| 国产午夜精品论理片| 高清日韩中文字幕在线| 99在线视频只有这里精品首页| 尾随美女入室| 丰满的人妻完整版| 国产人妻一区二区三区在| 99久久九九国产精品国产免费| 亚洲内射少妇av| 国产v大片淫在线免费观看| 成年女人看的毛片在线观看| 偷拍熟女少妇极品色| 国产综合懂色| 免费看a级黄色片| 亚洲精品影视一区二区三区av| 成人欧美大片| 日本一二三区视频观看| 欧美变态另类bdsm刘玥| eeuss影院久久| 午夜爱爱视频在线播放| 少妇高潮的动态图| 久久人人爽人人爽人人片va| 性欧美人与动物交配| 九九在线视频观看精品| 亚洲内射少妇av| 成熟少妇高潮喷水视频| 欧美精品一区二区大全| 69人妻影院| 99精品在免费线老司机午夜| eeuss影院久久| 久久九九热精品免费| 欧美日本视频| 91av网一区二区| 免费大片18禁| 高清毛片免费观看视频网站| 久久韩国三级中文字幕| 色播亚洲综合网| 最近手机中文字幕大全| 日韩高清综合在线| 国产成人影院久久av| 精品久久久久久久久久久久久| 成人毛片60女人毛片免费| 一夜夜www| 国产成人精品一,二区 | 久久精品国产亚洲网站| 欧美区成人在线视频| 欧美性猛交╳xxx乱大交人| 丝袜喷水一区| 亚洲国产欧美人成| 久久亚洲精品不卡| 国产精品久久久久久精品电影小说 | 亚洲美女视频黄频| 国产黄色视频一区二区在线观看 | 人妻制服诱惑在线中文字幕| 国产精品一区二区三区四区久久| 免费看av在线观看网站| 成年女人看的毛片在线观看| 一级毛片久久久久久久久女| 免费观看精品视频网站| 亚洲精品成人久久久久久| 给我免费播放毛片高清在线观看| 国内精品一区二区在线观看| 身体一侧抽搐| 一级二级三级毛片免费看| 久久精品久久久久久久性| 中国美白少妇内射xxxbb| 亚洲一级一片aⅴ在线观看| 中出人妻视频一区二区| 亚洲成人久久性| 好男人视频免费观看在线| 18+在线观看网站| 97超视频在线观看视频| 如何舔出高潮| 只有这里有精品99| 在线免费十八禁| 成人亚洲欧美一区二区av| 一级毛片久久久久久久久女| 免费人成在线观看视频色| av免费在线看不卡| 女人被狂操c到高潮| 亚洲成人久久性| 欧美成人一区二区免费高清观看| 熟妇人妻久久中文字幕3abv| 国国产精品蜜臀av免费| 精品久久久久久久久久免费视频| 久久精品国产自在天天线| 午夜激情欧美在线| 亚洲一区二区三区色噜噜| 丝袜美腿在线中文| 内射极品少妇av片p| 亚洲国产欧美人成| 日本一本二区三区精品| 久久久久久伊人网av| 欧美高清性xxxxhd video| 国产精品无大码| 美女 人体艺术 gogo| av在线蜜桃| 亚洲欧美清纯卡通| 国产精品1区2区在线观看.| 久久6这里有精品| 国产精品99久久久久久久久| 欧美一区二区国产精品久久精品| av在线老鸭窝| 草草在线视频免费看| 久久久国产成人精品二区| 神马国产精品三级电影在线观看| 人妻少妇偷人精品九色| 国产熟女欧美一区二区| 91在线精品国自产拍蜜月| 日韩成人伦理影院| 久久久精品94久久精品| 国产在线精品亚洲第一网站| 久久99热这里只有精品18| 晚上一个人看的免费电影| 国产三级在线视频| 在线免费观看不下载黄p国产| 两个人的视频大全免费| 99久国产av精品国产电影| 免费av不卡在线播放| 日本爱情动作片www.在线观看| 亚洲成av人片在线播放无| 国产成人福利小说| 精品国内亚洲2022精品成人| 亚洲av免费在线观看| 亚洲天堂国产精品一区在线| 久久这里只有精品中国| 午夜福利成人在线免费观看| 久久99蜜桃精品久久| 精品无人区乱码1区二区| 日韩欧美 国产精品| 日韩亚洲欧美综合| 国产极品精品免费视频能看的| 99riav亚洲国产免费| .国产精品久久| 毛片一级片免费看久久久久| 欧美激情久久久久久爽电影| 国产黄片视频在线免费观看| 亚洲性久久影院| 观看美女的网站| 久久精品国产亚洲av涩爱 | 99热精品在线国产| 亚洲国产欧洲综合997久久,| 亚洲人成网站高清观看| 亚州av有码| 国产成人a区在线观看| 国产精品,欧美在线| 级片在线观看| 草草在线视频免费看| 亚洲av电影不卡..在线观看| 国产 一区精品| 国产精品伦人一区二区| 国产 一区精品| 国产一级毛片在线| www日本黄色视频网| 午夜精品国产一区二区电影 | 欧美精品国产亚洲| 国产在线精品亚洲第一网站| 黑人高潮一二区| 国产探花在线观看一区二区| 少妇的逼好多水| 国产激情偷乱视频一区二区| 久久人人爽人人爽人人片va| 国产精品电影一区二区三区| 天堂√8在线中文| 青春草视频在线免费观看| 亚洲精品粉嫩美女一区| 国产伦精品一区二区三区四那| 久久久久久国产a免费观看| 国产午夜精品一二区理论片| 国国产精品蜜臀av免费| 日韩在线高清观看一区二区三区| 成年女人永久免费观看视频| 亚洲人成网站在线观看播放| a级毛片a级免费在线| 国产亚洲精品av在线| 一级二级三级毛片免费看| 亚洲无线观看免费| 91久久精品国产一区二区三区| 免费av观看视频| 美女xxoo啪啪120秒动态图| 国产午夜精品论理片| 亚洲av一区综合| 91精品一卡2卡3卡4卡| 天美传媒精品一区二区| 岛国毛片在线播放| 亚洲精品影视一区二区三区av| 免费av不卡在线播放| 美女cb高潮喷水在线观看| 国产女主播在线喷水免费视频网站 | 国产综合懂色| 亚洲自偷自拍三级| 免费人成在线观看视频色| 国产精品爽爽va在线观看网站| 亚洲性久久影院| 一个人观看的视频www高清免费观看| 国语自产精品视频在线第100页| 国模一区二区三区四区视频| 你懂的网址亚洲精品在线观看 | 精品一区二区三区视频在线| 日本免费一区二区三区高清不卡| 午夜激情欧美在线| 91aial.com中文字幕在线观看| 亚洲精品456在线播放app| 国产不卡一卡二| 老师上课跳d突然被开到最大视频| 国产成人精品婷婷| 日韩一区二区视频免费看| 国产一区二区三区在线臀色熟女| 日韩国内少妇激情av| 一边摸一边抽搐一进一小说|