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

    面向線性光譜混合分解的鄰域像元集螺線型構(gòu)建方法

    2017-12-05 05:43:45劉博宇邢華橋
    測繪學報 2017年11期
    關鍵詞:螺線鄰域分辨率

    劉博宇,陳 軍,邢華橋,3,武 昊,張 俊

    1.吉林大學地球探測科學與技術(shù)學院,吉林 長春 130026; 2.國家基礎地理信息中心,北京 100830; 3.中國礦業(yè)大學(北京)地球科學與測繪工程學院,北京 100083

    面向線性光譜混合分解的鄰域像元集螺線型構(gòu)建方法

    劉博宇1,2,陳 軍2,邢華橋2,3,武 昊2,張 俊2

    1.吉林大學地球探測科學與技術(shù)學院,吉林 長春 130026; 2.國家基礎地理信息中心,北京 100830; 3.中國礦業(yè)大學(北京)地球科學與測繪工程學院,北京 100083

    高時間分辨率遙感影像在地表景觀破碎區(qū)域易形成混合像元,難以發(fā)揮其高時間維度優(yōu)勢?,F(xiàn)有方式多是基于線性光譜混合模型,借助鄰域像元所構(gòu)成的像元集合組成線性方程組,求出組分光譜值的最小二乘解,提高其空間分辨率。然而,現(xiàn)有方法依賴窗口形式來構(gòu)建鄰域像元集合,在某些區(qū)域易造成方程組無解的欠定問題。本文在分析其問題原因的基礎上,引入阿基米德螺線代替?zhèn)鹘y(tǒng)的矩形窗口,對鄰域各像元依次遍歷,構(gòu)建空間鄰近、組分相近的鄰域像元集合來解決該問題。在GlobeLand 30數(shù)據(jù)上的試驗表明,螺線型構(gòu)建方法對5種混合尺度上多種類型地物均具有穩(wěn)定的精度,與傳統(tǒng)窗口構(gòu)建方法相比,可從構(gòu)建鄰域像元集合方面將總體理論精度提高2%,分解結(jié)果精度提高近1個數(shù)量級。

    線性光譜混合分解;混合像元;鄰域像元;螺線;窗口

    環(huán)境變化監(jiān)測、地表覆蓋數(shù)據(jù)更新、土地資源動態(tài)管理等迫切需要兼具高時間和高空間分辨率的遙感影像[1-3],但限于運行軌道、硬件等因素制約,遙感傳感器還難以同時滿足該高分辨率需求[4-5]。國內(nèi)外通常是借助高空間分辨率的地表分類圖,對多期低空間分辨率影像分別進行光譜分解,進而組成高空間分辨率影像序列[6-9]。設低空間分辨率像元P的光譜值為SP,對應高空間分辨率地物組分為Ci(i=1,2,…,n),其光譜值為SC,光譜分解是將SP看作由n個SC混合而成,基于表征兩者之間關系的光譜混合模型,由SP推算出對應的SC[10-12]。常用的線性光譜混合分解方式是基于地理學第一定律。假設鄰近地物有較強相關性,取當前和鄰近的m個像元,構(gòu)建其光譜值SP與組分光譜值SC之間的線性方程組L,在觀測數(shù)(像元數(shù))m大于等于未知數(shù)(組分數(shù))n的條件下,按最小二乘法求解出當前像元P內(nèi)的SC[13-14]。然而,當鄰近地物類別相差較大時,鄰域內(nèi)某些像元的組分明顯不同于當前像元,屬于無效觀測,這往往致使L在觀測數(shù)m大于未知數(shù)n時仍然引發(fā)導致SC無解的欠定問題。組成非欠定的光譜分解線性方程組,需確保所選擇的鄰域像元均為有效觀測,鄰域范圍像元集合的構(gòu)建是解決此問題的關鍵。

    通常,鄰域像元集主要采用規(guī)則窗口來構(gòu)建[15]。例如,文獻[16]直接使用以待分解像元為中心的窗口來收集鄰域像元;在文獻[16]的基礎上,文獻[17]通過結(jié)合空間距離權(quán)重,減小窗口中較遠處的像元對中心待分解像元組分值的影響;文獻[18]通過增加光譜權(quán)重,進一步改善了窗口內(nèi)部參與解算的像元種類,以此來提高分解精度;文獻[19]提出使用限制性最小二乘法來保障組分值位于合理區(qū)間;文獻[4]采用了一個移動窗口來收集周圍鄰域像元,來處理影像場景內(nèi)的不同地方的光譜可變問題,并通過減少解算的部分類數(shù),使得窗口中像元數(shù)大于組分數(shù),來保障順利分解;文獻[20]通過對比不同大小的窗體下的計算精度,先驗性地確定最佳窗口,構(gòu)建鄰域像元集。窗口型的鄰域像元集構(gòu)建方法原理簡潔而易于實現(xiàn),其效率較高,獲得了廣泛的應用。然而,這些在窗口基礎上作出的改進,以及在其后的分解中所采取的保障措施并不是構(gòu)建鄰域像元集的根本解決方案。究其原因,是由于窗口的覆蓋形狀難以順應地物分布走勢,在部分區(qū)域不能從鄰域內(nèi)界定出各組分相近像元分布的準確范圍,從而引發(fā)欠定問題。

    螺線是一種具有非規(guī)則空間形態(tài)的幾何曲線[22-23]。相較于窗口,螺線更加貼近鄰域周圍的非規(guī)則地表。沿其展開的軌跡,根據(jù)當前像元內(nèi)的組分分布,對其四周地物類別相近的像元范圍可進行最佳逼近,可構(gòu)建空間相鄰、組分相近的鄰域像元集合,有望在一定程度上解決欠定難題。本文第1節(jié)分析線性光譜分解理論將混合像元分解為組分值時面臨的問題,揭示傳統(tǒng)窗口鄰域像元集構(gòu)建方式的缺陷,提出基于螺線構(gòu)建非欠定的鄰域像元集合的解決思路;第2節(jié)介紹螺線在二維笛卡爾坐標系中的參數(shù)公式,推導出其在影像柵格空間上的離散化形式,進而提出螺線形式的鄰域像元集合構(gòu)建算法;第3節(jié)以GlobeLand 30地表覆蓋數(shù)據(jù)產(chǎn)品為基礎,對比不同像元混合尺度上螺線型方法與窗口方法的試驗結(jié)果;最后,總結(jié)并指出有待進一步研究的問題。

    1 相似地物的螺線逼近

    1.1 問題分析

    在線性光譜混合模型中,將混合像元的值R表達為該像元內(nèi)部各種地類(組分)值r與各地類所占面積比例(豐度)f的線性組合,如式(1)(n為組分總數(shù);ε為非線性殘差)?;谝延械牡乇矸诸悎D提取得到豐度后,采用線性光譜混合模型將混合像元分解為組分值,即求解以r為未知數(shù)的式(1)的過程。

    (1)

    按方程組求解的一般要求,方程數(shù)m應大于等于未知數(shù)的數(shù)量n(式(2)),此時屬于適定、超定問題,可用最小二乘方法求最優(yōu)解;當m小于n時,會發(fā)生無唯一解的欠定問題,求不出未知數(shù)。對線性光譜分解過程而言,即需列出滿足式(2)的多個線性方程,組成線性方程組(如式(3))。為此,通常是依據(jù)地理學第一定律的觀點,收集待分解像元空間附近的其他像元共同組成鄰域像元集合[13,20],以構(gòu)成待分解像元的解算線性方程組[24-26]。

    m≥n

    (2)

    (3)

    鄰域像元集的窗口構(gòu)建方式,就是收集一定矩形空間范圍內(nèi)的像元,組成鄰域像元集合。窗口的數(shù)學形式如式(4),大小為2s+1(s為自然數(shù),一般以像元為單位),l為其所能解算的最大組分數(shù)。已知地表分類圖的l為確定值,由式(4)即可計算窗口大小。由于地物類型多樣,分布形式復雜,鄰近地物的相關性難以一致,使得鄰域中的其他像元與待分解像元組分類型構(gòu)成并不完全相同。這些像元增大了中心像元方程組的方程數(shù)m,使得常出現(xiàn)方程數(shù)m遠多于組分數(shù)n,但仍然解不出組分值的情況。規(guī)則的窗口不能順應地物的分布形狀,造成其鄰域像元集合中常出現(xiàn)與待分解中心像元組分無關的其他像元,增加了解算方程組中無效方程的規(guī)模,給組分值的解算帶來困難;部分均質(zhì)像元因受到其他像元多種組分的影響,增加了分解其組分光譜值的復雜性。圖1示例列舉了窗口構(gòu)建鄰域像元集的4類典型情況:①像元A包含森林、草地、耕地3種組分類型,窗口代入了人造地表組分類型,其鄰域像元集組分與豐度分布如表1-1所示,以豐度為系數(shù)的解算方程組中混入了無關光譜,增加了方程組的分解難度;②像元B包含水體、濕地、耕地3種組分類型,其窗口代入了含有無關組分的像元d5,其對應的線性方程組無法分解d5內(nèi)的組分,如表1-2所示;③像元C包含5種組分類型,其窗口可收集6個像元(即m為6),其鄰域像元集總共包含6種組分類型(即n為6),但由于其周圍像元均不包含人造和裸地類型的光譜,造成對應方程組無法分解這兩類光譜值,如表1-3所示;④像元D僅包含耕地一種組分類型,其窗口引入了包括人造、草地其他無關組分的像元。

    (4)

    表1 圖1中鄰域像元集組分及豐度分布一覽表(子表1-1、1-2、1-3、1-4分別對應A、B、C、D的鄰域像元集)Tab.1 The endmembers and abundance of adjacent pixel sets in fig.1(sub-table 1-1,1-2,1-3,1-4 correspond to A,B,C,D) (%)

    1?31?4序號像元編碼耕地森林草地人造覆蓋裸地水體序號像元編碼耕地人造覆蓋草地1C6095141201D100002a85423600172g793073a10801830003g8771224b832253100124g993525b96912120085h798206b10871030006h982180

    1.2 螺線逼近

    阿基米德螺線(本文簡稱螺線)是一種沿中心向外逐漸旋轉(zhuǎn)展開而形成的曲線,呈現(xiàn)出非規(guī)則的幾何形態(tài)(圖2)。其形態(tài)更加貼近非規(guī)則的鄰域地表(圖3),沿其幾何軌跡可對鄰域內(nèi)周圍分布有多種地物的地表進行逼近,適于對各種空間尺度大小的對象進行局部遍歷,以螺線中心為起點,沿其展開軌跡,順次遍歷其周圍各空間位置的對象。這種由近及遠地遍歷周圍對象的方式,較符合地理學第一定律關于空間鄰近地物相關的假設對選擇鄰近地物的要求。螺線的極坐標形式如式(5),其中a為極徑,極角θ屬于[0,2π]區(qū)間。

    ρ=aθ,a>0

    (5)

    圖2 地表極坐標上的螺線Fig.2 Polar Archimedean spiral in land surface

    在與圖1相同背景下,圖3示例采用螺線對像元A、B、C、D周圍順次遍歷,在較近的區(qū)域選擇出組分類型相近的像元。結(jié)合表2可知,在像元A周圍,其構(gòu)建的鄰域像元集包含了均有耕地、林地、草地3種組分類型的A和a1、a2、b1、c1共5個像元,即其線性方程組的觀測方程均為有效觀測,觀測數(shù)m為5,組分值未知數(shù)n為3,較窗口方法構(gòu)建的線性方程組規(guī)模更小,易于解算;在像元B處,螺線構(gòu)建的像元集避開了像元d5的影響,利用該方程組可更加高效地分解出其內(nèi)的耕地、水體、濕地組分值;對像元C構(gòu)建的鄰域像元集還包含了該區(qū)域內(nèi)最鄰近的a10、b10、c10、c9共4個組分相關像元,其線性方程組能有效地解算出像元C的組分值;由于像元D內(nèi)部僅含有一種耕地組分類型,螺線構(gòu)建的鄰域像元集僅包含一個D像元。

    圖1 窗口構(gòu)建鄰域像元集示例(該窗口大小為3,A、B、C、D為待分解像元,背景為地表分類圖)Fig.1 A window-based construction example(window size is 3,A,B,C,D are low spatial pixel,background is a class map)

    圖3 螺線對地表各地物類別的逼近示例(沿螺線軌跡選定的其他像元用黑點標記)Fig.3 An example of fitting similar pixels with Archimedean spiral in land surface (the similar pixels labeled by black dot)

    表2 圖3中螺線逼近的鄰域像元集組分及豐度分布一覽表(A、B、C、D分別對應表2-1、2-2、2-3、2-4)Tab.2 The endmember and abundance of spiral-based adjacent pixel sets in Fig.3(sub-table 2-1,2-2,2-3,2-4 correspond to A,B,C,D) (%)

    2?32?4序號像元編碼耕地森林草地人造覆蓋裸地序號像元編碼耕地1C609514121D1002a1080183003b1087103004c108710575c9441915815

    1.3 螺線型構(gòu)建思路

    在異質(zhì)地表上,顧及地理學第一定律關于鄰近地物相關性強的假設,以及光譜分解線性方程組的解算要求,生成柵格化的螺線軌跡,在當前像元的周圍鄰域內(nèi)進行遍歷,并判別像元組分的相關性,發(fā)展螺線型的鄰域像元集構(gòu)建算法,以此來構(gòu)建空間鄰近、組分相近的鄰域像元集,是解決欠定問題的可取途徑。

    圖4 基于螺線的領域像元集構(gòu)建思路Fig.4 Strategy of constructing adjacent pixel sets with Archimedean spiral

    2 鄰域像元集的螺線構(gòu)建

    2.1 螺線軌跡的柵格化生成算法

    遙感影像是以矩陣組織像元的柵格格式,研究螺線在二維笛卡爾坐標系下的離散化形式形式,有助于螺線在影像處理中的應用。在二維笛卡爾坐標系下,螺線極坐標公式的參數(shù)形式為式(6),a=vt,θ=wt,v為線速度,w為角速度。

    (6)

    由式(6)可推導出螺線在矩陣空間上的離散化形式,如式(7),r為半徑,f1(θ)、f2(θ)為函數(shù),如式(8),θ為由0到2π周期性變化的角度,其定義域如式(8)和(9)所示,r是自然數(shù),k為整數(shù)。由該離散化的螺線公式,以一個像元為單位,以某像元中心為螺線的原點, 沿其幾何軌跡順次確定出周圍各像元的中心位置, 取其上離散化的整數(shù)為坐標值,周圍鄰域像元如圖5(b)所示。

    圖5 螺線在二維笛卡爾坐標系向矩陣空間轉(zhuǎn)換示例Fig.5 Archimedean spiral conversion from two dimensional Cartesian coordinate system to matrix system

    (7)

    (8)

    2.2 基于螺線軌跡的鄰域像元集生成

    在低空間分辨率影像上,采用螺線構(gòu)建鄰域像元集,需按低空間分辨率影像與該區(qū)域地表分類圖的空間位置,提取各像元范圍內(nèi)的組分類型C、組分數(shù)n和豐度f。螺線型算法的原理如圖6所示,設鄰域像元集中各像元豐度組成的矩陣為F,其秩記為rankF,其形式化描述如下:

    第1步,在影像柵格上,建立每個像元P的鄰域像元集SP,并將P加入SP。

    第2步,以P為螺線中心,沿離散化的螺線順次遍歷P鄰域內(nèi)的其他像元P′。

    第3步,判別P′內(nèi)的組分類型C′,若均為相似組分則將其加入SP,進入下一步驟;否則,繼續(xù)按當前的次序遍歷。

    第4步,判斷當前鄰域像元集的豐度矩陣的秩rankF是否等于待分解像元的組分數(shù)n,即rankF=n;若滿足,則完成構(gòu)建鄰域像元集SP過程;否則,按當前次序順次加1,繼續(xù)第2步過程。

    圖6 螺線構(gòu)建算法原理Fig.6 Archimedean spiral-based adjacent pixel sets construction approach

    3 試驗與分析

    3.1 試驗區(qū)與數(shù)據(jù)

    采用30 m空間分辨率的2010期GlobeLand 30地表覆蓋產(chǎn)品作為分類圖,按混合像元分別為4、6、8、10和12倍于分類像元的混合尺度關系,檢驗了螺線型方法對混合像元構(gòu)建領域像元集與分解時的精度及效率,并與窗口型方法相應結(jié)果對比。試驗區(qū)是面積約為1300 km2的北京地區(qū),其內(nèi)部廣泛分布有耕、林、草、灌叢等共8種地物類型(圖7)。各尺度的混合像元值按式(1)殘差為零,根據(jù)分類圖生成,值域為[0.1,0.9];由式(4)計算,窗口大小為3;螺線鄰域像元集像元數(shù)規(guī)模設為中心像元組分數(shù)的2倍。

    圖7 GlobeLand 30北京區(qū)域地表覆蓋圖Fig.7 GlobeLand 30 class map in Beijing area

    3.2 試驗結(jié)果與分析

    對試驗區(qū)主要地表覆蓋上的構(gòu)建統(tǒng)計結(jié)果如圖8所示??傮w上,人造地表、水體與濕地、耕地、林地與灌叢分布區(qū)域的欠定像元比例,隨混合尺度變大而逐漸增加;窗口型方法存在不同程度的欠定問題,且其比例隨著混合尺度呈線性增加趨勢;螺線型方法在各尺度上基本持平,均處于較低水平,說明其對多種尺度的混合像元均具有穩(wěn)定的解算像元集合構(gòu)建能力,可顯著減少欠定問題的發(fā)生。綜合5種混合尺度的欠定像元統(tǒng)計結(jié)果,窗口方法在構(gòu)建非欠定的線性解算方程組方面精度約為97%,而螺線方法達到99%以上。

    基于所構(gòu)建的鄰域像元集,對5種混合尺度的低空間分辨率像元值進行光譜分解,得出的部分高空間分辨率灰度圖像如圖9所示??梢?,隨著混合尺度增加,窗口型方法分解結(jié)果中越來越多像元值的誤差超出了辨識不同地物類別的界限(試驗中各地物類別最小光譜差值為0.1,絕對誤差超出其一半即意味著難以提供準確的光譜信息);螺線型方法分解效果較好。通過統(tǒng)計試驗區(qū)分解灰度圖像絕對誤差的絕對值,顯示螺線型方法均值處于[4.4×10-4,1.3×10-3],標準偏差處于[0.4,1.2];窗口型方法相應指標分別為[3.1×10-3,131.9]和[5.3×10-2,1.2×105],表明螺線型方法對分解結(jié)果精度提高作用顯著。水體、濕地類型廣泛地與其他類型地物交錯分布且數(shù)量較少,造成鄰域像元集構(gòu)建及分解難度較大;人造地表類型分布相對集中,具有較規(guī)則的形狀,對這些區(qū)域的混合像元構(gòu)建鄰域像元集及分解的難度較??;灌叢地類型零星分布于林地周圍,易造成欠定問題,但由于面積較大,其對整體精度的影響適中。通過對比兩種方法在5種混合尺度上構(gòu)建鄰域像元集及分解的精度可知,螺線型構(gòu)建方法在處理多種地表類型、復雜的地類分布上具有優(yōu)勢。

    圖8 試驗區(qū)主要地表覆蓋類型總體欠定像元比例Fig.8 The whole underdetermined pixels ratio among class map

    對每個混合像元及試驗區(qū)總體圖幅構(gòu)建鄰域像元集合與進行混合像元分解所需的時間如圖10所示。結(jié)果表明,螺線型方法構(gòu)建鄰域像元集合所需時間是窗口型方法的1.63倍,基于其構(gòu)建的鄰域像元集合進行解混,所需時間是窗口型方法的2.16倍,其構(gòu)建及分解效率均低于窗口型方法。

    圖9 部分分解灰度圖像(分解值的絕對偏差大于0.05的像元標記為黑白方塊)Fig.9 The grey unmixed images of the test area (unmixed pixels that absolute deviation of values greater than 0.05 are marked as black and white block)

    4 結(jié) 論

    利用線性光譜混合模型分解低空間分辨率的像元光譜值時,在部分地表區(qū)域會發(fā)生解算方程組欠定問題,解不出組分的光譜值。傳統(tǒng)基于窗口的鄰域像元集構(gòu)建方式難以解決該問題,其原因在于:一是規(guī)則的矩形區(qū)域難以應對地表上多種地物分布形態(tài),使該構(gòu)建方式難以準確獲得相似的像元;二是其構(gòu)建過程缺少對像元組分類型的判別,其鄰域像元集未達到線性方程解算要求。螺線作為一種非規(guī)則的幾何曲線,其形態(tài)更加貼近非規(guī)則的地表?;诼菥€柵格化形式的鄰域像元集構(gòu)建方法,在鄰域內(nèi)遍歷到足夠的相似像元時,可在異質(zhì)性不同的各地表區(qū)域構(gòu)建出空間鄰近、組分相似的鄰域像元集,即可組建非欠定的線性方程組。本質(zhì)上,螺線與窗口方式均是通過逼近混合像元的鄰域,來獲取其周圍像元組成鄰域像元集,其區(qū)別在于用來逼近鄰域地表的是面狀的窗口還是線狀的螺線。螺線以一種細粒度的方式對分布于鄰域地表的像元進行由內(nèi)而外,并顧及四周的順序遍歷方式,這是面狀的窗口所不具備的。鄰域像元多是以非規(guī)則的分布方式位于混合像元周圍,因此,細粒度的螺線較粗粒度的規(guī)則窗口在擬合鄰域地表方面更加靈活。通過結(jié)合對遍歷的各像元進行取舍,構(gòu)建的鄰域像元集可極大地避免欠定問題,進而提高分解精度。

    在北京區(qū)域GlobeLand 30地表覆蓋數(shù)據(jù)上的試驗表明,相較于傳統(tǒng)窗口構(gòu)建方法,螺線型構(gòu)建方法將欠定像元比例減少了近兩個數(shù)量級,從構(gòu)建鄰域像元集合方面可將總體理論精度提高2%;基于該方法的混合像元分解結(jié)果精度提高近一個數(shù)量級;其構(gòu)建鄰域像元集及進而分解所需時間約為原窗口構(gòu)建方法的2倍。在本次試驗中,所用螺線并未以不同旋轉(zhuǎn)方向逼近實際的地表分類圖,但應注意到,其在鄰域上的旋轉(zhuǎn)方式可能對構(gòu)建及分解結(jié)果產(chǎn)生影響。同時,螺線型構(gòu)建方法的作用主要是逼近地物在地表分布形式,當?shù)乇韮煞N以上的地類僅存在于極少稀疏分布的混合像元中時,該方法仍難以成功構(gòu)建鄰域像元集;在其他地表覆蓋數(shù)據(jù)產(chǎn)品和幾何光學、概率等其他光譜混合模型中的應用效果還不明確,這些問題仍需深入研究。

    [1] 陳軍,陳晉,廖安平,等.全球30 m地表覆蓋遙感制圖的總體技術(shù)[J].測繪學報,2014,43(6):551-557.DOI:10.13485/j.cnki.11-2089.2014.0089.

    CHEN Jun,CHEN Jin,LIAO Anping,et al.Concepts and Key Techniques for 30 m Global Land Cover Mapping[J].Acta Geodaetica et Cartographica Sinica,2014,43(6):551-557.DOI:10.13485/j.cnki.11-2089.2014.0089.

    [2] CHEN Xuehong,YANG Dedi,CHEN Jin,et al.An Improved Automated Land Cover Updating Approach by Integrating with Downscaled NDVI Time Series Data[J].Remote Sensing Letters,2015,6(1):29-38.

    [3] CHEN Jun,CHEN Jin,LIAO Anping,et al.Global Land cover Mapping at 30 m Resolution:A POK-based Operational Approach[J].ISPRS Journal of Photogrammetry and Remote Sensing,2015,103(1):7-27.

    [4] GEVAERT C M,GARCíA-HARO F J.A Comparison of STARFM and an Unmixing-based Algorithm for Landsat and MODIS Data Fusion[J].Remote Sensing of Environment,2015,156(1):34-44.

    [5] GAO Feng,MASEK J,SCHWALLER M,et al.On the Blending of the Landsat and MODIS Surface Reflectance:Predicting Daily Landsat Surface Reflectance[J].IEEE Transactions on Geoscience and Remote Sensing,2006,44(8):2207-2218.

    [6] PERERA K,TSUCHIYA K.Experiment for Mapping Land Cover and It’s Change in Southeastern Sri Lanka Utilizing 250 m Resolution MODIS Imageries[J].Advances in Space Research,2009,43(9):1349-1355.

    [7] SAVTCHENKO A,OUZOUNOV D,AHMAD S,et al.Terra and Aqua MODIS Products Available from NASA GES DAAC[J].Advances in Space Research,2004,34(4):710-714.

    [8] 吳波,熊助國.基于光譜最佳尺度分割特征的高光譜混合像元分解[J].測繪學報,2012,41(2):205-212.

    WU Bo,XIONG Zhuguo.Unmixing of Hyperspectral Mixture Pixels Based on Spectral Multiscale Segemented Features[J].Acta Geodaetica et Cartographica Sinica,2012,41(2):205-212.

    [9] ASNER G P,WESSMAN C A,PRIVETTE J L.Unmixing the Directional Reflectances of AVHRR Sub-pixel Landcovers[J].IEEE Transactions on Geoscience and Remote Sensing,1997,35(4):868-878.

    [10] 李慧,張金區(qū),曹陽,等.端元可變非線性混合像元分解模型[J].測繪學報,2016,45(1):80-86.DOI:10.11947/j.AGCS.2016.20140520.

    LI Hui,ZHANG Jinqu,CAO Yang,et al.Nonlinear Spectral Unmixing for Optimizing Per-pixel Endmember Sets[J].Acta Geodaetica et Cartographica Sinica,2016,45(1):80-86.DOI:10.11947/j.AGCS.2016.20140520.

    [11] 張良培,沈煥鋒.遙感數(shù)據(jù)融合的進展與前瞻[J].遙感學報,2016,20(5):1050-1061.

    ZHANG Liangpei,SHEN Huanfeng.Progress and Future of Remote Sensing Data Fusion[J].Journal of Remote Sensing,2016,20(5):1050-1061.

    [12] 鄔明權(quán),王潔,牛錚,等.融合MODIS與Landsat數(shù)據(jù)生成高時間分辨率Landsat數(shù)據(jù)[J].紅外與毫米波學報,2012,31(1):80-84.

    WU Mingquan,WANG Jie,NIU Zheng,et al.A Model for Spatial and Temporal Data Fusion[J].Journal of Infrared and Millimeter Waves,2012,31(1):80-84.

    [13] RAO Yuhan,ZHU Xiaolin,CHEN Jin,et al.An Improved Method for Producing High Spatial-resolution NDVI Time Series Datasets with Multi-temporal MODIS NDVI Data and Landsat TM/ETM+ Images[J].Remote Sensing,2015,7(6):7865-7891.

    [14] 陳晉,馬磊,陳學泓,等.混合像元分解技術(shù)及其進展[J].遙感學報,2016,20(5):1102-1109.

    CHEN Jin,MA Lei,CHEN Xuehong,et al.Research Progress of Spectral Mixture Analysis[J].Journal of Remote Sensing,2016,20(5):1102-1109.

    [15] 何青,熊楨.基于鄰域的混合象元解譯[J].北方交通大學學報,1999,23(4):118-121.

    HE Qing,XIONG Zhen.Neighbor-field-based Mixed Pixel Interpretation[J].Journal of Northern Jiaotong University,1999,23(4):118-121.

    [16] ZHUKOV B,OERTEL D,LANZL F,et al.Unmixing-based Multisensor Multiresolution Image Fusion[J].IEEE Transactions on Geoscience and Remote Sensing,1999,37(3):1212-1226.

    [17] MASELLI F.Definition of Spatially Variable Spectral Endmembers by Locally Calibrated Multivariate Regression Analyses[J].Remote Sensing of Environment,2001,75(1):29-38.

    [18] BUSETTO L,MERONI M,COLOMBO R.Combining Medium and Coarse Spatial Resolution Satellite Data to Improve the Estimation of Sub-pixel NDVI Time Series[J].Remote Sensing of Environment,2008,112(1):118-131.

    [19] ZURITA-MILLA R,CLEVERS J G P W,SCHAEPMAN M E.Unmixing-based Landsat TM and MERIS FR Data Fusion[J].IEEE Geoscience and Remote Sensing Letters,2008,5(3):453-457.

    [20] 謝登峰,張錦水,孫佩軍,等.結(jié)合像元分解和STARFM模型的遙感數(shù)據(jù)融合[J].遙感學報,2016,20(1):62-72.

    XIE Dengfeng,ZHANG Jinshui,SUN Peijun,et al.Remote Sensing Data Fusion by Combining STARFM and Downscaling Mixed Pixel Algorithm[J].Journal of Remote Sensing,2016,20(1):62-72.

    [21] 張建,李宗南,謝靜.亞像元定位之鄰域像元選取問題研究[J].華中師范大學學報(自然科學版),2013,47(3):404-409.

    ZHANG Jian,LI Zongnan,XIE Jing.Research on Neighbor Pixels Choosing Issue of the Sub-pixel Mapping[J].Journal of Huazhong Normal University (Natural Sciences),2013,47(3):404-409.

    [22] 王明華,楊繼緒.阿基米德螺線的性質(zhì)與應用[J].數(shù)學通報,1989(7):11-12.

    WANG Minghua,YANG Jixu.The Properties and Application of Archimedean Spiral[J].Bulletin des Sciences Mathematics,1989(7):11-12.

    [23] 姚建武.螺線與生物體上的擬螺線[J].科學,2004,56(4):61-63.

    YAO Jianwu.On Spiral and Organism Draft Spiral Research[J].Science,2004,56(4):61-63.

    [24] XU Mingming,ZHANG Liangpei,DU Bo.An Image-based Endmember Bundle Extraction Algorithm Using both Spatial and Spectral Information[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2015,8(6):2607-2617.

    [25] DU Bo,ZHANG Liangpei.A Discriminative Metric Learning Based Anomaly Detection Method[J].IEEE Transactions on Geoscience and Remote Sensing,2014,52(11):6844-6857.

    [26] DU Bo,ZHANG Liangpei.Target Detection Based on a Dynamic Subspace[J].Pattern Recognition,2014,47(1):344-358.

    A Spiral-based Construction of Adjacent Pixel Sets for Linear Spectral Unmixing

    LIU Boyu1,2,CHEN Jun2,XING Huaqiao2,3,WU Hao2,ZHANG Jun2

    1.College of Geo-exploration Science and Technology,Jilin University,Changchun 130026,China; 2.National Geomatics Center of China,Beijing 100830,China; 3.Geoscience and Surveying Engineering College,China University of Mining and Technology (Beijing),Beijing 100083,China

    The problem of spectral mixing in fragmented landscape prevents the potentiality of high temporal resolution image from land surface detail dynamic monitoring.The general way is unmixing the spectrum to components based on linear spectral mixing model,with the aid of neighborhood pixels and the least squares method.However,constructing the neighborhood pixel set by a window leads to the underdetermined problem in some areas.This paper analyses the cause of the problem,and introduces spiral to construct optimal neighborhood pixel set as a solution.Experiment on GlobeLand 30 shows that the spiral method has a good applicability for each land surface type.Compared with the traditional method,the overall accuracy can be enhanced by 2%,the accuracy of unmixed results can be improved by nearly one order of magnitude.

    linear spectral unmixing; mixed pixel; adjacent pixel; spiral; window

    The State Key Program of National Natural Science Foundation of China (No.41231172)

    LIU Boyu(1987—),male,PhD candidate,majors in spatio-temporal data fusion and land cover mapping.

    劉博宇,陳軍,邢華橋,等.面向線性光譜混合分解的鄰域像元集螺線型構(gòu)建方法[J].測繪學報,2017,46(11):1841-1849.

    10.11947/j.AGCS.2017.20170109.

    LIU Boyu,CHEN Jun,XING Huaqiao,et al.A Spiral-based Construction of Adjacent Pixel Sets for Linear Spectral Unmixing[J].Acta Geodaetica et Cartographica Sinica,2017,46(11):1841-1849.DOI:10.11947/j.AGCS.2017.20170109.

    P237

    A

    1001-1595(2017)11-1841-09

    國家自然科學基金重點項目(41231172)

    (責任編輯:叢樹平)

    2017-03-09

    修回日期:2017-07-12

    劉博宇(1987—),男,博士生,研究方向為遙感時空數(shù)據(jù)融合與地表覆蓋制圖。

    E-mail:liuby10@mails.jlu.edu.cn

    猜你喜歡
    螺線鄰域分辨率
    三維Minkowski空間中的k-型偽零螺線
    稀疏圖平方圖的染色數(shù)上界
    EM算法的參數(shù)分辨率
    走近等角螺線
    探秘等角螺線
    走近等角螺線
    原生VS最大那些混淆視聽的“分辨率”概念
    基于鄰域競賽的多目標優(yōu)化算法
    自動化學報(2018年7期)2018-08-20 02:59:04
    基于深度特征學習的圖像超分辨率重建
    自動化學報(2017年5期)2017-05-14 06:20:52
    一種改進的基于邊緣加強超分辨率算法
    欧美在线黄色| 老司机午夜十八禁免费视频| 校园春色视频在线观看| 国产伦在线观看视频一区| 国产高潮美女av| 九九久久精品国产亚洲av麻豆| 国内毛片毛片毛片毛片毛片| 久久精品国产99精品国产亚洲性色| 午夜视频国产福利| 日韩欧美 国产精品| 国产91精品成人一区二区三区| 亚洲av熟女| 99久久精品国产亚洲精品| 欧美一级a爱片免费观看看| 精品福利观看| 男人舔女人下体高潮全视频| 国产精品久久久久久人妻精品电影| 色视频www国产| 国产精品野战在线观看| 欧美不卡视频在线免费观看| 欧美绝顶高潮抽搐喷水| 欧美一区二区精品小视频在线| 99国产精品一区二区三区| 91狼人影院| 免费看光身美女| 波野结衣二区三区在线| 色哟哟哟哟哟哟| 日本 av在线| 亚洲内射少妇av| 日本 欧美在线| 日韩 亚洲 欧美在线| 校园春色视频在线观看| 天堂av国产一区二区熟女人妻| 亚洲va日本ⅴa欧美va伊人久久| 国产伦精品一区二区三区四那| 国产精品不卡视频一区二区 | 亚洲人与动物交配视频| 亚洲精品在线观看二区| 国产高潮美女av| 午夜福利在线观看免费完整高清在 | 一区二区三区四区激情视频 | 女人被狂操c到高潮| 波多野结衣巨乳人妻| 国内少妇人妻偷人精品xxx网站| 天堂动漫精品| 亚洲人与动物交配视频| 免费av观看视频| 日韩欧美 国产精品| 少妇的逼好多水| 国产淫片久久久久久久久 | 搡女人真爽免费视频火全软件 | 一区福利在线观看| 女生性感内裤真人,穿戴方法视频| 国产精品日韩av在线免费观看| 国产av不卡久久| 亚洲av成人av| 久久婷婷人人爽人人干人人爱| 亚洲精品在线观看二区| 成人毛片a级毛片在线播放| 好男人电影高清在线观看| 国产极品精品免费视频能看的| 国产亚洲精品久久久久久毛片| 99视频精品全部免费 在线| 国产免费一级a男人的天堂| 国产免费av片在线观看野外av| 男女视频在线观看网站免费| 极品教师在线视频| 精品久久久久久久久久久久久| 成年女人毛片免费观看观看9| 人妻久久中文字幕网| 欧美中文日本在线观看视频| 欧美黄色淫秽网站| 亚洲精品456在线播放app | 成年女人看的毛片在线观看| 亚洲av成人av| 麻豆av噜噜一区二区三区| 欧美高清成人免费视频www| 此物有八面人人有两片| 日本免费a在线| 女人被狂操c到高潮| 成人永久免费在线观看视频| 搡老岳熟女国产| 日本a在线网址| 极品教师在线视频| 久久精品91蜜桃| 国产精品一区二区三区四区久久| 少妇人妻精品综合一区二区 | 国产成年人精品一区二区| 亚洲成av人片免费观看| 色综合婷婷激情| av中文乱码字幕在线| 日韩欧美免费精品| 成人性生交大片免费视频hd| 亚洲成a人片在线一区二区| 51国产日韩欧美| 一边摸一边抽搐一进一小说| 精品一区二区三区视频在线| 精品久久久久久成人av| 脱女人内裤的视频| 在线观看舔阴道视频| 日日摸夜夜添夜夜添小说| 淫妇啪啪啪对白视频| 在线国产一区二区在线| 欧美日韩国产亚洲二区| 精品一区二区三区视频在线观看免费| 亚洲成人精品中文字幕电影| 此物有八面人人有两片| 国产一级毛片七仙女欲春2| 变态另类成人亚洲欧美熟女| 国产亚洲精品综合一区在线观看| 亚洲精品亚洲一区二区| 91在线精品国自产拍蜜月| 精品日产1卡2卡| 伦理电影大哥的女人| 国产精品久久久久久亚洲av鲁大| 色综合亚洲欧美另类图片| 波多野结衣高清作品| 美女xxoo啪啪120秒动态图 | 亚洲在线观看片| 免费在线观看日本一区| 99久久99久久久精品蜜桃| 免费观看的影片在线观看| 精品日产1卡2卡| 亚洲在线观看片| 国产熟女xx| 精品国产亚洲在线| 1000部很黄的大片| av中文乱码字幕在线| 欧美bdsm另类| 香蕉av资源在线| 免费高清视频大片| 精品久久久久久久末码| a级一级毛片免费在线观看| 国产精品影院久久| 亚洲国产精品成人综合色| 高清在线国产一区| 又粗又爽又猛毛片免费看| 制服丝袜大香蕉在线| 日韩欧美在线二视频| 国产精品亚洲一级av第二区| 99在线人妻在线中文字幕| 国产日本99.免费观看| 少妇熟女aⅴ在线视频| 2021天堂中文幕一二区在线观| 男女床上黄色一级片免费看| 久久久久久久精品吃奶| 国产麻豆成人av免费视频| 又黄又爽又刺激的免费视频.| 小说图片视频综合网站| 无遮挡黄片免费观看| 欧美精品国产亚洲| 91在线精品国自产拍蜜月| 一本一本综合久久| 久久伊人香网站| 麻豆久久精品国产亚洲av| 乱人视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 99视频精品全部免费 在线| 99在线人妻在线中文字幕| 国产亚洲欧美98| 欧美不卡视频在线免费观看| 久久精品综合一区二区三区| 男人狂女人下面高潮的视频| 九九在线视频观看精品| 舔av片在线| 日本三级黄在线观看| 我的老师免费观看完整版| 嫩草影院精品99| 99热这里只有精品一区| 成人国产一区最新在线观看| 中亚洲国语对白在线视频| 日韩大尺度精品在线看网址| 中文字幕人成人乱码亚洲影| 无人区码免费观看不卡| 精品不卡国产一区二区三区| 久久久久久大精品| 乱码一卡2卡4卡精品| 欧美+亚洲+日韩+国产| 天美传媒精品一区二区| 尤物成人国产欧美一区二区三区| 日本与韩国留学比较| 18美女黄网站色大片免费观看| 免费观看精品视频网站| 搡老岳熟女国产| 毛片女人毛片| 亚洲国产日韩欧美精品在线观看| 深夜精品福利| 亚洲avbb在线观看| 午夜福利在线在线| 久久性视频一级片| 成年版毛片免费区| 18禁在线播放成人免费| 丰满人妻熟妇乱又伦精品不卡| 看十八女毛片水多多多| 成人国产一区最新在线观看| 一级作爱视频免费观看| 亚洲欧美日韩高清在线视频| 在线观看66精品国产| 久久久久久久久大av| 亚洲aⅴ乱码一区二区在线播放| 能在线免费观看的黄片| 国产精品不卡视频一区二区 | 老鸭窝网址在线观看| 成人亚洲精品av一区二区| 久久精品国产亚洲av涩爱 | 男女做爰动态图高潮gif福利片| 日日干狠狠操夜夜爽| 国内久久婷婷六月综合欲色啪| 亚洲中文字幕日韩| 亚洲自偷自拍三级| 禁无遮挡网站| 观看免费一级毛片| 激情在线观看视频在线高清| 精品国产三级普通话版| 成人午夜高清在线视频| 亚洲不卡免费看| 午夜两性在线视频| www.999成人在线观看| 亚洲av免费高清在线观看| 精品久久久久久久久久免费视频| 国产精品99久久久久久久久| 国产三级在线视频| 欧美一区二区国产精品久久精品| 欧美激情在线99| 成年人黄色毛片网站| 日韩欧美精品v在线| 一个人免费在线观看的高清视频| 午夜激情福利司机影院| 亚洲欧美日韩高清在线视频| 亚洲精品一卡2卡三卡4卡5卡| www.色视频.com| 色精品久久人妻99蜜桃| 国产色爽女视频免费观看| 亚洲无线观看免费| 国产精品久久久久久久电影| 精品不卡国产一区二区三区| 人人妻人人澡欧美一区二区| x7x7x7水蜜桃| 精品久久国产蜜桃| 欧美黄色片欧美黄色片| 一进一出抽搐动态| 亚洲人成网站高清观看| 国内精品美女久久久久久| 久久久久久国产a免费观看| 欧美日韩福利视频一区二区| 搡老熟女国产l中国老女人| 日本 av在线| 1024手机看黄色片| 成人国产综合亚洲| 九九热线精品视视频播放| 精品午夜福利视频在线观看一区| .国产精品久久| 有码 亚洲区| 免费看日本二区| 久久精品91蜜桃| 国产精品精品国产色婷婷| 在线播放国产精品三级| 国产三级在线视频| 久久午夜亚洲精品久久| 欧美+日韩+精品| 久久久成人免费电影| 久久午夜亚洲精品久久| 波野结衣二区三区在线| 免费大片18禁| 99热这里只有是精品50| av欧美777| 亚洲成人久久爱视频| 国产白丝娇喘喷水9色精品| 国产精品,欧美在线| 国产视频内射| 悠悠久久av| 两个人视频免费观看高清| 桃色一区二区三区在线观看| 国产一区二区三区视频了| 免费大片18禁| 成人一区二区视频在线观看| 国产免费av片在线观看野外av| 亚洲熟妇中文字幕五十中出| 美女黄网站色视频| 亚洲成人精品中文字幕电影| 国产毛片a区久久久久| 床上黄色一级片| 日本 av在线| 一进一出好大好爽视频| 99热只有精品国产| 国产三级黄色录像| 日日干狠狠操夜夜爽| 亚洲五月婷婷丁香| 亚洲成av人片免费观看| x7x7x7水蜜桃| 人妻久久中文字幕网| 亚洲精品乱码久久久v下载方式| 欧洲精品卡2卡3卡4卡5卡区| 国产在视频线在精品| 久久国产精品人妻蜜桃| 亚洲av中文字字幕乱码综合| 日本熟妇午夜| 又紧又爽又黄一区二区| 草草在线视频免费看| 一本一本综合久久| av福利片在线观看| 搡老妇女老女人老熟妇| 99久久精品一区二区三区| 日本 av在线| 亚洲人成伊人成综合网2020| 国产日本99.免费观看| 美女高潮喷水抽搐中文字幕| 一本综合久久免费| 天堂av国产一区二区熟女人妻| 少妇高潮的动态图| 精品国产亚洲在线| 99视频精品全部免费 在线| 日本熟妇午夜| 欧美乱色亚洲激情| 人人妻人人看人人澡| 日日摸夜夜添夜夜添av毛片 | 老司机深夜福利视频在线观看| 欧美高清成人免费视频www| 日日摸夜夜添夜夜添小说| 欧美一级a爱片免费观看看| 成人一区二区视频在线观看| 中文字幕高清在线视频| 亚洲国产日韩欧美精品在线观看| 91九色精品人成在线观看| 国产午夜精品论理片| 成人国产综合亚洲| 成人特级黄色片久久久久久久| 亚洲第一区二区三区不卡| 悠悠久久av| av在线观看视频网站免费| eeuss影院久久| 亚洲激情在线av| 波多野结衣高清作品| 色综合站精品国产| 俄罗斯特黄特色一大片| 国产精品久久久久久人妻精品电影| 欧美成人一区二区免费高清观看| 俺也久久电影网| 国产美女午夜福利| 国产单亲对白刺激| 亚洲第一电影网av| 99热精品在线国产| 亚洲一区二区三区不卡视频| 人妻制服诱惑在线中文字幕| 亚洲经典国产精华液单 | 在线观看午夜福利视频| 久久6这里有精品| 最近最新免费中文字幕在线| 老鸭窝网址在线观看| 国产黄a三级三级三级人| 一级毛片久久久久久久久女| 国产成人影院久久av| 又爽又黄无遮挡网站| 啦啦啦韩国在线观看视频| 亚洲电影在线观看av| av天堂在线播放| www.999成人在线观看| 啪啪无遮挡十八禁网站| 91麻豆av在线| 国产伦一二天堂av在线观看| 久久久久久九九精品二区国产| 宅男免费午夜| av视频在线观看入口| 国产精华一区二区三区| 老熟妇乱子伦视频在线观看| 久久6这里有精品| 伊人久久精品亚洲午夜| 亚洲欧美日韩卡通动漫| 婷婷亚洲欧美| 色精品久久人妻99蜜桃| 免费av不卡在线播放| 久久6这里有精品| 757午夜福利合集在线观看| 亚洲av一区综合| 在线播放无遮挡| 国产亚洲欧美在线一区二区| a级一级毛片免费在线观看| 亚洲av成人av| 欧美潮喷喷水| 90打野战视频偷拍视频| 9191精品国产免费久久| 乱人视频在线观看| 国产成+人综合+亚洲专区| 我要看日韩黄色一级片| 日韩欧美 国产精品| 亚洲最大成人手机在线| 国产一区二区三区视频了| 成年人黄色毛片网站| 久久精品91蜜桃| 青草久久国产| 全区人妻精品视频| 欧美成狂野欧美在线观看| 成人特级黄色片久久久久久久| 国产成人影院久久av| 国产亚洲精品久久久久久毛片| 欧美xxxx黑人xx丫x性爽| av天堂中文字幕网| 桃红色精品国产亚洲av| 麻豆久久精品国产亚洲av| 免费大片18禁| 久久午夜亚洲精品久久| 国产在视频线在精品| 黄色配什么色好看| 99视频精品全部免费 在线| 国产精品综合久久久久久久免费| 无遮挡黄片免费观看| 97超视频在线观看视频| 国产精品精品国产色婷婷| 国产亚洲av嫩草精品影院| 亚洲精品456在线播放app | av在线老鸭窝| 熟女电影av网| 九九久久精品国产亚洲av麻豆| 噜噜噜噜噜久久久久久91| 三级毛片av免费| 韩国av一区二区三区四区| 国产三级在线视频| 久久天躁狠狠躁夜夜2o2o| 少妇高潮的动态图| 18禁裸乳无遮挡免费网站照片| 麻豆成人av在线观看| 国产精品久久视频播放| 亚洲精品在线美女| av女优亚洲男人天堂| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产色片| av中文乱码字幕在线| 国产成人欧美在线观看| 欧美bdsm另类| 高清在线国产一区| 亚洲,欧美精品.| 啦啦啦观看免费观看视频高清| 久久人人爽人人爽人人片va | 国产私拍福利视频在线观看| 怎么达到女性高潮| 亚洲av日韩精品久久久久久密| 中文字幕高清在线视频| 国产伦人伦偷精品视频| 少妇高潮的动态图| 亚洲五月婷婷丁香| 大型黄色视频在线免费观看| 日本黄色片子视频| 欧美在线黄色| 高清在线国产一区| 变态另类丝袜制服| 在线播放无遮挡| 亚洲av第一区精品v没综合| 欧美不卡视频在线免费观看| 中亚洲国语对白在线视频| 国产亚洲精品久久久久久毛片| 在现免费观看毛片| 亚洲av成人不卡在线观看播放网| 欧美丝袜亚洲另类 | 女生性感内裤真人,穿戴方法视频| 亚洲中文日韩欧美视频| 日本免费a在线| 在线十欧美十亚洲十日本专区| 国产精品久久久久久人妻精品电影| 日本成人三级电影网站| 麻豆国产av国片精品| 91九色精品人成在线观看| 搡老妇女老女人老熟妇| 欧美日韩综合久久久久久 | 亚洲一区高清亚洲精品| 亚洲av成人不卡在线观看播放网| 免费电影在线观看免费观看| 99久久无色码亚洲精品果冻| 欧美国产日韩亚洲一区| 中文亚洲av片在线观看爽| 欧美黑人欧美精品刺激| 真人一进一出gif抽搐免费| 国产伦精品一区二区三区四那| 人人妻,人人澡人人爽秒播| 免费观看人在逋| 亚洲欧美精品综合久久99| .国产精品久久| 国产精品一区二区三区四区久久| 国产精品美女特级片免费视频播放器| 亚洲中文日韩欧美视频| 午夜福利18| 成人无遮挡网站| 国产伦精品一区二区三区四那| 亚洲精品日韩av片在线观看| 欧美zozozo另类| 又黄又爽又免费观看的视频| 久久久久国内视频| 91久久精品电影网| 又爽又黄无遮挡网站| 九九久久精品国产亚洲av麻豆| 亚洲欧美日韩高清在线视频| 一个人免费在线观看的高清视频| 国内精品一区二区在线观看| 俺也久久电影网| 国产乱人视频| 亚洲欧美精品综合久久99| 91麻豆精品激情在线观看国产| 天堂网av新在线| 乱人视频在线观看| 亚洲专区国产一区二区| 国产真实乱freesex| 久久国产乱子伦精品免费另类| 国产精品永久免费网站| 久久精品影院6| 中文资源天堂在线| 夜夜爽天天搞| 老熟妇仑乱视频hdxx| 免费看日本二区| 色吧在线观看| 赤兔流量卡办理| 欧美激情在线99| 日本黄大片高清| 国内精品久久久久久久电影| 久久亚洲精品不卡| 欧美最黄视频在线播放免费| 日本免费a在线| 久久久精品大字幕| 亚洲精品一卡2卡三卡4卡5卡| 欧美日本视频| 男人舔奶头视频| 欧美性猛交╳xxx乱大交人| 精品一区二区三区视频在线| 成人亚洲精品av一区二区| 免费av不卡在线播放| 午夜福利免费观看在线| 精品午夜福利视频在线观看一区| 久久精品91蜜桃| 免费大片18禁| 亚洲第一电影网av| 亚洲成人精品中文字幕电影| 午夜福利在线观看吧| 一本精品99久久精品77| 一个人看视频在线观看www免费| 亚洲成av人片免费观看| 亚洲人成网站高清观看| 欧美日本亚洲视频在线播放| 国产精华一区二区三区| 国产精品亚洲一级av第二区| 女人十人毛片免费观看3o分钟| 国内精品一区二区在线观看| 狠狠狠狠99中文字幕| 亚洲国产精品合色在线| 色播亚洲综合网| 小说图片视频综合网站| a在线观看视频网站| 日本黄大片高清| 成人性生交大片免费视频hd| 久久精品人妻少妇| 最新在线观看一区二区三区| 国产三级在线视频| 91麻豆av在线| 久久这里只有精品中国| 国产久久久一区二区三区| 亚洲激情在线av| 欧美一级a爱片免费观看看| 俄罗斯特黄特色一大片| 免费看光身美女| 日本 欧美在线| 亚洲第一区二区三区不卡| 久久久国产成人免费| 男人舔奶头视频| 成人鲁丝片一二三区免费| 国产三级在线视频| 极品教师在线视频| 欧美午夜高清在线| 亚洲欧美日韩高清在线视频| 亚洲av电影在线进入| 欧美绝顶高潮抽搐喷水| 丝袜美腿在线中文| 身体一侧抽搐| 亚洲久久久久久中文字幕| 亚洲av日韩精品久久久久久密| 白带黄色成豆腐渣| 18禁黄网站禁片午夜丰满| 国产精品久久久久久精品电影| 中文字幕高清在线视频| 俺也久久电影网| 91av网一区二区| 一本精品99久久精品77| 变态另类丝袜制服| 亚洲av中文字字幕乱码综合| 日本黄色视频三级网站网址| 精品一区二区三区视频在线| 永久网站在线| av女优亚洲男人天堂| 亚洲真实伦在线观看| 成熟少妇高潮喷水视频| 国产精品,欧美在线| 中文在线观看免费www的网站| 男女床上黄色一级片免费看| 怎么达到女性高潮| 欧美精品国产亚洲| 少妇的逼水好多| 日本免费a在线| 国产毛片a区久久久久| 精品久久久久久久末码| 午夜福利在线观看免费完整高清在 | 麻豆久久精品国产亚洲av| 精品无人区乱码1区二区| 欧美在线黄色| 国产aⅴ精品一区二区三区波| 精品午夜福利在线看| 亚洲午夜理论影院| 丝袜美腿在线中文| 午夜福利在线观看免费完整高清在 | 俺也久久电影网| 国产高清视频在线播放一区| 久久久国产成人免费| 自拍偷自拍亚洲精品老妇| 色av中文字幕| 国产极品精品免费视频能看的| 哪里可以看免费的av片| 天堂网av新在线| 色吧在线观看| 国产色婷婷99| 长腿黑丝高跟|