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

    基于雙折射晶體的快拍穆勒矩陣成像測偏原理分析?

    2018-06-14 08:49:24曹奇志元昌安胡寶清任文藝趙銀軍張晶李建映鄧婷MingwuJin
    物理學(xué)報 2018年10期

    曹奇志元昌安胡寶清任文藝趙銀軍張晶李建映 鄧婷 Mingwu Jin

    1)(廣西師范學(xué)院物理與電子工程學(xué)院,南寧 530023)

    2)(廣西師范學(xué)院北部灣環(huán)境演變與資源利用教育部重點實驗室,廣西地表過程與智能模擬重點實驗室,南寧 530023)

    3)(西北農(nóng)林科技大學(xué)理學(xué)院,楊凌 712100)

    4)(Department of Physics,University of Texas at Arlington,Arlington,TX 76019,USA)

    5)(Department of Electrical and Computer Engineering,University of Delaware,Newark,DE 19716,USA)

    1 引 言

    偏振成像技術(shù)是利用光電成像器件獲取目標景物輻射的偏振信息,不僅可獲得目標光學(xué)輻射的光強度信息,而且可獲得目標的偏振信息,顯著增加被探測目標場景的信息量[1?3].目標的偏振信息與其自身的介電常數(shù)、粗糙度、組織結(jié)構(gòu)、含水量等有密切關(guān)系[4].穆勒矩陣是一種應(yīng)用廣泛的偏振表征量,能夠全面反映目標的偏振光學(xué)特性.它是一個4×4的矩陣,共16個陣元.穆勒矩陣成像測偏技術(shù)(Mueller matrix imaging polarimeter,MMIP)是一種能完全表征目標偏振屬性的偏振成像技術(shù)[5,6],是當今偏振成像領(lǐng)域研究的熱點之一[5?17].

    穆勒矩陣成像測偏技術(shù)具有光學(xué)成像非侵入、無輻射、無損傷等優(yōu)點,按照調(diào)制方式不同,可分為分時型和快拍型[10,18].分時型一般包含有機械轉(zhuǎn)動(如旋轉(zhuǎn)玻片等)[19?21]或相位延遲調(diào)制(如液晶等)[22?34]部件,其偏振態(tài)產(chǎn)生部分(PSG)和分析部分(PSA)都是沿著時間軸調(diào)制.分時型穆勒矩陣成像測偏儀雖然原理簡單,獲得的偏振圖像空間分辨率高,但系統(tǒng)中含有活動部件和需要多次測量,這限制了其測量精度和應(yīng)用范圍[25].

    快拍(snapshot,也有學(xué)者稱為同時或快照式)穆勒矩陣成像測偏技術(shù)(SMMIP)是采用不同的空間載頻將目標的全部偏振信息調(diào)制到一幀干涉圖中,通過一次曝光獲取目標的全部16個穆勒矩陣陣元圖像[12,25],其測量速度最快[12,25].SMMIP是在空間調(diào)制快拍斯托克斯(Stokes)成像測偏技術(shù)上發(fā)展起來的[35].空間調(diào)制快拍斯托克斯成像測偏技術(shù)是采用不同的載頻將目標的全部4個斯托克斯參量(S0,S1,S2和S3)調(diào)制到一幀干涉圖中,通過一次曝光獲取目標全部斯托克斯參量[35?40].它按核心調(diào)制器件可分為偏振光柵型[35]和雙折射型(包括楔形棱鏡[36]、薩瓦偏光鏡[37,38]和改進型薩瓦偏光鏡(modified savart polariscopes,MSP)[39?41]).2015年,Wang等[12]和Kudenov等[25,35]開展了以偏振光柵為核心調(diào)制器件的SMMIP研究,并展示了其在眼部疾病(青光眼)早期診斷方面的潛在應(yīng)用價值,但獲取目標的穆勒矩陣圖像空間分辨率較低;具有高消光比的雙折射晶體備受空間調(diào)制快拍斯托克斯成像測偏技術(shù)研究者們的青睞[26?41],并且已有研究表明以MSP為核心調(diào)制器件的技術(shù)方案在通道寬度、空間分辨率、信噪比和視場等方面具有顯著優(yōu)勢[39,40].本文主要對以MSP為核心調(diào)制器件的SMMIP的基本原理進行分析.

    2 光路設(shè)計與測量原理

    圖1是MSP-SMMIP光學(xué)設(shè)計示意圖,它包含一個偏振態(tài)產(chǎn)生部分和一個偏振態(tài)分析部分,主要包括光源(source),準直透鏡(f1),起偏器(P1),改進型薩瓦偏光鏡(MSP1),半波片(HWP1),MSP2,成像鏡(f2),被測樣品(sample),準直透鏡(f3),MSP3,HWP2,MSP4,檢偏器(P2),二次成像鏡(f4)及焦平面陣列(FPA).其中MSP-SMMIP系統(tǒng)中核心分光器件是MSP,可以將其看成兩塊薩瓦板(SP)夾一個半波片的“三明治”結(jié)構(gòu),其中MSP3和MSP4的SP厚度相等,它們是MSP1和MSP2的SP厚度的兩倍.準單色光源Source發(fā)出光束經(jīng)過f1準直后,入射到起偏器P1上,經(jīng)過P1的透射光可以表示成入射光的歸一化斯托克斯參量形式,Sin=(S0,in(x,y)/2)[1 0 1 0]T,上標T代表矩陣中的轉(zhuǎn)置,S0,in(x,y)是入射光的總光強.透射光經(jīng)過MSP1,HWP1和MSP2后分成四束光,這些光線經(jīng)過成像鏡f2后,形成干涉條紋定位在樣品上.這些條紋通過兩種載頻調(diào)制樣品的穆勒矩陣M(x,y),被調(diào)制后的光場通過f3準直進入到偏振態(tài)分析模塊中.從偏振態(tài)產(chǎn)生模塊中出來的四束光,經(jīng)過MSP3,HWP2和MSP4后分成16束光.最后,這16束光線經(jīng)過檢偏器P2和二次成像鏡f4干涉、成像于焦平面陣列FPA上.

    圖1 MSP-SMMIP光學(xué)設(shè)計示意圖 起偏器P1和檢偏器P2的偏振化方向為45?,半波片HWP1和HWP2的快軸方向為22.5?,改進型薩瓦偏光鏡MSP1和MSP3沿著y軸方向剪切,MSP2和MSP4沿著x軸方向剪切,MSP3和MSP4的SP單板厚度是MSP1和MSP2的兩倍Fig.1.Optical layout of MSP-SMMIP,P1and P2are linear polarizers at 45?,half wave-plates,HWP1and HWP2have fast axes oriented at 22.5?.MSP1and MSP3shear the beam along x while MSP2and MSP4shear along y.The single SP’s thickness of the generator’s and analyzer’s MSP are t and 2t,respectively.

    利用斯托克斯矢量-穆勒矩陣(Stokes vector-Mueller matrix)形式[6],可以方便地研究空間調(diào)制穩(wěn)態(tài)快拍穆勒矩陣成像測偏系統(tǒng)的探測原理.設(shè)入射光的Stokes參量為S0,in(x,y),用一個4×1的矩陣表示.出射光的Stokes參量S0,out(x,y)等于光學(xué)系統(tǒng)的Mueller矩陣M乘以入射光的Stokes參量S0,in(x,y),即S0,out(x,y)=MS0,in(x,y).基于本測量系統(tǒng)的基本結(jié)構(gòu)和探測原理,設(shè)f1=f2=f3=f4時,入射光經(jīng)過系統(tǒng)的矩陣傳輸方程可表示為

    其中,MP2(45?),M4,MH2(22.5?),M3,M(x,y),M2,MH1(22.5?),M1 和MP1(45?)分別表示檢偏器P2、改進型薩瓦偏光鏡MSP4、半波片HWP2、改進型薩瓦偏光鏡MSP3、樣品、改進型薩瓦偏光鏡MSP2、半波片HWP1、改進型薩瓦偏光鏡MSP1和起偏器P1的Mueller矩陣[6].

    通過計算(1)式,可得到焦平面FPA上光強I(x,y)=S0,out(x,y):

    其中k=2π?,?=?/(λf)是空間載頻,λ是入射光的波長,f是成像鏡的焦距,?是薩瓦偏光鏡單板橫向剪切量[39].為了重構(gòu)樣品的Mueller矩陣陣元mij(下標i和j是4×4的Mueller矩陣的行和列),需對干涉圖I(x,y)進行傅里葉變換將其轉(zhuǎn)換到頻域I′(x,y)(如圖2所示),可以看出該圖中有33個分離的峰(通道),這些通道復(fù)系數(shù)Sn分別包含樣品Mueller矩陣不同陣元信息(見表1).

    為解調(diào)出這些Mueller矩陣陣元,需要采用二維濾波器對相應(yīng)通道進行濾波.當某個k通道被濾波,再對該通道進行反傅里葉變換獲得Ck,它包含復(fù)系數(shù)Sn和調(diào)制相位因子.為了解調(diào)出Sn,需要一個或者多個已知Mueller矩陣Mr作為參考數(shù)據(jù)去定標每個通道的調(diào)制相位因子[6,8].

    (3)式中Csk和Crk分別是樣品和參考數(shù)據(jù)的第k通道數(shù)據(jù).此外,Cs0和Cr0是用來消去光源對樣品測試結(jié)果的影響.求解各Mueller矩陣分量的等式見表2.

    表1 每個通道復(fù)系數(shù)Sn(包含樣品被編碼到干涉圖中的Mueller矩陣陣元)Table 1.The complex coefficients Snof each channel(including the Mueller matrix elements encoded in the interferogram).

    表2 從頻域中求解Mueller矩陣分量Table 2.Resolved Muller matrix elements mijfrom the Fourier domain.

    圖2 干涉圖頻譜,包含33個通道Fig.2.Fourier spectra of interferogram.It includes 33 channels.

    3 光學(xué)指標

    以上對MSP-SMMIP的工作原理和特性進行了分析,下面將以特定的面陣探測器來分析設(shè)計系統(tǒng)的性能參數(shù).面陣探測器對空間調(diào)制快拍穆勒矩陣測偏儀的性能起著決定性作用,其幀頻決定快拍測量速度,其分辨率是決定測偏儀的空間分辨率的重要因素.由于MSP-SMMIP空間通道數(shù)達到33個,為了使得每個通道容納的信息盡可能多,必須充分利用面陣探測器的整個探測面陣,以獲得最大的空間分辨率.因此,應(yīng)以可獲得的面陣探測器來制定相關(guān)技術(shù)指標.本文以加拿大DALSA公司生產(chǎn)的Pantera TF 1M60為例制定指標.CCD的參數(shù):空間分辨率為1024×1024,像元大小為12μm×12μm,最高幀頻可達60幀/s.

    由于 CCD上沿y方向的一行像元數(shù)為N=1024,且相鄰像元間距為a=12μm.由圖3可知,在頻域中,沿著y方向需分成7個通道,每個通道頻域?qū)挾鹊扔诳臻g載頻?,為了盡可能地利用探測器陣列全部陣元,每個通道空域?qū)挾葹長=1/?=Na/7≈1.755 mm,則載頻?=1/L=0.57 mm?1. 在頻域中,如果采用矩形濾波器濾波,則濾波器截止頻率為?c=?/2=0.285 mm?1.系統(tǒng)的空間分辨率由系統(tǒng)的放大倍數(shù)和像元大小共同決定ε=7a/Q,其中a為像素尺寸,Q為光學(xué)系統(tǒng)的總放大倍率,乘以7是因為整個探測器陣面沿著y(和x)方向分成了7個通道.可見快拍穆勒矩陣成像測偏儀測量速度提高16倍是以空間分辨率降低7倍為代價的.若系統(tǒng)的總放大倍率為Q=1,空間分辨率就為ε=84μm.上述系統(tǒng)指標確定以后,由于?=?/λf,λ是入射光的波長,f是成像鏡的焦距,?是薩瓦偏光鏡單板橫向剪切量.因此可根據(jù)這些指標來設(shè)計關(guān)鍵光學(xué)元器件的參數(shù):如改進型薩瓦偏光鏡的厚度;物鏡、準直透鏡和成像鏡的焦距等.

    4 數(shù)值模擬

    參考上節(jié)光學(xué)指標,基于MATLAB軟件平臺,進行模擬實驗,除對模擬結(jié)果進行視覺定性評估,也采用相關(guān)系數(shù)(correlation coefficients,CC)和峰值信噪比(peak signal to noise ratio,PSNR)進行定量評估.

    4.1 模擬及其定性分析

    圖3(a)為輸入系統(tǒng)的Mueller矩陣目標圖像,該圖像十分有趣,圖像中心的空間頻率明顯高于邊緣,其中部分圖像尤為明顯(如m33等).圖4是基于(2)式模擬得到的干涉圖,從圖中可以看到清晰的干涉條紋.圖4的傅里葉變換譜見圖2,可以看出該圖中有33個分離的峰(通道),這表明模擬的圖像中包含了相應(yīng)的載頻.同時這33個峰(通道)分別包含著輸入Mueller矩陣16個陣元的信息,如果要解調(diào)出這些Mueller矩陣陣元,需要對相應(yīng)通道采用二維濾波器進行濾波,再進行二維反傅里葉變換.為了消去調(diào)制相位因子,需要對每個通道進行定標,被用于定標的參考目標必須確保每個通道的定標數(shù)據(jù)不為零.MSP-SMMIP光學(xué)系統(tǒng),只需要偏振化方向為22.5?的線性偏振片作為參考目標,一次即可將33個通道全部定標.相比文獻[25],采用多器件多角度的多次定標,MSP-SMMIP定標更加簡潔高效,兩者差異主要是因為兩種光學(xué)系統(tǒng)采用的核心分光器件不同,分光機理也不同,對偏振態(tài)的調(diào)制結(jié)果不同.圖3(b)為反演出的Mueller矩陣圖像,該圖很好地反演出了輸入Mueller矩陣圖像的低頻信息(如m12).圖3(b)部分分量(如m33)中心部分模糊,由第4節(jié)光學(xué)指標計算中可知,這主要是部分目標高頻信息超出了低通濾波器的截止頻率,而采用頻域反演技術(shù)時,低通濾波器會截斷超出截止頻域的高頻信息部分.這樣會導(dǎo)致高頻信息在濾波過程中丟失,而高頻信息對應(yīng)的是目標空間頻率高的部分,即目標的細節(jié)(如m33的中心部分),目標細節(jié)信息丟失會導(dǎo)致目標圖像中心高頻部分模糊.需要特別指出的是,圖3是偽彩色圖像,其顏色表示該數(shù)值大小,具體顏色代表那個數(shù)值可查看圖3下部的顏色條.

    圖3 (a)模擬輸入目標圖像;(b)測量結(jié)果圖像Fig.3.(a)Simulated input;(b)simulated measurement.

    圖4 MSP-SMMIP模擬干涉圖Fig.4.Simulated image of MSP-SMMIP.

    4.2 反演圖像質(zhì)量定量評價

    相關(guān)系數(shù)CC和PSNR可以較好地從圖像的相關(guān)性和絕對誤差兩方面評價圖像的重構(gòu)質(zhì)量.首先采用相似性評估的方法來評價反演圖像和輸入目標圖像的相似性.相似性測量主要是計算兩者的CC(其值范圍在0—1之間),CC越高表明這兩幅圖像越相似,也就是該反演圖像與輸入目標圖像越相似,重構(gòu)得越好.圖5所示為反演圖像與輸入目標圖像的CC,可以看出,CC均在0.82以上,表明反演圖像與輸入圖像之間存在強相關(guān),亦表明獲得了非常好的測量結(jié)果.

    圖5 反演圖像與輸入目標圖像的CCFig.5.The correlation coefficients between reconstructed image and input image.

    圖6 峰值信噪比Fig.6.Peak signal to noise ratio.

    PSNR經(jīng)常用作圖像重建質(zhì)量評價.PSNR表示信號最大可能功率和影響它的表示精度的破壞性噪聲功率的比值,由于許多信號都有非常寬的動態(tài)范圍,峰值信噪比常用對數(shù)分貝單位來表示.一般來說,PSNR值越大,重構(gòu)質(zhì)量越好,絕對誤差越小.圖6是反演穆勒矩陣每個陣元的PSNR,從圖中可知,單獨占有一個通道的穆勒矩陣陣元(如m12和m21)的PSNR比多個穆勒矩陣陣元(如m24和m42)共享一個通道高.此外,從圖中可知PSNR均在23以上,這表明測量結(jié)果絕對誤差較小.

    在此需特別指出,本文主要聚焦于介紹MSPSMMIP探測的新原理和系統(tǒng)的光學(xué)指標初步設(shè)計,給出基于理想情況下的理論分析和數(shù)值模擬結(jié)果.相應(yīng)的實驗室原理驗證實驗、定標實驗和定量分析等將在未來的工作中呈現(xiàn).

    5 結(jié) 論

    提出了一種以改進型薩瓦偏光鏡為核心器件的空間調(diào)制穩(wěn)態(tài)快拍穆勒矩陣成像測偏新技術(shù),它采用不同的空間載頻將16個穆勒矩陣陣元調(diào)制到一幅干涉圖中,能實現(xiàn)一次曝光獲取目標圖像和全部穆勒矩陣陣元圖像.采用斯托克斯矢量-穆勒矩陣形式闡明MSP-SMMIP的探測原理,研究了Mueller矩陣重構(gòu)和系列定標技術(shù),分析了系統(tǒng)的光學(xué)指標,采用計算機仿真實驗驗證了MSPSMMIP原理方案的正確性.對仿真結(jié)果的主觀和客觀評價表明模擬測量結(jié)果較好.與分時型相比,該技術(shù)的優(yōu)勢在于:1)測量系統(tǒng)無運動部件,結(jié)構(gòu)緊湊,體積小,系統(tǒng)的可靠性和穩(wěn)定性好;2)測量系統(tǒng)不受探測目標和系統(tǒng)自身相互運動以及外界環(huán)境擾動的影響,偏振探測精度高;3)探測系統(tǒng)通過一次測量即可獲取探測目標的全部偏振特性,探測速度快,可以用于靜、動態(tài)目標和動態(tài)場景的探測.所以,該技術(shù)可極大地滿足生物醫(yī)學(xué)(如眼部疾病)、遙感(如動態(tài)目標識別)和工業(yè)(如流水線質(zhì)量監(jiān)控)等領(lǐng)域?qū)铙w、運動目標或動態(tài)場景同時和實時快速檢測的實際需求.然而,現(xiàn)有頻域反演技術(shù)的局限性導(dǎo)致Mueller矩陣陣元重構(gòu)過程中目標高頻信息丟失,高頻信息對應(yīng)著目標的細節(jié),這導(dǎo)致該技術(shù)暫時不能應(yīng)用于高分辨領(lǐng)域,這個瓶頸急需突破.

    [1]Snik F,Craven-Jones J,Escuti M,Fineschid S,Harringtone D,Martinof A D,Mawetg D,Riedih J,Tyo J S『2014Proc.SPIE 9099,Polarization:Measurement,Analysis,and Remote Sensing XIBaltimore,Maryland,United States,March 24–28,2014 p90990B』

    [2]Tyo J,Goldstein L,Chenault B,Shaw A 2006Appl.Opt.45 5453

    [3]Li S J,Jiang H L,Zhu J P,Duan J,Fu Q,Fu Y G,Dong K Y 2013Chin.Opt.6 803(in Chinese)[李淑軍,姜會林,朱京平,段錦,付強,付躍剛,董科研 2013中國光學(xué) 6 803]

    [4]Bass M,Mahajan N 2010Handbook of Optics,3rd Edition,Volume I:Geometrical and Physical Optics,Polarized Light,Components and Instruments(New York:McGraw Hill)pp478–512

    [5]Wang Y,He H,Zeng N,Xie J,Liao R,Chang J,Sun M,Ma H 2015World J.Complex Med.1 74(in Chinese)[王曄,何宏輝,曾楠,謝軍,廖然,常金濤,孫明皓,馬輝 2015世界復(fù)合醫(yī)學(xué)1 74]

    [6]Goldstein D H 2010Polarized Light(3rd Ed.)(Boca Raton:CRC Press)pp628–631

    [7]Lu S Y,Chipman R 1996J.Opt.Soc.Am.A13 1106

    [8]Guo Y H 2014Ph.D.Dissertation(Shenzhen:Tsinghua University)pp10–30(in Chinese)[郭亦鴻2014博士學(xué)位論文 (深圳:清華大學(xué))第11—30頁]

    [9]He H,Zeng N,Du E,Guo Y,Li D,Liao R,Ma H 2013Photon.Lasers Med.2 129

    [10]Alali S,Vitkin A 2015J.Biomed.Opt.20 611041

    [11]Pierangelo A,Nazac A,Benali A,Validire P,Cohen H,Novikova T,Ibrahim B,Manhas S,Fallet C,Antonelli M,Martino A 2013Opt.Express21 14120

    [12]WangY,Kudenov M,KashaniA,KashaniA,Schwiegerling J,Escuti M『2015Proc.SPIE 9613,Polarization Science and Remote Sensing VIISan Diego,California,United States September 1,2015 p96130A』

    [13]Bueno M,Artal P 1999Opt.Lett.24 64

    [14]Li Y N,Sun X B,Qiao Y L,Zhang Q,Hong J 2010J.Atmos.Environ.Opt.5 203(in Chinese)[李雅男,孫曉兵,喬延利,張蕎,洪津2010大氣與環(huán)境光學(xué)學(xué)報5 203]

    [15]Qiao L F,Zhang Y M,Xie Q Y,Fang J,Wang J J,Zhang Q X 2009J.Combust.Sci.Technol.15 172(in Chinese)[喬利鋒,張永明,謝啟源,方俊,王進軍,張啟興2009燃燒科學(xué)與技術(shù)15 172]

    [16]Wang Y,Wang H,Chen S,Qi Y,Liu H 2011Proc.SPIE 8192,International Symposium on Photoelectronic Detection and Imaging 2011 Laser Sensing and Imaging;and Biological and Medical Applications of Photonics Sensing and ImagingBeijing,China,August 23,2011 p81923U1

    [17]Xiao S A,Li G H,Li J Z,Xue Q W 1991J.Qufu Normal Univ.18 44(in Chinese)[肖勝安,李國華,李繼仲,薛慶文1991曲阜師范大學(xué)學(xué)報18 44]

    [18]Wang X,Zhang M Y,Chen Z Y,Bai X F,Jin W Q 2013Infrared Laser Eng.42 2244(in Chinese)[王霞,張明陽,陳振躍,拜曉鋒,金偉其2013紅外與激光工程42 2244]

    [19]Azzam A 1978Opt.Lett.2 148

    [20]Collins W,Koh J 1999J.Opt.Soc.Am.A16 1997

    [21]Gerligand Y,Smith M,Chipman R 1999Opt.Express4 420

    [22]Laude-Boulesteix B,Martino A,Drévillon B,Schwartz L 2004Appl.Opt.43 2824

    [23]Benkelfat E,Horache H,Zou Q,Vinouze B 2003Opt.Commun.221 271

    [24]Jellison E,Modine A 1997Appl.Opt.36 8190

    [25]Kudenov M,Escuti M,Hagen N,Dereniak E,Oka K 2012Opt.Lett.37 1367

    [26]Goldstein D H 1992Appl.Opt.31 6676

    [27]Pezzaniti L,Chipman R 1995Opt.Eng.34 1558

    [28]Du E,He H,Zeng N,Sun M,Guo Y,Wu J,Ma H 2014J.Biomed.Opt.19 760131

    [29]Sun M,He H,Zeng N,Du E,Guo Y,Liu S,Wu J,He R,Ma H 2014Biomed.Opt.Express5 4223

    [30]Li J H,Zheng M,Zhang X B,Li Y Q 2016Laser Optoelectron.Prog.53 212021(in Chinese)[李建慧,鄭猛,張雪冰,李艷秋2016激光與光電子學(xué)進展53 212021]

    [31]Deng Y 2005Ph.D.Dissertation(Wuhan:Huazhong University of Science and Technology)pp83–97(in Chinese)[鄧勇2005博士學(xué)位論文 (武漢:華中科技大學(xué))第83–97頁]

    [32]Zhang X G,Jiang Y S,Lu X M 2008Acta Opt.Sin.28 1191(in Chinese)[張緒國,江月松,路小梅2008光學(xué)學(xué)報28 1191]

    [33]Chen X G,Liu S Y,Zhang C W,Wu Y P,Ma Z C,Sun T Y,Xu Z M 2014Acta Phys.Sin.63 180701(in Chinese)[陳修國,劉世元,張傳維,吳懿平,馬智超,孫堂友,徐智謀2014物理學(xué)報63 180701]

    [34]Li Y B,Li S B,Chen W J,Zeng Y X,Yang J Y 2010Opto-Electron.Eng.37 41(in Chinese)[李宇波,李世博,陳偉堅,曾宇驍,楊建義2010光電工程37 41]

    [35]Kudenov M,Escuti M,Dereniak E Oka K 2011Appl.Opt.50 2283

    [36]Oka K,Kaneko T 2003Opt.Express11 1510

    [37]Luo H,Oka K,DeHoog E,Schiewgerling J,Dereniak E L 2008Appl.Opt.47 4413

    [38]Luo H T 2008Ph.D.Dissertation(Arizona:Universityof Arizona)pp21–54

    [39]Cao Q Z,Zhang C M,DeHoog E 2012Appl.Opt.51 5791

    [40]Cao Q Z,Zhang J,DeHoog E,Zhang C M 2016Appl.Opt.55 954

    [41]Cao Q Z,Zhang J,DeHoog E,Lu Y,Hu B Q,Li W G,Li J Y,Fan D X,Deng T,Yan Y 2016Acta Phys.Sin.65 050702(in Chinese)[曹奇志,張晶,Edward DeHoog,盧遠,胡寶清,李武鋼,李建映,樊東鑫,鄧婷,閻妍2016物理學(xué)報65 050702]

    日韩精品青青久久久久久| 99re在线观看精品视频| 一二三四社区在线视频社区8| 51午夜福利影视在线观看| 亚洲人成电影免费在线| 精品一区二区三区av网在线观看| 一边摸一边抽搐一进一出视频| 女人高潮潮喷娇喘18禁视频| 久久久久久国产a免费观看| 国产精华一区二区三区| 国产日韩一区二区三区精品不卡| 又黄又爽又免费观看的视频| 国产在线观看jvid| 国产精品av久久久久免费| 国产精品美女特级片免费视频播放器 | 久久香蕉激情| 99精品在免费线老司机午夜| 欧美一区二区精品小视频在线| 在线免费观看的www视频| 亚洲精华国产精华精| 中出人妻视频一区二区| 亚洲国产精品合色在线| 国产主播在线观看一区二区| 人成视频在线观看免费观看| 亚洲avbb在线观看| 一级黄色大片毛片| av视频免费观看在线观看| 两个人视频免费观看高清| 国产成人一区二区三区免费视频网站| 亚洲 欧美一区二区三区| 亚洲专区国产一区二区| 国产不卡一卡二| 国产精品av久久久久免费| 日本vs欧美在线观看视频| 99国产精品99久久久久| 精品人妻1区二区| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲国产精品久久男人天堂| 亚洲欧美激情在线| 欧美黑人欧美精品刺激| 日韩有码中文字幕| 99国产精品99久久久久| 在线视频色国产色| 51午夜福利影视在线观看| 少妇的丰满在线观看| 国产亚洲精品av在线| 亚洲 国产 在线| 日韩三级视频一区二区三区| 99在线视频只有这里精品首页| 国产精品一区二区免费欧美| 成人特级黄色片久久久久久久| av视频免费观看在线观看| 中亚洲国语对白在线视频| 精品久久久久久成人av| 亚洲av五月六月丁香网| 午夜免费观看网址| 欧美色视频一区免费| 黑人欧美特级aaaaaa片| 老司机靠b影院| 一二三四社区在线视频社区8| 99在线视频只有这里精品首页| 中文字幕久久专区| 国产高清视频在线播放一区| 亚洲成av片中文字幕在线观看| 又大又爽又粗| xxx96com| 亚洲中文字幕一区二区三区有码在线看 | 免费看美女性在线毛片视频| 国产精品亚洲av一区麻豆| 50天的宝宝边吃奶边哭怎么回事| 亚洲成人久久性| 久久天堂一区二区三区四区| 麻豆国产av国片精品| √禁漫天堂资源中文www| 一级毛片精品| 国产精品亚洲美女久久久| 成人手机av| 可以在线观看的亚洲视频| 女人被狂操c到高潮| 国产高清videossex| 午夜日韩欧美国产| 亚洲熟女毛片儿| 色老头精品视频在线观看| 美女国产高潮福利片在线看| 国产亚洲精品av在线| 成人国产综合亚洲| 亚洲国产高清在线一区二区三 | 免费久久久久久久精品成人欧美视频| 国内久久婷婷六月综合欲色啪| 欧美一级a爱片免费观看看 | 久久午夜亚洲精品久久| 后天国语完整版免费观看| aaaaa片日本免费| 自拍欧美九色日韩亚洲蝌蚪91| 日韩国内少妇激情av| avwww免费| 亚洲五月天丁香| 级片在线观看| 国产99白浆流出| 国产午夜精品久久久久久| 国产精品美女特级片免费视频播放器 | 两个人免费观看高清视频| 久久久久国内视频| 国产成人免费无遮挡视频| www.自偷自拍.com| 免费看美女性在线毛片视频| 欧美大码av| 国产一级毛片七仙女欲春2 | 国产亚洲欧美98| 两性午夜刺激爽爽歪歪视频在线观看 | 在线十欧美十亚洲十日本专区| 国产一区在线观看成人免费| 午夜福利视频1000在线观看 | 免费观看人在逋| 老司机午夜福利在线观看视频| 9热在线视频观看99| 久热这里只有精品99| 国产亚洲av嫩草精品影院| 亚洲色图综合在线观看| 国产又爽黄色视频| 一级片免费观看大全| 国产精品永久免费网站| 欧美+亚洲+日韩+国产| videosex国产| 国产激情久久老熟女| 在线观看免费视频日本深夜| 免费高清视频大片| e午夜精品久久久久久久| 天天一区二区日本电影三级 | 麻豆一二三区av精品| 亚洲精品一区av在线观看| 日本免费a在线| 欧美另类亚洲清纯唯美| 亚洲国产欧美日韩在线播放| 国产伦人伦偷精品视频| xxx96com| 免费少妇av软件| 中亚洲国语对白在线视频| 亚洲欧美日韩高清在线视频| 久久久国产欧美日韩av| av超薄肉色丝袜交足视频| 欧美久久黑人一区二区| 久久国产乱子伦精品免费另类| av视频在线观看入口| 国产片内射在线| 999精品在线视频| 黑丝袜美女国产一区| 欧美一级毛片孕妇| 日韩中文字幕欧美一区二区| 大陆偷拍与自拍| 精品日产1卡2卡| 黄网站色视频无遮挡免费观看| 19禁男女啪啪无遮挡网站| 国产精品亚洲av一区麻豆| 久热爱精品视频在线9| 人人妻人人爽人人添夜夜欢视频| 国产97色在线日韩免费| 身体一侧抽搐| 黄频高清免费视频| 国内精品久久久久精免费| 日韩欧美国产一区二区入口| 夜夜爽天天搞| 国产精品一区二区免费欧美| 一级a爱片免费观看的视频| 欧美成人一区二区免费高清观看 | 性色av乱码一区二区三区2| 亚洲国产中文字幕在线视频| 精品电影一区二区在线| 这个男人来自地球电影免费观看| 真人做人爱边吃奶动态| 丰满的人妻完整版| 在线观看日韩欧美| 少妇熟女aⅴ在线视频| 国产一卡二卡三卡精品| 成年女人毛片免费观看观看9| 变态另类成人亚洲欧美熟女 | 欧美绝顶高潮抽搐喷水| 十八禁网站免费在线| 日本vs欧美在线观看视频| 国产精品1区2区在线观看.| 99国产精品免费福利视频| 99久久精品国产亚洲精品| 手机成人av网站| 99精品欧美一区二区三区四区| 久久这里只有精品19| 国产精品1区2区在线观看.| 亚洲中文av在线| 99久久综合精品五月天人人| 亚洲熟女毛片儿| 大陆偷拍与自拍| 亚洲欧洲精品一区二区精品久久久| 一级a爱视频在线免费观看| 欧美黄色片欧美黄色片| 亚洲成国产人片在线观看| 精品少妇一区二区三区视频日本电影| 黄色视频,在线免费观看| 天堂影院成人在线观看| 在线视频色国产色| 99久久国产精品久久久| 国产午夜福利久久久久久| 大陆偷拍与自拍| 国产精品一区二区在线不卡| 亚洲国产精品999在线| 成人手机av| 老司机深夜福利视频在线观看| 久久中文看片网| 国产欧美日韩综合在线一区二区| 亚洲午夜精品一区,二区,三区| 亚洲精品美女久久久久99蜜臀| 欧美黄色淫秽网站| 999久久久精品免费观看国产| 成在线人永久免费视频| 日韩大尺度精品在线看网址 | 亚洲情色 制服丝袜| 日本五十路高清| 久久天躁狠狠躁夜夜2o2o| 麻豆一二三区av精品| av视频免费观看在线观看| 精品无人区乱码1区二区| 九色亚洲精品在线播放| 国产亚洲欧美在线一区二区| 日本一区二区免费在线视频| 中文字幕高清在线视频| 中文亚洲av片在线观看爽| 丝袜在线中文字幕| 国产黄a三级三级三级人| 欧美色欧美亚洲另类二区 | 国产一区二区激情短视频| 桃红色精品国产亚洲av| 99久久综合精品五月天人人| 国产亚洲精品综合一区在线观看 | 亚洲一区高清亚洲精品| 亚洲av成人一区二区三| 韩国精品一区二区三区| 十八禁网站免费在线| 亚洲人成电影免费在线| 久久久久久久久中文| 在线观看www视频免费| 久久国产乱子伦精品免费另类| 日本免费a在线| 成熟少妇高潮喷水视频| 久久久久久久久中文| 久久久国产欧美日韩av| 色播亚洲综合网| 国产高清videossex| 国产三级在线视频| 久久性视频一级片| 欧美国产精品va在线观看不卡| 久久天堂一区二区三区四区| 校园春色视频在线观看| 精品久久久久久,| 国产色视频综合| 免费在线观看视频国产中文字幕亚洲| 国产精品自产拍在线观看55亚洲| 国产激情欧美一区二区| 国产一卡二卡三卡精品| svipshipincom国产片| 村上凉子中文字幕在线| 热99re8久久精品国产| 在线观看66精品国产| 中文字幕另类日韩欧美亚洲嫩草| 真人做人爱边吃奶动态| 99在线人妻在线中文字幕| 高清在线国产一区| 国产又爽黄色视频| 国产精品亚洲av一区麻豆| 12—13女人毛片做爰片一| 啪啪无遮挡十八禁网站| 亚洲一区二区三区不卡视频| videosex国产| 久久 成人 亚洲| 色综合亚洲欧美另类图片| 国产单亲对白刺激| 91成人精品电影| 亚洲成av人片免费观看| 国语自产精品视频在线第100页| 一个人免费在线观看的高清视频| 国产成人免费无遮挡视频| 国产私拍福利视频在线观看| 日韩欧美在线二视频| 色在线成人网| 国产精品美女特级片免费视频播放器 | 男女之事视频高清在线观看| 久久青草综合色| 久热这里只有精品99| 国产单亲对白刺激| 国产亚洲精品第一综合不卡| 日本vs欧美在线观看视频| 悠悠久久av| 亚洲精华国产精华精| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美+亚洲+日韩+国产| 午夜福利视频1000在线观看 | 在线观看免费午夜福利视频| 欧美日韩黄片免| 老鸭窝网址在线观看| 精品久久久久久,| 中文字幕精品免费在线观看视频| 神马国产精品三级电影在线观看 | 精品久久久久久成人av| 黄色成人免费大全| 国产精品久久久人人做人人爽| av超薄肉色丝袜交足视频| 九色国产91popny在线| 亚洲成人久久性| 亚洲人成77777在线视频| 免费在线观看亚洲国产| 欧美乱码精品一区二区三区| 久久狼人影院| 亚洲精品久久国产高清桃花| 国产av一区在线观看免费| 黄网站色视频无遮挡免费观看| 国产精品二区激情视频| 国产麻豆成人av免费视频| 韩国av一区二区三区四区| avwww免费| 99在线视频只有这里精品首页| 亚洲国产精品合色在线| 1024香蕉在线观看| 欧美人与性动交α欧美精品济南到| 黄色a级毛片大全视频| 欧美日韩福利视频一区二区| 午夜福利成人在线免费观看| 757午夜福利合集在线观看| 欧美性长视频在线观看| 丁香欧美五月| 日日干狠狠操夜夜爽| 精品日产1卡2卡| 国产成人欧美| 免费看美女性在线毛片视频| 国产精品影院久久| 国产三级黄色录像| 亚洲av电影不卡..在线观看| 亚洲精品一卡2卡三卡4卡5卡| 露出奶头的视频| 我的亚洲天堂| 亚洲精品久久成人aⅴ小说| 久久久久九九精品影院| 国产精品野战在线观看| 国产精品久久久人人做人人爽| 中文字幕色久视频| ponron亚洲| 电影成人av| 欧美 亚洲 国产 日韩一| 欧美激情极品国产一区二区三区| 一区在线观看完整版| 咕卡用的链子| 色综合亚洲欧美另类图片| 亚洲成国产人片在线观看| 日本免费一区二区三区高清不卡 | 国产日韩一区二区三区精品不卡| 久久久久亚洲av毛片大全| 18禁裸乳无遮挡免费网站照片 | av天堂久久9| 91字幕亚洲| 男女下面插进去视频免费观看| 中文字幕久久专区| 亚洲精品久久成人aⅴ小说| 国产熟女午夜一区二区三区| 午夜激情av网站| 91av网站免费观看| xxx96com| 国产一区在线观看成人免费| 国产成人av教育| 一二三四在线观看免费中文在| 国产伦人伦偷精品视频| 国产熟女午夜一区二区三区| 又黄又爽又免费观看的视频| 99国产精品免费福利视频| 国产97色在线日韩免费| 激情视频va一区二区三区| 免费高清在线观看日韩| 久久久久精品国产欧美久久久| 91成人精品电影| 国产麻豆成人av免费视频| 亚洲色图av天堂| 亚洲久久久国产精品| 婷婷六月久久综合丁香| 国产成人精品在线电影| 婷婷丁香在线五月| 手机成人av网站| 人人妻人人爽人人添夜夜欢视频| 免费不卡黄色视频| 欧美+亚洲+日韩+国产| 午夜精品在线福利| 一边摸一边做爽爽视频免费| 电影成人av| 91国产中文字幕| 热99re8久久精品国产| 欧美亚洲日本最大视频资源| av网站免费在线观看视频| 多毛熟女@视频| 国产真人三级小视频在线观看| 一进一出抽搐动态| 九色亚洲精品在线播放| 国产精品爽爽va在线观看网站 | av网站免费在线观看视频| 可以在线观看毛片的网站| 精品人妻在线不人妻| 麻豆国产av国片精品| 手机成人av网站| www.999成人在线观看| 国产麻豆成人av免费视频| 色老头精品视频在线观看| 长腿黑丝高跟| 成人18禁在线播放| 波多野结衣高清无吗| 丝袜人妻中文字幕| 香蕉久久夜色| 长腿黑丝高跟| 精品久久久久久成人av| 一进一出抽搐gif免费好疼| 91av网站免费观看| 亚洲色图综合在线观看| 国产高清激情床上av| 欧美色欧美亚洲另类二区 | 午夜福利视频1000在线观看 | av超薄肉色丝袜交足视频| 在线观看免费视频日本深夜| 亚洲一卡2卡3卡4卡5卡精品中文| 女人高潮潮喷娇喘18禁视频| 欧美日韩亚洲国产一区二区在线观看| 国产精品免费一区二区三区在线| 久久 成人 亚洲| 久久中文字幕人妻熟女| 久久久久国产精品人妻aⅴ院| 午夜福利在线观看吧| 亚洲欧美精品综合一区二区三区| 大香蕉久久成人网| 亚洲精品在线美女| 免费在线观看黄色视频的| 国产欧美日韩一区二区三区在线| 成人国语在线视频| 久久亚洲精品不卡| 午夜免费鲁丝| 久久久久久大精品| 国产精品乱码一区二三区的特点 | 每晚都被弄得嗷嗷叫到高潮| 国产亚洲欧美98| avwww免费| 欧美久久黑人一区二区| 国产成人精品在线电影| www日本在线高清视频| 1024视频免费在线观看| 国产精品综合久久久久久久免费 | 夜夜看夜夜爽夜夜摸| 亚洲人成电影免费在线| 欧美国产精品va在线观看不卡| 中文字幕另类日韩欧美亚洲嫩草| 变态另类成人亚洲欧美熟女 | 精品久久久久久,| 国产精品久久电影中文字幕| 亚洲av成人av| 又黄又粗又硬又大视频| av超薄肉色丝袜交足视频| 亚洲中文字幕一区二区三区有码在线看 | 天天一区二区日本电影三级 | 欧美+亚洲+日韩+国产| www国产在线视频色| 亚洲片人在线观看| 精品久久久久久久久久免费视频| 国产精品自产拍在线观看55亚洲| 欧美日韩精品网址| 欧美人与性动交α欧美精品济南到| 日本五十路高清| 亚洲男人的天堂狠狠| 99精品欧美一区二区三区四区| 看黄色毛片网站| 欧美成人性av电影在线观看| 亚洲人成网站在线播放欧美日韩| 热99re8久久精品国产| 欧美在线黄色| 最新美女视频免费是黄的| 欧美日本中文国产一区发布| 99国产精品一区二区三区| 亚洲五月天丁香| 国产高清视频在线播放一区| 搞女人的毛片| 视频在线观看一区二区三区| 大型黄色视频在线免费观看| 国产亚洲欧美精品永久| 亚洲男人天堂网一区| 在线国产一区二区在线| 亚洲国产欧美一区二区综合| 午夜视频精品福利| 欧美成人免费av一区二区三区| 99热只有精品国产| 午夜日韩欧美国产| 色精品久久人妻99蜜桃| 久久人人97超碰香蕉20202| 老司机午夜福利在线观看视频| 欧美老熟妇乱子伦牲交| 久久精品国产综合久久久| 母亲3免费完整高清在线观看| 老司机靠b影院| 国产成人欧美在线观看| 国产91精品成人一区二区三区| 成人欧美大片| 国产成人精品久久二区二区91| 亚洲欧美激情在线| av片东京热男人的天堂| 国产私拍福利视频在线观看| 91在线观看av| 日韩高清综合在线| 自线自在国产av| 成人三级做爰电影| 天堂动漫精品| 国产精品乱码一区二三区的特点 | 自线自在国产av| 亚洲av电影不卡..在线观看| 亚洲熟妇熟女久久| 色老头精品视频在线观看| 后天国语完整版免费观看| 波多野结衣av一区二区av| 高清在线国产一区| or卡值多少钱| av免费在线观看网站| 女警被强在线播放| 日韩 欧美 亚洲 中文字幕| 久久人人爽av亚洲精品天堂| 极品教师在线免费播放| 日韩中文字幕欧美一区二区| 在线av久久热| 国产91精品成人一区二区三区| 9191精品国产免费久久| 日韩高清综合在线| 男女午夜视频在线观看| 欧美久久黑人一区二区| 国产亚洲精品久久久久5区| 在线观看www视频免费| 久久久国产成人精品二区| 国产91精品成人一区二区三区| 黑人操中国人逼视频| a级毛片在线看网站| 国产欧美日韩一区二区三区在线| www.精华液| 亚洲精品粉嫩美女一区| 精品福利观看| 如日韩欧美国产精品一区二区三区| 亚洲国产精品sss在线观看| 欧美一级a爱片免费观看看 | 欧美日韩亚洲综合一区二区三区_| 国产亚洲精品久久久久久毛片| 操美女的视频在线观看| 88av欧美| www日本在线高清视频| 99精品久久久久人妻精品| 男女下面进入的视频免费午夜 | 曰老女人黄片| 人人妻人人澡人人看| 国产精品久久久久久精品电影 | 波多野结衣av一区二区av| 91老司机精品| 精品免费久久久久久久清纯| 激情在线观看视频在线高清| 欧美日韩乱码在线| 国产午夜福利久久久久久| 国产区一区二久久| 亚洲久久久国产精品| 欧美日韩瑟瑟在线播放| a在线观看视频网站| 国产午夜福利久久久久久| 国产欧美日韩综合在线一区二区| 一区福利在线观看| 在线国产一区二区在线| 69av精品久久久久久| 亚洲三区欧美一区| 一进一出抽搐动态| 国产精品久久久久久精品电影 | 搡老岳熟女国产| 日日夜夜操网爽| 亚洲精华国产精华精| 精品人妻1区二区| 久久人人爽av亚洲精品天堂| 两个人视频免费观看高清| 精品电影一区二区在线| 级片在线观看| 男女之事视频高清在线观看| 国产欧美日韩一区二区三区在线| 在线播放国产精品三级| 久久久精品欧美日韩精品| 国产午夜福利久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 色综合站精品国产| 亚洲中文字幕日韩| 精品午夜福利视频在线观看一区| 国产精品九九99| av免费在线观看网站| 夜夜夜夜夜久久久久| 美女高潮喷水抽搐中文字幕| 成人国产综合亚洲| 很黄的视频免费| 18禁国产床啪视频网站| 色婷婷久久久亚洲欧美| 精品国产一区二区久久| 国产一卡二卡三卡精品| 757午夜福利合集在线观看| 精品国产一区二区久久| 深夜精品福利| 欧美日韩精品网址| 女人精品久久久久毛片| 久久久水蜜桃国产精品网| 757午夜福利合集在线观看| 精品国产一区二区久久| 深夜精品福利| 午夜激情av网站| 亚洲一区高清亚洲精品| 女生性感内裤真人,穿戴方法视频| 在线观看免费视频日本深夜| 大码成人一级视频| 90打野战视频偷拍视频| 妹子高潮喷水视频| 亚洲av电影不卡..在线观看| 国产av精品麻豆| 亚洲精品中文字幕在线视频|