磁共振成像的圖像重建方法及圖像重建裝置的制作方法

            文檔序號:6444608閱讀:200來源:國知局
            專利名稱:磁共振成像的圖像重建方法及圖像重建裝置的制作方法
            技術領域
            本發明涉及磁共振成像技術領域,尤其涉及一種磁共振成像的圖像重建方法及圖
            像重建裝置。
            背景技術
            磁共振成像是利用磁共振現象進行成像的一種技術。磁共振現象的原理主要包括:包含單數質子的原子核,例如人體內廣泛存在的氫原子核,其質子具有自旋運動,猶如一個小磁體,并且這些小磁體的自旋軸無一定的規律,如果施加外在磁場,這些小磁體將按外在磁場的磁力線重新排列,具體為在平行于或反平行于外在磁場磁力線的兩個方向排列,將上述平行于外在磁場磁力線的方向稱為正縱向軸,將上述反平行于外在磁場磁力線的方向稱為負縱向軸,原子核只具有縱向磁化分量,該縱向磁化分量既具有方向又具有幅度。用特定頻率的射頻(RF, Radio Frequency)脈沖激發處于外在磁場中的原子核,使這些原子核的自旋軸偏離正縱向軸或負縱向軸,產生共振,這就是磁共振現象。上述被激發原子核的自旋軸偏離正縱向軸或負縱向軸之后,原子核具有了橫向磁化分量。停止發射射頻脈沖后,被激發的原子核發射回波信號,將吸收的能量逐步以電磁波的形式釋放出來,其相位和能級都恢復到激發前的狀態,將原子核發射的回波信號經過空間編碼等進一步處理即可重建圖像。現有的一種磁共振成像過程大致包括:首先將成像區域對應的厚層在片層(slice)方向劃分為多個子厚層(sub-slab, SS),每個子厚層包括多個片層。如圖1所示,圖1中以劃分8個子厚層的情況為例。之后對每個子厚層,利用3D快速自旋回波進行成像掃描,得到每個子厚層對應的圖像數據,對各個子厚層對應的圖像數據進行傅里葉變換重建出各個子厚層的臨時圖像,在圖像域對各個子厚層的臨時圖像進行合并,得到成像區域對應的磁共振圖像。其中,利用快速自旋回波進行成像掃描的具體方法可包括:在一個重復時間TR內發射不同角度的脈沖序列,同時在一定片層編碼梯度下改變相位編碼梯度,以便充滿一個片層編碼k空間;在另外一個重復時間TR內,保持射頻脈沖不變,改變片層編碼梯度,并在該片層編碼梯度下改變相位編碼梯度,得到另外一個片層編碼的k空間;依次類推;直到采集完整個k空間數據。在成像掃描過程中,為了克服片層變形,在相位編碼時,通常需要在子厚層片層方向的兩側按照事先確定的擴展因子進行擴展,以得到大于激發厚度(Excited Thickness)的片層編碼厚度(Slice Encoded Thickness),并對該編碼厚度對應的片層厚度進行編碼。如圖2所示,圖2的中間填充區域對應的層厚TH為當前子厚層的激發厚度,圖2的整個區域對應的層厚STH為當前子厚層的片層編碼厚度。其中,激發厚度通常等于所對應的子厚層的厚度。在圖像域對各個子厚層的臨時圖像進行合并時,通常是對成像區域的每個位置,將同一位置所對應的各子厚層的片層進行合并。如圖3所示,圖3為對成像區域S位置所對應的各子厚層的片層進行合并的示意圖。其中,對每個子厚層來說,除了激發厚度對應的片層(圖3中所示S位置對應的子厚層2的片層)主要包含圖像信息以外,其它擴展出來的片層中除個別片層中包含部分圖像信息外,大部分都是包含的噪聲。現有技術中,直接對成像區域同一位置對應的各子厚層的片層進行簡單疊加時,會對噪聲進行累加,從而降低磁共振圖像的信噪比(SNR, Signal-to-Noise Ratio),影響成像質量。

            發明內容
            有鑒于此,本發明一方面提出了一種磁共振成像的圖像重建方法,另一方面提出了一種磁共振成像的圖像重建裝置,用以提高磁共振圖像的信噪比,從而提高圖像質量。本發明提出的磁共振成像的圖像重建方法,包括:A、利用平均濾波器對圖像域的各個子厚層的片層進行濾波,得到對應各片層的蒙板;B、將各蒙板的像素值與一預先確定的噪聲閾值進行比較,對于小于所述噪聲閾值的像素值,將該蒙板所對應片層的相應位置像素值按照預設標準降低;C、對成像區域同一位置對應的各個子厚層的片層進行合并,得到成像區域的磁共振圖像。根據一個實施例,該方法進一步包括:al、分別從每個圖像域的子厚層中選取出至少一個邊緣片層;bl、分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域;Cl、計算每個矩形區域的灰度平均值及標準差,并將每個矩形區域的灰度平均值及標準差對應存儲為一個數據對;dl、將所述數據對按照灰度平均值的大小進行升序或降序的排序;el、將包含最小灰度平均值的數據對作為第一有效數據對,并執行步驟fl ;或者,將包含最小灰度平均值的數據對作為檢索數據對,計算該檢索數據對與相鄰的設定數目的數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述檢索數據對作為第一有效數據對,并執行步驟Π,否則,將相鄰的下一個數據對作為檢索數據對,并返回執行本步驟的差值計算過程;fl、將所述第一有效數據對作為當前有效數據對;gl、計算當前有效數據對與相鄰的下一個數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述相鄰的下一個數據對作為當前有效數據對,并返回執行本步驟;否則,執行下一步驟;h1、將所有有效數據對的灰度平均值的平均值作為所述噪聲閾值,或將所有有效數據對的灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。根據另一實施例,該方法進一步包括:a2、分別從每個圖像域的子厚層中選取出至少一個邊緣片層;b2、分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域;c2、計算每個矩形區域的灰度平均值;d2、將所述灰度平均值進行升序或降序的排序;e2、將最小灰度平均值作為第一有效灰度平均值,并執行步驟f2 ;或者,將最小灰度平均值作為檢索灰度平均值,計算該檢索灰度平均值與相鄰的設定數目的灰度平均值之間的差值,在所述差值小于第一設定閾值時,將所述檢索灰度平均值作為第一有效灰度平均值,并執行步驟f2,否則,將相鄰的下一個灰度平均值作為檢索灰度平均值,并返回執行本步驟的差值計算過程;f2、將所述第一有效灰度平均值作為當前有效灰度平均值;g2、計算當前有效灰度平均值與相鄰的下一個灰度平均值之間的差值,在所述差值小于所述第一設定閾值時,將所述相鄰的下一個灰度平均值作為當前有效灰度平均值,并返回執行本步驟;否則,執行下一步驟;h2、將所有有效灰度平均值的平均值作為所述噪聲閾值,或將所有有效灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。根據再一實施例,該方法進一步包括:a3、分別從每個圖像域的子厚層中選取出至少一個邊緣片層;b3、分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域;c3、計算每個矩形區域的標準差;d3、將所述標準差進行升序或降序的排序;e3、將最小標準差作為第一有效標準差,并執行步驟f3 ;或者,將最小標準差作為檢索標準差,計算該檢索標準差與相鄰的設定數目的標準差之間的差值,在所述差值小于第二設定閾值時,將所述檢索標準差作為第一有效標準差,并執行步驟f3,否則,將相鄰的下一個標準差作為檢索標準差,并返回執行本步驟的差值計算過程;f3、將所述第一有效標準差作為當前有效標準差;g3、計算當前有效標準差與相鄰的下一個標準差之間的差值,在所述差值小于所述第二設定閾值時,將所述相鄰的下一個標準差作為當前有效標準差,并返回執行本步驟;否則,執行下一步驟;h3、計算所有有效標準差對應的矩形區域的灰度平均值,將所述灰度平均值的平均值作為所述噪聲閾值,或將所述灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。可選地,所述將該蒙板所對應片層的相應位置像素值按照預設標準降低包括:將該蒙板所對應片層的相應位置像素值設置為O ;或者,將該蒙板所對應片層的相應位置像素值乘以小于I的系數。可選地,利用像素相加法或像素平方和法,對成像區域同一位置對應的各個子厚層的片層進行合并。本發明提出的磁共振成像的圖像重建裝置,包括:一個蒙板生成模塊,用于利用平均濾波器對圖像域的各個子厚層的片層進行濾波,得到對應各片層的蒙板;—個噪聲削弱模塊,用于將各蒙板的像素值與一預先確定的噪聲閾值進行比較,對于小于所述噪聲閾值的像素值,將該蒙板所對應片層的相應位置像素值按照預設標準降低;一個圖像合并模塊,用于對成像區域同一位置對應的各個子厚層的片層進行合并,得到成像區域的磁共振圖像。
            根據一個實施例,該裝置進一步包括:一個第一片層選取模塊,用于對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層;一個第一區域選取模塊,用于分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域;一個第一數據計算模塊,用于計算每個矩形區域的灰度平均值及標準差,并將每個矩形區域的灰度平均值及標準差對應存儲為一個數據對;一個第一排序模塊,用于將所述數據對按照灰度平均值的大小進行升序或降序的排序;一個第一初始值確定模塊,用于將包含最小灰度平均值的數據對作為第一有效數據對;或者從包含最小灰度平均值的數據對開始,將當前數據對作為檢索數據對,計算該檢索數據對與相鄰的設定數目的數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述檢索數據對作為第一有效數據對;一個第一有效值確定模塊,用于從所述第一有效數據對開始,將該第一有效數據對作為當前有效數據對;計算當前有效數據對與相鄰的下一個數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述相鄰的下一個數據對作為當前有效數據對,并繼續查找下一個有效數據對,否則停止有效數據對的查找;一個第一噪聲閾值確定模塊,用于將所有有效數據對的灰度平均值的平均值作為所述噪聲閾值,或將所有有效數據對的灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。根據另一個實施例,該裝置進一步包括:—個第二片層選取模塊,用于對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層;一個第二區域選取模塊,用于分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域;一個第二數據計算模塊,用于計算每個矩形區域的灰度平均值;一個第二排序模塊,用于將所述灰度平均值進行升序或降序的排序;一個第二初始值確定模塊,用于將最小灰度平均值作為第一有效灰度平均值;或者從最小灰度平均值開始,將當前灰度平均值作為檢索灰度平均值,計算該檢索灰度平均值與相鄰的設定數目的灰度平均值之間的差值,在所述差值小于第一設定閾值時,將所述檢索灰度平均值作為第一有效灰度平均值;—個第二有效值確定模塊,用于從所述第一有效灰度平均值開始,將該第一有效灰度平均值作為當前有效灰度平均值;計算當前有效灰度平均值與相鄰的下一個灰度平均值之間的差值,在所述差值小于所述第一設定閾值時,將所述相鄰的下一個灰度平均值作為當前有效灰度平均值,并繼續查找下一個有效灰度平均值,否則停止有效灰度平均值的查找;一個第二噪聲閾值確定模塊,用于將所有有效灰度平均值的平均值作為所述噪聲閾值,或將所有有效灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。根據再一實施例,該裝置進一步包括:
            一個第三片層選取模塊,用于對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層;一個第三區域選取模塊,用于分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域;一個第三數據計算模塊,用于計算每個矩形區域的標準差;一個第三排序模塊,用于將所述標準差進行升序或降序的排序;一個第三初始值確定模塊,用于將最小標準差作為第一有效標準差;或者從最小標準差開始,將當前標準差作為檢索標準差,計算該檢索標準差與相鄰的設定數目的標準差之間的差值,在所述差值小于第二設定閾值時,將所述檢索標準差作為第一有效標準差;一個第三有效值確定模塊,用于從所述第一有效標準差開始,將該第一有效標準差作為當前有效標準差;計算當前有效標準差與相鄰的下一個標準差之間的差值,在所述差值小于所述第二設定閾值時,將所述相鄰的下一個標準差作為當前有效標準差,并繼續查找下一個有效標準差,否則停止有效標準差的查找;一個第三噪聲閾值確定模塊,用于計算所有有效標準差對應的矩形區域的灰度平均值,將所述灰度平均值的平均值作為所述噪聲閾值,或將所述灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。從上述方案中可以看出,本發明實施例中首先獲取每個片層的蒙板,將該蒙板中的像素值逐一與預先確定的噪聲閾值進行比較,以確定對應片層中是否存在純噪聲,在確定其存在純噪聲時,降低該噪聲對應的像素值,之后再對降低了噪聲像素值的片層進行合并,從而削弱了噪聲的影響。進一步地,噪聲閾值可根據每個子厚層中最不可能存在圖像信息的區域的灰度平均值來確定,使得對于不同的成像條件有針對其的具體噪聲閾值,從而噪聲閾值更有針對性,也更準確。此外,將噪聲所對應的像素值降低時,可直接設為0,也可將其乘以一個小于I的系數,方式靈活多樣,可根據實際需要確定。


            下面將通過參照附圖詳細描述本發明的優選實施例,使本領域的普通技術人員更清楚本發明的上述及其它特征和優點,附圖中:圖1為子厚層劃分示意圖。圖2為每個子厚層的激發層厚與擴展層厚的關系示意圖。圖3為對成像區域S位置所對應的各子厚層的片層進行合并的示意圖。圖4為本發明實施例中磁共振成像的圖像重建方法的示例性流程圖。圖5為本發明實施例中噪聲閾值的一種確定方法的流程不意圖。圖6為本發明實施例中噪聲閾值的又一種確定方法的流程不意圖。圖7為本發明實施例中噪聲閾值的再一種確定方法的流程示意圖。圖8a和圖8b為本發明實施例中從圖像域的子厚層中選取噪聲區域的方法示意圖;其中,圖8a為從圖像域的子厚層中選取片層的示意圖;圖Sb為從片層中選取矩形區域的示意圖。圖9為本發明實施例中磁共振成像的圖像重建裝置的一個結構示意圖。圖10為本發明實施例中磁共振成像的圖像重建裝置的又一個結構示意圖。圖11為本發明實施例中磁共振成像的圖像重建裝置的再一個結構示意圖。圖12為本發明實施例中磁共振成像的圖像重建裝置的再一個結構示意圖。圖中,附圖標記如下:SD-片層方向TH-激發厚度STH-編碼厚度SSl-子厚層I SS2-子厚層2 SS3-子厚層3 SS4-子厚層4401-生成對應每個片層的蒙板 402-利用蒙板中的像素值與噪聲閾值的比較結果,確定并調低相應片層中的噪聲像素值403-對降低了噪聲像素值的片層進行合并501-從子厚層中選擇邊緣片層 502-從邊緣片層中選擇矩形區域503-計算矩形區域的由灰度平均值及標準差構成的數據對 504-將數據對按照灰度平均值的大小進行排序505-從數據對中確定第一個有效數據對506-從數據對中確定所有有效數據對507-根據所有有效數據對中的灰度平均值確定噪聲閾值601-從子厚層中選擇邊緣片層 602-從邊緣片層中選擇矩形區域 603-計算矩形區域的灰度平均值604-將灰度平均值進行排序605-確定第一個有效灰度平均值606-確定所有有效灰度平均值607-根據所有有效灰度平均值確定噪聲閾值701-從子厚層中選擇邊緣片層702-從邊緣片層中選擇矩形區域703-計算矩形區域的標準差704-將標準差進行排序705-確定第一個有效標準差706-確定所有有效標準差707-根據所有有效標準差對應的矩形區域的灰度平均值確定噪聲閾值901-蒙板生成模塊902-噪聲削弱模塊903-圖像合并模塊1001-第一片層選取模塊1002-第一區域選取模塊1003-第一數據計算模塊1004-第一排序模塊1005-第一初始值確定模塊1006-第一有效值確定模塊1007-第一噪聲閾值確定模塊1101-第二片層選取模塊1102-第二區域選取模塊1103-第二數據計算模塊1104-第二排序模塊1105-第二初始值確定模塊1106-第二有效值確定模塊1107-第
            二噪聲閾值確定模塊1201-第三片層選取模塊1202-第三區域選取模塊1203-第三數據計算模塊1204-第三排序模塊1205-第三初始值確定模塊1206-第三有效值確定模塊1207-第三噪聲閾值確定模塊
            具體實施例方式本發明實施例中,為了降低噪聲的累積效應,削弱各片層中純噪聲的影響,考慮在對各片層進行合并之前,降低各片層中的噪聲像素值。為此,本發明實施例中,首先確定一個用于標志噪聲水平的噪聲閾值,之后利用該噪聲閾值檢測各片層中的噪聲像素值,進而降低該噪聲像素值的值,最后再對降低了噪聲像素值的片層進行合并。為使本發明的目的、技術方案和優點更加清楚,以下舉實施例對本發明進一步詳細說明。 圖4為本發明實施例中磁共振成像的圖像重建方法的示例性流程圖。如圖4所示,該流程包括如下步驟:步驟401,對圖像域的各個子厚層,利用平均濾波器對各子厚層的片層進行濾波,得到對應各片層的蒙板。步驟402,對每個蒙板,將該蒙板的像素灰度值逐一與預先確定的噪聲閾值進行比較,對于小于所述噪聲閾值的像素灰度值,將該蒙板所對應片層的相應位置像素值按照預設標準降低。本實施例中,噪聲閾值用以標志噪聲水平,低于該噪聲閾值的像素值可認為是噪聲像素值,進而可降低其取值;高于該噪聲閾值的像素值可認為是有用的圖像像素值,可保持其取值不變。本步驟中,主要是根據蒙板中的像素值(即像素灰度值)與噪聲閾值的比較結果,確定相應片層中的噪聲像素值,并調低該噪聲像素值。其中,噪聲閾值可以根據經驗值或仿真值確定,也可以根據每個子厚層中最不可能存在圖像信息的區域的灰度平均值來確定。其中,對于根據每個子厚層中最不可能存在圖像信息的區域的灰度平均值來確定噪聲閾值的方法將在后面的圖5至圖7中詳細介紹。其中,將該蒙板所對應片層的相應位置像素值按照預設標準降低的方式可有很多。例如,可將該蒙板所對應片層的相應位置像素值設置為O ;或者也可將該蒙板所對應片層的相應位置像素值乘以小于I的系數,例如,0.1,0.2,0.3,0.4,0.5等。步驟403,對成像區域的每個位置,將同一位置對應的各個子厚層的片層進行合并,得到成像區域的磁共振圖像。本步驟中,是對降低了噪聲像素值的片層進行的合并。其中,對片層進行合并的方法可有多種,例如,可以是像素相加法,也可以是像素平方和法等。由于本實施例中是對降低了噪聲像素值的片層進行的合并,因此削弱了噪聲累積的影響,提高了磁共振圖像的信噪比,進而提高了圖像質量。對于步驟402中提到的噪聲閾值,在根據每個子厚層中最不可能存在圖像信息的區域的灰度平均值來確定時,可有多種實現方式。本實施例中,每個子厚層中最不可能存在圖像信息的區域可以認為是其邊緣片層的四個邊角的區域,當然也可以根據具體情況選擇其他區域。圖5至圖7中給出了三種噪聲閾值的確定過程。圖5為噪聲閾值的一種確定方法的流程示意圖。如圖5所示,該方法包括如下步驟:步驟501,對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層。本實施例中,可從每個子厚層中選取片層方向最外側的兩個片層,如圖8a所示,片層方向的上下兩個邊緣片層。實際應用中,也可只選取其中一個片層;或者,也可選取上面兩個片層和/或下面兩個片層等。步驟502,分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域。本實施例中,可如圖Sb所示,從每個邊緣片層的四個角中分別選取出一個設定大小的矩形區域。當然也可以從上述四個角的任一個或任兩個或任三個角中分別選取出一個設定大小的矩形區域。并且該矩形區域的大小可根據經驗值或仿真值預先設定,也可根據其它方法確定。
            步驟503,計算每個矩形區域的灰度平均值及標準差,并將每個矩形區域的灰度平均值及標準差對應存儲為一個數據對。具體實現時,可將各矩形區域的灰度平均值及標準差分別存儲在兩個隊列中,SP將各矩形區域的灰度平均值存儲在一個隊列中,將各矩形區域的標準差存儲在另一個隊列中,但二者的位置存在一一對應的關系。或者,也可將各矩形區域的灰度平均值及標準差對應存儲在一張表格中等。具體實現方式可根據實際需要確定,此處不對其進行限制。步驟504,將所述數據對按照灰度平均值的大小進行升序或降序的排序。步驟505,從包含最小灰度平均值的數據對開始,將該數據對作為檢索數據對,計算該檢索數據對與相鄰的設定數目的數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述檢索數據對作為第一有效數據對;否則,將相鄰的下一個數據對作為檢索數據對,并重復執行上述求差及比較過程,直到找到第一有效數據對。本步驟中,有效數據對指的是能表示純噪聲矩形區域的數據對。設定數目可以是根據經驗值或仿真值等確定的數值,如可以是相鄰的I個、2個、3個或4個等。本步驟中,灰度平均值之差對應的設定閾值和標準差之差對應的設定閾值均可根據經驗值或仿真值確定,也可根據其它方法確定。步驟506,從所述第一有效數據對開始,將該第一有效數據對作為當前有效數據對,計算當前有效數據對與相鄰的下一個數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述相鄰的下一個數據對作為當前有效數據對,并重復執行上述求差及比較過程,直到找出所有有效數據對。步驟507,將所有有效數據對的灰度平均值的平均值作為噪聲閾值,或將所有有效數據對的灰度平均值的平均值乘以預設系數后的值作為噪聲閾值。本步驟中,該預設系數可根據經驗值或仿真值確定,也可根據其它方法確定。實際應用中,也可不采用上述步驟505中描述的過程,而是直接將包含最小灰度平均值的數據對作為第一個有效數據對。圖5中所示的噪聲閾值確定方法,在確定表示純噪聲的有效數值時,是采用的灰度平均值和標準差的雙項判斷標準。實際應用中,也可采用二者之一進行單項的判斷標準。如圖6和圖7所示。圖6為噪聲閾值的又一種確定方法的流程示意圖。如圖6所示,該方法包括如下步驟:步驟601 步驟602的具體實現過程與步驟501 步驟502的具體實現過程一致。步驟603,計算每個矩形區域的灰度平均值。步驟604,將所述灰度平均值進行升序或降序的排序。步驟605,從最小灰度平均值開始,將其作為檢索灰度平均值,計算該檢索灰度平均值與相鄰的設定數目的灰度平均值之間的差值,在所述差值小于第一設定閾值時,將所述檢索灰度平均值作為第一有效灰度平均值;否則,將相鄰的下一個灰度平均值作為檢索灰度平均值,并重復執行上述求差及比較過程,直到找到第一有效灰度平均值。本步驟中,第一有效灰度平均值指定是能表示純噪聲矩形區域的灰度平均值。設定數目可以是根據經驗值或仿真值等確定的數值,如可以是相鄰的I個、2個、3個或4個坐寸O本步驟中,第一設定閾值可根據經驗值或仿真值確定,也可根據其它方法確定。步驟606,從所述第一有效灰度平均值開始,將該第一有效灰度平均值作為當前有效灰度平均值,計算當前有效灰度平均值與相鄰的下一個灰度平均值之間的差值,在所述差值小于所述第一設定閾值時,將所述相鄰的下一個灰度平均值作為當前有效灰度平均值,并重復執行上述求差及比較過程,直到找出所有有效標準差。步驟607,將所有有效灰度平均值的平均值作為噪聲閾值,或將所有有效灰度平均值的平均值乘以預設系數后的值作為噪聲閾值。本步驟中,該預設系數可根據經驗值或仿真值確定,也可根據其它方法確定。實際應用中,也可不采用上述步驟605描述的過程,而是直接將最小灰度平均值作為第一個有效灰度平均值。圖7為噪聲閾值的再一種確定方法的流程示意圖。如圖7所示,該方法包括如下步驟:步驟701 步驟702的具體實現過程與步驟501 步驟502的具體實現過程一致。步驟703,計算每個矩形區域的標準差。步驟704,將所述標準差進行升序或降序的排序。步驟705,從最小標準差開始,將其作為檢索標準差,計算該檢索標準差與相鄰的設定數目的標準差之間的差值,在所述差值小于第二設定閾值時,將所述檢索標準差作為第一有效標準差;否則,將相鄰的下一個標準差作為檢索標準差,并重復執行上述求差及比較過程,直到找到第一有效標準差。本步驟中,第一有效標準差指定是能表示純噪聲矩形區域的標準差。設定數目可以是根據經驗值或仿真值等確定的數值,如可以是相鄰的I個、2個、3個或4個等。本步驟中,第二設定閾值可根據經驗值或仿真值確定,也可根據其它方法確定。步驟706,從所述第一有效標準差開始,將該第一有效標準差作為當前有效標準差,計算當前有效標準差與相鄰的下一個標準差之間的差值,在所述差值小于所述第二設定閾值時,將所述相鄰的下一個標準差作為當前有效標準差,并重復執行上述求差及比較過程,直到找出所有有效標準差。步驟707,計算所有有效標準差對應的矩形區域的灰度平均值,將所述灰度平均值的平均值作為噪聲閾值,或將所述灰度平均值的平均值乘以預設系數后的值作為噪聲閾值。本步驟中,該預設系數可根據經驗值或仿真值確定,也可根據其它方法確定。實際應用中,也可不采用上述步驟705中描述的過程,而是直接將最小標準差作為第一個有效標準差。本發明實施例中,根據每個子厚層中最不可能存在圖像信息的區域的灰度平均值來確定噪聲閾值時,使得對于不同的成像條件有針對其的具體噪聲閾值,從而噪聲閾值更有針對性,也更準確。以上對本發明實施例中的磁共振成像的圖像重建方法進行了詳細描述,下面再對本發明實施例中磁共振成像的重建裝置進行詳細描述。圖9為本發明實施例中磁共振成像的圖像重建裝置的一個結構示意圖。如圖9所不,該裝置包括:一個蒙板生成模塊901、一個噪聲削弱模塊902和一個圖像合并模塊903。其中,蒙板生成模塊901用于對圖像域的各個子厚層,利用平均濾波器對成像區域同一位置對應的各個子厚層的片層進行濾波,得到對應各片層的蒙板。噪聲削弱模塊902用于對每個蒙板,將該蒙板的像素灰度值逐一與預先確定的噪聲閾值進行比較,對于小于所述噪聲閾值的像素灰度值,將該蒙板所對應片層的相應位置像素值按照預設標準降低。其中,噪聲削弱模塊902的具體操作過程可與圖4所示步驟402中的具體操作過程一致。例如,可將該蒙板所對應片層的相應位置像素值設置為O ;或者也可將該蒙板所對應片層的相應位置像素值乘以小于I的系數,例如,0.1,0.2,0.3,0.4,0.5等。噪聲閾值可以根據經驗值或仿真值確定,也可以根據每個子厚層中最不可能存在圖像信息的區域的灰度平均值來確定。圖像合并模塊903用于對成像區域同一位置對應的各個子厚層的片層進行合并,得到成像區域的磁共振圖像。圖像合并模塊903的具體操作過程可與圖4所示步驟403中的具體操作過程一致。例如,可以是像素相加法,也可以是像素平方和法等。本發明實施例中,由于是對降低了噪聲像素值的片層進行的合并,因此削弱了噪聲累積的影響,提高了磁共振圖像的信噪比,進而提高了圖像質量。在根據每個子厚層中最不可能存在圖像信息的區域的灰度平均值來確定時,可有多種實現方式。相應地,本實施例中的磁共振成像的圖像重建裝置可進一步包括如圖10至圖12中的任一附圖所示的功能模塊。圖10示出了本發明實施例中磁共振成像的圖像重建裝置的又一個結構示意圖。如圖10所示,該裝置除包括蒙板生成模塊901、噪聲削弱模塊902和圖像合并模塊903之夕卜,進一步還包括:一個第一片層選取模塊1001、一個第一區域選取模塊1002、一個第一數據計算模塊1003、一個第一排序模塊1004、一個第一初始值確定模塊1005、一個第一有效值確定模塊1006和一個第一噪聲閾值確定模塊1007。其中,第一片層選取模塊1001用于對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層。第一區域選取模塊1002用于分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域。第一數據計算模塊1003用于計算每個矩形區域的灰度平均值及標準差,并將每個矩形區域的灰度平均值及標準差對應存儲為一個數據對。第一排序模塊1004用于將所述數據對按照灰度平均值的大小進行升序或降序的排序。第一初始值確定模塊1005用于從包含最小灰度平均值的數據對開始,將該數據對作為檢索數據對,計算該檢索數據對與相鄰的設定數目的數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述檢索數據對作為第一有效數據對;否則,將相鄰的下一個數據對作為檢索數據對,并重復執行上述求差及比較過程,直到找到第一有效數據對。第一有效值確定模塊1006用于從所述第一有效數據對開始,將該第一有效數據對作為當前有效數據對;計算當前有效數據對與相鄰的下一個數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述相鄰的下一個數據對作為當前有效數據對,并繼續查找下一個有效數據對;否則,停止有效數據對的查找。第一噪聲閾值確定模塊1007用于將所有有效數據對的灰度平均值的平均值作為噪聲閾值,或將所有有效數據對的灰度平均值的平均值乘以預設系數后的值作為噪聲閾值。其中,各新增功能模塊的具體操作過程可與圖5所示各相應步驟中的具體實現過程一致,此處不再一一贅述。實際應用中,所述第一初始值確定模塊1005也可不米用上述描述的確定第一有效數據對的過程,而是直接將包含最小灰度平均值的數據對作為第一有效數據對。圖11示出了本發明實施例中磁共振成像的圖像重建裝置的再一個結構示意圖。如圖11所示,該裝置除包括蒙板生成模塊901、噪聲削弱模塊902和圖像合并模塊903之夕卜,進一步還包括:一個第二片層選取模塊1101、一個第二區域選取模塊1102、一個第二數據計算模塊1103、一個第二排序模塊1104、一個第二初始值確定模塊1105、一個第二有效值確定模塊1106和一個第二噪聲閾值確定模塊1107。其中,第二片層選取模塊1101用于對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層。第二區域選取模塊1102用于分別從每個邊緣片層的四個角中選取出設定大小的矩形區域。第二數據計算模塊1103用于計算每個矩形區域的灰度平均值。第二排序模塊1104用于將所述灰度平均值進行升序或降序的排序。第二初始值確定模塊1105用于從最小灰度平均值開始,將其作為檢索灰度平均值,計算該檢索灰度平均值與相鄰的設定數目的灰度平均值之間的差值,在所述差值小于第一設定閾值時,將所述檢索灰度平均值作為第一有效灰度平均值;否則,將相鄰的下一個灰度平均值作為檢索數據對,并重復執行上述求差及比較過程,直到找到第一有效灰度平均值。第二有效值確定模塊1106用于從所述第一有效灰度平均值開始,將該第一有效灰度平均值作為當前有效灰度平均值;計算當前有效灰度平均值與相鄰的下一個灰度平均值之間的差值,在所述差值小于所述第一設定閾值時,將所述相鄰的下一個灰度平均值作為當前有效灰度平均值,并繼續查找下一個有效灰度平均值;否則,停止有效灰度平均值的查找。第二噪聲閾值確定模塊1107用于將所有有效灰度平均值的平均值作為噪聲閾值,或將所有有效灰度平均值的平均值乘以預設系數后的值作為噪聲閾值。其中,各新增功能模塊的具體操作過程可與圖6所示各相應步驟中的具體實現過程一致,此處不再一一贅述。實際應用中,上述第二初始值確定模塊1105也可不米用上述描述的確定第一有效灰度平均值的過程,而是直接將最小灰度平均值作為第一有效灰度平均值。圖12示出了本發明實施例中磁共振成像的圖像重建裝置的再一個結構示意圖。如圖12所示,該裝置除包括蒙板生成模塊901、噪聲削弱模塊902和圖像合并模塊903之夕卜,進一步還包括:一個第三片層選取模塊1201、一個第三區域選取模塊1202、一個第三數據計算模塊1203、一個第三排序模塊1204、一個第三初始值確定模塊1205、一個第三有效值確定模塊1206和一個第三噪聲閾值確定模塊1207。第三片層選取模塊1201用于對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層。第三區域選取模塊1201用于分別從每個邊緣片層的四個角中選取出設定大小的矩形區域。第三數據計算模塊1203用于計算每個矩形區域的標準差。第三排序模塊1204用于將所述標準差進行升序或降序的排序。第三初始值確定模塊1205用于從最小標準差開始,將其作為檢索標準差,計算該檢索標準差與相鄰的設定數目的標準差之間的差值,在所述差值小于第二設定閾值時,將所述檢索標準差作為第一有效標準差;否則,將相鄰的下一個標準差作為檢索數據對,并重復執行上述求差及比較過程,直到找到第一有效標準差。第三有效值確定模塊1206用于從所述第一有效標準差開始,將該第一有效標準差作為當前有效標準差;計算當前有效標準差與相鄰的下一個標準差之間的差值,在所述差值小于所述第二設定閾值時,將所述相鄰的下一個標準差作為當前有效標準差,并繼續查找下一個有效標準差;否則,停止有效標準差的查找。第三噪聲閾值確定模塊1207用于計算所有有效標準差對應的矩形區域的灰度平均值,將所述灰度平均值的平均值作為噪聲閾值,或將所述灰度平均值的平均值乘以預設系數后的值作為噪聲閾值。其中,各新增功能模塊的具體操作過程可與圖7所示各相應步驟中的具體實現過程一致,此處不再一一贅述。實際應用中,上述第三初始值確定模塊1205也可不采用上述描述的確定第一有效標準差的過程,而是直接將最小標準差作為第一有效標準差。本發明實施例中,根據每個子厚層中最不可能存在圖像信息的區域的灰度平均值來確定噪聲閾值時,使得對于不同的成像條件有針對其的具體噪聲閾值,從而噪聲閾值更有針對性,也更準確。本領域技術人員可以理解附圖只是一個優選實施例的示意圖,附圖中的模塊或流程并不一定是實施本發明所必須的。本領域技術人員可以理解實施例中的裝置中的模塊可以按照實施例描述分布于實施例的裝置中,也可以進行相應變化位于不同于本實施例的一個或多個裝置中。上述實施例的模塊可以合并為一個模塊,也可以進一步拆分成多個子模塊。本發明實施例中的部分步驟,可以利用軟件實現,相應的軟件程序可以存儲在可讀取的存儲介質中,如光盤或硬盤等。 本發明涉及磁共振成像技術,并且公開了一種磁共振成像的圖像重建方法和一種磁共振成像的圖像重建裝置。其中,所述磁共振成像的圖像重建方法包括:利用平均濾波器對圖像域的各個子厚層的片層進行濾波,得到對應各片層的蒙板;將各蒙板的像素值與一預先確定的噪聲閾值進行比較,對于小于所述噪聲閾值的像素值,將該蒙板所對應片層的相應位置像素值按照預設標準降低;對成像區域同一位置對應的各個子厚層的片層進行合并,得到成像區域的磁共振圖像。本發明能夠削弱了噪聲的影響,提高磁共振圖像的信噪t匕,從而提高圖像質量。
            以上所述僅為本發明的較佳實施例而已,并不用以限制本發明,凡在本發明的精神和原則之內,所作的任何修改、等同替換、改進等,均應包含在本發明的保護范圍之內。
            權利要求
            1.一種磁共振成像的圖像重建方法,包括: 利用平均濾波器對圖像域的各個子厚層的片層進行濾波,得到對應各片層的蒙板;將各蒙板的像素值與一預先確定的噪聲閾值進行比較,對于小于所述噪聲閾值的像素值,將該蒙板所對應片層的相應位置像素值按照預設標準降低; 對成像區域同一位置對應的各個子厚層的片層進行合并,得到成像區域的磁共振圖像。
            2.根據權利要求1所述的方法,其特征在于,該方法進一步包括: al、分別從每個圖像域的子厚層中選取出至少一個邊緣片層; bl、分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域; Cl、計算每個矩形區域的灰度平均值及標準差,并將每個矩形區域的灰度平均值及標準差對應存儲為一個數據對; dl、將所述數據對按照灰度平均值的大小進行升序或降序的排序;el、將包含最小灰度平均值的數據對作為第一有效數據對,并執行步驟f I ;或者,將包含最小灰度平均值的數據對作為檢索數據對,計算該檢索數據對與相鄰的設定數目的數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述檢索數據對作為第一有效數據對,并執行步驟Π,否則,將相鄰的下一個數據對作為檢索數據對,并返回執行本步驟的差值計算過程; H、將所述第一有效數據對作為當前有效數據對; gl、計算當前有效數據對與相鄰的下一個數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述相鄰的下一個數據對作為當前有效數據對,并返回執行本步驟;否則,執行下一步驟; h1、將所有有效數據對的灰度平均值的平均值作為所述噪聲閾值,或將所有有效數據對的灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。
            3.根據權利要求1所述的方法,其特征在于,該方法進一步包括: a2、分別從每個圖像域的子厚層中選取出至少一個邊緣片層; b2、分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域; c2、計算每個矩形區域的灰度平均值; d2、將所述灰度平均值進行升序或降序的排序; e2、將最小灰度平均值作為第一有效灰度平均值,并執行步驟f2 ;或者,將最小灰度平均值作為檢索灰度平均值,計算該檢索灰度平均值與相鄰的設定數目的灰度平均值之間的差值,在所述差值小于第一設定閾值時,將所述檢索灰度平均值作為第一有效灰度平均值,并執行步驟f2,否則,將相鄰的下一個灰度平均值作為檢索灰度平均值,并返回執行本步驟的差值計算過程; f2、將所述第一有效灰度平均值作為當前有效灰度平均值; g2、計算當前有效灰度平均值與`相鄰的下一個灰度平均值之間的差值,在所述差值小于所述第一設定閾值時,將所述相鄰的下一個灰度平均值作為當前有效灰度平均值,并返回執行本步驟;否則,執行下一步驟; h2、將所有有效灰度平均值的平均值作為所述噪聲閾值,或將所有有效灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。
            4.根據權利要求1所述的方法,其特征在于,該方法進一步包括: a3、分別從每個圖像域的子厚層中選取出至少一個邊緣片層; b3、分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域; c3、計算每個矩形區域的標準差; d3、將所述標準差進行升序或降序的排序; e3、將最小標準差作為第一有效標準差,并執行步驟f3 ;或者,將最小標準差作為檢索標準差,計算該檢索標準差與相鄰的設定數目的標準差之間的差值,在所述差值小于第二設定閾值時,將所述檢索標準差作為第一有效標準差,并執行步驟f3,否則,將相鄰的下一個標準差作為檢索 標準差,并返回執行本步驟的差值計算過程;f3、將所述第一有效標準差作為當前有效標準差; g3、計算當前有效標準差與相鄰的下一個標準差之間的差值,在所述差值小于所述第二設定閾值時,將所述相鄰的下一個標準差作為當前有效標準差,并返回執行本步驟;否貝U,執行下一步驟; h3、計算所有有效標準差對應的矩形區域的灰度平均值,將所述灰度平均值的平均值作為所述噪聲閾值,或將所述灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。
            5.根據權利要求1至4中任一項所述的方法,其特征在于,所述將該蒙板所對應片層的相應位置像素值按照預設標準降低包括: 將該蒙板所對應片層的相應位置像素值設置為O ;或者, 將該蒙板所對應片層的相應位置像素值乘以小于I的系數。
            6.根據權利要求1至4中任一項所述的方法,其特征在于,利用像素相加法或像素平方和法,對成像區域同一位置對應的各個子厚層的片層進行合并。
            7.—種磁共振成像的圖像重建裝置,包括: 一個蒙板生成模塊,用于利用平均濾波器對圖像域的各個子厚層的片層進行濾波,得到對應各片層的蒙板; 一個噪聲削弱模塊,用于將各蒙板的像素值與一預先確定的噪聲閾值進行比較,對于小于所述噪聲閾值的像素值,將該蒙板所對應片層的相應位置像素值按照預設標準降低;一個圖像合并模塊,用于對成像區域同一位置對應的各個子厚層的片層進行合并,得到成像區域的磁共振圖像。
            8.根據權利要求7所述的裝置,其特征在于,該裝置進一步包括: 一個第一片層選取模塊,用于對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層; 一個第一區域選取模塊,用于分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域; 一個第一數據計算模塊,用于計算每個矩形區域的灰度平均值及標準差,并將每個矩形區域的灰度平均值及標準差對應存儲為一個數據對; 一個第一排序模塊,用于將所述數據對按照灰度平均值的大小進行升序或降序的排序; 一個第一初始值確定模塊,用于將包含最小灰度平均值的數據對作為第一有效數據對;或者從包含最小灰度平均值的數據對開始,將當前數據對作為檢索數據對,計算該檢索數據對與相鄰的設定數目的數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述檢索數據對作為第一有效數據對; 一個第一有效值確定模塊,用于從所述第一有效數據對開始,將該第一有效數據對作為當前有效數據對;計算當前有效數據對與相鄰的下一個數據對之間的灰度平均值之差及標準差之差,在兩個差值分別小于各自對應的設定閾值時,將所述相鄰的下一個數據對作為當前有效數據對,并繼續查找下一個有效數據對,否則停止有效數據對的查找; 一個第一噪聲閾值確定模塊,用于將所有有效數據對的灰度平均值的平均值作為所述噪聲閾值,或將所有有效數據對的灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。
            9.根據權利要求7所述的裝置,其特征在于,該裝置進一步包括: 一個第二片層選取模塊,用于對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層; 一個第二區域選取模塊,用于分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域; 一個第二數據計算模塊,用于計算每個矩形區域的灰度平均值; 一個第二排序模塊,用于將所述灰度平均值進行升序或降序的排序; 一個第二初始值確定模塊,用于將最小灰度平均值作為第一有效灰度平均值;或者從最小灰度平均值開始,將當前灰度平均值作為檢索灰度平均值,計算該檢索灰度平均值與相鄰的設定數目的灰度平均值之間的差值,在所述差值小于第一設定閾值時,將所述檢索灰度平均值作為第一有效灰度平均值; 一個第二有效值確定模塊,用于從所述第一有效灰度平均值開始,將該第一有效灰度平均值作為當前有效灰度平均值;計算當前有效灰度平均值與相鄰的下一個灰度平均值之間的差值,在所述差值小于 所述第一設定閾值時,將所述相鄰的下一個灰度平均值作為當前有效灰度平均值,并繼續查找下一個有效灰度平均值,否則停止有效灰度平均值的查找; 一個第二噪聲閾值確定模塊,用于將所有有效灰度平均值的平均值作為所述噪聲閾值,或將所有有效灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。
            10.根據權利要求7所述的裝置,其特征在于,該裝置進一步包括: 一個第三片層選取模塊,用于對圖像域的各個子厚層,分別從每個子厚層中選取出至少一個邊緣片層; 一個第三區域選取模塊,用于分別從每個邊緣片層的四個角中選取出至少一個設定大小的矩形區域; 一個第三數據計算模塊,用于計算每個矩形區域的標準差; 一個第三排序模塊,用于將所述標準差進行升序或降序的排序; 一個第三初始值確定模塊,用于將最小標準差作為第一有效標準差;或者從最小標準差開始,將當前標準差作為檢索標準差,計算該檢索標準差與相鄰的設定數目的標準差之間的差值,在所述差值小于第二設定閾值時,將所述檢索標準差作為第一有效標準差; 一個第三有效值確定模塊,用于從所述第一有效標準差開始,將該第一有效標準差作為當前有效標準差;計算當前有效標準差與相鄰的下一個標準差之間的差值,在所述差值小于所述第二設定閾值時,將所述相鄰的下一個標準差作為當前有效標準差,并繼續查找下一個有效標準差,否則停止有效標準差的查找; 一個第三噪聲閾值確定 模塊,用于計算所有有效標準差對應的矩形區域的灰度平均值,將所述灰度平均值的平均值作為所述噪聲閾值,或將所述灰度平均值的平均值乘以預設系數后作為所述噪聲閾值。
            全文摘要
            本發明涉及磁共振成像技術,并且公開了一種磁共振成像的圖像重建方法和一種磁共振成像的圖像重建裝置。其中,所述磁共振成像的圖像重建方法包括利用平均濾波器對圖像域的各個子厚層的片層進行濾波,得到對應各片層的蒙板;將各蒙板的像素值與一預先確定的噪聲閾值進行比較,對于小于所述噪聲閾值的像素值,將該蒙板所對應片層的相應位置像素值按照預設標準降低;對成像區域同一位置對應的各個子厚層的片層進行合并,得到成像區域的磁共振圖像。本發明能夠削弱了噪聲的影響,提高磁共振圖像的信噪比,從而提高圖像質量。
            文檔編號G06T5/00GK103186886SQ201110456518
            公開日2013年7月3日 申請日期2011年12月31日 優先權日2011年12月31日
            發明者李國斌 申請人:西門子(深圳)磁共振有限公司
            網友詢問留言 已有0條留言
            • 還沒有人留言評論。精彩留言會獲得點贊!
            1
            婷婷六月激情在线综合激情,亚洲国产大片,久久中文字幕综合婷婷,精品久久久久久中文字幕,亚洲一区二区三区高清不卡,99国产精品热久久久久久夜夜嗨 ,欧美日韩亚洲综合在线一区二区,99国产精品电影,伊人精品线视天天综合,精品伊人久久久大香线蕉欧美
            亚洲精品1区 国产成人一级 91精品国产欧美一区二区 亚洲精品乱码久久久久久下载 国产精品久久久久久久伊一 九色国产 国产精品九九视频 伊人久久成人爱综合网 欧美日韩亚洲区久久综合 欧美日本一道免费一区三区 夜夜爽一区二区三区精品 欧美日韩高清一区二区三区 国产成人av在线 国产精品对白交换绿帽视频 国产视频亚洲 国产在线欧美精品 国产精品综合网 国产日韩精品欧美一区色 国产日韩精品欧美一区喷 欧美日韩在线观看区一二 国产区精品 欧美视频日韩视频 中文字幕天天躁日日躁狠狠躁97 视频一二三区 欧美高清在线精品一区二区不卡 国产精品揄拍一区二区久久 99久久综合狠狠综合久久aⅴ 亚洲乱码视频在线观看 日韩在线第二页 亚洲精品无码专区在线播放 成人亚洲网站www在线观看 欧美三级一区二区 99久久精品免费看国产高清 91麻豆国产在线观看 最新日韩欧美不卡一二三区 成人在线观看不卡 日韩国产在线 在线亚洲精品 亚洲午夜久久久久中文字幕 国产精品成人久久久久久久 精品国产一区二区在线观看 欧美精品国产一区二区三区 中文在线播放 亚洲第一页在线视频 国产午夜精品福利久久 九色国产 精品国产九九 国产永久视频 久久精品人人做人人综合试看 国产一区二区三区免费观看 亚洲精品国产电影 9999热视频 国产精品资源在线 麻豆久久婷婷国产综合五月 国产精品免费一级在线观看 亚洲国产一区二区三区青草影视 中文在线播放 国产成人综合在线 国产在线观看色 国产亚洲三级 国产片一区二区三区 久久99精品久久久久久牛牛影视 亚洲欧美日韩国产 四虎永久免费网站 国产一毛片 国产精品视频在 九九热在线精品 99精品福利视频 色婷婷色99国产综合精品 97成人精品视频在线播放 精品久久久久久中文字幕 亚洲欧美一区二区三区孕妇 亚洲欧美成人网 日韩高清在线二区 国产尤物在线观看 在线不卡一区二区 91网站在线看 韩国精品福利一区二区 欧美日韩国产成人精品 99热精品久久 国产精品免费视频一区 高清视频一区 精品九九久久 欧美日韩在线观看免费 91欧美激情一区二区三区成人 99福利视频 亚洲国产精品91 久热国产在线 精品久久久久久中文字幕女 国产精品久久久久久久久99热 成人自拍视频网 国产精品视频久久久久久 久久影院国产 国产玖玖在线观看 99精品在线免费 亚洲欧美一区二区三区导航 久久久久久久综合 国产欧美日韩精品高清二区综合区 国产精品视频自拍 亚洲一级片免费 久久久久久九九 国产欧美自拍视频 视频一区二区在线观看 欧美日韩一区二区三区久久 中文在线亚洲 伊人热人久久中文字幕 日韩欧美亚洲国产一区二区三区 欧美亚洲国产成人高清在线 欧美日韩国产码高清综合人成 国产性大片免费播放网站 亚洲午夜综合网 91精品久久一区二区三区 国产无套在线播放 国产精品视频网站 国产成人亚洲精品老王 91在线网站 国产视频97 欧美黑人欧美精品刺激 国产一区二区三区免费在线视频 久久久国产精品免费看 99re6久精品国产首页 久久精品91 国产成人一级 国产成人精品曰本亚洲 日本福利在线观看 伊人成综合网 久久综合一本 国产综合久久久久久 久久精品成人免费看 久久福利 91精品国产91久久久久久麻豆 亚洲精品成人在线 亚洲伊人久久精品 欧美日本二区 国产永久视频 国产一区二 一区二区福利 国产一毛片 亚洲精品1区 毛片一区二区三区 伊人久久大香线蕉综合影 国产欧美在线观看一区 亚洲国产欧洲综合997久久 国产一区二区免费视频 国产91精品对白露脸全集观看 久久亚洲国产伦理 欧美成人伊人久久综合网 亚洲性久久久影院 久久99国产精一区二区三区! 91精品国产欧美一区二区 欧美日韩亚洲区久久综合 日韩精品一二三区 久久久夜色精品国产噜噜 国产在线精品福利91香蕉 久久久久久久亚洲精品 97se色综合一区二区二区 91国语精品自产拍在线观看性色 91久久国产综合精品女同我 日韩中文字幕a 国产成人亚洲日本精品 久久国产精品-国产精品 久久国产经典视频 久久国产精品伦理 亚洲第一页在线视频 国产精品久久久久三级 日韩毛片网 久久免费高清视频 麻豆国产在线观看一区二区 91麻豆国产福利在线观看 国产成人精品男人的天堂538 一区二区三区中文字幕 免费在线视频一区 欧美日韩国产成人精品 国产综合网站 国产资源免费观看 亚洲精品亚洲人成在线播放 精品久久久久久中文字幕专区 亚洲人成人毛片无遮挡 国产一起色一起爱 国产香蕉精品视频在 九九热免费观看 日韩亚洲欧美一区 九九热精品在线观看 精品久久久久久中文字幕专区 亚洲欧美自拍偷拍 国产精品每日更新 久久久久国产一级毛片高清板 久久天天躁狠狠躁夜夜中文字幕 久久精品片 日韩在线毛片 国产成人精品本亚洲 国产成人精品一区二区三区 九九热在线观看 国产r级在线观看 国产欧美日韩精品高清二区综合区 韩国电影一区二区 国产精品毛片va一区二区三区 五月婷婷伊人网 久久一区二区三区免费 一本色道久久综合狠狠躁篇 亚洲综合色站 国产尤物在线观看 亚洲一区亚洲二区 免费在线视频一区 欧洲精品视频在线观看 日韩中文字幕a 中文字幕日本在线mv视频精品 91精品在线免费视频 精品国产免费人成在线观看 精品a级片 中文字幕日本在线mv视频精品 日韩在线精品视频 婷婷丁香色 91精品国产高清久久久久 国产成人精品日本亚洲直接 五月综合视频 欧美日韩在线亚洲国产人 精液呈暗黄色 亚洲乱码一区 久久精品中文字幕不卡一二区 亚洲天堂精品在线 激情婷婷综合 国产免费久久精品久久久 国产精品亚洲二区在线 久久免费播放视频 五月婷婷丁香综合 在线亚洲欧美日韩 久久免费精品高清麻豆 精品久久久久久中文字幕 亚洲一区网站 国产精品福利社 日韩中文字幕免费 亚洲综合丝袜 91精品在线播放 国产精品18 亚洲日日夜夜 伊人久久大香线蕉综合影 亚洲精品中文字幕乱码影院 亚洲一区二区黄色 亚洲第一页在线视频 一区二区在线观看视频 国产成人福利精品视频 亚洲高清二区 国内成人免费视频 精品亚洲性xxx久久久 国产精品合集一区二区三区 97av免费视频 国产一起色一起爱 国产区久久 国产资源免费观看 99精品视频免费 国产成人一级 国产精品九九免费视频 欧美91精品久久久久网免费 99热国产免费 久久精品色 98精品国产综合久久 久久精品播放 中文字幕视频免费 国产欧美日韩一区二区三区在线 精品久久蜜桃 国产小视频精品 一本色道久久综合狠狠躁篇 91在线免费观看 亚洲精品区 伊人成综合网 伊人热人久久中文字幕 伊人黄色片 99国产精品热久久久久久夜夜嗨 久久免费精品视频 亚洲一区二区三区高清不卡 久久久久国产一级毛片高清板 国产片一区二区三区 久久狠狠干 99久久婷婷国产综合精品电影 国产99区 国产精品成人久久久久 久久狠狠干 青青国产在线观看 亚洲高清国产拍精品影院 国产精品一区二区av 九九热在线免费视频 伊人久久国产 国产精品久久久久久久久久一区 在线观看免费视频一区 国产精品自在在线午夜区app 国产精品综合色区在线观看 国产毛片久久久久久国产毛片 97国产免费全部免费观看 国产精品每日更新 国产尤物视频在线 九九视频这里只有精品99 一本一道久久a久久精品综合 久久综合给会久久狠狠狠 国产成人精品男人的天堂538 欧美一区二区高清 毛片一区二区三区 国产欧美日韩在线观看一区二区三区 在线国产二区 欧美不卡网 91在线精品中文字幕 在线国产福利 国内精品91久久久久 91亚洲福利 日韩欧美国产中文字幕 91久久精品国产性色也91久久 亚洲性久久久影院 欧美精品1区 国产热re99久久6国产精品 九九热免费观看 国产精品欧美日韩 久久久久国产一级毛片高清板 久久国产经典视频 日韩欧美亚洲国产一区二区三区 欧美亚洲综合另类在线观看 国产精品自在在线午夜区app 97中文字幕在线观看 视频一二三区 精品国产一区在线观看 国产欧美日韩在线一区二区不卡 欧美一区二三区 伊人成人在线观看 国内精品91久久久久 97在线亚洲 国产在线不卡一区 久久久全免费全集一级全黄片 国产精品v欧美精品∨日韩 亚洲毛片网站 在线不卡一区二区 99re热在线视频 久久激情网 国产毛片一区二区三区精品 久久亚洲综合色 中文字幕视频免费 国产视频亚洲 婷婷伊人久久 国产一区二区免费播放 久久99国产精品成人欧美 99国产在线视频 国产成人免费视频精品一区二区 国产不卡一区二区三区免费视 国产码欧美日韩高清综合一区 久久精品国产主播一区二区 国产一区电影 久久精品国产夜色 国产精品国产三级国产 日韩一区二区三区在线 久久97久久97精品免视看 久久国产免费一区二区三区 伊人久久大香线蕉综合电影网 99re6久精品国产首页 久久激情网 亚洲成人高清在线 国产精品网址 国产成人精品男人的天堂538 香蕉国产综合久久猫咪 国产专区中文字幕 91麻豆精品国产高清在线 久久国产经典视频 国产精品成人va在线观看 国产精品爱啪在线线免费观看 日本精品久久久久久久久免费 亚洲综合一区二区三区 久久五月网 精品国产网红福利在线观看 久久综合亚洲伊人色 亚洲国产精品久久久久久网站 在线日韩国产 99国产精品热久久久久久夜夜嗨 国产综合精品在线 国产区福利 精品亚洲综合久久中文字幕 国产制服丝袜在线 毛片在线播放网站 在线观看免费视频一区 国产精品久久久精品三级 亚洲国产电影在线观看 最新日韩欧美不卡一二三区 狠狠综合久久综合鬼色 日本精品1在线区 国产日韩一区二区三区在线播放 欧美日韩精品在线播放 亚洲欧美日韩国产一区二区三区精品 久久综合久久网 婷婷六月激情在线综合激情 亚洲乱码一区 国产专区91 97av视频在线观看 精品久久久久久中文字幕 久久五月视频 国产成人福利精品视频 国产精品网址 中文字幕视频在线 精品一区二区三区免费视频 伊人手机在线视频 亚洲精品中文字幕乱码 国产在线视频www色 色噜噜国产精品视频一区二区 精品亚洲成a人在线观看 国产香蕉尹人综合在线 成人免费一区二区三区在线观看 国产不卡一区二区三区免费视 欧美精品久久天天躁 国产专区中文字幕 久久精品国产免费中文 久久精品国产免费一区 久久无码精品一区二区三区 国产欧美另类久久久精品免费 欧美精品久久天天躁 亚洲精品在线视频 国产视频91在线 91精品福利一区二区三区野战 日韩中文字幕免费 国产精品99一区二区三区 欧美成人高清性色生活 国产精品系列在线观看 亚洲国产福利精品一区二区 国产成人在线小视频 国产精品久久久久免费 99re热在线视频 久久久久久久综合 一区二区国产在线播放 成人国产在线视频 亚洲精品乱码久久久久 欧美日韩一区二区综合 精品久久久久免费极品大片 中文字幕视频二区 激情粉嫩精品国产尤物 国产成人精品一区二区视频 久久精品中文字幕首页 亚洲高清在线 国产精品亚洲一区二区三区 伊人久久艹 中文在线亚洲 国产精品一区二区在线播放 国产精品九九免费视频 亚洲二区在线播放 亚洲狠狠婷婷综合久久久久网站 亚洲欧美日韩网站 日韩成人精品 亚洲国产一区二区三区青草影视 91精品国产福利在线观看 国产精品久久久久久久久99热 国产一区二区精品尤物 久碰香蕉精品视频在线观看 亚洲日日夜夜 在线不卡一区二区 国产午夜亚洲精品 九九热在线视频观看这里只有精品 伊人手机在线视频 91免费国产精品 日韩欧美中字 91精品国产91久久久久 国产全黄三级播放 视频一区二区三区免费观看 国产开裆丝袜高跟在线观看 国产成人欧美 激情综合丝袜美女一区二区 国产成人亚洲综合无 欧美精品一区二区三区免费观看 欧美亚洲国产日韩 日韩亚州 国产欧美日韩精品高清二区综合区 亚洲午夜国产片在线观看 精品久久久久久中文字幕 欧美精品1区 久久伊人久久亚洲综合 亚洲欧美日韩精品 国产成人精品久久亚洲高清不卡 久久福利影视 国产精品99精品久久免费 久久久久免费精品视频 国产日产亚洲精品 亚洲国产午夜电影在线入口 精品无码一区在线观看 午夜国产精品视频 亚洲一级片免费 伊人久久大香线蕉综合影 国产精品久久影院 久碰香蕉精品视频在线观看 www.欧美精品 在线小视频国产 亚洲国产天堂久久综合图区 欧美一区二区三区不卡 日韩美女福利视频 九九精品免视频国产成人 不卡国产00高中生在线视频 亚洲第一页在线视频 欧美日韩在线播放成人 99re视频这里只有精品 国产精品91在线 精品乱码一区二区三区在线 国产区久久 91麻豆精品国产自产在线观看一区 日韩精品成人在线 九九热在线观看 国产精品久久不卡日韩美女 欧美一区二区三区综合色视频 欧美精品免费一区欧美久久优播 国产精品网址 国产专区中文字幕 国产精品欧美亚洲韩国日本久久 日韩美香港a一级毛片 久久精品123 欧美一区二区三区免费看 99r在线视频 亚洲精品国产字幕久久vr 国产综合激情在线亚洲第一页 91免费国产精品 日韩免费小视频 亚洲国产精品综合一区在线 国产亚洲第一伦理第一区 在线亚洲精品 国产精品一区二区制服丝袜 国产在线成人精品 九九精品免视频国产成人 亚洲国产网 欧美日韩亚洲一区二区三区在线观看 在线亚洲精品 欧美一区二区三区高清视频 国产成人精品男人的天堂538 欧美日韩在线观看区一二 亚洲欧美一区二区久久 久久精品中文字幕首页 日本高清www午夜视频 久久精品国产免费 久久999精品 亚洲国产精品欧美综合 88国产精品视频一区二区三区 91久久偷偷做嫩草影院免费看 国产精品夜色视频一区二区 欧美日韩导航 国产成人啪精品午夜在线播放 一区二区视频在线免费观看 99久久精品国产自免费 精液呈暗黄色 久久99国产精品 日本精品久久久久久久久免费 精品国产97在线观看 99re视频这里只有精品 国产视频91在线 999av视频 亚洲美女视频一区二区三区 久久97久久97精品免视看 亚洲国产成人久久三区 99久久亚洲国产高清观看 日韩毛片在线视频 综合激情在线 91福利一区二区在线观看 一区二区视频在线免费观看 激情粉嫩精品国产尤物 国产成人精品曰本亚洲78 国产成人精品本亚洲 国产精品成人免费视频 国产成人啪精品视频免费软件 久久精品国产亚洲妲己影院 国产精品成人久久久久久久 久久大香线蕉综合爱 欧美一区二区三区高清视频 99热国产免费 在线观看欧美国产 91精品视频在线播放 国产精品福利社 欧美精品一区二区三区免费观看 国产一区二区免费视频 国产午夜精品一区二区 精品视频在线观看97 91精品福利久久久 国产一区福利 国产综合激情在线亚洲第一页 国产精品久久久久久久久久久不卡 九色国产 在线日韩国产 黄网在线观看 亚洲一区小说区中文字幕 中文字幕丝袜 日本二区在线观看 日本国产一区在线观看 欧美日韩一区二区三区久久 欧美精品亚洲精品日韩专 国产日产亚洲精品 久久综合九色综合欧美播 亚洲国产欧美无圣光一区 欧美视频区 亚洲乱码视频在线观看 久久无码精品一区二区三区 九九热精品免费视频 久久99精品久久久久久牛牛影视 国产精品成久久久久三级 国产一区福利 午夜国产精品视频 日本二区在线观看 99久久网站 国产亚洲天堂 精品国产一区二区三区不卡 亚洲国产日韩在线一区 国产成人综合在线观看网站 久久免费高清视频 欧美在线导航 午夜精品久久久久久99热7777 欧美久久综合网 国产小视频精品 国产尤物在线观看 亚洲国产精品综合一区在线 欧美一区二区三区不卡视频 欧美黑人欧美精品刺激 日本福利在线观看 久久国产偷 国产手机精品一区二区 国产热re99久久6国产精品 国产高清啪啪 欧美亚洲国产成人高清在线 国产在线第三页 亚洲综合一区二区三区 99r在线视频 99精品久久久久久久婷婷 国产精品乱码免费一区二区 国产在线精品福利91香蕉 国产尤物视频在线 五月婷婷亚洲 中文字幕久久综合伊人 亚洲精品一级毛片 99国产精品电影 在线视频第一页 久久99国产精品成人欧美 国产白白视频在线观看2 成人精品一区二区www 亚洲成人网在线观看 麻豆91在线视频 色综合合久久天天综合绕视看 久久精品国产免费高清 国产不卡一区二区三区免费视 欧美国产中文 99精品欧美 九九在线精品 国产中文字幕在线免费观看 国产一区中文字幕在线观看 国产成人一级 国产精品一区二区制服丝袜 国产一起色一起爱 亚洲精品成人在线 亚洲欧美精品在线 国产欧美自拍视频 99精品久久久久久久婷婷 久99视频 国产热re99久久6国产精品 视频一区亚洲 国产精品视频分类 国产精品成在线观看 99re6久精品国产首页 亚洲在成人网在线看 亚洲国产日韩在线一区 久久国产三级 日韩国产欧美 欧美在线一区二区三区 国产精品美女一级在线观看 成人午夜免费福利视频 亚洲天堂精品在线 91精品国产手机 欧美日韩视频在线播放 狠狠综合久久综合鬼色 九一色视频 青青视频国产 亚洲欧美自拍一区 中文字幕天天躁日日躁狠狠躁97 日韩免费大片 996热视频 伊人成综合网 亚洲天堂欧美 日韩精品亚洲人成在线观看 久久综合给会久久狠狠狠 日韩精品亚洲人成在线观看 日韩国产欧美 亚洲成aⅴ人片在线影院八 亚洲精品1区 99久久精品免费 国产精品高清在线观看 国产精品久久久免费视频 在线亚洲欧美日韩 91在线看视频 国产精品96久久久久久久 欧美日韩国产成人精品 91在线亚洲 热久久亚洲 国产精品美女免费视频观看 日韩在线毛片 亚洲永久免费视频 九九免费在线视频 亚洲一区网站 日本高清二区视频久二区 精品国产美女福利在线 伊人久久艹 国产精品久久久久三级 欧美成人精品第一区二区三区 99久久精品国产自免费 在线观看日韩一区 国产中文字幕一区 成人免费午夜视频 欧美日韩另类在线 久久99国产精品成人欧美 色婷婷中文网 久久天天躁夜夜躁狠狠躁2020 欧美成人伊人久久综合网 国产精品福利资源在线 国产伦精品一区二区三区高清 国产精品亚洲综合色区韩国 亚洲一区欧美日韩 色综合视频 国语自产精品视频在线区 国产高清a 成人国内精品久久久久影 国产在线精品香蕉综合网一区 国产不卡在线看 国产成人精品精品欧美 国产欧美日韩综合精品一区二区三区 韩国电影一区二区 国产在线视频www色 91中文字幕在线一区 国产人成午夜免视频网站 亚洲综合一区二区三区 色综合视频一区二区观看 久久五月网 九九热精品在线观看 国产一区二区三区国产精品 99久热re在线精品996热视频 亚洲国产网 在线视频亚洲一区 日韩字幕一中文在线综合 国产高清一级毛片在线不卡 精品国产色在线 国产高清视频一区二区 精品日本久久久久久久久久 亚洲国产午夜精品乱码 成人免费国产gav视频在线 日韩欧美一区二区在线观看 欧美曰批人成在线观看 韩国电影一区二区 99re这里只有精品6 日韩精品一区二区三区视频 99re6久精品国产首页 亚洲欧美一区二区三区导航 欧美色图一区二区三区 午夜精品视频在线观看 欧美激情在线观看一区二区三区 亚洲热在线 成人国产精品一区二区网站 亚洲一级毛片在线播放 亚洲一区小说区中文字幕 亚洲午夜久久久久影院 国产自产v一区二区三区c 国产精品视频免费 久久调教视频 国产成人91激情在线播放 国产精品欧美亚洲韩国日本久久 久久亚洲日本不卡一区二区 91中文字幕网 成人国产在线视频 国产视频91在线 欧美成人精品第一区二区三区 国产精品福利在线 久久综合九色综合精品 欧美一区二区三区精品 久久国产综合尤物免费观看 久久99青青久久99久久 日韩精品免费 久久国产精品999 91亚洲视频在线观看 国产精品igao视频 色综合区 在线亚洲欧国产精品专区 国产一区二区三区在线观看视频 亚洲精品成人在线 一区二区国产在线播放 中文在线亚洲 亚洲精品第一国产综合野 国产一区二区精品久久 一区二区三区四区精品视频 99热精品久久 中文字幕视频二区 国产成人精品男人的天堂538 99精品影视 美女福利视频一区二区 久久午夜夜伦伦鲁鲁片 综合久久久久久久综合网 国产精品国产欧美综合一区 国产99视频在线观看 国产亚洲女在线精品 婷婷影院在线综合免费视频 国产亚洲3p一区二区三区 91成人爽a毛片一区二区 亚洲一区二区高清 国产欧美亚洲精品第二区首页 欧美日韩导航 亚洲高清二区 欧美激情观看一区二区久久 日韩毛片在线播放 亚洲欧美日韩高清中文在线 亚洲日本在线播放 国产精品一区二区制服丝袜 精品国产一区二区三区不卡 国产不卡在线看 国产欧美网站 四虎永久在线观看视频精品 国产黄色片在线观看 夜夜综合 一本色道久久综合狠狠躁篇 欧美亚洲综合另类在线观看 国产91在线看 伊人久久国产 欧美一区二区在线观看免费网站 国产精品久久久久三级 久久福利 日韩中文字幕a 亚洲午夜久久久久影院 91在线高清视频 国产亚洲一区二区三区啪 久久人精品 国产精品亚洲午夜一区二区三区 综合久久久久久 久久伊人一区二区三区四区 国产综合久久久久久 日韩一区精品视频在线看 国产精品日韩欧美制服 日本精品1在线区 99re视频 无码av免费一区二区三区试看 国产视频1区 日韩欧美中文字幕一区 日本高清中文字幕一区二区三区a 亚洲国产欧美无圣光一区 国产在线视频一区二区三区 欧美国产第一页 在线亚洲欧美日韩 日韩中文字幕第一页 在线不卡一区二区 伊人久久青青 国产精品一区二区在线播放 www.五月婷婷 麻豆久久婷婷国产综合五月 亚洲精品区 久久国产欧美另类久久久 99在线视频免费 伊人久久中文字幕久久cm 久久精品成人免费看 久久这里只有精品首页 88国产精品视频一区二区三区 中文字幕日本在线mv视频精品 国产在线精品成人一区二区三区 伊人精品线视天天综合 亚洲一区二区黄色 国产尤物视频在线 亚洲精品99久久久久中文字幕 国产一区二区三区免费观看 伊人久久大香线蕉综合电影网 国产成人精品区在线观看 日本精品一区二区三区视频 日韩高清在线二区 久久免费播放视频 一区二区成人国产精品 国产精品免费精品自在线观看 亚洲精品视频二区 麻豆国产精品有码在线观看 精品日本一区二区 亚洲欧洲久久 久久中文字幕综合婷婷 中文字幕视频在线 国产成人精品综合在线观看 91精品国产91久久久久福利 精液呈暗黄色 香蕉国产综合久久猫咪 国产专区精品 亚洲精品无码不卡 国产永久视频 亚洲成a人片在线播放观看国产 一区二区国产在线播放 亚洲一区二区黄色 欧美日韩在线观看视频 亚洲精品另类 久久国产综合尤物免费观看 国产一区二区三区国产精品 高清视频一区 国产精品igao视频 国产精品资源在线 久久综合精品国产一区二区三区 www.五月婷婷 精品色综合 99热国产免费 麻豆福利影院 亚洲伊人久久大香线蕉苏妲己 久久电影院久久国产 久久精品伊人 在线日韩理论午夜中文电影 亚洲国产欧洲综合997久久 伊人国产精品 久草国产精品 欧美一区精品二区三区 亚洲成人高清在线 91免费国产精品 日韩精品福利在线 国产一线在线观看 国产不卡在线看 久久99青青久久99久久 亚洲精品亚洲人成在线播放 99久久免费看国产精品 国产日本在线观看 青草国产在线视频 麻豆久久婷婷国产综合五月 国产中文字幕一区 91久久精品国产性色也91久久 国产一区a 国产欧美日韩成人 国产亚洲女在线精品 一区二区美女 中文字幕在线2021一区 在线小视频国产 久久这里只有精品首页 国产在线第三页 欧美日韩中文字幕 在线亚洲+欧美+日本专区 精品国产一区二区三区不卡 久久这里精品 欧美在线va在线播放 精液呈暗黄色 91精品国产手机 91在线免费播放 欧美视频亚洲色图 欧美国产日韩精品 日韩高清不卡在线 精品视频免费观看 欧美日韩一区二区三区四区 国产欧美亚洲精品第二区首页 亚洲韩精品欧美一区二区三区 国产精品视频免费 在线精品小视频 久久午夜夜伦伦鲁鲁片 国产无套在线播放 久热这里只精品99re8久 欧美久久久久 久久香蕉国产线看观看精品蕉 国产成人精品男人的天堂538 亚洲人成网站色7799在线观看 日韩在线第二页 一本色道久久综合狠狠躁篇 国产一区二区三区不卡在线观看 亚洲乱码在线 在线观看欧美国产 久久福利青草精品资源站免费 国产玖玖在线观看 在线亚洲精品 亚洲成aⅴ人在线观看 精品91在线 欧美一区二三区 日韩中文字幕视频在线 日本成人一区二区 日韩免费专区 国内精品在线观看视频 久久国产综合尤物免费观看 国产精品系列在线观看 一本一道久久a久久精品综合 亚洲免费播放 久久精品国产免费 久久人精品 亚洲毛片网站 亚洲成a人一区二区三区 韩国福利一区二区三区高清视频 亚洲精品天堂在线 一区二区三区中文字幕 亚洲国产色婷婷精品综合在线观看 亚洲国产成人久久笫一页 999国产视频 国产精品香港三级在线电影 欧美日韩一区二区三区四区 日韩国产欧美 国产精品99一区二区三区 午夜国产精品理论片久久影院 亚洲精品中文字幕麻豆 亚洲国产高清视频 久久免费手机视频 日韩a在线观看 五月婷婷亚洲 亚洲精品中文字幕麻豆 中文字幕丝袜 www国产精品 亚洲天堂精品在线 亚洲乱码一区 国产日韩欧美三级 久久999精品 伊人热人久久中文字幕 久热国产在线视频 国产欧美日韩在线观看一区二区三区 国产一二三区在线 日韩国产欧美 91精品国产91久久久久 亚洲一区小说区中文字幕 精品一区二区免费视频 国产精品视频免费 国产精品亚洲综合色区韩国 亚洲国产精品成人午夜在线观看 欧美国产日韩精品 中文字幕精品一区二区精品