• <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
    一種改進的基于邊緣加強超分辨率算法
    狂野欧美白嫩少妇大欣赏| 免费观看的影片在线观看| 国产黄a三级三级三级人| 美女高潮喷水抽搐中文字幕| 精品久久久久久久久亚洲 | 成人午夜高清在线视频| 成人二区视频| 亚洲国产日韩欧美精品在线观看| 麻豆久久精品国产亚洲av| 亚洲美女视频黄频| 啪啪无遮挡十八禁网站| 一本一本综合久久| 日本黄色片子视频| 亚洲国产日韩欧美精品在线观看| 久久欧美精品欧美久久欧美| 久久久久久大精品| 国产麻豆成人av免费视频| 观看免费一级毛片| 九九在线视频观看精品| 日日撸夜夜添| 日韩欧美免费精品| 成年免费大片在线观看| 中文资源天堂在线| 最好的美女福利视频网| 偷拍熟女少妇极品色| av女优亚洲男人天堂| 欧洲精品卡2卡3卡4卡5卡区| 成人欧美大片| av.在线天堂| 色在线成人网| 日本黄大片高清| 国产精品久久久久久亚洲av鲁大| 久久久久久久午夜电影| 热99re8久久精品国产| 成人美女网站在线观看视频| 永久网站在线| 亚洲熟妇熟女久久| 少妇熟女aⅴ在线视频| 22中文网久久字幕| 国产精品日韩av在线免费观看| 桃红色精品国产亚洲av| 男女那种视频在线观看| 午夜福利在线观看免费完整高清在 | 亚洲av日韩精品久久久久久密| 桃红色精品国产亚洲av| 淫妇啪啪啪对白视频| АⅤ资源中文在线天堂| 一进一出好大好爽视频| 九九在线视频观看精品| 国产女主播在线喷水免费视频网站 | 亚洲人与动物交配视频| 最近在线观看免费完整版| 日本三级黄在线观看| 国产探花极品一区二区| 午夜福利欧美成人| 夜夜爽天天搞| 亚洲第一电影网av| 中文亚洲av片在线观看爽| 精品国产三级普通话版| 国产一区二区亚洲精品在线观看| 亚洲四区av| 亚洲人成网站在线播| 级片在线观看| 性欧美人与动物交配| 在线免费十八禁| 伊人久久精品亚洲午夜| 尾随美女入室| 日韩,欧美,国产一区二区三区 | 欧美一区二区精品小视频在线| 亚洲 国产 在线| 无遮挡黄片免费观看| 91麻豆精品激情在线观看国产| 波多野结衣高清作品| 色综合站精品国产| 欧美日韩乱码在线| 欧美黑人巨大hd| 精品日产1卡2卡| 免费av观看视频| 日韩人妻高清精品专区| 极品教师在线免费播放| 校园春色视频在线观看| 亚洲中文字幕日韩| 男女边吃奶边做爰视频| 很黄的视频免费| 赤兔流量卡办理| 精品人妻偷拍中文字幕| 欧美性猛交╳xxx乱大交人| 精品无人区乱码1区二区| 啦啦啦观看免费观看视频高清| 国产一级毛片七仙女欲春2| 美女黄网站色视频| 男人的好看免费观看在线视频| 成人二区视频| 天堂av国产一区二区熟女人妻| 人妻丰满熟妇av一区二区三区| 欧美潮喷喷水| 亚洲精华国产精华液的使用体验 | 国产亚洲av嫩草精品影院| 久久久国产成人免费| 成年免费大片在线观看| 亚洲人成网站在线播| 免费在线观看日本一区| 国产午夜精品论理片| 中文亚洲av片在线观看爽| 午夜影院日韩av| 精品午夜福利在线看| 一级毛片久久久久久久久女| 午夜福利视频1000在线观看| 一级黄片播放器| 亚洲图色成人| 十八禁国产超污无遮挡网站| 午夜免费男女啪啪视频观看 | 久久久久久久久久成人| 成人国产综合亚洲| 欧美一区二区精品小视频在线| 一区福利在线观看| 亚洲国产欧美人成| 性插视频无遮挡在线免费观看| 简卡轻食公司| 日韩欧美在线二视频| 真人做人爱边吃奶动态| 国产乱人伦免费视频| 久久精品91蜜桃| 俄罗斯特黄特色一大片| 国语自产精品视频在线第100页| 黄色一级大片看看| 国产精品乱码一区二三区的特点| av在线亚洲专区| 日本黄色视频三级网站网址| 亚洲欧美日韩无卡精品| 国产av不卡久久| 69av精品久久久久久| 午夜激情欧美在线| 天堂影院成人在线观看| 久久久久久久久久成人| 一进一出抽搐gif免费好疼| 草草在线视频免费看| 久久精品综合一区二区三区| 最近最新免费中文字幕在线| 午夜激情福利司机影院| 日韩欧美国产在线观看| 亚洲国产精品合色在线| 日韩欧美国产一区二区入口| 色噜噜av男人的天堂激情| 亚洲内射少妇av| 美女xxoo啪啪120秒动态图| 一本久久中文字幕| 男人狂女人下面高潮的视频| 大又大粗又爽又黄少妇毛片口| 伊人久久精品亚洲午夜| 又粗又爽又猛毛片免费看| 精品人妻视频免费看| 国国产精品蜜臀av免费| 嫁个100分男人电影在线观看| 欧美成人免费av一区二区三区| 美女被艹到高潮喷水动态| 欧美区成人在线视频| 国产一区二区在线观看日韩| 亚洲性夜色夜夜综合| 精品一区二区三区视频在线| 国产亚洲欧美98| 国产精品电影一区二区三区| 国产免费男女视频| 91久久精品电影网| 国产精品久久久久久精品电影| 一个人免费在线观看电影| 在现免费观看毛片| 精品久久久久久久久久久久久| 国产精品99久久久久久久久| 久久精品国产亚洲av涩爱 | 久久久国产成人免费| 无人区码免费观看不卡| 日本一本二区三区精品| 国内精品一区二区在线观看| 大型黄色视频在线免费观看| 黄片wwwwww| 久久人人爽人人爽人人片va| 亚洲欧美清纯卡通| 国产精品久久电影中文字幕| 窝窝影院91人妻| 白带黄色成豆腐渣| 国产男靠女视频免费网站| 88av欧美| 性色avwww在线观看| 国产精品久久久久久久久免| 人人妻人人看人人澡| 国产亚洲精品久久久久久毛片| 精品久久久久久久久亚洲 | 色综合婷婷激情| 国产女主播在线喷水免费视频网站 | 别揉我奶头~嗯~啊~动态视频| 亚洲av电影不卡..在线观看| 国产精品人妻久久久久久| 精品福利观看| 日本黄色视频三级网站网址| 无遮挡黄片免费观看| 九九热线精品视视频播放| 午夜精品在线福利| 日本一本二区三区精品| 亚洲电影在线观看av| 久久久久久伊人网av| 国产伦人伦偷精品视频| av女优亚洲男人天堂| 观看免费一级毛片| 国产白丝娇喘喷水9色精品| 成人三级黄色视频| 少妇裸体淫交视频免费看高清| 一级黄片播放器| netflix在线观看网站| 91午夜精品亚洲一区二区三区 | 精品99又大又爽又粗少妇毛片 | 亚洲性夜色夜夜综合| 少妇人妻精品综合一区二区 | 久久精品夜夜夜夜夜久久蜜豆| 国产高清有码在线观看视频| 嫁个100分男人电影在线观看| 美女高潮的动态| 精品国内亚洲2022精品成人| 女同久久另类99精品国产91| 一区福利在线观看| 观看美女的网站| 国产亚洲欧美98| 看黄色毛片网站| 日韩一本色道免费dvd| 久久精品影院6| 国内精品宾馆在线| 亚洲av成人精品一区久久| av天堂中文字幕网| eeuss影院久久| 成人午夜高清在线视频| 亚洲av电影不卡..在线观看| 成人高潮视频无遮挡免费网站| 亚洲乱码一区二区免费版| 国产精品精品国产色婷婷| 美女高潮的动态| 欧美+日韩+精品| 国产 一区精品| 日本a在线网址| 内地一区二区视频在线| 日日摸夜夜添夜夜添av毛片 | 性插视频无遮挡在线免费观看| 色尼玛亚洲综合影院| 人妻制服诱惑在线中文字幕| 午夜福利18| 免费看美女性在线毛片视频| 女人被狂操c到高潮| 两性午夜刺激爽爽歪歪视频在线观看| 国产午夜精品久久久久久一区二区三区 | 中文字幕人妻熟人妻熟丝袜美| 国产老妇女一区| 久久中文看片网| 性插视频无遮挡在线免费观看| 日韩欧美在线二视频| 国产精品98久久久久久宅男小说| 99热网站在线观看| 给我免费播放毛片高清在线观看| 亚洲av成人精品一区久久| 小说图片视频综合网站| 国产蜜桃级精品一区二区三区| 精品久久久久久久末码| 中文字幕人妻熟人妻熟丝袜美| 午夜福利在线观看吧| 国产日本99.免费观看| 婷婷精品国产亚洲av| a级毛片a级免费在线| 久久久久精品国产欧美久久久| av天堂在线播放| 国产精品免费一区二区三区在线| 亚洲欧美激情综合另类| 麻豆一二三区av精品| 久久精品国产亚洲av香蕉五月| 无人区码免费观看不卡| 亚洲成人久久性| 啪啪无遮挡十八禁网站| 日韩国内少妇激情av| 国产av不卡久久| 日本三级黄在线观看| 变态另类成人亚洲欧美熟女| 超碰av人人做人人爽久久| 午夜福利18| 国产高清有码在线观看视频| 如何舔出高潮| 在线观看66精品国产| 色吧在线观看| а√天堂www在线а√下载| 国产蜜桃级精品一区二区三区| 啪啪无遮挡十八禁网站| 真实男女啪啪啪动态图| 国产极品精品免费视频能看的| 91久久精品国产一区二区三区| 久久久国产成人精品二区| 亚洲色图av天堂| 亚洲av美国av| 麻豆精品久久久久久蜜桃| 精品福利观看| 国产精品国产三级国产av玫瑰| 日韩欧美精品免费久久| 久久久久久久亚洲中文字幕| 免费在线观看影片大全网站| 亚洲精品粉嫩美女一区| 在线免费十八禁| 日韩在线高清观看一区二区三区 | 在线免费观看的www视频| 别揉我奶头~嗯~啊~动态视频| 亚洲午夜理论影院| 久久国产乱子免费精品| 国产精品久久久久久av不卡| av在线老鸭窝| 国产三级中文精品| 日本a在线网址| 级片在线观看| 人人妻人人澡欧美一区二区| 一进一出抽搐动态| 国产女主播在线喷水免费视频网站 | 国产精品电影一区二区三区| 91av网一区二区| 俺也久久电影网| 欧美日韩国产亚洲二区| 最近视频中文字幕2019在线8| 淫妇啪啪啪对白视频| 久久99热这里只有精品18| 亚洲一区高清亚洲精品| 亚洲美女黄片视频| 成年免费大片在线观看| 亚洲精品在线观看二区| 99热6这里只有精品| 亚洲专区国产一区二区| 国产精品日韩av在线免费观看| 男人狂女人下面高潮的视频| 欧美日韩综合久久久久久 | 成人av在线播放网站| 成人高潮视频无遮挡免费网站| 国产私拍福利视频在线观看| 国产欧美日韩一区二区精品| 精品久久久久久久人妻蜜臀av| 国产主播在线观看一区二区| 精品午夜福利视频在线观看一区| 亚洲自拍偷在线| 精品无人区乱码1区二区| 悠悠久久av| 国产亚洲精品久久久com| 变态另类成人亚洲欧美熟女| 欧美xxxx黑人xx丫x性爽| 天天一区二区日本电影三级| 91在线精品国自产拍蜜月| 亚洲午夜理论影院| 黄色视频,在线免费观看| 真人做人爱边吃奶动态| 亚洲av日韩精品久久久久久密| 亚洲三级黄色毛片| 亚洲熟妇熟女久久| 日韩一区二区视频免费看| 少妇丰满av| 精品免费久久久久久久清纯| 日韩 亚洲 欧美在线| 午夜福利在线观看免费完整高清在 | 日本 欧美在线| 欧美激情久久久久久爽电影| 亚洲av一区综合| 欧美高清成人免费视频www| 国产精品日韩av在线免费观看| 老女人水多毛片| 色在线成人网| xxxwww97欧美| 亚洲自拍偷在线| 亚洲成a人片在线一区二区| 中文字幕av成人在线电影| 欧美色视频一区免费| 亚洲中文字幕日韩| 美女cb高潮喷水在线观看| 女同久久另类99精品国产91| 99riav亚洲国产免费| 小说图片视频综合网站| 日本一二三区视频观看| 亚洲一区二区三区色噜噜| 禁无遮挡网站| 人妻夜夜爽99麻豆av| 亚洲最大成人手机在线| 国产女主播在线喷水免费视频网站 | 亚洲国产欧洲综合997久久,| 九色国产91popny在线| 国产精品久久久久久精品电影| 直男gayav资源| 久久婷婷人人爽人人干人人爱| 欧美高清性xxxxhd video| 五月伊人婷婷丁香| 亚洲国产精品久久男人天堂| 久久精品国产清高在天天线| 国产aⅴ精品一区二区三区波| 精品一区二区免费观看| 国产精品永久免费网站| 日韩中字成人| 久久精品国产亚洲网站| 人人妻,人人澡人人爽秒播| av女优亚洲男人天堂| 波野结衣二区三区在线| 搡老妇女老女人老熟妇| 不卡视频在线观看欧美| 久久亚洲真实| 国产老妇女一区| 桃红色精品国产亚洲av| 1000部很黄的大片| 国产美女午夜福利| av中文乱码字幕在线| 色视频www国产| 深夜a级毛片| 99久国产av精品| 国产不卡一卡二| 亚洲av电影不卡..在线观看| 久久久成人免费电影| 欧美国产日韩亚洲一区| 美女高潮喷水抽搐中文字幕| 成人二区视频| 午夜福利在线在线| 国产亚洲av嫩草精品影院| 国产乱人伦免费视频| 欧美黑人巨大hd| 一级黄色大片毛片| 久久午夜亚洲精品久久| 国产午夜福利久久久久久| 亚洲国产精品成人综合色| 两性午夜刺激爽爽歪歪视频在线观看| 啪啪无遮挡十八禁网站| 夜夜看夜夜爽夜夜摸| av在线亚洲专区| 免费观看人在逋| 黄色配什么色好看| 中亚洲国语对白在线视频| 看免费成人av毛片| 一边摸一边抽搐一进一小说| 不卡一级毛片| 一进一出抽搐动态| 午夜福利18| 女的被弄到高潮叫床怎么办 | 国产在视频线在精品| 国产精品,欧美在线| 日本在线视频免费播放| 五月玫瑰六月丁香| 中文字幕av在线有码专区| 中文字幕熟女人妻在线| 国产精品乱码一区二三区的特点| 色精品久久人妻99蜜桃| 动漫黄色视频在线观看| 欧美色视频一区免费| 欧美最新免费一区二区三区| 亚洲狠狠婷婷综合久久图片| 精品久久久久久久久久免费视频| 99在线人妻在线中文字幕| 中文字幕高清在线视频| 日本爱情动作片www.在线观看 | 成人高潮视频无遮挡免费网站| 少妇猛男粗大的猛烈进出视频 | 欧美xxxx性猛交bbbb| 神马国产精品三级电影在线观看| 如何舔出高潮| 色吧在线观看| 精品久久久久久久人妻蜜臀av| 欧美性感艳星| 亚洲av二区三区四区| 禁无遮挡网站| 神马国产精品三级电影在线观看| 婷婷丁香在线五月| 免费电影在线观看免费观看| 国产精品野战在线观看| 我要搜黄色片| 亚洲av二区三区四区| 国产亚洲精品久久久com| 精品人妻偷拍中文字幕| 国产精品伦人一区二区| 高清在线国产一区| 不卡一级毛片| 嫁个100分男人电影在线观看| 黄色日韩在线| 国产伦在线观看视频一区| 久久久久久久久大av| 一级av片app| 久久欧美精品欧美久久欧美| 国产欧美日韩一区二区精品| 成人高潮视频无遮挡免费网站| 国产精品自产拍在线观看55亚洲| 美女被艹到高潮喷水动态| 99热只有精品国产| 国产久久久一区二区三区| 婷婷色综合大香蕉| 免费看av在线观看网站| 亚洲精品一区av在线观看| 国产日本99.免费观看| 神马国产精品三级电影在线观看| 国产中年淑女户外野战色| 日韩欧美国产一区二区入口| 春色校园在线视频观看| 亚洲经典国产精华液单| 一个人观看的视频www高清免费观看| 欧美日韩国产亚洲二区| 国产精品自产拍在线观看55亚洲| netflix在线观看网站| 亚洲狠狠婷婷综合久久图片| 日本一本二区三区精品| 91av网一区二区| 韩国av在线不卡| 亚洲成人久久性| 男插女下体视频免费在线播放| 99久久九九国产精品国产免费| 深爱激情五月婷婷| 窝窝影院91人妻| 88av欧美| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产精品sss在线观看| 亚洲欧美激情综合另类| 亚洲第一区二区三区不卡| 亚洲五月天丁香| 欧美日韩瑟瑟在线播放| 我要搜黄色片| 99热精品在线国产| 88av欧美| 国产高清三级在线| 亚洲综合色惰| 亚洲图色成人| 久久精品国产自在天天线| 三级毛片av免费| 我的老师免费观看完整版| 久久国内精品自在自线图片| 亚洲专区国产一区二区| 99久久无色码亚洲精品果冻| 国产白丝娇喘喷水9色精品| 夜夜爽天天搞| 看十八女毛片水多多多| 久久6这里有精品| 国产免费男女视频| 国产精品女同一区二区软件 | 深夜a级毛片| 亚洲天堂国产精品一区在线| 日韩强制内射视频| 免费av不卡在线播放| 熟女人妻精品中文字幕| 又粗又爽又猛毛片免费看| xxxwww97欧美| 一a级毛片在线观看| 亚洲在线自拍视频| 午夜亚洲福利在线播放| 亚洲电影在线观看av| 日本黄色片子视频| 精品乱码久久久久久99久播| 亚洲性久久影院| 桃红色精品国产亚洲av| 午夜影院日韩av| 老女人水多毛片| 俄罗斯特黄特色一大片| 熟女电影av网| 国产69精品久久久久777片| 高清日韩中文字幕在线| 91狼人影院| 俺也久久电影网| 亚洲天堂国产精品一区在线| 我要看日韩黄色一级片| 亚洲人成网站在线播放欧美日韩| 亚洲欧美激情综合另类| 我的女老师完整版在线观看| 亚洲午夜理论影院| 高清毛片免费观看视频网站| 久久人人精品亚洲av| 天堂影院成人在线观看| 天堂动漫精品| 天堂av国产一区二区熟女人妻| 嫩草影院精品99| 日本熟妇午夜| 嫩草影院入口| 在线天堂最新版资源| 精品欧美国产一区二区三| 97碰自拍视频| 久久精品国产99精品国产亚洲性色| 色噜噜av男人的天堂激情| 一本精品99久久精品77| 美女高潮喷水抽搐中文字幕| 国产伦精品一区二区三区四那| 91久久精品国产一区二区三区| 日韩在线高清观看一区二区三区 | 无人区码免费观看不卡| 亚洲av中文av极速乱 | 亚洲 国产 在线| 夜夜夜夜夜久久久久| 五月伊人婷婷丁香| 蜜桃久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 99热精品在线国产| 我要看日韩黄色一级片| 日韩国内少妇激情av| 亚洲av熟女| 一级黄色大片毛片| 一本一本综合久久| 亚洲 国产 在线| 欧美3d第一页| 成人国产麻豆网| 十八禁国产超污无遮挡网站| 亚洲精品粉嫩美女一区| 在线播放无遮挡| 亚洲欧美日韩卡通动漫| 欧美日韩综合久久久久久 | 可以在线观看毛片的网站| 亚洲经典国产精华液单| av女优亚洲男人天堂| 热99在线观看视频| 91在线观看av| 一本一本综合久久| 午夜a级毛片| netflix在线观看网站| 亚洲最大成人手机在线| 亚洲美女黄片视频| 国内久久婷婷六月综合欲色啪| 哪里可以看免费的av片| 亚洲最大成人中文| 国产成人av教育| 免费不卡的大黄色大毛片视频在线观看 | 中文字幕av在线有码专区| 一本久久中文字幕|