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

    城鎮(zhèn)陰影區(qū)植物分布信息提取與亮度修復(fù)的方法

    2011-12-20 00:56:40周堅(jiān)華周一凡郭曉華
    關(guān)鍵詞:植被指數(shù)樹冠陰影

    周堅(jiān)華, 周一凡, 郭曉華, 任 珍

    (1.華東師范大學(xué) 地理信息科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,上海 200062;2.斯坦福大學(xué) 地球科學(xué)學(xué)院能源工程系,CA94305-2220,USA)

    城鎮(zhèn)陰影區(qū)植物分布信息提取與亮度修復(fù)的方法

    周堅(jiān)華1, 周一凡2, 郭曉華1, 任 珍1

    (1.華東師范大學(xué) 地理信息科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,上海 200062;2.斯坦福大學(xué) 地球科學(xué)學(xué)院能源工程系,CA94305-2220,USA)

    在建筑林立的城鎮(zhèn)建成區(qū),其綠化地通常有約20%的植物處于陰影區(qū),為了實(shí)現(xiàn)綠化植物群碳捕獲模型的全遙感驅(qū)動,必須解決的基礎(chǔ)性難題之一是分離與提取陰影區(qū)植物.為此提出了歸一化陰影指數(shù)(Normalized Difference Umbra Index,NDUI)的概念.當(dāng)采用歸一化植被指數(shù)(Normalized Difference Vegetation Index,NDVI)確定植被范圍后,可以通過NDUI進(jìn)一步分離和提取陰影區(qū)植物.為了改善識別效果,還可以在此基礎(chǔ)上對陰影區(qū)植被進(jìn)行亮度修復(fù).經(jīng)過仿真測試,證明該算法具有很好的穩(wěn)健性.

    城鎮(zhèn)綠地; 碳捕獲; 陰影指數(shù); 亮度修復(fù)

    0 引 言

    提取遙感圖像上的陰影信息和進(jìn)行亮度補(bǔ)償是遙感圖像處理和信息提取研究領(lǐng)域的一個(gè)較早引起關(guān)注的困難課題.近年來,國內(nèi)外學(xué)者就此提出過一些算法[1-5].然而,至今尚沒有一種通用有效的方法.這主要是因?yàn)椴煌膱鼍昂腿蝿?wù)對陰影提取和補(bǔ)償?shù)囊笥休^大區(qū)別.比如一些研究將陰影考慮為低頻信號,如較大山體背陰坡形成的陰影等,此時(shí)可以通過低通濾波的方法提取陰影[6];而對于那些周期性重復(fù)、邊緣灰度梯度明顯的陰影,比如城鎮(zhèn)居住小區(qū)建筑物形成的陰影等,則適合采用一些對影像紋理敏感的描述符表征它們[7-9].另有一些研究根據(jù)陰影區(qū)亮度在各個(gè)波段均體現(xiàn)為低值的特性,來區(qū)分陰影與暗色調(diào)地物[10].由于遙感影像上每個(gè)像素的亮度都是輻射源輻射通量作用于地物反射(包括部分吸收和散射)特性的復(fù)雜函數(shù),在光照條件和地物的某些信息(比如太陽高度角、地形坡度、坡向、建筑物高度等)已知的情況下,也有一些研究考慮逐像素直接計(jì)算陰影改正值[11-13].

    我們的研究小組一直致力于城鎮(zhèn)綠化植物群碳捕獲的遙感驅(qū)動模型研究.模型的主要驅(qū)動參數(shù)可以從遙感圖像上自動識別和模擬計(jì)算得到,比如樹種、冠徑、株數(shù)、郁閉度、陰影遮蔽率等.在這些參數(shù)的提取中,一個(gè)必須解決的基礎(chǔ)性難題是陰影區(qū)植物信息的分離與提取.

    在景觀空間異質(zhì)度極高的城鎮(zhèn)建成區(qū),通常有約20%的植物處于陰影區(qū),即那些對植物形成光線遮擋的建筑物陰影區(qū)、高大植物對低矮植物形成的光線遮擋區(qū),以及樹冠的本影區(qū);與此相對的是對植被基本沒有遮擋的亮區(qū).無論采用何種遙感圖像數(shù)據(jù),陰影區(qū)和亮區(qū)植被的反射波譜差異都很大;不僅如此,在陰影區(qū),植物的樹種、植株大小和分布等信息均大幅減弱.因此,解決陰影區(qū)植物信息提取問題,不僅要求提取植被覆蓋范圍,而且要求區(qū)分植株,即提取喬、灌木的輪廓.本文的研究和實(shí)驗(yàn)正是圍繞這一要求展開的.

    1 方法和數(shù)據(jù)

    這部分將從分析城鎮(zhèn)綠化植物群的波譜特征入手,選擇適當(dāng)?shù)闹脖恢笖?shù),盡可能完整地提取植被覆蓋范圍.為了區(qū)分亮區(qū)和陰影區(qū)植被,這里提出了歸一化陰影指數(shù)NDUI的概念,并通過NDUI和歸一化植被指數(shù)NDVI聯(lián)合提取陰影區(qū)植被,以及對陰影區(qū)植被做亮度修復(fù).下面展開介紹它們的理論和方法.所有測試樣本取自上海市1988~2003年1∶5 000~1∶10 000的彩紅外航片.

    1.1 歸一化植被指數(shù)和提取植被總覆蓋

    為了分析城市建成區(qū)綠化植物群與背景下墊面的反射波譜特征差異,這里采集植被和其他主要城鎮(zhèn)下墊面類型樣本的RGB值(分別大致對應(yīng)景物波譜的IR,R,G值).眾所周知,含葉綠素的健康植物在IR波段具有高反射特性,而在其他兩個(gè)波段具有強(qiáng)吸收的特點(diǎn).以下分別用R,G,B,和IR來表示相應(yīng)的RGB和IR值.樹冠無論在亮區(qū)和陰影區(qū)都有IR?R和IR?G;而其他下墊面類型的波譜沒有這個(gè)特征,通常表現(xiàn)為IR、R、G比較接近.

    因此,對于彩紅外圖像IM(IR,R,G),植被像素集V可以由(1)式獲得.

    其中,V是圖像IM 滿足條件IM(IR-G)>T的子集,該條件也可以寫作IM(IR-R)>T.T為閾值,根據(jù)對具體圖像的波譜分析給出.

    各種植被指數(shù)一般都依據(jù)(1)式的原理給出.但為了使用方便,通常要求植被指數(shù)存在一個(gè)比較確定的分割閾值,該閾值與圖像類型和獲取時(shí)的輻射照度水平等環(huán)境條件的差異無關(guān).一些歸一化植被指數(shù)符合這一要求.本文比較了以下三種植被指數(shù)提取城鎮(zhèn)綠化植物覆蓋范圍的效果(見圖1).

    圖1 以代表性植被指數(shù)提取城鎮(zhèn)綠化植物群覆蓋范圍的效果比較Fig.1 Effective comparison of extracting urban landscape vegetation by three typical vegetation indexes

    由圖1可見,使用NDVI提取的植被覆蓋范圍(包括陰影區(qū)植被)比較準(zhǔn)確完整.故本文選用NDVI提取城鎮(zhèn)綠化植被總覆蓋.

    對于區(qū)分植被覆蓋與背景區(qū)域,NDVI具有比較確定的分割閾值;此外,它對圖像輻照度水平差異不太敏感.圖2顯示了亮區(qū)、陰影區(qū)植被和城鎮(zhèn)其它主要下墊面類型的NDVI(樣本沿橫軸排列,樣本序號省略未注出,圖4同).由圖可見,無論在亮區(qū)還是陰影區(qū),植被的NDVI顯著大于其它下墊面類型;根據(jù)試驗(yàn),以NDVI>0.18為條件,可以準(zhǔn)確提取植被范圍.因?yàn)閷τ谄渌聣|面類型的樣本,在亮區(qū)和陰影區(qū)都有-0.36<NDVI<1.8.

    1.2 歸一化陰影指數(shù)和提取陰影區(qū)植被

    據(jù)前述,NDVI用于提取植被總覆蓋非常合理和有效.由于歸一化,NDVI對亮區(qū)和陰影區(qū)的照度水平差不敏感,這有利于從多種下墊面類型中提取植被,卻不利于在植被覆蓋區(qū)域中區(qū)分亮區(qū)和陰影區(qū).此時(shí)需要一個(gè)新指數(shù)用于陰影區(qū)植被分離和提取.

    在陰影區(qū)植被提取試驗(yàn)中,飽和度S的特殊表現(xiàn)引人注意.S是HSI色彩系統(tǒng)中的分量.該系統(tǒng)以色調(diào)H、飽和度S和亮度I表達(dá)色彩空間.當(dāng)以HSI系統(tǒng)分析城鎮(zhèn)下墊面的色彩特征時(shí),不難發(fā)現(xiàn)在城鎮(zhèn)地區(qū),植被覆蓋區(qū)具有較高的飽和度,而大多數(shù)其它下墊面類型,如人工建材表面,表現(xiàn)為消色體,飽和度很低.遙感圖像波譜特征分析也能給出類似結(jié)論.當(dāng)將彩紅外RGB圖像轉(zhuǎn)換為HSI圖像,無論在亮區(qū)或陰影區(qū),樹冠影像具有較高的飽和度;而水泥、瀝青等人工建材表面多為消色體,無論在亮區(qū)或陰影區(qū),都具有低飽和度特征.與明度I不同,飽和度S在陰影區(qū)并不顯著下降,有些甚至高于亮區(qū).所以,可以用(S-I)>T1(T1是以S-I提取陰影區(qū)植被時(shí)的下限)來判斷是否陰影區(qū)植被(見圖3).

    圖2 亮區(qū)、陰影區(qū)植被以及典型城鎮(zhèn)下墊面類型的NDVI比較Fig.2 Comparing NDVI of landscape vegetation and other surface species sampling at darken and bright areas

    圖3 以S-I提取陰影區(qū)植被Fig.3 Extracting plants at darken area of urban landscape vegetation referring the difference between saturation and brightness

    與植被指數(shù)歸一化的原理相似,為了增大亮區(qū)、陰影區(qū)S-I的差異,以減少對陰影區(qū)細(xì)小亮度差異的敏感性,本文提出了歸一化陰影指數(shù)(NDUI)的概念,即

    歸一化的實(shí)質(zhì)是用可變系數(shù)1/(S+I(xiàn))對S-I值做非線性拉伸.這里通過亮區(qū)和陰影區(qū)共140個(gè)植被樣本,分析了這種拉伸對陰影區(qū)樹冠提取效果的影響(見圖4).由圖可見,亮區(qū)和陰影區(qū)植被的S-I值差異不顯著,陰影區(qū)值僅略高于亮區(qū)值,一些亮區(qū)、陰影區(qū)樣本甚至出現(xiàn)粘結(jié)和相互穿插的情況.所以,直接以S-I來分離陰影區(qū)植被效果不佳.但當(dāng)使用NDUI時(shí),情況有了很大改善.由圖4可見,植被的NDUI與S-I在亮區(qū)相差不大,但在陰影區(qū),NDUI與S-I的差異顯著增大,即亮區(qū)與陰影區(qū)植被NDUI的差異顯著增大.通過分析具體樣本,NDUI跳高較多的是陰影區(qū)內(nèi)具有較高飽和度的樹冠.這就是使用NDUI可以較好提取陰影區(qū)樹冠覆蓋范圍的原因.不僅如此,就像經(jīng)過歸一化的NDVI對亮區(qū)和陰影區(qū)的照度水平差不敏感一樣,經(jīng)過歸一化的NDUI對陽高角變化等全局性的照度變化也不太敏感,這也是使用NDUI的優(yōu)點(diǎn)之一.

    根據(jù)試驗(yàn),可以用NDUI>0.4來提取陰影區(qū)中的樹冠.對于圖像IM,陰影區(qū)樹冠范圍DV可寫作:

    1.3 陰影區(qū)植被亮度修復(fù)

    在本文的試驗(yàn)中,陰影區(qū)植被亮度修復(fù)在前述陰影區(qū)植被提取的基礎(chǔ)上進(jìn)行,亮度修復(fù)值的計(jì)算采用直方圖均衡化(Histogram Equalization,HE)方法.

    HE是圖像增強(qiáng)中普遍使用的算法.它屬于基于累計(jì)分布函數(shù)(Cumulative Distribution Function,CDF)的直方圖修正法,其結(jié)果是使圖像灰度級分布具有均勻概率密度的特點(diǎn).這里采用HE方法,是為了將陰影區(qū)植被的亮度調(diào)整到與亮區(qū)植被盡可能接近.

    本文使用一種稱為CLAHE(Contrast-Limited Adaptive Histogram Equalization)的HE算法[17].CLAHE屬于自適應(yīng)HE,它采用一定尺寸的滑動窗口對每個(gè)像素進(jìn)行局部HE,使每個(gè)窗口的直方圖接近平坦的目標(biāo)直方圖.這里采用CLAHE,是為了避免由于陰影區(qū)植被像素在全圖像總像素中的占比較少,而受到壓抑.窗口尺寸由陰影區(qū)圖斑的平均尺寸確定,在本文的實(shí)驗(yàn)中,窗口尺寸設(shè)為12×12(像素).

    在對整個(gè)植被覆蓋區(qū)域做CLAHE后,以其結(jié)果替換陰影區(qū)內(nèi)的相應(yīng)像素.圖5是一些陰影區(qū)植被亮度修復(fù)的例子,由圖可見,修復(fù)運(yùn)算使陰影區(qū)植被獲得了比較自然和接近亮區(qū)植被的色彩.

    圖4 亮區(qū)、陰影區(qū)植被S-I與NDUI的比較Fig.4 Taking a comparison between S-Iand NDUI of landscape vegetation sampling at both darken and bright areas

    2 分析和討論

    這里著重討論本文算法的理論依據(jù)、使用效果和一些可能引起爭議的問題.

    2.1 為何NDVI對亮區(qū)、陰影區(qū)照度水平的差異不敏感

    NDVI對植被處于亮區(qū)、陰影區(qū)不敏感的主要原因,是NDVI對IR-G采用了非線性拉伸.對于陰影區(qū),分子IR-G相對于亮區(qū)減小,但同時(shí)分母IR+G也相應(yīng)減??;而對于亮區(qū),分子相對于陰影區(qū)增大,但同時(shí)分母也相應(yīng)增大.所以亮區(qū)NDVI并未明顯高于陰影區(qū).

    區(qū)別植被與背景的主要標(biāo)志是植被的IR-G大大高于其它下墊面類型.但直接使用IRG提取植被效果并不理想.這是因?yàn)镮R-G對亮區(qū)、陰影區(qū)的照度水平變化過于敏感.根據(jù)對圖2中樣本的統(tǒng)計(jì),亮區(qū)植被的IR-G均值為75.88,相比陰影區(qū)的這一值僅為33.16,二者相差56.27%,如果直接使用(IR-G)>T(T為以IR-G提取植被時(shí)的下限)提取植被,可能造成對亮區(qū)植被的過提取(植被中含有其它下墊面信息)和陰影區(qū)植被的欠提取.

    而使用NDVI提取植被的情況很不相同.由于NDVI使用可變系數(shù)1/(IR+G)非線性拉伸IR-G值,對于圖2所示的樣本,在陰影區(qū),植被的IR-G值比亮區(qū)平均降低56.27%的同時(shí),IR+G值平均比亮區(qū)也降低了45.83%.它意味著,二者的比值NDVI僅降低了21.94%,即這種拉伸降低了NDVI對亮區(qū)和陰影區(qū)照度水平差的敏感性.所以使用NDVI比直接使用IR-G提取的植被覆蓋范圍更準(zhǔn)確和完整,為進(jìn)一步分離和提取陰影區(qū)植被奠定了基礎(chǔ).5

    圖 陰影區(qū)植被提取和亮度修復(fù)實(shí)例Fig.5 Samples of extracting plants at darken area and repairing their brightness

    2.2 分離陰影區(qū)植被中使用NDUI和明度I的差別

    為了盡可能完整地提取陰影區(qū)植被,在使用NDVI提取植被覆蓋范圍時(shí),可以適當(dāng)放寬下限,比如放寬到NDVI>0.16.此時(shí)植被覆蓋范圍中,可能包含了部分陰影區(qū)背景,如裸土地、空禿草地等.如果使用明度I<T2(T2為陰影區(qū)明度上限)來提取陰影區(qū)植被,提取的植被里會包含這些背景信息;而使用NDUI則不會出現(xiàn)這個(gè)問題,因?yàn)樵诓始t外圖像上具有足夠高飽和度的只能是樹冠.圖6顯示了當(dāng)取NDVI>0.16時(shí),以這兩種方法提取和修復(fù)陰影區(qū)植物亮度的差別.

    因?yàn)镮直接反映輻射照度水平,當(dāng)遙感圖像的類型和獲取的環(huán)境條件不同時(shí),T2值都需要做出相應(yīng)調(diào)整.其次,以一個(gè)明度閾值來劃分陰影區(qū)常常不夠準(zhǔn)確.試驗(yàn)證明,以I方法檢測的邊界做陰影區(qū)植被亮度修復(fù),修復(fù)范圍會包括樹冠在草地和裸土地上的投影,甚至背光的混凝土墻面和建筑物在混凝土地面的投影等(如圖6(c)中箭頭標(biāo)示之處).而以NDUI>0.4作為陰影區(qū)植被提取條件時(shí)(簡稱“NDUI方法”),只提取了陰影區(qū)樹冠,并對其進(jìn)行了亮度修復(fù).在圖6(d)中,以紅、綠色細(xì)線分別標(biāo)示以I方法和NDUI方法提取的植被陰影區(qū).由圖可見,后者并非前者的子集,以I方法提取的陰影區(qū)邊界,有非植被覆蓋區(qū)被提取,也有陰影區(qū)植被被遺漏;而NDUI方法提取的邊界則比較合理,它比較完整地提取了陰影區(qū)樹冠.

    圖6 以I方法和NDUI方法進(jìn)行陰影區(qū)植被提取和亮度修復(fù)效果的比較Fig.6 A case of extracting plants at darken area referring both Iand NDUI and repairing brightness of these plants

    2.3 算法的穩(wěn)健性和實(shí)用性分析

    在陰影區(qū)植被提取算法中使用了兩個(gè)閾值:0.18(以NDVI>0.18提取植被總覆蓋)和0.4(以NDUI>0.4提取陰影區(qū)植被樹冠,該閾值根據(jù)實(shí)驗(yàn)數(shù)據(jù)歸納得到).下面討論這兩個(gè)閾值的穩(wěn)定性.

    NDVI對輻射照度水平的差異不太敏感,所以提取植被覆蓋范圍的閾值具有較好的穩(wěn)定性.比如圖2的樣本取自兩個(gè)年份不同季節(jié)的4張彩紅外航片,證明該閾值對時(shí)間和季節(jié)變化具有較好的適應(yīng)性.從相關(guān)報(bào)道來看,該閾值也與遙感圖像的類型無關(guān),如張友水使用IKONOS多光譜圖像,以NDVI>0.18提取城鎮(zhèn)綠地植被覆蓋效果很好[18],說明該閾值對不同類型的遙感圖像也具有穩(wěn)定性.

    在HSI彩色系統(tǒng)中,飽和度S是對輻射照度水平最不敏感的分量.而明度I卻正相反;所以用差值S-I來提取陰影區(qū)植被非常合理;但S-I對陰影區(qū)亮度水平的細(xì)小差異依然敏感.而NDUI通過歸一化進(jìn)一步加大了亮區(qū)和陰影區(qū)植被的差異,而使NDUI對陰影區(qū)亮度水平的細(xì)小差異不敏感,同時(shí)也使閾值0.4對陰影區(qū)亮度水平的細(xì)小差異不敏感.圖4中的140個(gè)樣本取自兩個(gè)年份不同季節(jié)的6張彩紅外航片,證明該閾值對時(shí)間和季節(jié)變化,以及陰影區(qū)亮度水平的細(xì)小差異具有較好的穩(wěn)定性.

    本文算法的穩(wěn)健性也從陰影區(qū)植被的提取實(shí)例中得到了證實(shí)(參見圖5和圖6).

    從實(shí)用性來看,本文的算法除了能為城鎮(zhèn)綠化植物群碳捕獲模型的遙感驅(qū)動提供必須的基礎(chǔ)數(shù)據(jù)外,還有可能服務(wù)于需要分離、提取和修復(fù)陰影區(qū)植被亮度的任何應(yīng)用領(lǐng)域.在實(shí)際應(yīng)用中,可以結(jié)合使用I方法和NDUI方法分離和提取陰影區(qū)植物.可以用I方法提取整個(gè)陰影區(qū)綠地范圍,以幫助定量一些植物生長的環(huán)境壓力參數(shù),如遮蔽率等.同時(shí)用NDUI方法提取陰影區(qū)樹冠邊界,并對其進(jìn)行亮度修復(fù),以幫助樹冠的定性和定量信息提取,如統(tǒng)計(jì)株數(shù)、計(jì)算冠徑、識別樹種和計(jì)算樹冠郁閉飽和度等.

    3 結(jié) 論

    本文提出了歸一化陰影指數(shù)NDUI的概念以及城鎮(zhèn)地區(qū)陰影區(qū)植被提取和亮度修復(fù)的算法.經(jīng)過樣本統(tǒng)計(jì)分析和遙感圖像實(shí)例測試,證明該算法具有很好的穩(wěn)健性.該算法以歸一化植被指數(shù)NDVI確定植被總覆蓋,和以NDUI分離和提取陰影區(qū)樹冠.在陰影區(qū)植被提取中,通過“I方法”與“NDUI方法”的比較,證明前者對提取陰影區(qū)植被總覆蓋有幫助;后者對于陰影區(qū)樹冠提取性能優(yōu)異.所以,可以根據(jù)需要選用適當(dāng)?shù)姆椒?當(dāng)強(qiáng)調(diào)提取整個(gè)陰影區(qū)綠地范圍時(shí),比如自動計(jì)算陰影區(qū)遮蔽率等,使用I方法比較有效,但需要針對遙感圖像類型和獲取的環(huán)境條件等做統(tǒng)計(jì)分析來確定I的閾值,還需要防止陰影區(qū)過提取問題.當(dāng)側(cè)重于單植株分析時(shí),比如要求自動統(tǒng)計(jì)株數(shù)、計(jì)算冠徑、識別樹種和計(jì)算樹冠郁閉飽和度等,則NDUI方法有無可替代的優(yōu)勢.本文還通過亮度修復(fù)試驗(yàn),證明對于陰影區(qū)植被采用自適應(yīng)直方圖均衡化亮度修復(fù),能獲得比較自然和接近亮區(qū)植被的色彩.同時(shí)這個(gè)試驗(yàn)也直觀地顯示了I方法和NDUI方法在提取對象和范圍上的顯著差別.

    [1] TSAI V J D.A comparative study on shadow compensation of color aerial images in invariant color models[J].IEEE Transactions Geoscience and Remote Sensins,2006,44(6):1661-1667.

    [2] YOUNG S D,KAKARLAPUDI S,UIJT DE HAAG M.A shadow detection and extraction algorithm using digital elevation models and x-band weather radar measurements[J].International Journal of Remote Sensing,2005,26:1531-1549.

    [3] SARABANDI P,YAMAZAKI F,MATSUOKAV M,et al.Shadow detection and radiometric restoration in satellite high resolution images[C]//Proceedings of the IEEE International Geoscience and Remote Sensing Symposium.Anchorage,AK:[s.n.],2004(6):3733-3747.

    [4] FINLAYSON G D,HORDLEY S D,DREW M S.Removing Shadows from Images[M].Berlin:Springer Verlag,2002.

    [5] AKIRA S,AKIO S,HIROYUKI A,et al.Dynamic shadow compensation of aerial images based on color and spatial analysis[C]//The International Conference on Pattern Recognition(ICPR’00).Barcelona:[s.n.],2000(1):317-320.

    [6] 王樹根,郭澤金,李德仁.彩色航空影像上陰影區(qū)域信息補(bǔ)償?shù)姆椒ǎ跩].武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2003,28(5):514-516.

    WANG S G,GUO Z J,LI D R.Shadow compensation of color aerial images[J].Geomatics and Information Sci-ence of Wuhan University,2003,28(5):514-516.

    [7] 何凱,趙紅穎,劉晶晶,等.基于分形及紋理分析的城市遙感影像陰影去除[J].天津大學(xué)學(xué)報(bào),2008,41(7):800-804.

    HE K,ZHAO H Y,LIU J J,et al.Shadow removal of city’s remote sensing image based on fractal and texture analysis[J].Journal of Tianjin University,2008,41(7):800-804.

    [8] 許妙忠,余志惠.高分辨衛(wèi)星影像中陰影的自動提取與處理[J].測繪信息與工程,2003,28(1):20-22.

    XU M Z,YU Z H.Automated extraction of shadow in very-high resolution spatial satellite images[J].Journal of Geomatrics,2003,28(1):20-22.

    [9] DARE P M.Shadow analysis in high-resolution satellite imagery of urban areas[J].Photogrammetric Engineering and Remote Sensing,2005,71:169-177.

    [10] 李存軍,劉良云,王紀(jì)華,等.基于主成分融合信息失真的城市IKONOS影像陰影自動提取研究[J].武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2008,33(9):947-950.

    LI C J,LIU L Y,WANG J H,et al.IKONOS shadow extraction in urban region based on principal component fusion information distortion[J].Geomatrics and Information Science of Wuhan University,2008,33(9):947-950.

    [11] 李先華,劉順喜,黃薇,等.月球表面遙感圖像陰影消除及其信息恢復(fù)研究[J].應(yīng)用光學(xué),2009,30(3):423-426.

    LI X H,LIU S X,HUANG W,et al.Shadow removal and information recovery for remote sensing images of lunar surface[J].Journal of Applied Optics,2009,30(3):423-426.

    [12] 高龍華.基于 Quickbird影像的陰影校正模型研究[J].重慶郵電學(xué)院學(xué)報(bào),2005,17(2):226-229.

    GAO L H.Research into the shadow emendation model based on Quickbird image[J].Journal of Chongqing University of Posts and Telecommunications,2005,17(2):226-229.

    [13] 李先華.遙感信息的地形影響與改正[J].測繪學(xué)報(bào),1986,15(2):102-109.

    LI X H.A radiometric correction of relief disturbances in the remote sensed data[J].Acta Geodaetica et Cartographica Sinica,1986,15(2):102-109.

    [14] LIU H Q,HUETE A R.A feedback based modification of the NDVI to minimize soil and atmospheric noise[J].IEEE Transactions Geoscience Remote Sensing,1995,33:457-465.

    [15] HUETE A R.A soil-adjusted vegetation index(SAVI)[J].Remote Sensing Environment,1988,25:295-309.

    [16] GAO X,HUETE A R,NI W.Optical-biophysical relationships of vegetation spectra without background contamination[J].Remote Sensing Environment,2000,74(3):609-620.

    [17] STEFANOYIANNIS A P,COSTARIDOU L,SKIADOPOULOS S.A digital equalization technique improving visualization of dense mammary gland and breast periphery in mammography[J].European Journal of Radiology,2003,45:139-149.

    [18] 張友水,馮學(xué)智,都金康,等.IKONOS影像在城市綠地提取中的應(yīng)用[J].地理研究,2004,23(2):275-280.

    ZHANG Y S,F(xiàn)ENG X Z,DU J K,et al.The applications of IKONOS images in the extraction of urban landscape vegetation[J].Geographical Research,2004,23(2):275-280.

    Methods of extracting distribution information of plants at urban darken areas and repairing their brightness

    ZHOU Jian-hua1, ZHOU Yi-fan2, GUO Xiao-hua1, REN Zhen1

    (1.Key Laboratory of Geographical Information Science,Ministry of Education,East China Normal University,Shanghai 200062,China;2.Department of Energy Resources Engineering,School of Earth Science,Stanford University,CA94305-2220,USA)

    In the urban built-up areas with densely high-rising buildings,about 20%of landscape vegetation is growing at darken areas shaded by buildings or other higher plants.In order to drive a model of measuring captured carbon by remote sensing,one of the puzzles to be solved firstly is to separate and extract those plants at darken areas.So a new conception of NDUI(Normalized Difference Umbra Index)has been put forward here.With NDUI,the total area of landscape vegetation is extracted from an infrared remote sensing image by means of NDVI(Normalized Difference Vegetation Index),and then the plants at darken areas are separated from the total area through another segmentation threshold referring to NDUI.These separated plants at darken area can be enhanced with an algorithm of repairing brightness.It is demonstrated by experiments that the algorithm proposed here can work robustly to meet the requirement.

    urban landscape vegetation; captured carbon; NDUI(Normalized Difference Umbra Index); repairing brightness

    TP79;TP751.1

    A

    10.3969/j.issn.1000-5641.2011.06.001

    1000-5641(2011)06-0001-09

    2010-08

    國家自然科學(xué)基金(40671177,41071275)

    周堅(jiān)華,女,副教授,研究方面為圖像識別、生態(tài)遙感.E-mail:jhzhou@geo.ecnu.edu.cn.

    猜你喜歡
    植被指數(shù)樹冠陰影
    你來了,草就沒有了陰影
    文苑(2020年11期)2020-11-19 11:45:11
    樹冠羞避是什么原理?
    中國詩歌(2019年6期)2019-11-15 00:26:47
    榕樹
    樹冠
    文學(xué)港(2019年5期)2019-05-24 14:19:42
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    讓光“驅(qū)走”陰影
    陰影魔怪
    主要植被指數(shù)在生態(tài)環(huán)評中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29
    国产精品久久久久久人妻精品电影 | 精品卡一卡二卡四卡免费| 久久精品人人爽人人爽视色| 久久久久网色| av片东京热男人的天堂| 夫妻性生交免费视频一级片| a级片在线免费高清观看视频| 久久久国产欧美日韩av| 欧美精品亚洲一区二区| 精品卡一卡二卡四卡免费| 高清av免费在线| 亚洲一区二区三区欧美精品| 黄片小视频在线播放| 免费黄色在线免费观看| 18禁裸乳无遮挡动漫免费视频| av在线观看视频网站免费| 99re6热这里在线精品视频| 丝袜在线中文字幕| 最近手机中文字幕大全| 免费久久久久久久精品成人欧美视频| 中文字幕另类日韩欧美亚洲嫩草| 成年动漫av网址| 51午夜福利影视在线观看| 国产精品久久久久久精品电影小说| 成人影院久久| 久久这里只有精品19| 精品福利永久在线观看| 成年动漫av网址| 国产精品久久久久久精品电影小说| 国产片特级美女逼逼视频| 在线观看免费高清a一片| 啦啦啦中文免费视频观看日本| 免费不卡黄色视频| 在线观看免费日韩欧美大片| 欧美日本中文国产一区发布| 国产xxxxx性猛交| 伊人久久大香线蕉亚洲五| 国产日韩欧美视频二区| 亚洲自偷自拍图片 自拍| 可以免费在线观看a视频的电影网站 | 国产成人一区二区在线| 亚洲少妇的诱惑av| 青春草国产在线视频| 精品少妇一区二区三区视频日本电影 | 亚洲七黄色美女视频| 国产福利在线免费观看视频| 国产一区二区激情短视频 | 国产一区二区三区av在线| 国产伦人伦偷精品视频| 久久久久久久久久久免费av| 日韩 欧美 亚洲 中文字幕| 国产日韩欧美亚洲二区| 国产爽快片一区二区三区| 黄色视频在线播放观看不卡| 精品少妇内射三级| 国产成人啪精品午夜网站| 只有这里有精品99| 欧美xxⅹ黑人| 午夜福利在线免费观看网站| 亚洲中文av在线| 中文字幕人妻丝袜一区二区 | 国产在线视频一区二区| 欧美黄色片欧美黄色片| 日韩精品免费视频一区二区三区| 国语对白做爰xxxⅹ性视频网站| 日韩中文字幕视频在线看片| 男女边摸边吃奶| 亚洲欧美精品自产自拍| 欧美在线黄色| 国产成人系列免费观看| 亚洲国产精品国产精品| 亚洲专区中文字幕在线 | 高清视频免费观看一区二区| 丝瓜视频免费看黄片| 狂野欧美激情性bbbbbb| 国产免费一区二区三区四区乱码| 亚洲国产中文字幕在线视频| 日本av免费视频播放| 男女下面插进去视频免费观看| 你懂的网址亚洲精品在线观看| 亚洲男人天堂网一区| 美国免费a级毛片| 国产精品免费大片| 亚洲av中文av极速乱| 成年人免费黄色播放视频| 亚洲国产精品国产精品| 老司机亚洲免费影院| 在线观看免费高清a一片| 成人18禁高潮啪啪吃奶动态图| 国产成人免费观看mmmm| 欧美日韩视频精品一区| 黄频高清免费视频| 亚洲国产精品999| 国产日韩欧美亚洲二区| 久久久久国产精品人妻一区二区| 啦啦啦在线免费观看视频4| 国产精品一国产av| 日韩av不卡免费在线播放| 亚洲视频免费观看视频| 熟妇人妻不卡中文字幕| 国产精品久久久久久精品电影小说| 亚洲精品av麻豆狂野| 一级片'在线观看视频| 日日撸夜夜添| 国产女主播在线喷水免费视频网站| 亚洲第一青青草原| kizo精华| 欧美97在线视频| 亚洲在久久综合| 久久精品aⅴ一区二区三区四区| 日本91视频免费播放| 欧美成人午夜精品| 婷婷色综合www| 毛片一级片免费看久久久久| 91老司机精品| 在线观看国产h片| 精品卡一卡二卡四卡免费| 中文字幕高清在线视频| 婷婷色综合www| 超碰成人久久| 国产精品99久久99久久久不卡 | 天天躁狠狠躁夜夜躁狠狠躁| 国产精品蜜桃在线观看| 黄色 视频免费看| 国产成人欧美| 激情视频va一区二区三区| 亚洲精品成人av观看孕妇| 日韩制服丝袜自拍偷拍| 日韩电影二区| 亚洲在久久综合| 国产伦人伦偷精品视频| 国产av精品麻豆| 久久国产精品大桥未久av| 国产免费现黄频在线看| 久久久久久久久久久免费av| 免费在线观看完整版高清| 另类亚洲欧美激情| 巨乳人妻的诱惑在线观看| 久久精品国产a三级三级三级| 99久久人妻综合| 亚洲五月色婷婷综合| 国产老妇伦熟女老妇高清| 最近手机中文字幕大全| 高清欧美精品videossex| 亚洲国产最新在线播放| 久久99热这里只频精品6学生| netflix在线观看网站| 美国免费a级毛片| 午夜91福利影院| av卡一久久| 9191精品国产免费久久| 日日摸夜夜添夜夜爱| 国产99久久九九免费精品| 少妇猛男粗大的猛烈进出视频| 国语对白做爰xxxⅹ性视频网站| 一本大道久久a久久精品| 精品卡一卡二卡四卡免费| av卡一久久| 国产片特级美女逼逼视频| 18禁观看日本| 日韩电影二区| 90打野战视频偷拍视频| 午夜福利网站1000一区二区三区| 亚洲精品第二区| 一级a爱视频在线免费观看| 久久久久久人人人人人| tube8黄色片| 亚洲男人天堂网一区| av有码第一页| a级毛片在线看网站| 美女主播在线视频| 午夜激情久久久久久久| 丁香六月欧美| 侵犯人妻中文字幕一二三四区| 久久久精品区二区三区| 新久久久久国产一级毛片| 亚洲天堂av无毛| 免费观看a级毛片全部| 成人国产av品久久久| 午夜精品国产一区二区电影| 亚洲精品国产区一区二| 97在线人人人人妻| h视频一区二区三区| 国产一区有黄有色的免费视频| 婷婷色av中文字幕| 国产乱人偷精品视频| 在线观看一区二区三区激情| 在线观看免费视频网站a站| 精品人妻一区二区三区麻豆| 男人舔女人的私密视频| 亚洲欧美色中文字幕在线| 亚洲成人国产一区在线观看 | 丝袜喷水一区| 免费看不卡的av| videosex国产| 色精品久久人妻99蜜桃| 国产精品成人在线| www.精华液| 精品卡一卡二卡四卡免费| 欧美日韩精品网址| 国产亚洲av片在线观看秒播厂| 国产在线视频一区二区| 中文天堂在线官网| 日本vs欧美在线观看视频| 久久影院123| 啦啦啦视频在线资源免费观看| 成年人午夜在线观看视频| 欧美日韩亚洲高清精品| 另类亚洲欧美激情| 色网站视频免费| 十八禁人妻一区二区| 高清av免费在线| 无限看片的www在线观看| 纯流量卡能插随身wifi吗| 亚洲国产精品999| 日韩av在线免费看完整版不卡| 一级片'在线观看视频| 一本大道久久a久久精品| 国产成人免费观看mmmm| 午夜日韩欧美国产| 亚洲av在线观看美女高潮| 成人亚洲欧美一区二区av| 视频区图区小说| 国产高清国产精品国产三级| 国产精品久久久久久精品电影小说| 在线观看一区二区三区激情| 国产成人精品久久二区二区91 | 99精国产麻豆久久婷婷| 久久性视频一级片| 欧美国产精品一级二级三级| av电影中文网址| 国产亚洲av片在线观看秒播厂| av有码第一页| 天天操日日干夜夜撸| 一本色道久久久久久精品综合| 日日爽夜夜爽网站| 老司机靠b影院| 尾随美女入室| a级毛片在线看网站| 精品久久久久久电影网| 国产人伦9x9x在线观看| 久久久久网色| 在线观看人妻少妇| av.在线天堂| 亚洲国产中文字幕在线视频| 超碰成人久久| 国产免费又黄又爽又色| 免费观看人在逋| 亚洲av成人精品一二三区| 精品少妇久久久久久888优播| 欧美亚洲日本最大视频资源| 人人妻人人爽人人添夜夜欢视频| 久久鲁丝午夜福利片| 日韩精品有码人妻一区| 久久综合国产亚洲精品| 日韩 亚洲 欧美在线| 久久久精品免费免费高清| 一级片免费观看大全| 久久99热这里只频精品6学生| 99久久人妻综合| 一级片免费观看大全| 999精品在线视频| 人人妻人人爽人人添夜夜欢视频| 在线观看三级黄色| 999久久久国产精品视频| www.av在线官网国产| 色吧在线观看| 一区二区三区激情视频| 天天添夜夜摸| 久久午夜综合久久蜜桃| 精品人妻一区二区三区麻豆| 日韩大码丰满熟妇| 9191精品国产免费久久| 久久久精品免费免费高清| 在线观看免费午夜福利视频| 夜夜骑夜夜射夜夜干| avwww免费| 精品一区二区免费观看| 蜜桃国产av成人99| 中文天堂在线官网| 欧美日韩国产mv在线观看视频| 考比视频在线观看| 亚洲成国产人片在线观看| 在线亚洲精品国产二区图片欧美| 亚洲国产欧美一区二区综合| 王馨瑶露胸无遮挡在线观看| 亚洲男人天堂网一区| 日韩大片免费观看网站| 日韩 亚洲 欧美在线| 国产精品香港三级国产av潘金莲 | 亚洲伊人久久精品综合| 在线观看www视频免费| 中文字幕人妻丝袜制服| 国产免费视频播放在线视频| 宅男免费午夜| 国产视频首页在线观看| 日日啪夜夜爽| 一区二区三区精品91| 婷婷色综合www| 久久精品国产亚洲av涩爱| 亚洲国产欧美在线一区| 女的被弄到高潮叫床怎么办| 午夜av观看不卡| 老司机影院成人| 老汉色av国产亚洲站长工具| 一级毛片电影观看| 搡老岳熟女国产| 欧美日韩视频精品一区| 亚洲一区中文字幕在线| 国产成人精品在线电影| 久久久久精品国产欧美久久久 | 一区在线观看完整版| 老司机亚洲免费影院| 国产欧美亚洲国产| 亚洲av电影在线观看一区二区三区| 97在线人人人人妻| 精品国产一区二区久久| 最近的中文字幕免费完整| 国产免费现黄频在线看| 青春草亚洲视频在线观看| 老汉色∧v一级毛片| 国产精品.久久久| 欧美国产精品一级二级三级| av在线老鸭窝| 日韩制服骚丝袜av| 免费看av在线观看网站| 男人操女人黄网站| 曰老女人黄片| 国产精品一区二区在线观看99| 在线观看国产h片| 午夜影院在线不卡| 桃花免费在线播放| 亚洲成人免费av在线播放| 日本色播在线视频| 亚洲欧美日韩另类电影网站| 亚洲国产精品一区二区三区在线| 精品卡一卡二卡四卡免费| 1024香蕉在线观看| av网站在线播放免费| 国产在线一区二区三区精| 大香蕉久久成人网| 一级黄片播放器| 制服丝袜香蕉在线| 热re99久久国产66热| av有码第一页| 国产人伦9x9x在线观看| 搡老岳熟女国产| 人人妻人人澡人人看| 天堂俺去俺来也www色官网| 满18在线观看网站| 一级毛片我不卡| 波野结衣二区三区在线| xxxhd国产人妻xxx| 国产不卡av网站在线观看| 日韩中文字幕视频在线看片| 国产深夜福利视频在线观看| 婷婷色综合大香蕉| 如何舔出高潮| 超色免费av| 五月天丁香电影| 亚洲成人av在线免费| 亚洲精品久久午夜乱码| 久久狼人影院| 黑人猛操日本美女一级片| 波多野结衣av一区二区av| bbb黄色大片| 久久久久久免费高清国产稀缺| 美女大奶头黄色视频| 日韩欧美精品免费久久| 你懂的网址亚洲精品在线观看| 久久韩国三级中文字幕| 久久精品国产亚洲av涩爱| 国产亚洲最大av| 精品国产一区二区三区四区第35| e午夜精品久久久久久久| 日日啪夜夜爽| 人人妻人人澡人人看| 日韩中文字幕欧美一区二区 | 男女下面插进去视频免费观看| 大片免费播放器 马上看| 精品少妇黑人巨大在线播放| 国产精品女同一区二区软件| 午夜免费男女啪啪视频观看| 777米奇影视久久| 青春草国产在线视频| 国产精品国产三级专区第一集| 亚洲欧美中文字幕日韩二区| 欧美xxⅹ黑人| 亚洲四区av| 蜜桃在线观看..| 一区二区av电影网| 9色porny在线观看| 激情五月婷婷亚洲| 韩国av在线不卡| 大香蕉久久网| 一本—道久久a久久精品蜜桃钙片| 国产 精品1| 亚洲男人天堂网一区| 婷婷色综合大香蕉| 女性被躁到高潮视频| 精品酒店卫生间| 2018国产大陆天天弄谢| 亚洲图色成人| 日韩大码丰满熟妇| 国产精品熟女久久久久浪| 欧美激情极品国产一区二区三区| 欧美少妇被猛烈插入视频| 新久久久久国产一级毛片| 欧美精品一区二区免费开放| 亚洲欧洲精品一区二区精品久久久 | 一区二区av电影网| 97精品久久久久久久久久精品| 国产不卡av网站在线观看| 这个男人来自地球电影免费观看 | 免费少妇av软件| 欧美亚洲 丝袜 人妻 在线| 日本爱情动作片www.在线观看| 国产精品麻豆人妻色哟哟久久| 免费高清在线观看视频在线观看| 少妇猛男粗大的猛烈进出视频| 国产淫语在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲三区欧美一区| 亚洲av国产av综合av卡| 欧美成人午夜精品| 久久精品亚洲av国产电影网| 日韩中文字幕欧美一区二区 | 亚洲一码二码三码区别大吗| 国产精品三级大全| 国产男人的电影天堂91| 欧美成人午夜精品| 亚洲精品美女久久av网站| 七月丁香在线播放| 水蜜桃什么品种好| 精品人妻熟女毛片av久久网站| 亚洲av福利一区| 老司机靠b影院| 亚洲av欧美aⅴ国产| 中文字幕高清在线视频| 国产不卡av网站在线观看| 少妇人妻 视频| 在线免费观看不下载黄p国产| 女性生殖器流出的白浆| 男女边吃奶边做爰视频| 黑人巨大精品欧美一区二区蜜桃| 午夜日本视频在线| 亚洲欧美一区二区三区久久| 亚洲精品一二三| 女人精品久久久久毛片| 亚洲人成77777在线视频| 黑人猛操日本美女一级片| 丝瓜视频免费看黄片| 亚洲激情五月婷婷啪啪| 国产一区二区 视频在线| 伊人久久大香线蕉亚洲五| 国产精品 国内视频| 大片免费播放器 马上看| 热99久久久久精品小说推荐| 国产一区二区在线观看av| 国产精品国产三级国产专区5o| 捣出白浆h1v1| 在线观看免费午夜福利视频| 一本大道久久a久久精品| 久久精品国产亚洲av涩爱| xxx大片免费视频| 久久久久久久国产电影| 日韩av在线免费看完整版不卡| 人妻人人澡人人爽人人| 老汉色av国产亚洲站长工具| 少妇人妻久久综合中文| 国产黄频视频在线观看| 欧美人与性动交α欧美软件| 十八禁高潮呻吟视频| 亚洲av中文av极速乱| 一级片'在线观看视频| 999精品在线视频| 精品国产露脸久久av麻豆| 丰满乱子伦码专区| 国产一级毛片在线| 香蕉丝袜av| 中文字幕最新亚洲高清| 无遮挡黄片免费观看| 免费高清在线观看视频在线观看| 亚洲精品国产av成人精品| 老汉色av国产亚洲站长工具| 夫妻午夜视频| 女人高潮潮喷娇喘18禁视频| 亚洲美女黄色视频免费看| 2021少妇久久久久久久久久久| 国产欧美亚洲国产| 国产视频首页在线观看| 午夜影院在线不卡| 国产国语露脸激情在线看| 99热全是精品| 天堂中文最新版在线下载| 国产成人av激情在线播放| 国产又爽黄色视频| 日韩不卡一区二区三区视频在线| 精品少妇一区二区三区视频日本电影 | 少妇的丰满在线观看| 99九九在线精品视频| 精品一区在线观看国产| 丝袜在线中文字幕| 免费高清在线观看日韩| 91国产中文字幕| av不卡在线播放| 精品人妻熟女毛片av久久网站| 男女边吃奶边做爰视频| 亚洲精品国产区一区二| 欧美在线一区亚洲| 99久久综合免费| 国产一区二区三区综合在线观看| www.精华液| 国产亚洲精品第一综合不卡| 别揉我奶头~嗯~啊~动态视频 | 97精品久久久久久久久久精品| 大香蕉久久网| 亚洲综合色网址| 成年av动漫网址| 亚洲激情五月婷婷啪啪| 欧美黑人欧美精品刺激| 日韩av不卡免费在线播放| 999久久久国产精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲,欧美,日韩| 中文字幕精品免费在线观看视频| 亚洲精品,欧美精品| 丁香六月天网| 一级片'在线观看视频| 精品少妇黑人巨大在线播放| 国产精品熟女久久久久浪| 亚洲美女黄色视频免费看| av在线老鸭窝| 午夜免费鲁丝| 免费女性裸体啪啪无遮挡网站| 麻豆乱淫一区二区| 最近最新中文字幕免费大全7| 欧美激情 高清一区二区三区| 国产精品 欧美亚洲| 国产av一区二区精品久久| 街头女战士在线观看网站| 免费在线观看黄色视频的| 秋霞在线观看毛片| 欧美精品人与动牲交sv欧美| 五月天丁香电影| 在线天堂最新版资源| 夜夜骑夜夜射夜夜干| 精品免费久久久久久久清纯 | 国产淫语在线视频| 国产高清国产精品国产三级| 亚洲欧美成人精品一区二区| 中文字幕最新亚洲高清| 在线精品无人区一区二区三| 精品第一国产精品| 欧美日本中文国产一区发布| 91精品国产国语对白视频| 国产精品一区二区在线观看99| 美国免费a级毛片| 欧美在线一区亚洲| 在线免费观看不下载黄p国产| 欧美日韩视频精品一区| 自线自在国产av| 在线观看免费午夜福利视频| 亚洲少妇的诱惑av| 国产精品亚洲av一区麻豆 | 一本色道久久久久久精品综合| 国产xxxxx性猛交| 伦理电影免费视频| 91精品三级在线观看| 国产1区2区3区精品| 在线天堂最新版资源| 国产黄色视频一区二区在线观看| 亚洲精品久久久久久婷婷小说| 另类精品久久| 久久精品久久久久久噜噜老黄| 80岁老熟妇乱子伦牲交| 午夜91福利影院| 只有这里有精品99| 日韩视频在线欧美| 国产成人a∨麻豆精品| 又大又爽又粗| 波多野结衣av一区二区av| 亚洲精品久久午夜乱码| 电影成人av| 国产日韩欧美亚洲二区| 国产亚洲精品第一综合不卡| 精品少妇久久久久久888优播| 交换朋友夫妻互换小说| 精品亚洲乱码少妇综合久久| 成人18禁高潮啪啪吃奶动态图| 欧美日韩亚洲综合一区二区三区_| 超碰97精品在线观看| 美女扒开内裤让男人捅视频| 久久久国产精品麻豆| 韩国av在线不卡| xxx大片免费视频| 久久久国产精品麻豆| 亚洲熟女毛片儿| 大片电影免费在线观看免费| 欧美最新免费一区二区三区| 久久精品久久久久久久性| 亚洲国产欧美在线一区| 90打野战视频偷拍视频| 热re99久久精品国产66热6| 日韩av免费高清视频| 老司机影院毛片| 亚洲av综合色区一区| 日本猛色少妇xxxxx猛交久久| 最近中文字幕2019免费版| 国产一级毛片在线| av一本久久久久| 99精国产麻豆久久婷婷| 老司机靠b影院| 国产亚洲最大av| 在线看a的网站|