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

    一種雙重判斷機制的音頻篡改盲檢測算法*

    2014-08-16 08:00:24呂志勝胡永健李晗劉琲貝
    關(guān)鍵詞:子塊偏移量斜率

    呂志勝 胡永健 李晗 劉琲貝

    (華南理工大學 電子與信息學院,廣東 廣州 510640)

    如果數(shù)字錄音設備由市電供電,則在錄音過程中電網(wǎng)頻率(ENF)會在錄制的內(nèi)容中留下印跡,此印跡可用于音頻篡改檢測[1-2].文獻[3-6]中通過比較從待測語音信號中提取的ENF 與實時記錄的市電ENF 的連續(xù)性和一致性來判斷是否發(fā)生篡改.文獻[7-8]中指出日光燈的閃爍頻率是ENF 的兩倍,并將此類方法推廣到對室內(nèi)照明環(huán)境中錄制的視頻真實性進行檢測.這類方法的不足是需要大量參考ENF信號,且靈敏度不高,只能對較長的音頻(持續(xù)時間為幾分鐘甚至幾個小時)進行檢測.文獻[9-10]中利用不同電網(wǎng)的ENF 存在差異這一特點估計音頻錄制的區(qū)域,但不涉及對音頻篡改的討論.目前利用ENF 實現(xiàn)音頻篡改盲檢測的代表性工作有文獻[11-13].文獻[11]中利用頻譜距離和ENF 相位的突變來檢測音頻真實性,但對頻率的分辨率較敏感,只有分辨率足夠大時才能準確估計相位,否則誤差較大.文獻[12]中采用高精度傅里葉分析方法改進了文獻[11]中的方法,但計算量較大.筆者曾在文獻[13]中通過引入一個理想的正弦信號作為參考信號,計算待測音頻中ENF 信號各子塊與參考信號取得最大相關(guān)時的偏移量,通過最大相關(guān)偏移的變化情況來檢測音頻篡改.

    文中從兩個方面改進文獻[13]中的方法:一是提出一種不需要額外參考信號來計算各塊ENF 信號最大相關(guān)偏移的方法;二是利用最大相關(guān)偏移的變化及其極值點斜率變化聯(lián)合判斷篡改區(qū)域.所提出的方法更加便捷和準確.

    1 文獻回顧與問題的提出

    設讀取的待測音頻信號為s(n),其采樣頻率為fs,其中n 為采樣時刻.為減小運算量,對音頻信號進行下采樣得到x(n)=s(nM),其中M=fs/fd,fd為下采樣頻率.設電網(wǎng)固有頻率為f0,為了防止頻率混疊,取fd≥2f0.將x(n)通過一個中心頻率為f0的窄帶零相位帶通濾波器濾波,得到ENF 信號y(n)=xENF(n)+d(n),其中xENF(n)為理想的ENF 信號,d(n)為窄帶噪聲,n=1,2,…,L,其中L 為音頻數(shù)據(jù)的長度.

    為了利用所抽取的ENF 信號檢測音頻是否發(fā)生篡改,可將y(n)劃分為相互重疊的子塊,各子塊長度為ENF 固有周期的整數(shù)倍.設子塊的采樣點數(shù)為NW,則NW=C0N0,這里N0=fdf0,為一個ENF 周期內(nèi)的采樣點數(shù).通常相鄰子塊之間取(C0-1)N0的數(shù)據(jù)重疊.由于理想的xENF(n)為正弦信號,故第i 個子塊的信號可表示為

    式中,k=(i-1)N0+l,l=1,2,3,…,NW,i=1,2,…,NBLOCK,NBLOCK為子塊數(shù).θi為第i 個子塊的初始相位.A0是ENF 信號的幅度.具體估算初始相位的方法詳見文獻[12].

    若待測音頻沒有被改動過,則在不考慮干擾的情況下ENF 信號是一個連續(xù)光滑的正弦波,這個特性反映到各個子塊上,則表現(xiàn)為各個子塊的初始相位相等.文獻[12]中通過評估各子塊ENF 信號初始相位是否相等來確定所測音頻是否發(fā)生插入或刪除篡改.當處理真實世界的待測音頻時,必須考慮電網(wǎng)的波動.受電網(wǎng)波動干擾的ENF 信號是一個有畸變的正弦信號,此時各子塊的初始相位不再嚴格相等,而是有小幅度的波動,但在篡改區(qū)域,相鄰塊的相位會出現(xiàn)大幅度的跳變.下面用一個例子來說明.圖1 中的待測音頻有部分語音片段(第9.502 至第13.887s)被刪除,文中用文獻[12]中的方法來檢測篡改.先從待測音頻中抽取ENF 信號,對應的信號見圖1(a),然后分塊并計算各子塊的初始相位.圖1(b)中虛點線為正常未篡改時的各子塊相位的正常波動曲線,星型曲線為篡改后各子塊相位的波動曲線.可以看到大約在第500 塊處(即對應篡改處)相鄰塊的初始相位出現(xiàn)了很大的變化,這說明文獻[12]中的方法可以檢測出發(fā)生了篡改.但也注意到由于其跳變不陡峭,且跳變的起始點和終止點不清晰,故較難精確確定篡改發(fā)生的位置.

    圖1 文獻[12]和[13]中刪除篡改檢測的結(jié)果Fig.1 Detection results of deletion for Ref.[12]and[13]

    文中作者的早期工作于文獻[13]對此提出了一種改進算法,先采取同樣的分塊策略劃分ENF 信號,考慮到未篡改的音頻中理想的ENF 信號是一條正弦曲線,故引入一個標準的正弦信號作為參考信號,對其也進行相應分塊,然后計算出ENF 信號各子塊和參考信號對應子塊的相關(guān)性,將出現(xiàn)相關(guān)系數(shù)最大處的偏移稱為該子塊的最大相關(guān)偏移量(MOCC).若未發(fā)生篡改,各個子塊的最大相關(guān)偏移在理想情況下應該相等,而對來自真實世界的音頻信號,MOCC 值也應該相差不大.基于MOCC 進行篡改檢測的原理和技術(shù)細節(jié)詳見文獻[13],其檢測效果見圖1(c)和(d).由圖可見在篡改發(fā)生處相鄰子塊的MOCC 值發(fā)生了很大跳變,且跳變的邊緣較陡峭.文獻[13]中的不足是在未發(fā)生篡改的區(qū)域MOCC 曲線有城墻一樣的齒,容易引起誤判(即虛警),且緊鄰跳變處也有齒,不方便精確定位.文中針對這兩個問題進行進一步研究,提出利用斜率來解決上述問題,同時還提出一種免除引入?yún)⒖夹盘栔苯佑嬎鉓OCC 值的方法.

    2 文中算法

    文中直接以y(n)的起始子塊作為基準信號,計算各子塊與它的最大相關(guān)偏移.設基準信號y1為

    式中,偏移量=0,1,2,…,N0-1,l=1,2,3,…,NW,θ1為起始子塊的初始相位.當 取不同值時,分別計算各子塊與基準信號的相關(guān)系數(shù),對于第i 個子塊有:

    式中,η 為噪聲之間以及噪聲與ENF 信號之間的互相關(guān)之和.式(3)的具體推導和文獻[13]中計算子塊與參考塊之間的互相關(guān)類似.由第1 項可見,當-2 /N0+θi-θ1=0 時,Ri)取得最大值

    類似于文獻[13],將對應式(4)的偏移量 稱為最大相關(guān)偏移MOCC,記作max(i).對于信噪比較高的音頻信號有2?η,此時Ri()受噪聲的影響很小.

    針對前述MOCC 曲線的城墻齒,文中提出一種雙處理機制從MOCC 曲線得到一條新的曲線,它在未篡改區(qū)域平坦而在遭受篡改的地方突變.第1 層機制根據(jù)偏移量變化的子塊間隔時間和偏移量變化的累積確定篡改區(qū)域.首先計算相鄰子塊的最大相關(guān)偏移之差

    式中,i=1,2,3,…,NBLOCK-1.將δmax(i)≠0 的子塊進行標記,若有P 個此類子塊,則記為i1,i2,i3,…,iP.顯然,這類塊是出現(xiàn)篡改的潛在塊.然后計算子塊im和im-1之間的時間間隔Dt(m)以及偏移量變化的累積Dos(m):

    式中,m=2,3,…,P.文中通過設定相鄰偏移量變化累積的門限Tos來減小因噪聲干擾而引起的虛警.如果Dos(m)≥Tos,則 max的差異可能是由篡改引起.否則,max的差異被認為是由噪聲引起的波動.與此同時,文中還設定時間間隔門限Tt來減小因ENF 偏差引起的虛警.若Dt(m)<Tt,則表明在短時間內(nèi)有突變,應該是由篡改導致;否則,非短時間內(nèi)突變,即由非篡改引起的變化.文中設定Tt=αNT0,其中NT0= NBLOCK/(P +1)為 max 發(fā)生變化的平均時間間隔,α 為時間因子.通常,篡改引起的變化時間間隔會遠小于NT0,而非篡改引起的則大于NT0.將同時滿足Dos(m)≥Tos和Dt(m)≤Tt的im值記為 max突變的子塊序號.假設共有U 個這樣的子塊,用集合K 記為K=[k1k2k3… kU].如果K 不為空,則該音頻信號被篡改過;否則未被篡改.

    第2 層機制根據(jù)MOCC 的斜率變化確定篡改位置.首先計算δmax(i)≠0 的P 個子塊之間MOCC 的斜率,并令其余子塊的斜率為0,則全體子塊的斜率G1可表示為

    式中,m=2,3,…P,i=1,2,…,NBLOCK.對于未篡改的原始音頻,即使ENF 存在一定的偏差,其斜率也基本保持不變;相反,被篡改過的音頻在篡改邊界的斜率則會出現(xiàn)較大的突變.據(jù)此可設定斜率判別門限TG=βG0,其中,是斜率不為0 的P 個子塊的斜率絕對值均值,β 為幅度因子.通常,篡改處的斜率會大于G0,非篡改處的則小于G0.把的子塊認為是可能發(fā)生篡改的位置.設共有Q 個這樣的子塊,用集合J 記錄,則有J=[j1j2j3… jQ].

    將上述兩種判決機制聯(lián)合使用,取兩者的交集為最終篡改判決結(jié)果,即Z=J∩K.若Z 不為空,則音頻信號被篡改過;否則為原始信號.根據(jù)集合Z對斜率G1進行處理,保留篡改子塊的斜率,而其余子塊的斜率用所有子塊斜率的均值mG1代替,則可得到G2:

    在G2所對應的曲線上,Z 的第一個子塊和最后一個子塊之間即為發(fā)生篡改的部分.結(jié)合ENF 的時間周期(例如1/50 Hz=0.02 s),進一步定位到音頻信號中的篡改區(qū)域.定位精度為一個子塊的時間長度.根據(jù)定位的篡改區(qū)域為靜音部分還是語音部分,判斷篡改操作是刪除還是插入.

    3 實驗結(jié)果

    文中實驗所用測試音頻均在市電頻率為50 Hz時錄制,音頻的保存格式為WAV,下采樣頻率fd=1000 Hz.用零相位FIR 窄帶濾波器提取的ENF 信號中會有噪聲干擾以及音頻的泄漏頻譜,盡管減小濾波器帶寬可降低這兩方面的影響,但這種做法會使濾波的過渡過程變長,從而使篡改邊界變平坦,影響定位準確性.事實上,當帶寬小于0.6 Hz 時甚至會出現(xiàn)無法檢測到篡改點的情況.經(jīng)驗上帶寬在0.6~1.4 Hz 之間取值效果較好,文中取1.2 Hz.另一方面,數(shù)據(jù)子塊越長,得到的MOCC 曲線越平滑,但篡改定位誤差也會越大.通常子塊長度可在60~200之間取值,即對應于C0=3~10.文中取C0=3.此外,文中直接以待測音頻的起始子塊作為基準信號替代文獻[13]中的標準正弦信號來計算各塊的MOCC 值.

    3.1 刪除篡改檢測

    仍以圖1 中的音頻為例,圖2(a)是用最大相關(guān)偏移法計算的MOCC,圖2(b)顯示各段的斜率.圖2(c)是第1 層和第2 層機制聯(lián)合作用的結(jié)果.可以看到,除了篡改區(qū)間斜率不為0 外,其余區(qū)間的斜率均為0.圖2(d)給出在對應音頻上定位的篡改區(qū)域為第9.37~9.77 s,這與實際音頻刪除位置第9.502 s接近.定位存在微小誤差主要有兩個方面的原因:一是前述的窄帶濾波后信號突變位置存在過度過程[14];二是文中檢測篡改是以一個數(shù)據(jù)塊為最小單位,這對定位精度也有一定的影響.文中兩層機制中所使用的參數(shù)Tos=2,α=0.7,β=1.3.

    圖2 刪除篡改檢測Fig.2 Detection of deletion forgery

    3.2 插入篡改檢測

    圖3(a)是從一個有插入的音頻提取的ENF 信號,插入發(fā)生在第8.921~13.568 s.相位法和最大相關(guān)偏移法的結(jié)果分別見圖3(b)和(c).相位法難以準確檢測篡改邊界,而MOCC 法則難以檢測插入段的后邊界.圖3(d)則是根據(jù)MOCC 曲線得到的斜率,圖3(e)是雙機制聯(lián)合作用的結(jié)果.由圖3(f)可看到文中方法估計的插入?yún)^(qū)域為第8.91~13.71 s,與實際插入?yún)^(qū)域的起點僅有0.011 s 的誤差,而與終點有0.142 s 的誤差,明顯好于文獻[12-13].

    為了考察不同斜率門限TG對檢測結(jié)果的影響,分別取不同的幅度因子β 對該插入音頻進行實驗,結(jié)果如表1 所示.通過多個音頻實驗可發(fā)現(xiàn)一般β 取值1.1~1.5 時結(jié)果較好.作為折中,文中取1.3.

    圖3 插入篡改檢測Fig.3 Detection of insertion forgery

    表1 取不同斜率門限TG 的定位誤差Table 1 Estimated error with different slope threshold TG

    3.3 機制2 的應用推廣

    如前所述,文獻[12]無法對篡改區(qū)域精確定位.文中將機制2 推廣應用到文獻[12]中解決定位問題,用圖3 中的插入篡改為例來說明具體步驟.首先尋找相位的極值點,計算各相鄰極值點之間的斜率,然后用機制2 的準則排除虛警和干擾,對篡改區(qū)域定位,結(jié)果如圖4 所示.

    圖4 (a)為篡改音頻ENF 信號各子塊的相位及相位的極值點,圖4(b)為各相鄰極值點之間的斜率,圖4(c)為處理后的結(jié)果.圖4(d)為在音頻上對應的篡改定位,其起始和終止位置分別為第9.48 和第14.32 s,雖然存在一定的誤差,但也估計出了篡改的大致位置.

    圖4 機制2 應用到文獻[12]的例子Fig.4 An example of applying the second mechanism to themethod in[12]

    同理將機制2 推廣到文獻[13]中,結(jié)果如圖5所示.直接用文獻[13]的方法定位出的篡改區(qū)域為第8.71~14.01 s,而聯(lián)合機制2 定位的篡改區(qū)域為第8.91~13.71 s,精度明顯提高.

    圖5 機制2 應用到文獻[13]的例子Fig.5 An example of applying the second mechanism to the method in[13]

    3.4 魯棒性檢測場景一

    重采樣是音頻中常見的操作,首先考察文中方法能否抗重采樣.利用音頻編輯軟件CoolEdit 將圖3中原始采樣率為16 kHz 的音頻分別作上、下采樣處理后進行檢測,結(jié)果如圖6 所示,可發(fā)現(xiàn)其與未重采樣時的結(jié)果圖3(e)和(f)沒有明顯差異.其中圖6(a)和(b)是采樣率為44kHz 的結(jié)果,圖6(c)和(d)是采樣率為8 kHz 的結(jié)果.

    圖6 采樣率為44 kHz 和8 kHz 的檢測結(jié)果Fig.6 Detection results of the re-sampled signal with frequency 44 kHz and 8 kHz

    文獻[12]中的語音庫[15]包含100 個原始音頻文件,對其進行編輯得到130 個篡改文件,其中刪除操作和插入操作的文件各占65 個.利用文中方法對其進行篡改檢測,得到正檢率93.08%,虛警率8%.其中,正檢率是指將篡改文件正確判斷為篡改文件的概率,虛警率是指將原始文件誤判為篡改文件的概率.將語音庫中所有文件分別以44、32、11、8 和6 kHz重采樣后,再用文中方法對其進行檢測,發(fā)現(xiàn)在不同重采樣頻率下正檢率仍為93.08%,虛警率仍為8%.說明重采樣對文中方法的檢測效果無影響.

    3.5 魯棒性檢測場景二

    壓縮是另一種常見的音頻操作.仍以圖3 中的音頻為例,利用壓縮軟件Lame 3.98 將其壓縮為MP3 格式文件,采樣率仍為16 kHz,比特率為32 kb/s,文中方法檢測結(jié)果如圖7(a)、(b)所示.由于在WAV 文件壓縮為MP3 的過程中存在頻域掩蔽效應,從而損失了ENF 信號的部分信息,使得篡改定位誤差與未壓縮時比較,起始點差異不大,但終止點差異較大,說明文中方法不能完全抵御壓縮操作,但有一定的抵御能力.

    圖7 MP3 的檢測結(jié)果(采樣率為16 kHz 和8 kHz)Fig.7 Detection results of the MP3 with frequency 16 kHz and 8 kHz

    圖7(c)、(d)是將圖3 中音頻壓縮為采樣率8 kHz,比特率32 kb/s 的MP3 的檢測結(jié)果,這相當于對音頻同時進行了下采樣和壓縮操作,此時MP3 丟棄的信息更多,故對ENF 信號造成更大的損失.此時出現(xiàn)明顯誤判,說明文中方法對下采樣和較大壓縮比(例如,32 kb/s)的聯(lián)合操作魯棒性較差.

    將3.4 節(jié)所用的語音庫分別以采樣率16、8 kHz壓縮成比特率為32 kb/s 的MP3,再用文中方法進行篡改檢測.采樣率為16 kHz 時,正檢率為93.08%,虛警率為8%,故算法性能未受影響.當采樣率為8 kHz時,正檢率為62.31%,虛警率為37%,算法性能明顯下降.

    3.6 魯棒性檢測場景三

    向音頻中添加噪聲也是掩蓋篡改痕跡的常用手段.為了檢測文中算法抗噪能力,對圖3 中的音頻添加不同功率的高斯白噪聲后,再用文中方法進行檢測.圖8(a)和(b)是信噪比為30 dB 時的結(jié)果,圖8(c)和(d)是信噪比為20dB 時的結(jié)果.前者與未加噪時的結(jié)果(見圖3(e)和(f))相比,無明顯差異;而后者出現(xiàn)虛警,且對篡改區(qū)域的定位也出現(xiàn)較大誤差.不過在信噪比為20 dB 時噪聲已使人耳感到明顯的不適.一般而言,信噪比在30 dB 以上時文中方法能取得較好的檢測效果.

    圖8 信噪比為30 dB 和20 dB 時的檢測結(jié)果Fig.8 Detection results of the noisy signal with SNR 30 dB and 20 dB

    對語音庫中音頻添加高斯白噪聲后再用文中方法進行篡改檢測.信噪比為30 dB 時,正檢率為85.38%,虛警率為15%,與未加噪時的正檢率93.08%和虛警率8%相比,算法性能有所下降.當信噪比為20 dB 時,正檢率為66.92%,虛警率為33%,可見算法性能進一步下降.

    4 結(jié)論

    針對現(xiàn)有基于ENF 信號的音頻篡改盲檢測方法對插入和刪除定位精度不足的問題進行改進,提出了一種使用雙重判斷機制聯(lián)合進行篡改定位的檢測算法.為簡化算法實現(xiàn),在計算子塊間最大相關(guān)偏移量時,還提出一種無需引入額外ENF 參考信號的計算方法.與現(xiàn)有文獻中兩種代表性算法相比,文中方法的篡改定位精度明顯提高.在抗重采樣和MP3壓縮以及抗噪方面,也有一定的魯棒性.文中的篡改定位機理也可容易地推廣到以前的算法.未來將在ENF 信號的抽取質(zhì)量改善、同時存在插入和刪除的音頻篡改檢測、自適應檢測門限的設置以及算法的魯棒性等幾個方面進行進一步研究.

    [1]Daéid N N,Houck M M.Interpol's forensic science review[M].Lyon:CRC Press,2010:379-380.

    [2]Brixen E.ENF quantification of the magnetic field[C]∥AES 33rd International Conference on Audio Forensics,Theory and Practice.Denver,Colorado:AES,2008.

    [3]Grigoras C.Digital audio recording analysis the electricnetwork frequency criterion[J].International Journal of Speech Language and the Law,2005,12(1):63-76.

    [4]Grigoras C.Applications of ENF criterion in forensic audio,video,computer and telecommunication analysis[J].Forensic Science International,2007,167(2):136-145.

    [5]Cooper A J.The electric network frequency (ENF)as an aid to authenticating forensic digital audio recordings-An automated approach[C]∥AES 33rd International Conference on Audio Forensics,Theory and Practice.Denver,Colorado:AES,2008.

    [6]Huijbregtse M,Geradts Z.Using the ENF criterion for determining the time of recording of short digital audio recordings[C]∥3rd International Workshop on Computational Forensics.Berlin:Springer-Verlag,2009:116-124.

    [7]Garg R,Varna A L,Wu M.Seeing ENF:natural time stamp for digital video via optical sensing and signal processing[C]∥19th ACM International Conference on Multimedia.Scottsdale:ACM,2011:23-32.

    [8]Garg R,Varna A L,Hajj-Ahmad A,et al.“Seeing”ENF:power-signature-based timestamp for digital multimedia via optical sensing and signal processing [J].IEEE Transactions on Information Forensics and Security,2013,8(9):1417-1432.

    [9]Hajj-Ahmad A,Garg R,Wu M.Instantaneous frequency estimation and localization for ENF signals[C]∥2012 Asia-Pacific Signal & Information Processing Association Annual Summit and Conference(APSIPA ASC).Hollywood,California:IEEE,2012:1-10.

    [10]Hajj-Ahmad A,Garg R,Wu M.ENF based location classification of sensor recordings[C]∥2013 IEEE International Workshop on Information Forensics and Security(WIFS).Guangzhou,China:IEEE,2013:138-143.

    [11]Nicolalde D P,Apolinario J A.Evaluating digital audio authenticity with spectral distances and ENF phase change[C]∥IEEE International Conference on Acoustics,Speech and Signal Processing.Taipei:IEEE,2009:1417-1420.

    [12]Nicolalde D P,Apolinario J A,Biscainho L W P.Audio authenticity:detecting ENF discontinuity with high precision phase analysis[J].IEEE Transactions on Information Forensics and Security,2010,5(3):534-543.

    [13]Hu Yong-jian,Li Chang-Tsun,Lü Zhi-sheng,et al.Audio forgery detection based on max offsets for cross correlation between ENF and reference signal[C]∥11th International Workshop on Digital-Forensics and Watermarking.Shanghai:Springer,2013:253-266.

    [14]常廣,鄢素云,王毅.零相位數(shù)字濾波器在非平穩(wěn)信號處理中的應用[J].北京交通大學學報,2011,35(6):49-56.Chang Guang,Yan Su-yun,Wang Yi.Application of zero-phase digital filter on non-stationary signal processing[J].Journal of Beijing Jiao Tong University,2011,35(6):49-56.

    [15]Ortega-Garcia J,Gonzalez-Rodriguez J,Marrero-Aguiar V.Ahumada:a large speech corpus in Spanish for speaker characterization and identification[J].Speech Communication,2000,31(2):255-264.

    猜你喜歡
    子塊偏移量斜率
    基于八叉樹的地震數(shù)據(jù)多級緩存方法
    基于八叉樹的地震數(shù)據(jù)分布式存儲方法研究
    基于格網(wǎng)坐標轉(zhuǎn)換法的矢量數(shù)據(jù)脫密方法研究
    基于特征值算法的圖像Copy-Move篡改的被動取證方案
    物理圖像斜率的變化探討
    物理之友(2020年12期)2020-07-16 05:39:16
    基于波浪式矩陣置換的稀疏度均衡分塊壓縮感知算法
    攪拌針不同偏移量對6082-T6鋁合金接頭勞性能的影響
    基于最小二乘平差的全極化SAR配準偏移量估計方法
    測繪工程(2017年3期)2017-12-22 03:24:50
    求斜率型分式的取值范圍
    基于子孔徑斜率離散采樣的波前重構(gòu)
    久久国内精品自在自线图片| 成人亚洲精品av一区二区| 午夜精品国产一区二区电影 | 久久精品久久久久久噜噜老黄| 高清av免费在线| 午夜免费鲁丝| 综合色丁香网| 热re99久久精品国产66热6| 大香蕉97超碰在线| 亚洲在线观看片| 九九在线视频观看精品| 免费不卡的大黄色大毛片视频在线观看| 最近中文字幕高清免费大全6| 亚洲精品第二区| 成人一区二区视频在线观看| 日韩人妻高清精品专区| 中文精品一卡2卡3卡4更新| 亚洲国产最新在线播放| 亚洲天堂国产精品一区在线| 成年版毛片免费区| 久久久久久伊人网av| 卡戴珊不雅视频在线播放| 美女视频免费永久观看网站| 在线精品无人区一区二区三 | 18禁在线播放成人免费| 免费看日本二区| 日日摸夜夜添夜夜添av毛片| 亚洲欧美精品专区久久| 亚洲色图综合在线观看| 亚洲国产精品国产精品| 精品国产露脸久久av麻豆| 国产伦精品一区二区三区视频9| 久久97久久精品| 嘟嘟电影网在线观看| 欧美激情国产日韩精品一区| 免费观看a级毛片全部| 在线观看美女被高潮喷水网站| 久久精品国产亚洲网站| 伊人久久精品亚洲午夜| 看非洲黑人一级黄片| 黄色怎么调成土黄色| 国产视频内射| 久久鲁丝午夜福利片| 亚洲成人精品中文字幕电影| 偷拍熟女少妇极品色| 天堂俺去俺来也www色官网| 日韩一区二区视频免费看| 两个人的视频大全免费| 综合色av麻豆| 香蕉精品网在线| 日韩国内少妇激情av| 少妇人妻一区二区三区视频| a级毛片免费高清观看在线播放| 国产一区二区亚洲精品在线观看| 在线播放无遮挡| 赤兔流量卡办理| 日韩成人伦理影院| 国产成人精品福利久久| .国产精品久久| 欧美精品人与动牲交sv欧美| 黄片无遮挡物在线观看| 亚洲伊人久久精品综合| 国产一区二区亚洲精品在线观看| 天美传媒精品一区二区| 男人爽女人下面视频在线观看| 一级黄片播放器| 亚洲美女搞黄在线观看| 国产淫语在线视频| 亚洲在久久综合| 亚洲欧美精品自产自拍| 成人毛片60女人毛片免费| 免费看a级黄色片| 亚洲国产精品专区欧美| 国产精品精品国产色婷婷| 日韩免费高清中文字幕av| 国产精品久久久久久精品古装| 干丝袜人妻中文字幕| 国产淫语在线视频| 欧美xxⅹ黑人| 免费播放大片免费观看视频在线观看| a级毛片免费高清观看在线播放| 天天躁日日操中文字幕| 热99国产精品久久久久久7| 亚洲精品456在线播放app| 日韩在线高清观看一区二区三区| av网站免费在线观看视频| 中国三级夫妇交换| 亚洲国产最新在线播放| 久久精品国产亚洲av天美| 亚洲丝袜综合中文字幕| 免费少妇av软件| 免费看光身美女| 神马国产精品三级电影在线观看| 男女啪啪激烈高潮av片| 久久99蜜桃精品久久| 亚洲怡红院男人天堂| 成人特级av手机在线观看| 免费观看无遮挡的男女| 水蜜桃什么品种好| 国产久久久一区二区三区| 人妻 亚洲 视频| 白带黄色成豆腐渣| 男的添女的下面高潮视频| 18禁裸乳无遮挡动漫免费视频 | 久久精品国产鲁丝片午夜精品| 免费看不卡的av| 成人亚洲欧美一区二区av| 97人妻精品一区二区三区麻豆| 欧美激情在线99| 亚洲精品一区蜜桃| 波野结衣二区三区在线| 五月玫瑰六月丁香| 日本黄色片子视频| 国产高清不卡午夜福利| 久久久精品94久久精品| 精品99又大又爽又粗少妇毛片| 国产成人精品福利久久| 七月丁香在线播放| 熟女人妻精品中文字幕| 尾随美女入室| 在线观看一区二区三区| 国产淫片久久久久久久久| 国产午夜精品一二区理论片| 亚洲熟女精品中文字幕| 亚洲欧美清纯卡通| 联通29元200g的流量卡| 精品久久久久久久久亚洲| 国产成人a∨麻豆精品| 最近最新中文字幕大全电影3| 身体一侧抽搐| 校园人妻丝袜中文字幕| 日本一本二区三区精品| 欧美日韩在线观看h| 一个人观看的视频www高清免费观看| 涩涩av久久男人的天堂| 国产精品国产三级专区第一集| 精品少妇久久久久久888优播| 蜜桃久久精品国产亚洲av| 国产av码专区亚洲av| 有码 亚洲区| 性色av一级| 色综合色国产| 国产精品国产三级国产av玫瑰| 亚洲成人久久爱视频| 欧美+日韩+精品| 一区二区av电影网| 大香蕉97超碰在线| 天天一区二区日本电影三级| 啦啦啦在线观看免费高清www| 国产中年淑女户外野战色| 看免费成人av毛片| 自拍欧美九色日韩亚洲蝌蚪91 | 男人和女人高潮做爰伦理| 亚洲人成网站在线观看播放| 九九久久精品国产亚洲av麻豆| 国产精品嫩草影院av在线观看| 91精品一卡2卡3卡4卡| eeuss影院久久| 美女cb高潮喷水在线观看| 午夜视频国产福利| 免费播放大片免费观看视频在线观看| 精品一区二区三区视频在线| 少妇高潮的动态图| 久久久久久伊人网av| 水蜜桃什么品种好| 色网站视频免费| 国产成人精品婷婷| 国产一级毛片在线| 肉色欧美久久久久久久蜜桃 | 久久女婷五月综合色啪小说 | 中文字幕人妻熟人妻熟丝袜美| 成人毛片60女人毛片免费| 欧美一区二区亚洲| 青春草视频在线免费观看| 国产精品偷伦视频观看了| 欧美zozozo另类| 一级爰片在线观看| 欧美日韩精品成人综合77777| 18禁裸乳无遮挡免费网站照片| 国产伦精品一区二区三区四那| 国产 一区精品| 一本久久精品| 免费看光身美女| 亚洲精品色激情综合| 久久精品国产亚洲网站| 大话2 男鬼变身卡| 欧美成人a在线观看| 日韩一区二区视频免费看| 欧美日韩精品成人综合77777| 最近手机中文字幕大全| 亚洲欧美日韩东京热| 亚洲精品色激情综合| 久久精品久久久久久久性| 国产有黄有色有爽视频| 久久99热这里只频精品6学生| 国产精品麻豆人妻色哟哟久久| 久久99精品国语久久久| 在线 av 中文字幕| 成人免费观看视频高清| 日本猛色少妇xxxxx猛交久久| 国产一级毛片在线| 国产永久视频网站| 18禁裸乳无遮挡免费网站照片| 亚洲精品一二三| 亚洲怡红院男人天堂| 欧美xxxx性猛交bbbb| 成人免费观看视频高清| 久久久久久久午夜电影| 亚洲美女搞黄在线观看| 深夜a级毛片| 久久国内精品自在自线图片| 成人亚洲精品一区在线观看 | 在线观看免费高清a一片| 久久精品熟女亚洲av麻豆精品| 精华霜和精华液先用哪个| 亚洲精品日本国产第一区| 特级一级黄色大片| 啦啦啦在线观看免费高清www| 极品教师在线视频| 久久久午夜欧美精品| 嫩草影院入口| 特级一级黄色大片| 啦啦啦在线观看免费高清www| 成人亚洲精品av一区二区| 欧美精品人与动牲交sv欧美| 91精品伊人久久大香线蕉| 日韩欧美一区视频在线观看 | 亚洲自拍偷在线| 久久国内精品自在自线图片| 亚洲国产日韩一区二区| 在线观看国产h片| 真实男女啪啪啪动态图| 精品国产三级普通话版| 免费看日本二区| 1000部很黄的大片| 国产在视频线精品| 三级经典国产精品| 观看美女的网站| 国产日韩欧美在线精品| 99热网站在线观看| 中文字幕亚洲精品专区| 亚洲自拍偷在线| 成年人午夜在线观看视频| 国产精品一区二区三区四区免费观看| 搡老乐熟女国产| 日本一本二区三区精品| 少妇人妻久久综合中文| 99热这里只有精品一区| 亚洲av中文av极速乱| 亚洲va在线va天堂va国产| 欧美+日韩+精品| 国产男女内射视频| 久久精品国产a三级三级三级| 国产亚洲5aaaaa淫片| 禁无遮挡网站| 国产高清不卡午夜福利| 视频中文字幕在线观看| 熟女av电影| 精品午夜福利在线看| 熟妇人妻不卡中文字幕| 日产精品乱码卡一卡2卡三| 毛片女人毛片| 成人免费观看视频高清| 亚洲,一卡二卡三卡| 国产成人午夜福利电影在线观看| 2018国产大陆天天弄谢| 国产日韩欧美在线精品| 国产精品国产三级国产av玫瑰| 欧美日韩亚洲高清精品| 久久国产乱子免费精品| 国产精品国产三级国产专区5o| 日本黄大片高清| 色播亚洲综合网| 色5月婷婷丁香| 国产成人福利小说| 日本一本二区三区精品| 色吧在线观看| 精品亚洲乱码少妇综合久久| 亚洲av成人精品一二三区| 国产精品国产三级国产专区5o| 国产精品一区二区性色av| 久久午夜福利片| 亚洲欧美日韩东京热| 噜噜噜噜噜久久久久久91| 黄色欧美视频在线观看| 国产一区亚洲一区在线观看| 真实男女啪啪啪动态图| 成人国产av品久久久| 五月伊人婷婷丁香| 美女高潮的动态| 日韩人妻高清精品专区| 国产亚洲av嫩草精品影院| 白带黄色成豆腐渣| 嫩草影院精品99| 亚洲国产精品国产精品| 亚洲经典国产精华液单| 高清毛片免费看| 黄色一级大片看看| 九九久久精品国产亚洲av麻豆| 婷婷色综合www| 国产又色又爽无遮挡免| 国产精品成人在线| 国精品久久久久久国模美| 极品教师在线视频| 欧美日韩综合久久久久久| 综合色av麻豆| 国产高清三级在线| 22中文网久久字幕| 欧美xxxx黑人xx丫x性爽| 国产高清三级在线| 久久久亚洲精品成人影院| 亚洲精品国产色婷婷电影| 男人狂女人下面高潮的视频| 九草在线视频观看| 亚洲综合色惰| 国产成人aa在线观看| 亚洲在久久综合| 我的女老师完整版在线观看| 51国产日韩欧美| 岛国毛片在线播放| av又黄又爽大尺度在线免费看| 99热全是精品| 波野结衣二区三区在线| 国产亚洲午夜精品一区二区久久 | 欧美xxⅹ黑人| 在线观看美女被高潮喷水网站| 亚洲国产精品国产精品| 亚洲精品日韩在线中文字幕| 午夜福利视频精品| 午夜免费观看性视频| 少妇丰满av| 国产淫片久久久久久久久| 黄色配什么色好看| 久久影院123| 小蜜桃在线观看免费完整版高清| 久久人人爽人人爽人人片va| 一级爰片在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 身体一侧抽搐| 男插女下体视频免费在线播放| 三级男女做爰猛烈吃奶摸视频| 男男h啪啪无遮挡| 黄色一级大片看看| 久久97久久精品| 天天一区二区日本电影三级| 国产精品女同一区二区软件| 人人妻人人看人人澡| 麻豆国产97在线/欧美| 午夜精品一区二区三区免费看| 爱豆传媒免费全集在线观看| 又爽又黄a免费视频| 亚洲美女视频黄频| 中文字幕免费在线视频6| 黄片wwwwww| 国产v大片淫在线免费观看| 少妇人妻一区二区三区视频| 超碰97精品在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲精品第二区| 亚洲av国产av综合av卡| 97精品久久久久久久久久精品| 国产爱豆传媒在线观看| 丝袜脚勾引网站| 久久久久久九九精品二区国产| 下体分泌物呈黄色| 亚洲精品aⅴ在线观看| 久久久久网色| 人妻夜夜爽99麻豆av| 精品一区二区三区视频在线| 亚洲精品aⅴ在线观看| 久久久久精品性色| 国内少妇人妻偷人精品xxx网站| 精品人妻熟女av久视频| 精品人妻偷拍中文字幕| 欧美xxxx性猛交bbbb| 亚洲精品自拍成人| 欧美xxxx性猛交bbbb| 亚洲一区二区三区欧美精品 | 久久精品久久久久久噜噜老黄| 午夜精品国产一区二区电影 | 欧美3d第一页| 国产 一区精品| 特大巨黑吊av在线直播| 男女下面进入的视频免费午夜| 欧美成人午夜免费资源| 最新中文字幕久久久久| 亚洲精品乱码久久久v下载方式| 99热国产这里只有精品6| 内地一区二区视频在线| 国产精品秋霞免费鲁丝片| 国产乱人偷精品视频| 日日撸夜夜添| 又黄又爽又刺激的免费视频.| 亚州av有码| 大又大粗又爽又黄少妇毛片口| 午夜激情久久久久久久| 久久这里有精品视频免费| 国产一区二区亚洲精品在线观看| 赤兔流量卡办理| 欧美性猛交╳xxx乱大交人| 91精品一卡2卡3卡4卡| 在线观看美女被高潮喷水网站| 日本熟妇午夜| 亚洲国产精品成人综合色| 欧美最新免费一区二区三区| 亚洲va在线va天堂va国产| 国产午夜福利久久久久久| 国产淫片久久久久久久久| 久久久亚洲精品成人影院| 久久精品国产亚洲网站| 91aial.com中文字幕在线观看| 欧美 日韩 精品 国产| 精品国产三级普通话版| 新久久久久国产一级毛片| 亚洲精品乱久久久久久| 久久精品国产亚洲av天美| 亚洲精品日韩av片在线观看| 国产美女午夜福利| 日韩国内少妇激情av| 国产精品av视频在线免费观看| 欧美丝袜亚洲另类| 欧美xxⅹ黑人| 国产免费视频播放在线视频| 久久久久久久久久成人| 大片电影免费在线观看免费| 老师上课跳d突然被开到最大视频| 久久精品国产自在天天线| 久久久久性生活片| 成人毛片a级毛片在线播放| 国产综合精华液| 欧美一级a爱片免费观看看| 少妇被粗大猛烈的视频| 97精品久久久久久久久久精品| 日本免费在线观看一区| 女人十人毛片免费观看3o分钟| 91久久精品国产一区二区三区| 亚洲天堂国产精品一区在线| eeuss影院久久| 建设人人有责人人尽责人人享有的 | 视频区图区小说| 女的被弄到高潮叫床怎么办| 日产精品乱码卡一卡2卡三| 男人爽女人下面视频在线观看| 精品国产露脸久久av麻豆| 一个人观看的视频www高清免费观看| 在线精品无人区一区二区三 | 一级毛片我不卡| 亚洲av二区三区四区| 色视频在线一区二区三区| 一级毛片 在线播放| 免费看日本二区| 国产精品av视频在线免费观看| 男女那种视频在线观看| 国产视频首页在线观看| 久久久久国产网址| 日本一二三区视频观看| 亚洲精品乱码久久久久久按摩| 欧美丝袜亚洲另类| 又大又黄又爽视频免费| 日本一本二区三区精品| 久久久精品欧美日韩精品| 亚洲欧美日韩东京热| 亚洲国产最新在线播放| 天天躁夜夜躁狠狠久久av| 女人被狂操c到高潮| 成人国产麻豆网| 成人美女网站在线观看视频| 嫩草影院精品99| 久久亚洲国产成人精品v| 亚洲精品国产色婷婷电影| 青青草视频在线视频观看| 丝袜喷水一区| 亚洲美女视频黄频| 国产精品久久久久久久久免| 日产精品乱码卡一卡2卡三| 国产精品熟女久久久久浪| 国产av码专区亚洲av| 在线观看免费高清a一片| 国产成人freesex在线| 欧美一区二区亚洲| 中文资源天堂在线| 日韩强制内射视频| 国产伦精品一区二区三区四那| 国产成人精品婷婷| 天堂网av新在线| 亚洲国产欧美在线一区| 男女那种视频在线观看| 99久久中文字幕三级久久日本| 国产日韩欧美亚洲二区| www.色视频.com| 久久精品国产a三级三级三级| 亚洲欧洲国产日韩| 内地一区二区视频在线| 少妇丰满av| 亚洲精品乱码久久久v下载方式| 精品国产三级普通话版| 亚洲人成网站在线播| 日韩伦理黄色片| 免费黄网站久久成人精品| 男女下面进入的视频免费午夜| 老女人水多毛片| 国产亚洲最大av| 伦理电影大哥的女人| 一本久久精品| 日本与韩国留学比较| 极品少妇高潮喷水抽搐| 久久影院123| 成人亚洲欧美一区二区av| 欧美xxxx性猛交bbbb| 日韩精品有码人妻一区| 国产探花在线观看一区二区| 综合色丁香网| 大话2 男鬼变身卡| 极品教师在线视频| 高清av免费在线| 日韩,欧美,国产一区二区三区| 国内精品宾馆在线| av免费观看日本| 国产色婷婷99| 亚洲,一卡二卡三卡| 成人毛片60女人毛片免费| 久热久热在线精品观看| 一级片'在线观看视频| 天天一区二区日本电影三级| 国产老妇伦熟女老妇高清| 亚洲成人中文字幕在线播放| 国产色婷婷99| 国产av国产精品国产| 三级国产精品欧美在线观看| 亚洲精品456在线播放app| xxx大片免费视频| 免费看光身美女| 中文乱码字字幕精品一区二区三区| 欧美日韩视频精品一区| 新久久久久国产一级毛片| 亚洲av电影在线观看一区二区三区 | 国产伦精品一区二区三区视频9| 久久久精品免费免费高清| 日日摸夜夜添夜夜爱| 久久久久久九九精品二区国产| 亚洲久久久久久中文字幕| 九色成人免费人妻av| 亚洲最大成人中文| 丝袜脚勾引网站| 精品国产三级普通话版| 国产成人午夜福利电影在线观看| 午夜老司机福利剧场| 国产亚洲5aaaaa淫片| 真实男女啪啪啪动态图| 亚洲精品久久午夜乱码| 永久免费av网站大全| 亚洲人成网站在线播| 久久久精品94久久精品| 欧美区成人在线视频| 啦啦啦啦在线视频资源| 亚洲av福利一区| 大香蕉久久网| 国产中年淑女户外野战色| 七月丁香在线播放| 亚洲,欧美,日韩| 热99国产精品久久久久久7| 日韩中字成人| 久久这里有精品视频免费| 亚洲婷婷狠狠爱综合网| 2022亚洲国产成人精品| 国产综合懂色| 国产黄片视频在线免费观看| av网站免费在线观看视频| 国产久久久一区二区三区| 舔av片在线| 午夜福利在线观看免费完整高清在| 晚上一个人看的免费电影| 亚洲国产欧美人成| 人体艺术视频欧美日本| 女人被狂操c到高潮| 日韩av免费高清视频| 网址你懂的国产日韩在线| 日韩一区二区三区影片| 亚洲精品色激情综合| 岛国毛片在线播放| 秋霞在线观看毛片| 亚洲精品色激情综合| 国产毛片a区久久久久| 亚洲精品第二区| 爱豆传媒免费全集在线观看| 岛国毛片在线播放| 韩国av在线不卡| 亚洲av免费在线观看| 久久久久国产网址| 亚洲精品第二区| 亚洲精品乱久久久久久| 国产成人精品福利久久| 免费在线观看成人毛片| 国产一区亚洲一区在线观看| 成年版毛片免费区| 天堂俺去俺来也www色官网| 欧美精品国产亚洲| 亚洲成色77777| 成人鲁丝片一二三区免费| 有码 亚洲区| 九九在线视频观看精品| 18禁动态无遮挡网站| 91久久精品国产一区二区三区| 免费观看av网站的网址| 欧美日韩亚洲高清精品| 亚洲最大成人av| 国产精品一区二区性色av| 99久久人妻综合| 综合色丁香网| 中文字幕av成人在线电影| 男女那种视频在线观看| 午夜福利视频精品| 91aial.com中文字幕在线观看| 久热这里只有精品99| 高清欧美精品videossex|