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

    南京大學(xué)65 cm反射式望遠(yuǎn)鏡CCD相機(jī)的性能研究

    2023-02-02 04:11:20何秋會(huì)楊天寧
    天文研究與技術(shù) 2023年1期
    關(guān)鍵詞:測(cè)光天文望遠(yuǎn)鏡

    何秋會(huì),楊天寧,王 旭

    (1. 南京大學(xué)天文與空間科學(xué)學(xué)院,江蘇 南京 210046;2. 南京大學(xué)匡亞明學(xué)院,江蘇 南京 210046)

    天文學(xué)是一門觀測(cè)與理論研究緊密聯(lián)系的學(xué)科,也是現(xiàn)代科學(xué)與高新技術(shù)完美結(jié)合的典范。天文學(xué)的發(fā)展離不開天文望遠(yuǎn)鏡,現(xiàn)代天文望遠(yuǎn)鏡離不開探測(cè)器。1609年伽利略發(fā)明了天文望遠(yuǎn)鏡[1],之后的200多年,天文學(xué)家都是依靠望遠(yuǎn)鏡和肉眼進(jìn)行天文觀測(cè),并通過觀察或素描記錄他們的發(fā)現(xiàn)[2]。1839年美國(guó)化學(xué)家約翰的·雷柏首次用照相機(jī)拍攝了月亮的照片,從此照相技術(shù)開始用于天文觀測(cè),天文測(cè)光技術(shù)也逐漸建立。20世紀(jì)后,隨著電子技術(shù)的發(fā)展,人們開始把光信號(hào)轉(zhuǎn)化成電信號(hào),天文觀測(cè)探測(cè)器逐步出現(xiàn)了光電器件,如光電倍增管,天文觀測(cè)進(jìn)入了光電測(cè)光時(shí)代[3]。底片照相的優(yōu)點(diǎn)是探測(cè)面積大,可以長(zhǎng)時(shí)間曝光,缺點(diǎn)是響應(yīng)非線性。光電測(cè)光的優(yōu)點(diǎn)是高靈敏度而且線性,但探測(cè)面積小。1969年,美國(guó)貝爾實(shí)驗(yàn)室發(fā)明了CCD,它兼有底片照相和光電測(cè)光的優(yōu)點(diǎn)。1975年科學(xué)級(jí)CCD相機(jī)應(yīng)用于天文望遠(yuǎn)鏡,很快廣泛應(yīng)用于光學(xué)波段觀測(cè),從而徹底改變了天文學(xué)的觀測(cè)方法。20世紀(jì)90年代,光學(xué)天文探測(cè)器基本完全由科學(xué)級(jí)CCD探測(cè)器代替。

    科學(xué)級(jí)CCD相機(jī)在天文上的廣泛應(yīng)用,不僅因?yàn)镃CD相機(jī)可以得到數(shù)字化的圖像[4],還因?yàn)樗哂懈哽`敏、低噪聲、優(yōu)良的線性和大動(dòng)態(tài)范圍等良好性能[3]。這些性能使科學(xué)級(jí)CCD相機(jī)能探測(cè)到更暗弱的天體??茖W(xué)級(jí)CCD相機(jī)性能雖然良好,但在實(shí)際使用時(shí)會(huì)隨著環(huán)境改變或使用年限等因素變化而變化。當(dāng)購置CCD相機(jī)時(shí),生產(chǎn)廠家提供部分性能參數(shù),但按照慣例,生產(chǎn)廠家并不提供CCD相機(jī)各個(gè)像元的快門或線性等性能參數(shù)。在觀測(cè)亮天體時(shí),CCD相機(jī)短時(shí)間曝光進(jìn)行圖像采集,需考慮CCD相機(jī)快門性能對(duì)測(cè)光精度的影響,所以天文觀測(cè)工作者必須對(duì)CCD相機(jī)的性能參數(shù)進(jìn)行測(cè)定[5-6]。本文對(duì)南京大學(xué)65 cm天文望遠(yuǎn)鏡使用的探測(cè)器PIXIS 2048B CCD相機(jī)進(jìn)行了一些特性測(cè)定,分析它的快門函數(shù)和線性,以便于做精細(xì)測(cè)光時(shí)進(jìn)行必要修正。

    1 CCD的線性和快門函數(shù)

    1.1 快門函數(shù)

    天文CCD 相機(jī)系統(tǒng)要求有快門功能。由于相機(jī)快門的動(dòng)作獨(dú)立于光學(xué)系統(tǒng),所以在光學(xué)系統(tǒng)對(duì)光通量控制的基礎(chǔ)上,相機(jī)快門進(jìn)一步增強(qiáng)了對(duì)信號(hào)電荷量控制的靈活性[7]。 天文CCD相機(jī)的快門有光瞳式和簾幕式兩種類型,PIXIS 2048B CCD相機(jī)使用光瞳式。光瞳式CCD相機(jī)在曝光開始后,相機(jī)快門在CCD靶面中心最先打開,曝光結(jié)束時(shí)相機(jī)快門在CCD靶面中心最后關(guān)閉, 因此CCD靶面中心像元的實(shí)際曝光時(shí)長(zhǎng)與邊角的實(shí)際曝光時(shí)長(zhǎng)不同。其次相機(jī)快門的開關(guān)速度有限,當(dāng)使用CCD相機(jī)做短時(shí)間曝光進(jìn)行圖像采集時(shí),快門效應(yīng)很明顯。對(duì)于觀測(cè)比較亮的天體目標(biāo),圖像采集時(shí)一般只能進(jìn)行短時(shí)間曝光,因此必須對(duì)使用的CCD相機(jī)的快門進(jìn)行測(cè)定,了解其響應(yīng)大小,然后進(jìn)行必要的修正。

    測(cè)定快門修正的方法有文[8]提出的圓頂平場(chǎng)法和晨昏蒙影天光平場(chǎng)法。其中圓頂平場(chǎng)法不要求定平場(chǎng)函數(shù)而只要求定快門函數(shù),方法簡(jiǎn)單準(zhǔn)確,只需要光源穩(wěn)定即可[9],因此本文使用該方法進(jìn)行快門函數(shù)的測(cè)定。

    快門函數(shù)的修正首先需要測(cè)定快門函數(shù),快門函數(shù)的獲得類似于平場(chǎng)修正:假設(shè)CCD第i列第j行的像元響應(yīng)系數(shù)為a(i,j),穩(wěn)定的圓頂平場(chǎng)燈光照射到此像元的光強(qiáng)為L(zhǎng)(i,j),取CCD相機(jī)長(zhǎng)短不同的兩個(gè)積分時(shí)間t1和t2,測(cè)量該像元的實(shí)際輸出為I(i,j)t1和I(i,j)t2,則有

    I1≡I(i,j)t1=a(i,j)L(i,j)[t1+Δ(i,j)],

    (1)

    I2≡I(i,j)t2=a(i,j)L(i,j)[t2+Δ(i,j)],

    (2)

    其中,快門函數(shù)Δ(i,j)(對(duì)PIXIS 2048B CCD,i=1, 2, 3, ..., 2 048;j=1, 2, 3, ..., 2 048,(i,j)為CCD第i列第j行像元)為相機(jī)快門影響導(dǎo)致的延遲時(shí)間,即拍攝程序設(shè)置的積分時(shí)間t與真實(shí)的積分時(shí)間T(i,j)的差值。因?yàn)閍(i,j)和L(i,j)不變化,可消除它們的影響,(1)式和(2)式相除有

    I1/I2=[t1+Δ(i,j)]/[t2+Δ(i,j)] ,

    (3)

    整理可得快門函數(shù)

    Δ(i,j)=(I2t1-I1t2)/(I1-I2),

    (4)

    所以第i列第j行像元真實(shí)的曝光時(shí)間為

    T(i,j)=t+Δ(i,j).

    (5)

    1.2 線 性

    CCD的線性是動(dòng)態(tài)范圍內(nèi)入射光強(qiáng)變化與CCD響應(yīng)變化之間的線性程度。一般用代表CCD相機(jī)產(chǎn)生的電信號(hào)與入射光強(qiáng)之間的線性關(guān)系的線性度表征線性的好壞。在典型的視頻或電視應(yīng)用中,人們并不關(guān)心CCD相機(jī)的線性度,甚至允許高達(dá)20%的偏差[10]。然而,在成像科學(xué)研究領(lǐng)域,響應(yīng)的線性度是必須考慮的重要參數(shù)。理想情況是每一個(gè)像元可以作為一個(gè)高保真的光子測(cè)量器件,它產(chǎn)生的電信號(hào)與入射光的強(qiáng)度成精確的正比。因此,保證 CCD相機(jī)的線性度對(duì)于成像科學(xué)在定量研究中的應(yīng)用十分重要。一個(gè)高性能的科學(xué)級(jí)CCD相機(jī)的非線性度在整個(gè)動(dòng)態(tài)范圍內(nèi)甚至小于0.1%[8],同時(shí)像元之間具有良好的響應(yīng)均勻性。

    (6)

    根據(jù)不同數(shù)據(jù)的處理對(duì)精度的不同要求,可以在NL大于一定值的區(qū)域范圍,才認(rèn)為是非線性區(qū)域。

    2 觀測(cè)情況及觀測(cè)結(jié)果

    2.1 觀測(cè)情況

    南京大學(xué)天文與空間科學(xué)學(xué)院本科教學(xué)中心實(shí)驗(yàn)室在左滌江天文臺(tái)(東經(jīng)118°57′34″,北緯32°07′45″)擁有一架65 cm口徑的天文光學(xué)望遠(yuǎn)鏡,該望遠(yuǎn)鏡使用的探測(cè)器是型號(hào)為PIXIS 2048B的一款CCD相機(jī),這款CCD相機(jī)主要用于南京大學(xué)本科實(shí)驗(yàn)課程教學(xué)的實(shí)驗(yàn)觀測(cè)。實(shí)驗(yàn)觀測(cè)項(xiàng)目有大氣消光、像元比例尺、變星較差測(cè)光等。實(shí)驗(yàn)觀測(cè)項(xiàng)目的主要目的是訓(xùn)練天文本科學(xué)生的觀測(cè)能力以及數(shù)據(jù)處理能力,最初對(duì)實(shí)驗(yàn)結(jié)果的精度要求不高,CCD相機(jī)快門和線性對(duì)實(shí)驗(yàn)的影響歸算在誤差分析中。但是,隨著實(shí)驗(yàn)教學(xué)質(zhì)量提升計(jì)劃的實(shí)施,研究生參與觀測(cè)項(xiàng)目的增多以及和外單位合作的聯(lián)合觀測(cè)項(xiàng)目的開展,實(shí)驗(yàn)觀測(cè)項(xiàng)目對(duì)觀測(cè)精度的要求逐步提高。為了在觀測(cè)時(shí)更好地了解CCD相機(jī)性能,確定在精準(zhǔn)測(cè)光時(shí)合適的曝光時(shí)間,必須對(duì)CCD相機(jī)的參數(shù)有詳細(xì)的了解。

    2021年5~6月,我們通過使用望遠(yuǎn)鏡室內(nèi)的平場(chǎng)燈先對(duì)快門函數(shù)進(jìn)行測(cè)定。測(cè)定時(shí),選擇CCD最靈敏的I濾光片,曝光時(shí)間為一短一長(zhǎng),短的為1 s,長(zhǎng)的為10 s。CCD采集的圖像最大值為65 535 adu,1 s曝光的圖像減去本底后的值為5 000 adu左右,10 s曝光的圖像減去本底后的值為50 000 adu左右。測(cè)定快門函數(shù)的拍攝順序?yàn)? s曝光連續(xù)拍攝5幅,10 s曝光連續(xù)拍攝2幅,1 s曝光再連續(xù)拍攝5幅。在快門函數(shù)測(cè)定拍攝前后各拍攝10幅本底圖像,用作處理曝光拍攝的圖像和檢查CCD狀態(tài)是否穩(wěn)定。在整個(gè)測(cè)定過程中,保持CCD相機(jī)制冷后溫度和平場(chǎng)燈照度穩(wěn)定。

    同時(shí),我們對(duì)CCD相機(jī)的線性進(jìn)行了測(cè)定。使用經(jīng)過光學(xué)設(shè)計(jì)的出射光較為均勻的平場(chǎng)燈做平場(chǎng)拍攝。為了等待光源趨于穩(wěn)定,在開燈半小時(shí)后才開始拍攝,在剛開始拍攝時(shí),每隔幾分鐘拍1幅平場(chǎng)圖像,計(jì)算圖像的平均值,數(shù)值穩(wěn)定后正式開始采集。在進(jìn)行平場(chǎng)圖像采集時(shí),CCD相機(jī)的曝光序列分為圖像采集序列和監(jiān)測(cè)采集序列,如表1,監(jiān)測(cè)采集序列穿插在圖像采集序列中。根據(jù)測(cè)得的快門函數(shù),本文作者觀測(cè)時(shí)選擇的最低曝光時(shí)間從5 s開始,直到180 s飽和。CCD相機(jī)曝光序列的曝光時(shí)間為10 s, 5 s, 10 s, 15 s, 10 s, 20 s, 10 s, 25 s, 10 s, 30 s, 10 s, 35 s, 10 s, 40 s, 10 s, 45 s, ...., 10 s, 180 s,其中穿插的10 s曝光為監(jiān)測(cè)曝光,所以測(cè)定線性的圖像采集序列的時(shí)間為5 s, 10 s, 15 s, 20 s, 25 s, 30 s, 35 s, 40 s, 45 s, ...., 180 s。穿插10 s曝光是因?yàn)槿绻庠戳炼炔环€(wěn)定,光源亮度會(huì)慢慢漂移,為了減小光源不穩(wěn)定帶來的影響,中間不斷加拍圖像以監(jiān)測(cè)亮度變化。

    表1 CCD圖像采集的曝光時(shí)間順序Table 1 Exposure time sequence of CCD imaging

    通過對(duì)觀測(cè)數(shù)據(jù)的處理,我們得到快門函數(shù)、線性、增益和讀出噪聲。

    2.2 觀測(cè)結(jié)果

    2.2.1 快門函數(shù)

    我們使用通用的IRAF軟件[11]對(duì)觀測(cè)得到的數(shù)據(jù)進(jìn)行處理,得到相機(jī)的快門函數(shù)SHUTTER.fits(圖1)和快門函數(shù)的等光圖如圖2(其中不同的閉合曲線代表不同的快門修正值)。根據(jù)處理得到的CCD相機(jī)快門函數(shù)可知,快門中心處極大值約0.04 s,邊緣極小值約0.006 s,CCD相機(jī)的中心與邊緣曝光時(shí)間相差約0.04 s。因此,根據(jù)測(cè)量的快門函數(shù)結(jié)果,若需要整個(gè)圖像的測(cè)光精度差異小于1%,設(shè)置的曝光時(shí)間需大于4 s;若需要整個(gè)圖像的測(cè)光精度差異小于0.1%,設(shè)置的曝光時(shí)間則需大于40 s。

    圖1快門函數(shù)示意圖和圖2快門函數(shù)等光圖都只是示意圖,它們可以用來半定量地判斷CCD上有關(guān)區(qū)域的快門影響大小。為了定量修正快門效應(yīng),還需要做圖像處理[12]。在實(shí)際觀測(cè)的修正中, ‘SHUTTER’ 表示快門函數(shù),對(duì)于觀測(cè)目標(biāo)曝光t時(shí)間得到的星像圖像,先減去CCD本底 ‘BIAS’,再做平場(chǎng)處理后,得到的星像圖像表示為 ‘IMAGE’,使用任何圖像處理軟件例如IRAF,進(jìn)行圖像運(yùn)算處理,

    IMAGE0=[IMAGE/(SHUTTER+t)]t

    (7)

    這樣就得到經(jīng)過修正快門效應(yīng)的所需圖像 ‘IMAGE0’。

    快門函數(shù)測(cè)定后,可以使用一段時(shí)間,但長(zhǎng)時(shí)間多次反復(fù)開關(guān)快門,可能引起快門的機(jī)件老化,需要定期重新檢測(cè)。

    圖1 相機(jī)快門函數(shù)示意圖Fig.1 Diagram of the shutter function

    圖2 相機(jī)快門函數(shù)的等光圖Fig.2 Contour of the shutter function

    2.2.2 線 性

    通過處理觀測(cè)數(shù)據(jù)得到CCD的線性情況,如圖3和圖4。圖3是CCD的線性實(shí)測(cè)值和擬合曲線,圖4是測(cè)量值偏離線性的百分比。從處理的數(shù)據(jù)可以清晰地得到CCD的線性。實(shí)驗(yàn)使用的PIXIS 2048B CCD相機(jī)的動(dòng)態(tài)范圍是0到65 535 adu,測(cè)試結(jié)果顯示一直到61 900 adu,CCD的非線性都小于1%。需要說明的是,圖3和圖4是對(duì)面光源(圓頂平場(chǎng))的線性測(cè)定,雖然直到61 900 adu還保持線性,但是對(duì)于點(diǎn)光源(恒星)的拍攝,情況有所不同。在2021年的1月到5月,我們多次對(duì)巨蟹座中的疏散星團(tuán)M67進(jìn)行觀測(cè),以M67中的一顆亮星HD75700(08:51:49 +11:53:38)為例,曝光5 s得到圖像,經(jīng)過IRAF處理得這顆亮星的面亮度分布如圖5,星的徑向輪廓如圖6。由圖5和圖6可見,在約38 000 adu時(shí),星像已經(jīng)飽和,所以曝光后圖像數(shù)值大于38 000 adu的圖像部分是不能使用的。

    3 討 論

    過早飽和(溢出)是本CCD相機(jī)的弱點(diǎn),這限制了點(diǎn)光源測(cè)光的動(dòng)態(tài)范圍,以及進(jìn)一步使用點(diǎn)擴(kuò)散函數(shù)(Point Spread Function, PSF)分析。但如果控制曝光時(shí)間,避免亮星過度曝光,仍然可以使用本CCD相機(jī)對(duì)點(diǎn)光源做高精度測(cè)光[13]。

    圖3 CCD的線性實(shí)測(cè)值和擬合曲線Fig.3 Experimental value and fit curve of the CCD linearity

    圖4 測(cè)量值偏離線性的百分比Fig.4 Linearity deviation measument of the CCD

    圖5 拍攝的HD75700的面亮度分布圖Fig.5 Surface brightness profile of the observed HD75700

    圖6 拍攝的HD75700的徑向輪廓圖Fig.6 Radial profile of the observed HD75700

    本文在使用CCD相機(jī)進(jìn)行圖像采集時(shí),使用基特峰天文臺(tái)的MaxImDL控制程序。為探究使用MaxImDL程序控制CCD相機(jī)進(jìn)行圖像采集出現(xiàn)過早溢出的原因,我們聯(lián)系了云南天文臺(tái)1 m望遠(yuǎn)鏡的負(fù)責(zé)人和國(guó)家天文臺(tái)2.16 m望遠(yuǎn)鏡的觀測(cè)專家進(jìn)行了細(xì)致討論。云南天文臺(tái)1 m望遠(yuǎn)鏡使用不同型號(hào)的CCD相機(jī)溢出情況不同,PI公司大小為1 024 × 1 024的CCD相機(jī)沒有過早溢出,但是,后來購進(jìn)的ANDOR公司的CCD相機(jī)卻有過早溢出,經(jīng)過各種嘗試,也沒有從根本上解決過早溢出問題。國(guó)家天文臺(tái)2.16 m望遠(yuǎn)鏡使用MaxImDL控制程序,當(dāng)有過早溢出時(shí),改變bin參數(shù),過早溢出會(huì)得到改善。本文也曾嘗試改變參數(shù),但是情況沒有改變。在文獻(xiàn)中得知當(dāng)Massey等在基特峰天文臺(tái)處理老的RCA CCD相機(jī)資料時(shí),也遇到類似情況,他們把這歸之于 “細(xì)微的電荷傳遞問題”[13]。過早溢出究竟是電荷傳遞問題還是CCD芯片問題或控制程序問題,還需要進(jìn)一步分析探究。

    4 總 結(jié)

    對(duì)于在實(shí)驗(yàn)中使用的CCD等儀器設(shè)備,盡管廠家提供了參考的性能參數(shù),但為了更細(xì)致深入地了解該儀器設(shè)備,仍需要使用者測(cè)定相關(guān)參數(shù)。本文的研究測(cè)定方法對(duì)于新購CCD的使用者具有參考性和指導(dǎo)性。本文對(duì)南京大學(xué)65 cm光學(xué)望遠(yuǎn)鏡的后端探測(cè)器PIXIS 2048B CCD相機(jī)性能進(jìn)行了研究和實(shí)際測(cè)定,得到了CCD相機(jī)的快門函數(shù)、線性、增益和噪聲。測(cè)得的快門函數(shù)已存儲(chǔ)于觀測(cè)用的計(jì)算機(jī)中,以后的觀測(cè)中,必要時(shí)把觀測(cè)數(shù)據(jù)和快門修正函數(shù)一同給與數(shù)據(jù)處理者,方便進(jìn)行修正。這款CCD相機(jī)的讀出值從0到61 900 adu,非線性都小于1%,CCD相機(jī)的増益約1.02e-/adu,讀出噪聲約13e-。在后續(xù)的實(shí)驗(yàn)教學(xué)觀測(cè)以及研究項(xiàng)目觀測(cè)中,這些測(cè)得的性能參數(shù)將指導(dǎo)研究者更好地選擇合適的觀測(cè)目標(biāo),設(shè)置觀測(cè)參數(shù),以及后期處理數(shù)據(jù)的校正,有利于更好地發(fā)揮設(shè)備性能。

    致謝:非常感謝上海天文臺(tái)姚保安老師給予的指導(dǎo)和幫助。

    猜你喜歡
    測(cè)光天文望遠(yuǎn)鏡
    天文動(dòng)手做/中國(guó)古天文系列 三星在天、在隅、在戶
    軍事文摘(2024年4期)2024-01-09 09:08:34
    神奇的“千里眼”——望遠(yuǎn)鏡
    天文篇
    神奇的千里眼——望遠(yuǎn)鏡
    打造最大望遠(yuǎn)鏡
    本月話題機(jī)內(nèi)測(cè)光系統(tǒng)
    影像視覺(2017年10期)2018-01-31 23:20:26
    超級(jí)望遠(yuǎn)鏡頭
    不同種手機(jī)測(cè)光類軟件在林業(yè)中的試用探析
    一群“玩”天文的癡迷少年
    相機(jī)探秘測(cè)光指南
    日本黄色视频三级网站网址| 亚洲va日本ⅴa欧美va伊人久久| 少妇 在线观看| 亚洲国产毛片av蜜桃av| 少妇裸体淫交视频免费看高清 | 高清欧美精品videossex| 亚洲全国av大片| 一进一出好大好爽视频| 亚洲一区二区三区欧美精品| 国产精品永久免费网站| 九色亚洲精品在线播放| 国产又爽黄色视频| 免费人成视频x8x8入口观看| 老司机福利观看| 91九色精品人成在线观看| 麻豆久久精品国产亚洲av | 欧美最黄视频在线播放免费 | 亚洲欧美一区二区三区黑人| 久久草成人影院| 韩国av一区二区三区四区| 午夜精品国产一区二区电影| 久久香蕉国产精品| x7x7x7水蜜桃| 国产伦一二天堂av在线观看| 国产成人系列免费观看| 亚洲av熟女| 久久精品国产99精品国产亚洲性色 | 麻豆av在线久日| 国产亚洲欧美精品永久| 水蜜桃什么品种好| 精品国产一区二区三区四区第35| 亚洲国产精品一区二区三区在线| 高清av免费在线| 欧美性长视频在线观看| 无遮挡黄片免费观看| 国产av精品麻豆| 三级毛片av免费| 国产精品野战在线观看 | 国产97色在线日韩免费| 久久精品亚洲熟妇少妇任你| 日韩欧美一区视频在线观看| 亚洲精品中文字幕一二三四区| 97超级碰碰碰精品色视频在线观看| 丁香六月欧美| 亚洲专区中文字幕在线| 成在线人永久免费视频| 久久久国产成人精品二区 | 色精品久久人妻99蜜桃| 两个人看的免费小视频| 一进一出抽搐gif免费好疼 | 少妇 在线观看| 男女做爰动态图高潮gif福利片 | 亚洲精华国产精华精| 一边摸一边做爽爽视频免费| 欧美乱色亚洲激情| 女人被躁到高潮嗷嗷叫费观| 国产精品香港三级国产av潘金莲| 女性生殖器流出的白浆| 黄片大片在线免费观看| 国产人伦9x9x在线观看| 97人妻天天添夜夜摸| 国产精品影院久久| 国产精品香港三级国产av潘金莲| 长腿黑丝高跟| 久久中文字幕人妻熟女| 色综合欧美亚洲国产小说| 一本大道久久a久久精品| 色在线成人网| 在线看a的网站| 精品熟女少妇八av免费久了| 欧美人与性动交α欧美精品济南到| 欧美日本中文国产一区发布| 午夜福利一区二区在线看| 18禁国产床啪视频网站| 五月开心婷婷网| 男人的好看免费观看在线视频 | 十分钟在线观看高清视频www| 亚洲情色 制服丝袜| 美国免费a级毛片| 日韩免费av在线播放| 欧美黑人精品巨大| 免费高清在线观看日韩| 一区二区三区精品91| 国产区一区二久久| 正在播放国产对白刺激| 高清欧美精品videossex| 欧美精品亚洲一区二区| 嫁个100分男人电影在线观看| 99精品久久久久人妻精品| 国产精品一区二区在线不卡| 777久久人妻少妇嫩草av网站| 国产精品久久久av美女十八| 黄色视频,在线免费观看| 成人手机av| tocl精华| 亚洲一区高清亚洲精品| 亚洲欧洲精品一区二区精品久久久| 制服人妻中文乱码| 精品国内亚洲2022精品成人| 日韩精品中文字幕看吧| www.自偷自拍.com| av视频免费观看在线观看| 极品教师在线免费播放| 99香蕉大伊视频| 亚洲精品国产色婷婷电影| 国产免费男女视频| 成人av一区二区三区在线看| 丰满的人妻完整版| 黄色视频,在线免费观看| 十八禁人妻一区二区| 久久久久久久精品吃奶| av超薄肉色丝袜交足视频| 国产色视频综合| 久久久久国产一级毛片高清牌| 法律面前人人平等表现在哪些方面| 午夜福利影视在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 午夜福利一区二区在线看| 亚洲视频免费观看视频| 免费女性裸体啪啪无遮挡网站| 国产麻豆69| 自拍欧美九色日韩亚洲蝌蚪91| 午夜免费激情av| 老司机午夜福利在线观看视频| svipshipincom国产片| 脱女人内裤的视频| 精品福利观看| 久久人人97超碰香蕉20202| 国产在线观看jvid| 国产深夜福利视频在线观看| 久久精品国产清高在天天线| 亚洲av第一区精品v没综合| 日韩有码中文字幕| 亚洲精品国产色婷婷电影| 日韩国内少妇激情av| 大型黄色视频在线免费观看| 欧美黑人欧美精品刺激| 亚洲国产欧美一区二区综合| 成人国语在线视频| 99精品在免费线老司机午夜| 亚洲五月色婷婷综合| 波多野结衣av一区二区av| 国产三级黄色录像| 99riav亚洲国产免费| 高清av免费在线| 变态另类成人亚洲欧美熟女 | 国产一区二区激情短视频| 亚洲免费av在线视频| 色老头精品视频在线观看| 欧美老熟妇乱子伦牲交| 精品无人区乱码1区二区| 国产成+人综合+亚洲专区| 久久精品aⅴ一区二区三区四区| 超色免费av| 免费一级毛片在线播放高清视频 | 51午夜福利影视在线观看| 国产精品永久免费网站| av天堂在线播放| 国产精品久久久久久人妻精品电影| 欧美大码av| 精品电影一区二区在线| 电影成人av| 日韩大尺度精品在线看网址 | av在线天堂中文字幕 | 精品电影一区二区在线| 色精品久久人妻99蜜桃| 午夜91福利影院| 国产高清激情床上av| 香蕉国产在线看| 91九色精品人成在线观看| 日韩欧美一区视频在线观看| 精品国产一区二区久久| 又紧又爽又黄一区二区| av在线天堂中文字幕 | 最好的美女福利视频网| 国产黄色免费在线视频| 色播在线永久视频| 国产一卡二卡三卡精品| 夜夜夜夜夜久久久久| 老司机亚洲免费影院| 啦啦啦免费观看视频1| 99精国产麻豆久久婷婷| www.www免费av| 别揉我奶头~嗯~啊~动态视频| 女人被狂操c到高潮| 高清在线国产一区| 亚洲一区二区三区不卡视频| 欧美日韩国产mv在线观看视频| 美女国产高潮福利片在线看| 亚洲 国产 在线| 老汉色av国产亚洲站长工具| av福利片在线| 久久青草综合色| 身体一侧抽搐| 99国产精品免费福利视频| 国产成人影院久久av| 高清av免费在线| 国产成人av激情在线播放| 香蕉久久夜色| 色哟哟哟哟哟哟| 一区福利在线观看| 日本精品一区二区三区蜜桃| 久久久精品欧美日韩精品| 久久久久九九精品影院| 日本vs欧美在线观看视频| 国产成人av激情在线播放| 一级毛片精品| 一边摸一边做爽爽视频免费| 香蕉久久夜色| 日本免费一区二区三区高清不卡 | 水蜜桃什么品种好| 韩国精品一区二区三区| 国产一区二区三区综合在线观看| 欧美成狂野欧美在线观看| av天堂久久9| 黄色成人免费大全| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品sss在线观看 | 国产精品久久久久久人妻精品电影| 国产成人欧美在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产免费男女视频| 国产精品久久久av美女十八| 亚洲成av片中文字幕在线观看| 91国产中文字幕| 精品一区二区三区av网在线观看| 国产精品1区2区在线观看.| 婷婷六月久久综合丁香| 夜夜爽天天搞| 免费久久久久久久精品成人欧美视频| 国产精品综合久久久久久久免费 | 欧美日韩亚洲综合一区二区三区_| 伦理电影免费视频| 国产野战对白在线观看| 久久香蕉激情| 一区二区三区激情视频| 老鸭窝网址在线观看| 看黄色毛片网站| 1024香蕉在线观看| 老司机福利观看| 亚洲专区中文字幕在线| 成年人黄色毛片网站| 老司机亚洲免费影院| 日韩三级视频一区二区三区| 国产精品永久免费网站| 精品一区二区三卡| 精品一区二区三区四区五区乱码| 亚洲成人久久性| 电影成人av| 亚洲av成人一区二区三| 在线天堂中文资源库| 精品久久久久久成人av| 久久人人精品亚洲av| 欧美日韩视频精品一区| 亚洲熟女毛片儿| 精品国产乱码久久久久久男人| 50天的宝宝边吃奶边哭怎么回事| 丰满的人妻完整版| 国产亚洲欧美精品永久| 亚洲精品在线观看二区| 午夜精品久久久久久毛片777| 成年版毛片免费区| 天堂中文最新版在线下载| 日日夜夜操网爽| 久久精品aⅴ一区二区三区四区| 性少妇av在线| 色哟哟哟哟哟哟| 日韩人妻精品一区2区三区| 日韩欧美在线二视频| 波多野结衣高清无吗| 成人手机av| 欧美日韩瑟瑟在线播放| 俄罗斯特黄特色一大片| 国产成人精品久久二区二区免费| 午夜福利,免费看| 国产亚洲欧美精品永久| 美女大奶头视频| 久99久视频精品免费| 一级a爱视频在线免费观看| 亚洲九九香蕉| 91在线观看av| 91成年电影在线观看| 精品一区二区三区av网在线观看| 18禁国产床啪视频网站| 久久久久精品国产欧美久久久| 91成人精品电影| 精品熟女少妇八av免费久了| 国产无遮挡羞羞视频在线观看| 99国产精品一区二区蜜桃av| 日韩国内少妇激情av| 午夜福利免费观看在线| 国产在线精品亚洲第一网站| 1024香蕉在线观看| www.精华液| 色在线成人网| 亚洲一区高清亚洲精品| 国产亚洲精品第一综合不卡| 亚洲成a人片在线一区二区| 老熟妇仑乱视频hdxx| 男人的好看免费观看在线视频 | tocl精华| 欧美日韩黄片免| 亚洲国产欧美一区二区综合| 欧美日韩精品网址| 亚洲一码二码三码区别大吗| 欧美国产精品va在线观看不卡| 午夜福利免费观看在线| 热re99久久国产66热| 中文字幕人妻丝袜一区二区| 国产日韩一区二区三区精品不卡| 久久精品aⅴ一区二区三区四区| 国产成人欧美在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 首页视频小说图片口味搜索| 日韩成人在线观看一区二区三区| 99久久精品国产亚洲精品| 99久久久亚洲精品蜜臀av| 99riav亚洲国产免费| 精品福利观看| 国产精品 国内视频| 日韩国内少妇激情av| 日韩视频一区二区在线观看| 亚洲午夜理论影院| 亚洲性夜色夜夜综合| 可以在线观看毛片的网站| 黄片小视频在线播放| 丰满人妻熟妇乱又伦精品不卡| 在线观看一区二区三区激情| x7x7x7水蜜桃| 乱人伦中国视频| 99精品久久久久人妻精品| 精品国产超薄肉色丝袜足j| 中国美女看黄片| 777久久人妻少妇嫩草av网站| 欧美成人性av电影在线观看| 美女国产高潮福利片在线看| 久久久久久久久久久久大奶| 在线观看免费日韩欧美大片| 女同久久另类99精品国产91| 午夜免费激情av| 一级黄色大片毛片| 十八禁人妻一区二区| 少妇 在线观看| 久久久久久久久久久久大奶| 变态另类成人亚洲欧美熟女 | 老司机午夜福利在线观看视频| 国产av在哪里看| 很黄的视频免费| 天天躁夜夜躁狠狠躁躁| 电影成人av| 91九色精品人成在线观看| 久久中文看片网| a在线观看视频网站| 麻豆一二三区av精品| 国产男靠女视频免费网站| 高清黄色对白视频在线免费看| 在线av久久热| 女人爽到高潮嗷嗷叫在线视频| 国产精品国产av在线观看| 亚洲人成电影免费在线| 91在线观看av| 1024香蕉在线观看| 97超级碰碰碰精品色视频在线观看| 一级a爱视频在线免费观看| 久久精品人人爽人人爽视色| 男女高潮啪啪啪动态图| 亚洲欧美日韩无卡精品| 欧美成狂野欧美在线观看| 国产精品久久久久成人av| 99精品欧美一区二区三区四区| 国产单亲对白刺激| 国产成人欧美| 国产精品久久久久久人妻精品电影| 亚洲精品av麻豆狂野| 婷婷六月久久综合丁香| 无人区码免费观看不卡| 在线观看免费视频日本深夜| 黄色视频,在线免费观看| 久久热在线av| 免费少妇av软件| 99热国产这里只有精品6| 黄网站色视频无遮挡免费观看| 男女之事视频高清在线观看| 欧美激情久久久久久爽电影 | 国产av又大| 在线观看66精品国产| 欧美激情极品国产一区二区三区| 中出人妻视频一区二区| 亚洲 欧美 日韩 在线 免费| 国产av精品麻豆| 日日夜夜操网爽| 欧美激情 高清一区二区三区| 精品久久蜜臀av无| 欧美+亚洲+日韩+国产| 久久亚洲真实| 精品国产国语对白av| 国产黄a三级三级三级人| 波多野结衣一区麻豆| 国产精品秋霞免费鲁丝片| 女同久久另类99精品国产91| 高清黄色对白视频在线免费看| 亚洲精品久久午夜乱码| 欧美黑人精品巨大| e午夜精品久久久久久久| 午夜两性在线视频| 满18在线观看网站| 免费在线观看日本一区| 欧美在线黄色| 一级作爱视频免费观看| 成年版毛片免费区| 午夜福利在线免费观看网站| 黄片播放在线免费| 国产不卡一卡二| 国产亚洲欧美98| 三上悠亚av全集在线观看| 1024视频免费在线观看| 亚洲欧美一区二区三区久久| 美女高潮喷水抽搐中文字幕| 性欧美人与动物交配| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看完整版高清| 亚洲成国产人片在线观看| 超碰97精品在线观看| 51午夜福利影视在线观看| 夜夜躁狠狠躁天天躁| 久久亚洲精品不卡| aaaaa片日本免费| 热re99久久精品国产66热6| 嫩草影院精品99| 国产成人av激情在线播放| 日韩 欧美 亚洲 中文字幕| 久久香蕉激情| 精品国产亚洲在线| 这个男人来自地球电影免费观看| 国产精品国产高清国产av| 俄罗斯特黄特色一大片| 久久国产亚洲av麻豆专区| av中文乱码字幕在线| 国产成人av教育| 午夜视频精品福利| 日日干狠狠操夜夜爽| 亚洲自偷自拍图片 自拍| 国产精品二区激情视频| 女人精品久久久久毛片| 日本免费一区二区三区高清不卡 | 他把我摸到了高潮在线观看| 久久精品国产99精品国产亚洲性色 | 日本黄色日本黄色录像| 热99re8久久精品国产| 成人av一区二区三区在线看| 久久精品亚洲av国产电影网| 一级毛片高清免费大全| 看片在线看免费视频| 国产成年人精品一区二区 | 自线自在国产av| 69精品国产乱码久久久| 老汉色∧v一级毛片| 69av精品久久久久久| 女人爽到高潮嗷嗷叫在线视频| 欧美日本中文国产一区发布| 亚洲色图综合在线观看| 岛国视频午夜一区免费看| 国产成人啪精品午夜网站| 村上凉子中文字幕在线| 中文字幕人妻熟女乱码| 日韩中文字幕欧美一区二区| 最近最新中文字幕大全免费视频| 精品卡一卡二卡四卡免费| 久热这里只有精品99| 国产av又大| 美国免费a级毛片| 99精品欧美一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 男女下面进入的视频免费午夜 | 韩国精品一区二区三区| 国产精品久久视频播放| 亚洲中文日韩欧美视频| 极品教师在线免费播放| 色播在线永久视频| 女性生殖器流出的白浆| 人妻丰满熟妇av一区二区三区| 亚洲久久久国产精品| 精品第一国产精品| 99re在线观看精品视频| 色婷婷av一区二区三区视频| 大型黄色视频在线免费观看| 性少妇av在线| 乱人伦中国视频| 国产麻豆69| 国产精品综合久久久久久久免费 | 日韩欧美免费精品| 欧美乱码精品一区二区三区| 国产午夜精品久久久久久| 高清在线国产一区| 日韩欧美国产一区二区入口| 一边摸一边抽搐一进一小说| av电影中文网址| 国产有黄有色有爽视频| 成人精品一区二区免费| 久久香蕉国产精品| 女人被狂操c到高潮| 热re99久久精品国产66热6| 又大又爽又粗| 欧美乱码精品一区二区三区| 中文字幕人妻丝袜制服| 丝袜美足系列| 一个人免费在线观看的高清视频| 中文字幕av电影在线播放| 亚洲精品粉嫩美女一区| 怎么达到女性高潮| 免费在线观看完整版高清| 狂野欧美激情性xxxx| 久久人人97超碰香蕉20202| 婷婷六月久久综合丁香| 宅男免费午夜| 久久青草综合色| 久久精品人人爽人人爽视色| 高清欧美精品videossex| 人成视频在线观看免费观看| 亚洲五月色婷婷综合| 欧美精品一区二区免费开放| 国产欧美日韩一区二区三区在线| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩无卡精品| 精品午夜福利视频在线观看一区| 国产不卡一卡二| 国产精品久久电影中文字幕| 18禁国产床啪视频网站| 日韩大尺度精品在线看网址 | 天堂动漫精品| 成人精品一区二区免费| 免费女性裸体啪啪无遮挡网站| av免费在线观看网站| 大香蕉久久成人网| 大型黄色视频在线免费观看| 欧美中文日本在线观看视频| 日本 av在线| 在线播放国产精品三级| 国产精品二区激情视频| 午夜日韩欧美国产| 欧美激情高清一区二区三区| 水蜜桃什么品种好| 亚洲色图 男人天堂 中文字幕| 亚洲av成人不卡在线观看播放网| 99久久精品国产亚洲精品| 国产又色又爽无遮挡免费看| 精品第一国产精品| 91精品三级在线观看| 欧美乱码精品一区二区三区| 黑丝袜美女国产一区| 岛国在线观看网站| 高清毛片免费观看视频网站 | 免费高清在线观看日韩| 日韩免费高清中文字幕av| 一夜夜www| 国产乱人伦免费视频| 十八禁人妻一区二区| 成人av一区二区三区在线看| 国产精品秋霞免费鲁丝片| 天天影视国产精品| 在线播放国产精品三级| 久久天堂一区二区三区四区| 国产日韩一区二区三区精品不卡| 女生性感内裤真人,穿戴方法视频| 一本综合久久免费| a级毛片在线看网站| 国产精品一区二区精品视频观看| 一级片免费观看大全| 亚洲欧美激情在线| 男人操女人黄网站| 韩国精品一区二区三区| 首页视频小说图片口味搜索| 老熟妇仑乱视频hdxx| x7x7x7水蜜桃| 亚洲国产精品999在线| 亚洲国产中文字幕在线视频| 精品国产美女av久久久久小说| 亚洲欧美日韩高清在线视频| videosex国产| 欧美国产精品va在线观看不卡| 日韩人妻精品一区2区三区| 不卡一级毛片| 极品人妻少妇av视频| 精品国产一区二区三区四区第35| 亚洲成a人片在线一区二区| 午夜a级毛片| av网站在线播放免费| 欧美激情 高清一区二区三区| 啦啦啦在线免费观看视频4| 国内毛片毛片毛片毛片毛片| 97人妻天天添夜夜摸| 欧美色视频一区免费| 国产精品影院久久| 美女扒开内裤让男人捅视频| 亚洲男人天堂网一区| 91精品国产国语对白视频| 成人三级黄色视频| 最近最新中文字幕大全免费视频| 久久这里只有精品19| 国产蜜桃级精品一区二区三区| 怎么达到女性高潮| 中文字幕另类日韩欧美亚洲嫩草| av免费在线观看网站| 制服诱惑二区| 亚洲国产毛片av蜜桃av| 黑人欧美特级aaaaaa片| 国产av精品麻豆| 久久香蕉精品热| 桃色一区二区三区在线观看| 男女做爰动态图高潮gif福利片 | 日韩三级视频一区二区三区| 精品久久久精品久久久| 亚洲五月天丁香| 99久久久亚洲精品蜜臀av| 97超级碰碰碰精品色视频在线观看|