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

    基于Landsat8 OLI/TIRS數據的重慶市不透水面與城市熱環(huán)境關系研究

    2017-02-09 05:52:42胡日查王曉婷

    鄧 睿,胡日查,劉 亮,王曉婷

    (1.重慶交通大學 建筑與城市規(guī)劃學院,重慶 400074;2.重慶交通大學 水利水運工程教育部重點實驗室,重慶 400074)

    基于Landsat8 OLI/TIRS數據的重慶市不透水面與城市熱環(huán)境關系研究

    鄧 睿1,2,胡日查1,劉 亮1,王曉婷1

    (1.重慶交通大學 建筑與城市規(guī)劃學院,重慶 400074;2.重慶交通大學 水利水運工程教育部重點實驗室,重慶 400074)

    以2014年Landsat8 OLI/TIRS影像為研究數據,結合遙感與GIS技術,利用Carlson方法提取不透水面,采用Jiménez-Muoz等的劈窗算法反演地表溫度,分析了重慶市主城區(qū)不透水面與熱島效應的分布格局,并對二者進行對比分析,研究了不透水面與熱島效應之間的關系。研究結果表明:重慶市主城區(qū)不透水面與熱島區(qū)域分布較集中,大部分集中在建成區(qū),不透水面與地表溫度呈正相關關系,對城市熱島具有顯著的影響。

    環(huán)境工程;都市熱環(huán)境;不透水面;地表溫度;Landsat8 OLI/TIRS

    0 引 言

    城市熱島效應是指城市中心氣溫明顯高于其周邊郊區(qū)環(huán)境氣溫的一種現象[1]。城市熱島效應不單單意味著城市氣溫升高,更重要的是城市中心會形成一個低壓旋渦,導致人們生產生活所產生的污染氣體聚集在這個漩渦中,使城市氣候更加惡化,影響人類的健康生活[2]。城市不透水面是影響城市熱環(huán)境的重要因素,它通過改變城市表面和邊界層的顯熱和潛熱通量影響城市地表溫度[3]。熱環(huán)境是城市生態(tài)循環(huán)的重要環(huán)節(jié),溫度異常會打亂城市的正常運作,影響人們的正常生活,阻礙經濟發(fā)展。因此,對城市不透水面的研究是尋求緩解熱島效應的有效措施之一。

    目前,不透水面和熱島效應已受到國內外學者的廣泛關注。M.K.RIDD[4]于1995年提出城市地表覆蓋的V-I-S模型,認為城市的地表覆蓋由植被、不透水面、水體和土壤組成;S.R.PHINN等[5]運用約束光譜分離方法成功估算了不透水面的空間分布;YUAN Fei等[6]研究了4個季節(jié)不透水面、植被與地表溫度的關系,結果表明地表溫度和不透水面面積在任何季節(jié)都存在明顯的正相關關系。國內學者岳文澤等[7]從像元尺度與街鎮(zhèn)尺度兩方面,對上海市地表溫度與不透水面覆蓋度的關系進行了研究,結果表明,街鎮(zhèn)尺度上二者具有明確的線性關系;林云杉等[8]利用不透水面與植被覆蓋度二者間的負相關關系,提取了泉州市的不透水面信息,進而對不透水面與城市熱島效應的關系進行了研究;邱建壯等[9]利用分類回歸樹方法對城市不透水面覆蓋度進行了估算,對不透水面覆蓋度(ISP)與地表溫度的正相關關系進行了驗證。

    筆者針對Landsat8 OLI/TIRS影像數據的特點,選擇較佳的提取與反演算法,研究重慶市主城區(qū)不透水面與熱島效應的分布格局,并利用數學方法分析了二者之間的關系。針對目前研究較少的Landsat8 OLI/TIRS影像的地表溫度反演方面,筆者參考了多種算法,采用了劈窗算法,利用TIRS 10,TIRS 11兩個波段消除大氣及時相差別導致的誤差,反演的結果合理,可用于熱島效應研究。

    1 研究區(qū)域與研究數據

    1.1 研究區(qū)域概況

    筆者以重慶主城9區(qū)為研究對象,包括渝中區(qū)、南岸區(qū)、江北區(qū)、九龍坡區(qū)、大渡口區(qū)、沙坪壩區(qū)、渝北區(qū)、巴南區(qū)和北碚區(qū)。其位置如圖1。

    圖1 重慶市主城區(qū)位置Fig.1 Location map of the main city of Chongqing

    重慶主城區(qū)面積有5 473 km2;據2014年的統(tǒng)計資料顯示,主城區(qū)常住人口為818.98×104人。主城區(qū)是人口密集區(qū)域,也是城市的經濟中心。重慶市工業(yè)以化工業(yè)和重工業(yè)為主,大型企業(yè)和工廠相對集中,工業(yè)污染排放使主城區(qū)氣候受到嚴重影響,加之都市區(qū)相對郊區(qū),不透水面比率高,熱島效應較明顯。

    1.2 研究數據

    Landsat 8衛(wèi)星于2013年2月11號發(fā)射,其上攜帶有兩個主要載荷:OLI(全稱:operational land imager,陸地成像儀)和TIRS(全稱:thermal infrared sensor,熱紅外傳感器)。OLI陸地成像儀包括9個波段,空間分辨率為30 m。TIRS傳感器包含兩個分辨率為100 m的熱紅外波段。

    筆者以2014年7月30日與8月6日的四景Landsat8 OLI/TIRS影像為研究數據。重慶主城區(qū)范圍Landsat8 OLI/TIRS影像由四景拼接而成。數據信息如表1。

    表1 遙感影像信息

    2 研究方法

    采用T.N.CARLSON等[10]所提出的方法對重慶主城區(qū)的不透水面進行了提取。該方法主要思想是利用植被覆蓋度與不透水面率的負相關關系,先求得植被覆蓋度;再取相反,得到不透水面率。地表溫度反演受參數影響比較敏感,筆者采取J.C.JIMéNEZ-MUOZ等[11]的劈窗算法進行地表溫度反演,需要確定亮度溫度、地表比輻射率、大氣水汽含量幾個參數。最后利用線性回歸方法,定量分析不透水面率與地表溫度、植被覆蓋度與地表溫度之間的關系。

    2.1 不透水面提取

    根據“建成區(qū)內,不透水面覆蓋度與植被覆蓋度呈顯著負相關”的理論,T.N.CARLSON等[10]提出了一種提取不透水面方法。該方法考慮了兩種地類的不同特性,根據植被覆蓋度求取不透水面覆蓋度,精度較高(圖2)。

    圖2 不透水面提取流程Fig.2 Flow chart of extraction of impervious surface

    由于水體反射率較低,且水中多含有泥沙,在提取不透水面時,會對提取效果有一定的干擾。所以,在進行不透水面提取時,首先要做的就是水體掩膜。

    徐涵秋[12]提出一種改進的歸一化差異水體指數MNDWI(λMNDWI),其精度較高,波段運算如式(1):

    λMNDWI=(λGreen-λMIR)/(λGreen+λMIR)

    (1)

    式中:λGreen代表綠光波段;λMIR代表中紅外波段。

    筆者根據λMNDWI指數進行水體掩膜。

    植被覆蓋度可定義為單位面積內植被的垂直投影面積,它是衡量地表植被狀況的一個最重要的指標。T.N.CARLSON等[10]通過研究,得出了植被覆蓋度Fr的計算公式:

    (2)

    (3)

    歸一化植被指數NDVI(λNDVI)是檢測植被生長狀態(tài)、計算植被覆蓋度和消除部分輻射誤差的一個指數,它能反映出植物冠層的背景影響。Landsat TIRS影像數據的NDVI值利用4,5波段計算[13],即紅光波段(λR)與近紅外波段(λNIR)。式中:λNDVI-veg和λNDVI-soil分別為完全植被覆蓋型像元的NDVI值和完全裸土覆蓋類型像元的NDVI值,即純植被像元與純裸圖像元的NDVI值,二者分別代表完全被植被覆蓋的區(qū)域和完全裸露的區(qū)域,一般情況下,可利用NDVI值的最大值與最小值代替λNDVI-veg和λNDVI-soil的取值。

    在計算得到植被覆蓋度Fr后,利用Carlson方法計算不透水面率ISA(λISA)。具體計算方法如式(4):

    λISA=(1-Fr)dev

    (4)

    式中:λISA為不透水面率;Fr為植被覆蓋度;dev代表只適用于建成區(qū)。

    得到不透水面率后,采用多種閾值分離影像,并同監(jiān)督分類的結果作比較,確定不透水面的最佳分離閾值,進而分離出不透水面。

    2.2 地表溫度反演

    地面除了反射掉太陽光,會吸收一部分熱量,而地面的實測溫度就是地表溫度。地表溫度反演方法有很多種,針對Landsat8 OLI/TIRS影像,可利用USGS官網提供的光譜輻射值和定標參數等波段參數進行地表溫度反演(圖3)。徐涵秋[14]曾對比研究了基于Landsat8 OLI/TIRS影像多種算法,結果表明利用單通道算法要優(yōu)于其他算法。但筆者經過多次試驗發(fā)現,由于四景影像在時間上稍有差異,加之單通道算法只對一個波段進行運算,導致單通道算法反演出的溫度對兩期影像的效果明顯不同:7月30日的兩景影像的溫度明顯偏低,且對比起來,同一地物的溫度存在較大差異。而Landsat OLI有兩個熱紅外波段,非常適合利用劈窗算法,可以利用兩個熱紅外波段的不同吸收作用差值消除大氣影響。

    圖3 地表溫度反演流程Fig.3 Flow chart of inverse calculation of surface temperature

    輻射定標是將記錄的原始DN值轉換為大氣外層表面反射率,目的是消除傳感器本身產生的誤差。這樣可以統(tǒng)一對不同傳感器、不同日期影像進行定量比較。其輻射定標公式如式(5):

    Lλ=MLQcal+AL

    (5)

    式中:Lλ為λ波段的大氣頂部光譜輻射值;ML為λ波段的調整因子,可從影像數據的頭文件中得到,如10波段的調整因子值為頭文件中語句“RADIANCE_MULT_BAND_10”后的值,具體值為3.342 0E-04;AL為λ波段的調整參數,也可從頭文件中得到,如10波段的調整參數為語句“RADIANCE_ADD_BAND_ 10”后的值,具體值為0.100 00;Qcal為影像的灰度值即DN值。

    亮度溫度是遙感器所觀測到的熱輻射強度所對應的溫度,該值結果受大氣以及地表對熱輻射傳導的影響,雖不能代表地表溫度值,但它是溫度反演的重要參數。其計算量溫如式(6)。

    (6)

    式中:Lλ是輻射亮度;K1和K2是元數據熱轉換常數,可從影像頭文件中得到,為語句“K1_CONSTANT_BAND_x”和“K2_CONSTANT _BA ND_x”后的數值,x表示波段。對10波段,K1=774.89,K2=1 321.08;對11波段,K1=480.89,K2=1 201.14。K1單位為W/(m2·srad·μm),K2單位為K。

    在劈窗算法中,大氣水汽含量是一個非常重要的參數。其反演公式考慮到大氣對輻射的影響,利用當時的大氣水汽含量進行校正。大氣水汽含量值可以利用當時的氣溫和相對濕度進行推算,但這樣推算的結果比較粗糙,劈窗算法對參數的準確性比較敏感。楊槐[15]在研究中利用MODTRAD模型模擬了大氣水汽含量與大氣透過率的關系,并得到了二者的線性轉換公式。而NASA官網提供了Landsat OLI影像10波段的實時大氣透過率,根據其推算出來的值比較準確。所以筆者利用此方法計算了四景影像的大氣水汽含量。大氣水汽含量與10波段大氣透過率的轉換公式為

    w=-(τ10-1.040 2)/0.106 7

    (7)

    式中:τ10為10波段的大氣透過率;w為大氣水汽含量,g/cm2。

    將四景影像的成像時間和中心經緯度輸入網站的輸入框中,即可查詢當時10波段對應的大氣透過率。127/39,127/40,128/39,128/40四景影像的查詢結果分別是:0.49,0.56,0.44,0.47;計算得到大氣水汽含量分別是:5.18,4.53,5.65,5.37。

    覃志豪等[16]提出的地表比輻射率的經驗計算公式是基于NDVI指數計算各項參數的,植被覆蓋度計算公式為

    (8)

    利用Carlson方法提取不透水面時,根據不透水面與植被覆蓋區(qū)的負相關性,在λNDVI-veg和λNDVI-soil的取值上大多取最大值和最小值兩個極端值;而在利用覃志豪等[16]經驗公式時,有些研究者通常取經驗值0.7和0.05。而筆者考慮到四景影像的不同時相問題,分別對四景影像算取了植被覆蓋度,一概取經驗值會造成一定的誤差,所以在λNDVI-veg和λNDVI-soil二者的取值上,結合目視判斷,取NDVI累計百分比的95%和5%所在的值,當λNDVI≥λNDVI-veg時,植被覆蓋度Pv=1,即視為完全植被覆蓋區(qū)域;λNDVI≤λNDVI-soil時,植被覆蓋度Pv=0,視為完全裸露區(qū)域。

    計算地表比輻射率時需要植被溫度比Rv、裸土溫度比Rs和熱輻射校正項dε這3個參數。這3項參數均由植被覆蓋度計算而得,如式(9)、式(10)。

    Rv=0.933 2+0.058 5Pv

    (9)

    Rs=0.990 2+0.106 8Pv

    (10)

    式中:Rv為植被溫度比;Rs為裸土溫度比;dε為熱輻射校正項,根據植被覆蓋度的不同取值范圍賦值為:當Pv=0或1時,dε=0;當0

    根據以上各項參數,即可利用覃志豪等[16]提出的地表比輻射率的經驗計算公式計算地表比輻射率,計算如式(11):

    εi=PvRvεi,v+(1-Pv)RSεi,s+dε

    (11)

    式中:εi為i波段的地表比輻射率;εi, v為混合像元中植被的i波段地表比輻射率;εi,s為混合像元中裸土的i波段地表比輻射率。

    在利用覃志豪等[16]經驗公式計算比輻射率時,根據常用地物比輻射率光譜以及Landsat8OLI/TIRS數據特點,對10,11波段的取值:ε10,v=0.986 72,ε10,s=0.967 67,ε11,v=0.989 90,ε10,s=0.977 90[17]。另外,計算結果要在這二者的取值之間,即εi≥εi, v時,令εi=εi, v,εi≤εi, s時,令εi=εi, s。

    將以上計算的各項參數,帶入劈窗算法中,即可求出地表溫度,具體計算如式(12):

    TS=T10+c1(T10-T11)+c2(T10-T11)2+cO+(c3+c4w)(1-ε)+(c5+c6w)Δε

    (12)

    式中:Ts為地表溫度,WK;T10和T11為10,11波段的亮度溫度;w為大氣水汽含量,g/cm2;ε和Δε是平均比輻射率和比輻射率差值;常數參數c0,c1,c2,c3,c4,c5,c6的值分別是:-0.268,1.378,0.183,54.3,-2.238,-129.2,16.4。

    3 結果與分析

    3.1 不透水面分布格局及其特征

    為從不透水面率數據中得到不透水面分離閾值,筆者利用監(jiān)督分類方法將研究區(qū)影像分為“不透水面”與“其他”兩類,將分類結果作為對比數據確定閾值。監(jiān)督分類即利用確定的樣本類型識別未知類別的像元過程[18]。選取25個樣本區(qū),利用最大似然法進行監(jiān)督分類,在經多次分類和ROI調整后,選取精度較高的一次作為分類結果。監(jiān)督分類影像只是作為確定不透水面閾值的參考,所以只需要將建成區(qū)直觀的不透水面,即可直接判斷的建筑物、道路等明顯的不透水面區(qū)域進行分離即可。表2為主城區(qū)范圍內直觀不透水面監(jiān)督分類精度。

    表2 監(jiān)督分類精度

    不透水面的提取需要根據不透水面率確定閾值,進而提取區(qū)域,一般建成區(qū)內的經驗值為70%。筆者考慮到能夠準確提取研究區(qū)域內的不透水面信息,以精度評價的方法,在經驗值附近選取兩個值,以70%,72.5%,75%分別為閾值,將提取的不透水面與不透水面監(jiān)督分類圖進行對比。以3幅不同分類閾值的影像為待檢驗數據,以監(jiān)督分類影像為檢驗數據,進行精度評價(表3)。

    表3 不同閾值分類精度對比Table 3 Comparison of different thresholds classification accuracy /%

    綜合考慮表3中3個閾值分類的各項精度評價指數,最終選取72.5%為分離閾值,從所得的不透水面率影像中,提取出研究區(qū)域不透水面,如表4和圖4。

    表4 主城區(qū)不透水面面積及百分比

    圖4 不透水面率Fig.4 Impervious surfaces percentage

    由表4和圖4可以看出:主城區(qū)不透水面在空間分布上大概分為東、中、西這3個部分。東部大部分不透水面集中在沙坪壩區(qū)南部和九龍坡區(qū)北部;中部不透水面面積最大,也最集中,主要分布在渝北區(qū)南部、渝中區(qū)、大渡口區(qū)北部,巴南區(qū)和南岸區(qū)分布較少;東部大部分集中在江北區(qū)。渝北、沙坪壩、九龍坡等區(qū)不透水面面積較高,渝北區(qū)不透水面面積最高,為221.80 km2,而不透水面所占面積百分比最高的區(qū)是渝中區(qū),高達65.17%。因為面積大,植被覆蓋度高,巴南區(qū)不透水面所占百分比最低,只占總面積的5.42%。

    3.2 熱島分布格局及其特征

    檢索中國氣象數據網中國地面國際交換站氣候資料日值數據得知,2014年8月6日重慶沙坪壩氣象臺站(站臺號:57516)監(jiān)測地面0 cm溫度數據顯示:當日最高地表溫度68.7 ℃,最低27.7 ℃,平均地表溫度41.7 ℃。查詢沙坪壩氣象站位置反演數據為36.9 ℃,衛(wèi)星于凌晨03:27分過境,此時地表溫度尚未達到較高值,反演結果在監(jiān)測數據最高值與最低值之間,并接近平均值。

    水體表現為最低溫在25~30 ℃;植被溫度略高于水體,絕大部分城區(qū)溫度分布在35~45 ℃;兩江溫度有較小差異,長江溫度較嘉陵江低1~2 ℃;最高溫出現在江北、渝北、沙坪壩和九龍坡等區(qū)(表5)。兩期影像溫度差異在1~3 ℃不等,為使兩期影像溫度數值能盡量相同,筆者將四景影像的地表溫度標準化,并拉伸到25~45 ℃范圍內。

    表5 主城區(qū)平均地溫統(tǒng)計

    (13)

    式中:Tmin為地表溫度最小值;Tmax為最大值。

    計算結果即歸一化并拉伸后的結果,對結果進行影像鑲嵌和掩膜工作,得到最終的地表溫度影像(圖5)。

    圖5 重慶市地表溫度Fig.5 Surface temperature of Chongqing

    將地表溫度的歸一化值均分為6個等級,從高到低依次為高溫區(qū)、亞高溫區(qū)、中溫區(qū)、弱低溫區(qū)、亞低溫區(qū)、低溫區(qū),并統(tǒng)計各區(qū)各溫度級別的面積百分比,如表6。

    表6 主城區(qū)各溫度級別面積所占百分比統(tǒng)計

    九龍坡區(qū)高溫區(qū)面積最大,為6.2 km2;渝中區(qū)高溫區(qū)面積及其所占百分比最低,為0。但渝中區(qū)平均溫度較高,因為其低溫區(qū)百分比較低,溫度多集中在中溫區(qū),導致其平均溫度較高。巴南區(qū)弱低溫區(qū)與亞低溫區(qū)面積較大,所以其平均溫度最低;而大渡口區(qū)中溫區(qū)與亞高溫區(qū)面積百分比較高,所以,其平均溫度最高。中心城區(qū)并非溫度最高的區(qū)域,經過查詢地圖得知,高溫區(qū)大多出現在在大型工業(yè)園區(qū),工業(yè)園區(qū)內的生產設施散熱及工業(yè)排放均會造成高溫,而城市區(qū)溫度處于中溫區(qū)以上等級,植被覆蓋區(qū)與水域均在較低等級。

    地表溫度作為城市熱環(huán)境監(jiān)測與研究的重要參數之一,可以利用其有效的模擬城市熱島的格局分布,并對熱島效應的分布特征進行分析。熱島強度是反應城市溫度差異的一個指數,以地表溫度計算出熱島強度指數,并對其進行閾值分類,可以得到不同強度城市熱島的分布[19],如式(14)。

    THI=(T-Tmean)/Tmean

    (14)

    式中:T為地表溫度;Tmean為研究區(qū)域地表溫度平均值;THI為熱島強度。

    利用ENVI軟件中的Density Slice(密度分割)工具將熱島強度值等分6個等級,計算得到的熱島強度:強綠島、中綠島、弱綠島、弱熱島、中熱島、強熱島(圖6)。

    圖6 熱島強度分級Fig.6 Heat island intensity levels

    主城區(qū)范圍內,茂密植被覆蓋區(qū)域為高等級綠島,城市居住區(qū)為中級熱島??傮w看來,熱島分布比較集中,絕大部分集中在建成區(qū)內,向外熱島等級逐漸降低,南北部最低,若以長江為分界線,長江以北熱島面積較高,長江以南部分熱島面積非常低,而且,熱島多在以江岸向外延伸。

    熱島總面積為1 567.60 km2,占研究區(qū)域總面積的29.16%,強熱島面積為131.07 km2,占研究區(qū)域總面積的2.35%(表7)。渝中區(qū)與巴南區(qū)的熱島百分比分別為最高與最低值為79.55%和13.46%。渝中區(qū)位于主城九區(qū)的中心,在兩江交匯處,其面積較小,植被覆蓋度較低,多為建筑區(qū)域,所以其熱島百分比最高;而巴南區(qū)南部大多為植被覆蓋區(qū)域,只有北部沿江熱島集中,加之區(qū)域面積很大,所以熱島百分比較低。總體而言,沿江的幾個城區(qū)熱島百分比較高,南北幾個區(qū)由于其植被覆蓋面積較大,熱島百分比較低(表8)。

    表7 各強度熱島面積及所占百分比統(tǒng)計

    表8 主城各區(qū)各熱島級別面積所占百分比統(tǒng)計

    3.3 熱島與地表參數的關系分析

    不透水面提高了城市地表溫度,同時影響著熱量散發(fā),對熱環(huán)境有較大負面影響;而植被平衡著城市熱環(huán)境,隨著城市擴張與經濟發(fā)展,城市植被覆蓋逐漸減少,對熱環(huán)境造成了一定的負面影響。為研究不透水面與熱島效應的關系,筆者對不透水面率與地表溫度以及植被覆蓋度與地表溫度進行了相關性分析與回歸分析(圖7)。

    二者的線性回歸方程為:

    y=0.102x1+26.902,R2=0.436

    (15)

    y=-0.055x2+35.658,R2=0.421

    (16)

    式中:y為不透水面率;x1為地表溫度;x2為植被覆蓋度。

    不透水面率與地表溫度呈正相關關系,決定系數為0.436,相關系數為0.661,相關性顯著。由式(15)可知:不透水面率越高,地表溫度越高,不透水面率每升高10%,地表溫度升高1.02 ℃;植被覆蓋度與地表溫度呈負相關關系,決定系數為0.421,相關系數為0.649,相關性較強。由式(16)可知,植被覆蓋度越高,地表溫度越低,植被覆蓋度每升高10%,地表溫度下降0.55 ℃。

    圖7 不透水面率和植被覆蓋度與地表溫度散點圖Fig.7 Scatter diagram of impervious surfaces percentage and vegetation coverage with surface temperature

    4 結 論

    1)重慶主城區(qū)熱島效應比較明顯,尤其是高人口密度與高建筑物密度區(qū)域與工業(yè)區(qū)尤為明顯,人類活動、城市擴張和工業(yè)排熱是造成熱島效應的主要原因。不透水面與城市熱島區(qū)域在空間分布上也具有顯著的一致性,可見,不透水面是影響城市熱環(huán)境的重要因素之一。而植被與水域溫度普遍較低,對于緩解熱島效應具有明顯的作用。

    2)重慶主城區(qū)不透水面主要分布在建成區(qū),其分布與城市擴張緊密相關。從各區(qū)地表溫度來看,地表溫度與下墊面類型關系密切,不透水面區(qū)域地表溫度明顯較高,且不透水面率與地表溫度呈正相關關系,植被覆蓋度與地表溫度呈負相關關系。

    3)筆者針對Landsat8 OLI/TIRS影像數據的特點以JIMéNEZ-MUOZ等的劈窗算法,通過計算亮度溫度、比輻射率,大氣水汽含量反演地表溫度,與監(jiān)測站監(jiān)測數據對比,反演精度可分析熱島效應,對重慶市規(guī)劃建設有一定參考意義。

    [1] 郭儀南,張宏超,王健,等.城市熱島效應對瀝青路面溫度場及其力學性能的影響[J].重慶交通大學學報(自然科學版),2010,29(4):548-551. GUO Yinan, ZHANG Hongchao, WANG Jian, et al. Effect of urban heat island on temperature fields and performance of asphalt concrete pavement[J].JournalofChongqingJiaotongUniversity(NaturalScience), 2010,29(4):548-551.

    [2] 貢璐,呂光輝.綠洲城市熱島效應研究[M].烏魯木齊:新疆人民出版社,2010. GONG Lu, LV Guanghui.TheUrbanHeatIslandStudyofOasis[M]. Urrmqi: Xinjiang People’s Publishing House, 2010.

    [3] 孫立軍.基于城市地表覆被類型變化的城市熱島效應研究[D].杭州:浙江農林大學,2011. SUN Lijun.AnalysisonUrbanHeatIslandEffectBasedontheDynamicalChangeofUrbanSurfaceBiophysicalDescriptors[D].Hangzhou: Zhejiang A & F University, 2011.

    [4] RIDD M K. Exploring a V-I-S(Vegetation-Impervious Surface-Soil) model for urban ecosystem analysis through remote sensing: comparative anatomy for cities[J].InternationalJournalofRemoteSensing, 1995,16(12):2165-2185.

    [5] PHINN S R, STANFORD M,SCARTH P F, et al. Monitoring the composition and form of urban environments based on the vegetation-impervious surface-soil (VIS) model by sub-pixel analysis techniques[J].InternationalJournalofRemoteSensing, 2002,23(20):4131-4153.

    [6] YUAN Fei, BAUER M E. Comparison of impervious surface area and normalized difference vegetation index as indcators of surface urban heat island effects in Landsat imagery[J].RemoteSensingofEnvironment, 2007,106(3):375-386.

    [7] 岳文澤,吳次芳.基于混合光譜分解的城市不透水面分布估算[J].遙感學報,2007,11(6):914-922. YUE Wenze, WU Cifang. Urban impervious surface distribution estimation by spectral mixture analysis[J].JournalofRemoteSensing, 2007,11(6):914-922.

    [8] 林云杉,徐涵秋,周榕.城市不透水面及其與城市熱島的關系研究——以泉州市區(qū)為例[J].遙感技術與應用,2007,22(1):14-19. LIN Yunshan, XU Hanqiu, ZHOU Rong. A study on urban impervious surface area and its relation with urban heat island: Quanzhou City, China[J].RemoteSensingTechnologyandApplication, 2007,22(1):14-19.

    [9] 邱健壯,桑峰勇,高志宏.城市不透水面覆蓋度與地面溫度遙感估算與分析[J].測繪科學,2011,36(4):211-213. QIU Jianzhuang, SANG Fengyong, GAO Zhihong. RS estimating and analysis of urban impervious surface percentage and land surface temperature[J].ScienceofSurveyingandMapping, 2011,36(4):211-213.

    [10] CARLSON T N, ARTHUR S T. The impact of land use—land cover changes due to urbanization on surface microclimate and hydrology: a satellite perspective[J].GlobalandPlanetaryChange, 2000,25(1/2):49-65

    [11] JIMéNEZ-MUOZ J C, SOBRINO J A,SKOKOVIC D,et al. Land surface temperature retrieval methods from landsat-8 thermal infrared sensor data[J].IEEEGeoscience&RemoteSensingLetters, 2014,11(10):1840-1843.

    [12] 徐涵秋.利用改進的歸一化差異水體指數(MNDWI)提取水體信息的研究[J].遙感學報,2005,9(5):589-595. XU Hanqiu. A study on information extraction of water body with the modified normalized difference water index(MNDWI)[J].JournalofRemoteSensing, 2005,9(5):589-595.

    [13] 張風霖,緱變彩,李靖琳.Landsat7 ETM+與Landsat8 OLI植被和非植被定量研究[J].山西建筑,2014,40(11):240-241. ZHANG Fenglin, GOU Biancai, LI Jinglin. Quantitative study of vegetation and non-vegetation Landsat7 ETM+ with Landsat8 OLI[J].ShanxiArchitecture, 2014,40(11):240-241.

    [14] 徐涵秋.新型Landsat 8衛(wèi)星影像的反射率和地表溫度反演[J].地球物理學報,2015,58(3):741-747. XU Hanqiu. Retrieval of the reflectance and land surface temperature of the newly-launched Landsat 8 satellite[J].ChineseJournalofGeophysics, 2015,58(3):741-747.

    [15] 楊槐.從Landsat 8影像反演地表溫度的劈窗算法研究[J].測繪地理信息,2014,39(4):73-77. YANG Huai. Research of split-window algorithm for retrieval of land surface temperature from Landsat8[J].JournalofGeomatics, 2014,39(4):73-77.

    [16] 覃志豪,李文娟,徐斌,等.陸地衛(wèi)星TM6波段范圍內地表比輻射率的估計[J].國土資源遙感,2004(3):28-32. QIN Zhihao, LI Wenjuan, XU Bin, et al. The estimation of land surface emissivity for Landsat TM6[J].RemoteSensingforLand&Resources,2004(3):28-32.

    [17] 李瑤,潘竟虎.基于Landsat8劈窗算法與混合光譜分解的城市熱島效應空間格局分析——以蘭州中心城區(qū)為例[J].干旱區(qū)地理,2015,38(1):111-119. LI Yao, PAN Jinghu. Spatial pattern on urban heat environment using split window algorithm and spectral mixture analysis based on Landsat 8 images: a case of Lanzhou City[J].AridLandGeography, 2015,38(1):111-119.

    [18] 牟鳳云,羅丹,官冬杰,等.面向對象的土地覆蓋信息提取方法研究及應用[J].重慶交通大學學報(自然科學版),2014,33(6):104-108. MU Fengyun, LUO Dan, GUAN Dongjie, et al. Research on object-oriented land cover information extraction[J].JournalofChongqingJiaotongUniversity(NaturalScience), 2014,33(6):104-108.

    [19] 沈德才,楊燕瓊,吳振彪,等.基于Landsat 8的東莞市熱島效應研究[J].廣東林業(yè)科技,2014,30(6):20-24. SHEN Decai, YANG Yanqiong, WU Zhenbiao, et al. Research on urban heat island effect of Dongguan based on Landsat 8[J].GuangdongForestryScienceandTechnology, 2014,30(6):20-24.

    (責任編輯 劉 韜)

    Study on the Relationship between the Impervious Surface and the Urban Thermal Environment in Chongqing Based on Landsat8 OLI/TIRS Image

    DENG Rui1,2,HU Richa1,LIU Liang1,WANG Xiaoting1

    (1. College of Architecture and Urban Planning, Chongqing Jiaotong University, Chongqing 400074,P.R.China;2.Key Laboratory of Waterway Engineering, Chongqing Jiaotong University, Chongqing 400074,P.R.China)

    In this paper, On in combinatiom with remote sensing methods and GIS technology and used the Carlson methods Jiménez-Muoz split window algorithms wasused to extract impervious surface and inverse surface temperature based on Landsat-8 image data in 2014. Then analyzed the distribution pattern of impervious surface and thermal island and the correlation of them was analyzed to study the relationship between impervious surface and heat island effect on the main city of Chongqing. The results show that the impervious surface area and the heat island in main city of Chongqing? are more concentrated and mostly in built-up areas. Impervious surface and surface temperature has a positive correlation and it has a significant effect on the urban heat island.

    environment engineering; urban thermal environment; impervious surfaces; surface temperature; Landsat8 OLI/TIRS

    10.3969/j.issn.1674-0696.2017.01.13

    2015-10-14;

    2015-12-22

    國家自然科學基金項目(51208531);重慶市教委科學技術研究項目(KJ 120405);國土資源部地學空間信息技術重點實驗室開放基金項目(KLGSIT 2015-07);重慶交通大學國家內河航道整治工程技術研究中心暨水利水運工程教育部重點實驗室開放基金項目(SLK2014B04)

    鄧 睿(1983—),女,重慶人,講師,博士,主要從事地理信息系統(tǒng)、遙感技術等方面的研究。E-mail: trdeng@sina.com。

    胡日查(1993—),男(蒙古族),內蒙古興安人,本科生,主要從事地理信息系統(tǒng)方面的學習。E-mail:1412943469@qq.com。

    X87

    A

    1674-0696(2017)01-068-09

    日韩精品免费视频一区二区三区| 久久午夜亚洲精品久久| av国产免费在线观看| 亚洲自偷自拍图片 自拍| 国产精品久久电影中文字幕| 青草久久国产| 长腿黑丝高跟| av免费在线观看网站| 亚洲成av人片免费观看| 亚洲熟妇中文字幕五十中出| 久久国产乱子伦精品免费另类| 欧美精品亚洲一区二区| 黄色丝袜av网址大全| 亚洲人成伊人成综合网2020| 午夜福利免费观看在线| 国产高清激情床上av| 亚洲熟女毛片儿| 久久婷婷人人爽人人干人人爱| 好看av亚洲va欧美ⅴa在| 欧美中文综合在线视频| 色综合站精品国产| 成人av一区二区三区在线看| 搡老熟女国产l中国老女人| 久久精品国产99精品国产亚洲性色| 欧美黑人巨大hd| 亚洲性夜色夜夜综合| 免费电影在线观看免费观看| 在线a可以看的网站| 国产精品电影一区二区三区| 两人在一起打扑克的视频| 日韩av在线大香蕉| 哪里可以看免费的av片| 一区二区三区高清视频在线| 91老司机精品| 亚洲人成电影免费在线| 一个人免费在线观看电影 | 精品国产乱码久久久久久男人| 日韩欧美 国产精品| 怎么达到女性高潮| 日日干狠狠操夜夜爽| 亚洲av美国av| 啦啦啦免费观看视频1| 日本免费一区二区三区高清不卡| 国内精品一区二区在线观看| 亚洲男人的天堂狠狠| 麻豆av在线久日| 在线观看一区二区三区| 在线国产一区二区在线| 精品第一国产精品| 免费在线观看黄色视频的| 欧美日韩中文字幕国产精品一区二区三区| 18禁观看日本| 国产精品一区二区三区四区免费观看 | 欧美日本亚洲视频在线播放| 脱女人内裤的视频| 精品熟女少妇八av免费久了| 国模一区二区三区四区视频 | 亚洲18禁久久av| 久久久久久九九精品二区国产 | 国产伦在线观看视频一区| 非洲黑人性xxxx精品又粗又长| 老汉色av国产亚洲站长工具| av天堂在线播放| 亚洲人成网站在线播放欧美日韩| 午夜日韩欧美国产| 黄色a级毛片大全视频| 婷婷丁香在线五月| 狂野欧美白嫩少妇大欣赏| 免费人成视频x8x8入口观看| 老汉色av国产亚洲站长工具| 美女午夜性视频免费| 色哟哟哟哟哟哟| 欧美高清成人免费视频www| 亚洲国产精品久久男人天堂| 亚洲色图 男人天堂 中文字幕| 黄频高清免费视频| 黄色视频不卡| 男女下面进入的视频免费午夜| 国产av在哪里看| 欧美三级亚洲精品| 最近最新中文字幕大全电影3| 99精品久久久久人妻精品| 十八禁网站免费在线| 亚洲色图av天堂| 国产99久久九九免费精品| 操出白浆在线播放| 免费看十八禁软件| 久久久国产精品麻豆| 无遮挡黄片免费观看| 精品久久久久久久久久久久久| 午夜免费激情av| 露出奶头的视频| 精品福利观看| 制服人妻中文乱码| 91国产中文字幕| 午夜亚洲福利在线播放| 久久精品夜夜夜夜夜久久蜜豆 | 免费电影在线观看免费观看| 国产精品久久久久久精品电影| 色在线成人网| 男女午夜视频在线观看| 精品一区二区三区视频在线观看免费| 精品福利观看| 国产激情偷乱视频一区二区| 中文字幕av在线有码专区| 99久久99久久久精品蜜桃| 91老司机精品| 亚洲欧美日韩高清在线视频| 国产在线精品亚洲第一网站| 制服人妻中文乱码| 蜜桃久久精品国产亚洲av| 可以免费在线观看a视频的电影网站| 久久这里只有精品19| 别揉我奶头~嗯~啊~动态视频| 毛片女人毛片| 麻豆国产av国片精品| a在线观看视频网站| 中国美女看黄片| 夜夜爽天天搞| 啪啪无遮挡十八禁网站| 最近最新免费中文字幕在线| 亚洲精华国产精华精| 亚洲国产欧美一区二区综合| 亚洲国产精品成人综合色| 啪啪无遮挡十八禁网站| 国产亚洲精品久久久久5区| 听说在线观看完整版免费高清| 亚洲av日韩精品久久久久久密| 琪琪午夜伦伦电影理论片6080| 十八禁网站免费在线| 午夜免费观看网址| 深夜精品福利| 亚洲精品国产精品久久久不卡| 亚洲熟妇中文字幕五十中出| 国产三级黄色录像| 久久久水蜜桃国产精品网| 亚洲,欧美精品.| 国产精品久久久久久精品电影| 午夜福利免费观看在线| 亚洲九九香蕉| 日韩三级视频一区二区三区| 国产成人系列免费观看| 久久欧美精品欧美久久欧美| 久久精品成人免费网站| 国产亚洲欧美在线一区二区| 国产一区二区激情短视频| 天堂动漫精品| 长腿黑丝高跟| 国产91精品成人一区二区三区| 草草在线视频免费看| 亚洲国产看品久久| 黄色 视频免费看| 最新在线观看一区二区三区| 非洲黑人性xxxx精品又粗又长| 非洲黑人性xxxx精品又粗又长| aaaaa片日本免费| 欧美丝袜亚洲另类 | 在线观看66精品国产| 国产精品免费视频内射| 亚洲va日本ⅴa欧美va伊人久久| netflix在线观看网站| 中文字幕高清在线视频| 国产午夜精品论理片| 亚洲成av人片在线播放无| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av电影不卡..在线观看| 国产麻豆成人av免费视频| 99精品久久久久人妻精品| 制服诱惑二区| 人成视频在线观看免费观看| 亚洲18禁久久av| 十八禁网站免费在线| 天堂√8在线中文| 亚洲男人天堂网一区| 老司机深夜福利视频在线观看| 国产伦一二天堂av在线观看| 精品高清国产在线一区| 欧美激情久久久久久爽电影| 中文字幕久久专区| 给我免费播放毛片高清在线观看| 免费人成视频x8x8入口观看| 这个男人来自地球电影免费观看| 国产亚洲精品久久久久5区| 国产精品免费一区二区三区在线| 又粗又爽又猛毛片免费看| 亚洲人成网站高清观看| 很黄的视频免费| 成年人黄色毛片网站| 国产亚洲欧美在线一区二区| 国产一区二区在线av高清观看| 亚洲,欧美精品.| 琪琪午夜伦伦电影理论片6080| 搡老妇女老女人老熟妇| 久久精品成人免费网站| 最近最新免费中文字幕在线| 18禁裸乳无遮挡免费网站照片| 亚洲欧美日韩高清专用| 美女 人体艺术 gogo| 午夜成年电影在线免费观看| 亚洲欧美精品综合久久99| 日本三级黄在线观看| 日日爽夜夜爽网站| 成人精品一区二区免费| 一级毛片女人18水好多| 黑人操中国人逼视频| 婷婷亚洲欧美| 久久精品91蜜桃| 好男人电影高清在线观看| 三级国产精品欧美在线观看 | 成人一区二区视频在线观看| 亚洲人成电影免费在线| 在线看三级毛片| 亚洲国产中文字幕在线视频| 99热这里只有是精品50| 亚洲 欧美一区二区三区| 国产成人精品久久二区二区免费| 成人高潮视频无遮挡免费网站| 亚洲欧美精品综合一区二区三区| 精品一区二区三区视频在线观看免费| 此物有八面人人有两片| 国产男靠女视频免费网站| 国产高清激情床上av| 88av欧美| 日本免费a在线| 国产亚洲av嫩草精品影院| 精品欧美一区二区三区在线| 黑人欧美特级aaaaaa片| 中文资源天堂在线| 色老头精品视频在线观看| 国产亚洲精品久久久久久毛片| 天天躁夜夜躁狠狠躁躁| 啦啦啦免费观看视频1| 亚洲精品美女久久av网站| 一个人免费在线观看的高清视频| 中国美女看黄片| 无人区码免费观看不卡| 日韩高清综合在线| 在线观看免费午夜福利视频| 精品午夜福利视频在线观看一区| 999久久久国产精品视频| 国产精品影院久久| 在线观看一区二区三区| 亚洲精品色激情综合| 黄频高清免费视频| 欧美成狂野欧美在线观看| 蜜桃久久精品国产亚洲av| 成人一区二区视频在线观看| 韩国av一区二区三区四区| 国产亚洲av高清不卡| 国产精品久久久久久亚洲av鲁大| 久久婷婷人人爽人人干人人爱| 日本精品一区二区三区蜜桃| 国产亚洲欧美在线一区二区| 日韩大码丰满熟妇| 嫁个100分男人电影在线观看| 欧美人与性动交α欧美精品济南到| 国产69精品久久久久777片 | a级毛片a级免费在线| 欧美日韩亚洲综合一区二区三区_| 亚洲五月天丁香| 亚洲人成电影免费在线| 亚洲第一电影网av| 9191精品国产免费久久| 黄色成人免费大全| 伊人久久大香线蕉亚洲五| 两个人免费观看高清视频| 亚洲无线在线观看| 我要搜黄色片| 久久久久九九精品影院| 欧美一级毛片孕妇| 啦啦啦韩国在线观看视频| 日本成人三级电影网站| 又粗又爽又猛毛片免费看| 黄片小视频在线播放| 久久久久久久精品吃奶| 一进一出抽搐动态| 欧洲精品卡2卡3卡4卡5卡区| 窝窝影院91人妻| 老司机在亚洲福利影院| 亚洲av成人不卡在线观看播放网| 在线视频色国产色| 国产成年人精品一区二区| 久久中文字幕人妻熟女| 国产精品 国内视频| 亚洲人成伊人成综合网2020| 正在播放国产对白刺激| 国内精品一区二区在线观看| 国产高清视频在线播放一区| 久久亚洲精品不卡| 日韩欧美免费精品| 免费在线观看亚洲国产| 制服人妻中文乱码| 不卡av一区二区三区| 欧美一区二区精品小视频在线| 国产精品一及| а√天堂www在线а√下载| 久久午夜亚洲精品久久| √禁漫天堂资源中文www| 亚洲av五月六月丁香网| 亚洲va日本ⅴa欧美va伊人久久| av免费在线观看网站| 午夜福利18| 99在线人妻在线中文字幕| 国产私拍福利视频在线观看| 91av网站免费观看| 19禁男女啪啪无遮挡网站| 免费搜索国产男女视频| 三级毛片av免费| 日韩欧美国产在线观看| 亚洲av片天天在线观看| 妹子高潮喷水视频| 久久婷婷人人爽人人干人人爱| 大型av网站在线播放| 给我免费播放毛片高清在线观看| 久久久精品国产亚洲av高清涩受| 亚洲在线自拍视频| 国产99白浆流出| 亚洲美女视频黄频| 久久人妻福利社区极品人妻图片| 久9热在线精品视频| www.精华液| 成人国产综合亚洲| 亚洲真实伦在线观看| 女生性感内裤真人,穿戴方法视频| cao死你这个sao货| 午夜老司机福利片| 久久国产精品人妻蜜桃| 亚洲午夜精品一区,二区,三区| 999久久久国产精品视频| 亚洲激情在线av| 啪啪无遮挡十八禁网站| 老司机午夜十八禁免费视频| 日日爽夜夜爽网站| 国产亚洲精品久久久久久毛片| 久久亚洲精品不卡| 制服丝袜大香蕉在线| 香蕉国产在线看| 亚洲天堂国产精品一区在线| 搞女人的毛片| 婷婷精品国产亚洲av在线| 亚洲熟妇中文字幕五十中出| 夜夜爽天天搞| 久久久久国产一级毛片高清牌| 国产精品电影一区二区三区| 亚洲人成网站在线播放欧美日韩| 午夜a级毛片| 搡老妇女老女人老熟妇| 久久这里只有精品19| 日日爽夜夜爽网站| 国产一区二区三区视频了| 两个人免费观看高清视频| 色综合站精品国产| 久久亚洲真实| 久久午夜亚洲精品久久| 午夜久久久久精精品| 在线观看免费视频日本深夜| 欧美高清成人免费视频www| 国产99白浆流出| 久久国产乱子伦精品免费另类| 久久精品综合一区二区三区| 久久性视频一级片| 成人18禁在线播放| 国语自产精品视频在线第100页| 男女下面进入的视频免费午夜| 性色av乱码一区二区三区2| 久久久久久亚洲精品国产蜜桃av| 欧美乱码精品一区二区三区| 日韩中文字幕欧美一区二区| 黄片小视频在线播放| 亚洲色图 男人天堂 中文字幕| 成人一区二区视频在线观看| www.自偷自拍.com| 亚洲精品粉嫩美女一区| 最近在线观看免费完整版| 两个人看的免费小视频| 亚洲中文字幕一区二区三区有码在线看 | 精品第一国产精品| 日韩欧美 国产精品| 俄罗斯特黄特色一大片| 最新在线观看一区二区三区| 国产免费男女视频| 午夜福利18| 精品免费久久久久久久清纯| 色综合欧美亚洲国产小说| 日韩国内少妇激情av| 一a级毛片在线观看| 久久中文字幕一级| 国产真实乱freesex| 丁香六月欧美| 夜夜看夜夜爽夜夜摸| avwww免费| 日本 欧美在线| 91大片在线观看| 一边摸一边做爽爽视频免费| avwww免费| 成人手机av| 国产激情偷乱视频一区二区| 亚洲专区中文字幕在线| 久久久精品欧美日韩精品| 熟女少妇亚洲综合色aaa.| 一二三四在线观看免费中文在| 99精品欧美一区二区三区四区| 久久久久九九精品影院| 一级毛片女人18水好多| 少妇粗大呻吟视频| 少妇裸体淫交视频免费看高清 | 窝窝影院91人妻| 国产三级在线视频| 啪啪无遮挡十八禁网站| 亚洲人成网站高清观看| 一区二区三区高清视频在线| 国产亚洲欧美在线一区二区| 国产69精品久久久久777片 | 欧美性长视频在线观看| 国产精品av视频在线免费观看| 亚洲人成77777在线视频| 女警被强在线播放| 级片在线观看| 成人永久免费在线观看视频| 中文字幕人成人乱码亚洲影| 这个男人来自地球电影免费观看| 国产精品电影一区二区三区| 18美女黄网站色大片免费观看| 久久国产精品人妻蜜桃| 成人国产一区最新在线观看| 久久精品国产亚洲av高清一级| 国产又黄又爽又无遮挡在线| 国产爱豆传媒在线观看 | 免费一级毛片在线播放高清视频| 午夜影院日韩av| www.www免费av| 一级毛片精品| 欧美在线黄色| 香蕉国产在线看| 特级一级黄色大片| 88av欧美| 亚洲成a人片在线一区二区| 手机成人av网站| 天堂动漫精品| 婷婷丁香在线五月| 女同久久另类99精品国产91| 久久久久久久午夜电影| 婷婷丁香在线五月| 亚洲人成网站在线播放欧美日韩| 中文字幕熟女人妻在线| 欧美中文综合在线视频| 一二三四在线观看免费中文在| 变态另类丝袜制服| 少妇的丰满在线观看| 手机成人av网站| 长腿黑丝高跟| 又爽又黄无遮挡网站| 国产一区二区三区视频了| 日本a在线网址| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品sss在线观看| 嫁个100分男人电影在线观看| 日日爽夜夜爽网站| 国内揄拍国产精品人妻在线| 好男人在线观看高清免费视频| 成人永久免费在线观看视频| 美女高潮喷水抽搐中文字幕| 亚洲一区高清亚洲精品| 少妇熟女aⅴ在线视频| 亚洲男人天堂网一区| 国产成人av教育| 国产精品一区二区三区四区免费观看 | 成年人黄色毛片网站| 国产精品永久免费网站| 黄色 视频免费看| 国产伦人伦偷精品视频| 日本一本二区三区精品| 国产精品久久久av美女十八| 久久久久九九精品影院| 国产午夜精品论理片| 99热只有精品国产| 欧美黄色淫秽网站| 亚洲中文字幕日韩| 91成年电影在线观看| √禁漫天堂资源中文www| 国产精品久久久久久亚洲av鲁大| 国产精品 国内视频| 欧美最黄视频在线播放免费| 搡老岳熟女国产| 欧美成狂野欧美在线观看| 久久久久国产精品人妻aⅴ院| 久久精品影院6| 午夜成年电影在线免费观看| 欧美日韩黄片免| 亚洲av第一区精品v没综合| 脱女人内裤的视频| 一a级毛片在线观看| 制服诱惑二区| 99国产精品一区二区三区| 久久精品影院6| svipshipincom国产片| 999久久久国产精品视频| 人成视频在线观看免费观看| 日日爽夜夜爽网站| 国产精品久久久人人做人人爽| av片东京热男人的天堂| 身体一侧抽搐| 国产精品 国内视频| 欧美日韩国产亚洲二区| 88av欧美| 亚洲国产中文字幕在线视频| 色在线成人网| 久9热在线精品视频| 国产精品免费视频内射| 琪琪午夜伦伦电影理论片6080| 中文字幕最新亚洲高清| 亚洲乱码一区二区免费版| 国产av又大| 手机成人av网站| 精品国产乱码久久久久久男人| 国产探花在线观看一区二区| 国产一区在线观看成人免费| 欧美日韩福利视频一区二区| 国产精品1区2区在线观看.| 国产亚洲精品久久久久5区| 久久 成人 亚洲| 久久中文看片网| 法律面前人人平等表现在哪些方面| 女人高潮潮喷娇喘18禁视频| 制服诱惑二区| av免费在线观看网站| 一个人免费在线观看电影 | 99国产精品99久久久久| 久久这里只有精品中国| 亚洲成a人片在线一区二区| 老汉色∧v一级毛片| 18美女黄网站色大片免费观看| 成人高潮视频无遮挡免费网站| 中文字幕熟女人妻在线| 50天的宝宝边吃奶边哭怎么回事| 巨乳人妻的诱惑在线观看| 亚洲成a人片在线一区二区| 久久婷婷成人综合色麻豆| 舔av片在线| 亚洲在线自拍视频| 黄色女人牲交| 亚洲第一欧美日韩一区二区三区| 国产爱豆传媒在线观看 | 在线观看免费视频日本深夜| 一夜夜www| 看黄色毛片网站| 九色成人免费人妻av| 午夜久久久久精精品| avwww免费| bbb黄色大片| 精品国产亚洲在线| 久久精品成人免费网站| 成人国产综合亚洲| 国产精品自产拍在线观看55亚洲| 欧美最黄视频在线播放免费| 国产成年人精品一区二区| 久久热在线av| 国产高清激情床上av| 精品久久蜜臀av无| 男人舔女人的私密视频| 久久午夜亚洲精品久久| 欧美日韩瑟瑟在线播放| 成在线人永久免费视频| 亚洲av中文字字幕乱码综合| АⅤ资源中文在线天堂| 91大片在线观看| 国产免费av片在线观看野外av| 婷婷丁香在线五月| 天天躁狠狠躁夜夜躁狠狠躁| 日本撒尿小便嘘嘘汇集6| 999精品在线视频| 两性夫妻黄色片| 一级黄色大片毛片| 久久亚洲真实| 色老头精品视频在线观看| 国内毛片毛片毛片毛片毛片| 精品高清国产在线一区| 久久精品影院6| 50天的宝宝边吃奶边哭怎么回事| 丰满的人妻完整版| 校园春色视频在线观看| 国产精品久久久av美女十八| 精品一区二区三区四区五区乱码| 悠悠久久av| 欧美zozozo另类| 亚洲va日本ⅴa欧美va伊人久久| 精品国产乱码久久久久久男人| 日韩欧美在线乱码| 亚洲av日韩精品久久久久久密| 亚洲精品av麻豆狂野| АⅤ资源中文在线天堂| 国产探花在线观看一区二区| www.999成人在线观看| xxx96com| 免费av毛片视频| 高清在线国产一区| 两个人免费观看高清视频| 日本三级黄在线观看| av欧美777| 久久伊人香网站| 国产高清激情床上av| 欧洲精品卡2卡3卡4卡5卡区| 男人舔女人下体高潮全视频| 欧美另类亚洲清纯唯美| 欧美在线黄色| 亚洲第一欧美日韩一区二区三区| 国产成人系列免费观看| 国产区一区二久久| 99国产极品粉嫩在线观看| 麻豆成人午夜福利视频| 一进一出抽搐gif免费好疼| 好男人在线观看高清免费视频| 一进一出抽搐gif免费好疼| 狂野欧美白嫩少妇大欣赏| 亚洲中文字幕一区二区三区有码在线看 |