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

    基于多源遙感數(shù)據(jù)的人口空間網(wǎng)格化研究
    ——以福州市為例

    2023-08-01 07:22:56余盛華
    科技和產(chǎn)業(yè) 2023年12期

    余盛華

    (澳門城市大學(xué) 商學(xué)院, 澳門 999078)

    人口分布的空間化研究是人地關(guān)系研究的基礎(chǔ),在城市化進(jìn)程和城市空間擴(kuò)張、城市基礎(chǔ)設(shè)施的建設(shè)、城市經(jīng)濟(jì)產(chǎn)業(yè)的擴(kuò)張和集聚、城市生態(tài)環(huán)境的治理及公共服務(wù)設(shè)施的合理配置等方面的意義不容小覷,對于合理規(guī)劃區(qū)域人口、社會(huì)發(fā)展、經(jīng)濟(jì)建設(shè)和生態(tài)環(huán)境保護(hù)具有重要的參考價(jià)值[1]。

    傳統(tǒng)的以行政單元為基礎(chǔ)的人口統(tǒng)計(jì)數(shù)值并未結(jié)合空間數(shù)據(jù)和空間信息,使得人口數(shù)據(jù)無法同空間地理研究相結(jié)合,無法同社會(huì)經(jīng)濟(jì)要素、自然要素相結(jié)合,阻礙了部門間、行業(yè)間數(shù)據(jù)的共享和交流。同時(shí)作為地理信息存儲(chǔ)的自然地理數(shù)據(jù)(如地形、地勢、河流、氣候、植被等)及人文地理要素(如行政區(qū)、面積、土地利用類型等)一般都為柵格數(shù)據(jù),所以對區(qū)域人口分布進(jìn)行人口空間網(wǎng)格化研究的意義斐然,它不僅可以實(shí)現(xiàn)人口統(tǒng)計(jì)數(shù)據(jù)的空間化,還可以促進(jìn)人口數(shù)據(jù)與自然要素?cái)?shù)據(jù)在人地關(guān)系系統(tǒng)中的融合和分析[2]。

    隨著遙感技術(shù)(RS)及地理信息系統(tǒng)(GIS)的迅猛發(fā)展,采用遙感估算法對人口空間網(wǎng)格化的研究日趨成熟與完善[3-6]。目前運(yùn)用較多的是美國國家防務(wù)氣象資源衛(wèi)星(Defense Meteorological Satellite Program,DMSP),在這顆衛(wèi)星上搭載了名為Operational Linescan System(OLS)的傳感器。該傳感器對于燈光的敏感度極高,不僅能監(jiān)測到一般的較亮的光源,甚至能探測到火光、城市小規(guī)模居民地等發(fā)出的低強(qiáng)度燈光,所以它可以用來探究人口夜間的動(dòng)態(tài)分布,這樣就能在較少數(shù)據(jù)源,但又能保證一定準(zhǔn)確度的情況下實(shí)現(xiàn)人口的空間網(wǎng)格化分布,因此它具有其獨(dú)特的實(shí)用性和價(jià)值性,國內(nèi)外的學(xué)者們都對其有較高的期望[7-10]。梁友嘉等[11]采用夜間燈光數(shù)據(jù)輻射值以及LUCC(土地利用/土地覆蓋變化)指數(shù)對甘肅省張掖市甘州區(qū)的人口數(shù)據(jù)進(jìn)行了空間網(wǎng)格分布的構(gòu)建;卓莉等[12]采用夜間燈光總強(qiáng)度信息進(jìn)行了基于單元網(wǎng)格上的中國人口密度分布模擬計(jì)算;馬鈺琪等[13]利用多地理因子與夜間燈光數(shù)據(jù)相結(jié)合的方法對遼寧省人口數(shù)據(jù)進(jìn)行了人口空間網(wǎng)格化的探究。但國外學(xué)者在實(shí)際應(yīng)用中發(fā)現(xiàn)夜間燈光數(shù)據(jù)存在像元溢出現(xiàn)象(overglow)和像元過飽和(saturation)的問題[14-17],為了解決這一問題需要輔助以其他遙感數(shù)據(jù)對夜間燈光遙感數(shù)據(jù)進(jìn)行修正和適當(dāng)?shù)奶幚?。故本研究使用夜間燈光遙感數(shù)據(jù)的同時(shí)結(jié)合歸一化植被覆蓋指數(shù)(NDVI)的遙感影像數(shù)據(jù)來減少夜間燈光數(shù)據(jù)產(chǎn)生的像元溢出現(xiàn)象和像元過飽和的問題[18],并構(gòu)建相應(yīng)模型,通過模型計(jì)算人居指數(shù),最終完成福州市人口數(shù)據(jù)信息在1 km×1 km的柵格網(wǎng)上的時(shí)空分布。實(shí)現(xiàn)人口數(shù)據(jù)空間網(wǎng)格化不僅對促進(jìn)區(qū)域?qū)用嫔先丝凇h(huán)境和資源的優(yōu)化發(fā)展有重要的參考意義,也能為災(zāi)害風(fēng)險(xiǎn)預(yù)防、抗災(zāi)減災(zāi)提供較為精細(xì)的數(shù)據(jù)支持和保障[19-21]。

    1 研究區(qū)概況及基礎(chǔ)數(shù)據(jù)預(yù)處理

    1.1 研究區(qū)概況

    福州市是福建省省會(huì)城市,位于福建省東部沿海地區(qū),臺(tái)灣海峽西側(cè)、閩江下游。福州市下轄5區(qū)2縣級市6縣,包括鼓樓區(qū)、臺(tái)江區(qū)、倉山區(qū)、晉安區(qū)、馬尾區(qū)、長樂市、福清市、閩侯縣、連江縣、羅源縣、閩清縣、永泰縣、平潭縣,總體土地面積約為1.2×104km2,福州市行政區(qū)劃如圖1所示。根據(jù)2021年底福州市統(tǒng)計(jì)報(bào)告,其中實(shí)現(xiàn)GDP 11 300億元,總?cè)丝谶_(dá)到832萬人,城鎮(zhèn)人口603.2萬人,城鎮(zhèn)化率達(dá)72.5%[22]。福州市不僅是福建省政治和文化中心,同時(shí)也是海峽西岸經(jīng)濟(jì)區(qū)中心城市之一。近年來福州市正致力打造為“21世紀(jì)海上絲綢之路”的戰(zhàn)略核心城市。

    審圖號:閩S(2018)37號圖1 福州市行政區(qū)劃圖

    1.2 基礎(chǔ)數(shù)據(jù)預(yù)處理

    基礎(chǔ)數(shù)據(jù)預(yù)處理:①在地理信息系統(tǒng)中,提取基礎(chǔ)地理信息數(shù)據(jù)中的工作地域;②2013年夜間燈光遙感影像DMSP/OLS,獲取于美國國家環(huán)境信息研究中心網(wǎng)站所提供的Average Visible夜間燈光掃描影像,該影像的空間分辨率約為1 km,在地理信息系統(tǒng)中改變其投影方式由原來的地理經(jīng)緯度投影Geographic(LAT/LON)轉(zhuǎn)化為雙標(biāo)準(zhǔn)經(jīng)緯等積圓錐投影Albers以便后續(xù)空間化處理,緊接著裁模、提取、配準(zhǔn)等一系列地理標(biāo)準(zhǔn)化處理以提取福州市夜間燈光數(shù)據(jù);③將歸一化植被覆蓋指數(shù)(NDVI)的遙感影像也由原地理經(jīng)緯度投影Geographic(LAT/LON)轉(zhuǎn)化為雙標(biāo)準(zhǔn)經(jīng)緯等積圓錐投影Albers,這樣有利于不同遙感影像數(shù)據(jù)之間的空間疊加與分析;④在地理信息系統(tǒng)中建立福州市各轄區(qū)、市、縣的數(shù)據(jù)庫,并對數(shù)據(jù)庫進(jìn)行相應(yīng)字段的建立與編碼;⑤將福州市高程DEM數(shù)據(jù)的投影轉(zhuǎn)換為雙標(biāo)準(zhǔn)經(jīng)緯等積圓錐投影Albers,并對投影轉(zhuǎn)化后的數(shù)據(jù)進(jìn)行重采樣處理,最終生成空間分辨率為1 km的柵格數(shù)據(jù)。

    1.3 研究技術(shù)路線

    技術(shù)路線:①制作人口空間網(wǎng)格化區(qū)域的工作底圖;②對DMSP/OLS夜間燈光遙感影像與NDVI(植被覆蓋指數(shù))遙感影像進(jìn)行空間疊加;③基于遙感影像疊加后的關(guān)系構(gòu)建人居指數(shù)模型;④建立人居指數(shù)模型與統(tǒng)計(jì)人口之間的線性回歸方程;⑤觀察人居指數(shù)模型與統(tǒng)計(jì)人口的線性關(guān)系,并進(jìn)行人口模擬與模擬結(jié)果誤差的測算;⑥疊加DEM高程數(shù)據(jù)對人居指數(shù)模型進(jìn)行修正;⑦生成人口數(shù)據(jù)在空間網(wǎng)格上的展布。具體技術(shù)路線如圖2所示。

    2 研究方法

    2.1 確定適宜研究區(qū)與不適宜研究區(qū)

    由于福州市轄區(qū)內(nèi)有較多的山地、丘陵、巖體及水體,這些地區(qū)都不適宜人類居住與人口分布,故采用2013年SPOT 10 d合成的最大的NDVI,得到的NDVI大于0.9的區(qū)域基本都為植被覆蓋較多的森林,而NDVI小于0.1的區(qū)域基本都為植被覆蓋較少的巖體或水體,這些不適宜作為研究的區(qū)域都可以先進(jìn)行剔除以免對后續(xù)的研究產(chǎn)生干擾。

    2.2 使用夜間燈光數(shù)據(jù)和NDVI建立人居指數(shù)

    提取2013年福州市DMSP/OLS夜間燈光數(shù)據(jù)、NDVI遙感影像的空間分布圖,并對相應(yīng)的指數(shù)值進(jìn)行不同深度顏色上的渲染,比便于更好地對結(jié)果進(jìn)行觀察,如圖3、圖4所示。從兩幅圖中不難發(fā)現(xiàn),夜間燈光數(shù)據(jù)的空間分布與植被覆蓋指數(shù)的空間分布有較強(qiáng)的空間上的對應(yīng)關(guān)系,具體表現(xiàn)為夜間燈光數(shù)值較高的地方,植被覆蓋指數(shù)反而較低,反之亦然。其次,對比兩張圖,可以明顯看出NDVI在柵格像元中的數(shù)據(jù)精細(xì)度要高于DMSP/OLS數(shù)據(jù),這主要源于NDVI所使用傳感器的初始分辨率要高于DMSP/OLS數(shù)據(jù)的傳感器,同時(shí)福州市中心城區(qū)夜間燈光數(shù)據(jù)存在亮度過飽和而導(dǎo)致的像元溢出等問題。所以,本文將這兩種遙感影像數(shù)據(jù)進(jìn)行空間疊加分析不僅可以提高像元柵格的精細(xì)度還可以很大程度上削減夜間燈光像元溢出所帶來的問題。為確定遙感數(shù)據(jù)與人口數(shù)據(jù)的擬合關(guān)系,利用夜間燈光數(shù)據(jù)和植被覆蓋指數(shù)構(gòu)建人居指數(shù)模型,借鑒Lu等[23]的研究,利用式(1)構(gòu)建人居指數(shù)(human settlement index,HSI):

    審圖號:閩S(2018)37號圖3 2013年DMSP/OLS穩(wěn)定夜間燈光

    審圖號:閩S(2018)37號圖4 2013年最大NDVI

    式中:NDVImax為SPOTNDVI在2013年的最大值;OLSnor為2013年經(jīng)過標(biāo)準(zhǔn)化的夜間燈光指數(shù)。

    2.3 人口模擬的誤差值檢驗(yàn)

    人口模擬誤差檢驗(yàn)的計(jì)算公式為

    式中:MPE為平均相對誤差;RE為相對誤差;POPm為市(區(qū))縣模擬人口數(shù)量;POPa為市(區(qū))縣統(tǒng)計(jì)人口數(shù)量;n為福州市下轄市(區(qū))縣的數(shù)量。

    3 結(jié)果與分析

    3.1 利用人居指數(shù)對福州市人口進(jìn)行初步模擬

    通過融合的夜間燈光數(shù)據(jù)與植被覆蓋指數(shù)構(gòu)建出人居指數(shù)模型,同時(shí)建立人居指數(shù)與福州市實(shí)際人口數(shù)量之間的線性回歸關(guān)系,如圖5所示,從圖5中不難看出基于DMSP/OLS夜間燈光數(shù)據(jù)和NDVI植被覆蓋指數(shù)所建立的人居指數(shù)與實(shí)際統(tǒng)計(jì)人口數(shù)量之間的線性回歸關(guān)系相關(guān)性較弱(R2=0.352 1)。進(jìn)一步通過線性回歸方程計(jì)算模擬人口數(shù)量,其中長樂區(qū)672 316人、閩侯縣900 661人、福清市890 429人、連江縣551 546人、閩清縣487 763人、平潭縣405 195人、倉山區(qū)559 141人、永泰縣 525 615人、馬尾區(qū)445 383人、羅源縣491 092人、晉安區(qū)489 687人、臺(tái)江區(qū)340 549人、鼓樓區(qū)355 927人。 人口模擬相對誤差百分比絕對值為長樂區(qū)1.51%、閩侯縣36.03%、福清市27.89%、連江縣1.77%、閩清縣105.25%、平潭縣13.26%、倉山區(qū)26.69%、永泰縣110.71%、馬尾區(qū)92.03%、羅源縣136.47%、晉安區(qū)38.21%、臺(tái)江區(qū)23.80%、鼓樓區(qū)48.24%,平均相對誤差為50.91%,如圖6所示。從數(shù)據(jù)中可以看出其誤差百分比極大、導(dǎo)致線性回歸弱相關(guān)的區(qū)域基本都為地勢相對較高的福州市西部和北部地區(qū),如閩清縣、永泰縣、羅源縣等。進(jìn)一步探究造成這種問題的原因,分析判斷是由于此前基于夜間燈光數(shù)據(jù)與植被覆蓋指數(shù)構(gòu)建的人居指數(shù)模型并沒有考慮到海拔因素對人口分布的影響作用,導(dǎo)致人口模擬結(jié)果誤差顯著偏大,因此需進(jìn)一步疊加反映海拔的DEM高程數(shù)據(jù)對人居指數(shù)模型進(jìn)行數(shù)據(jù)修正。

    圖5 福州市各市(區(qū))縣人居指數(shù)累計(jì)值與統(tǒng)計(jì)人口責(zé)任編輯:數(shù)量的回歸關(guān)系

    審圖號:閩S(2018)37號圖6 福州市各市(區(qū))縣人口模擬誤差分布

    3.2 對初步模擬結(jié)果進(jìn)行海拔修正

    首先對海拔對于人口數(shù)量分布的影響進(jìn)行研究,福州市DEM高程影像顯示的海拔分布情況如圖7所示,從圖7中可以看出,福州市地勢西高東低,由西北向東南呈遞減趨勢,地勢最高處約為1 654 m,最低處約為-7 m。其次建立福州市各市(區(qū))縣平均人口密度與平均海拔高度的線性關(guān)系,使用指數(shù)函數(shù)進(jìn)行回歸分析,如圖8所示,得到的線性回歸方程的決定系數(shù)為R2=0.764,表明海拔高度與人口密度之間有比較強(qiáng)的負(fù)相關(guān)性,即可作為影響人口分布模擬的因素之一。

    審圖號:閩S(2018)37號圖7 福州市各市(區(qū))縣海拔分布

    圖8 福州市各市(區(qū))縣平均海拔與平均人口密度的關(guān)系

    通過研究證明海拔高度與人口分布密度有較強(qiáng)的相關(guān)性,所以可用DEM高程數(shù)據(jù)對之前基于DMSP/OLS夜間燈光遙感數(shù)據(jù)疊加NDVI植被覆蓋指數(shù)遙感數(shù)據(jù)后所建立的人居指數(shù)模型進(jìn)行數(shù)據(jù)修正,經(jīng)過DEM高程數(shù)據(jù)訂正后的人居指數(shù)計(jì)算公式為

    HSInew=

    式中:NDVImax為SPOTNDVI在2013年的最大值;OLSnor為2013年經(jīng)過標(biāo)準(zhǔn)化的夜間燈光指數(shù);DEM為高程。

    圖9給出了經(jīng)過DEM高程修正后的人居指數(shù)與統(tǒng)計(jì)人口數(shù)量之間的線性回歸關(guān)系。結(jié)果顯示,將DEM數(shù)據(jù)與人居指數(shù)進(jìn)行融合后所構(gòu)建的線性回歸方程的決定系數(shù)較之前的(R2=0.352 1)有了明顯的提高(R2=0.534 9),但依舊偏低仍不能表明人居指數(shù)與人口統(tǒng)計(jì)數(shù)據(jù)之間的高相關(guān)性。通過DEM高程數(shù)據(jù)修正后的線性回歸方程計(jì)算各市(區(qū))縣的模擬人口數(shù)量,其中長樂區(qū)798 637人、閩侯縣733 921人、福清市1 053 981人、連江縣464 660人、閩清縣279 922人、平潭縣290 932人、倉山區(qū)599 643人、永泰縣308 865人、馬尾區(qū)316 383人、羅源縣310 220人、晉安區(qū)303 438人、臺(tái)江區(qū)169 034人、鼓樓區(qū)197 330人。人口模擬相對誤差百分比絕對值為長樂區(qū)16.99%、閩侯縣10.84%、福清市14.65%、連江縣17.25%、閩清縣17.79%、平潭縣18.68%、倉山區(qū)21.38%、永泰縣23.82%、馬尾區(qū)36.41%、羅源縣49.38%、晉安區(qū)61.71%、臺(tái)江區(qū)62.18%、鼓樓區(qū)71.31%,平均相對誤差為32.49%。從數(shù)據(jù)中可以看出,原來人口模擬誤差值大于100%的地區(qū),經(jīng)過DEM高程修正后都有明顯的減小,總體數(shù)據(jù)也較之前更加趨于合理,如圖10所示。但仍有3個(gè)區(qū)域的模擬誤差值大于50%,分別是鼓樓區(qū)、臺(tái)江區(qū)、晉安區(qū)。這3個(gè)區(qū)域是構(gòu)成福州市中心城區(qū)的主要組成區(qū)域,而作為中心城區(qū)同時(shí)也是夜間燈光極其飽和的地域,故燈光溢出、像元過飽和的現(xiàn)象十分嚴(yán)重,對人口模擬數(shù)值產(chǎn)生了較大的干擾。

    圖9 福州市各市(區(qū))縣經(jīng)過DEM高程修正后的人居指數(shù)責(zé)任編輯:累計(jì)值與統(tǒng)計(jì)人口數(shù)量的回歸關(guān)系

    審圖號:閩S(2018)37號圖10 福州市各市(區(qū))縣人口模擬誤差分布

    3.3 對中心城區(qū)的人居指數(shù)暫時(shí)剔除

    通過再次對福州市下轄的各市(區(qū))縣模擬誤差的分析發(fā)現(xiàn),福州市區(qū)的鼓樓區(qū)、臺(tái)江區(qū)、晉安區(qū)在線性回歸方程中距離回歸直線離散的較遠(yuǎn)、誤差較大,這是由于這3個(gè)區(qū)人口密度較大,夜間燈光強(qiáng)度較大,導(dǎo)致像元異常飽和及夜間燈光嚴(yán)重溢出導(dǎo)致的,因此決定先將這3個(gè)區(qū)域的人居指數(shù)進(jìn)行暫時(shí)的數(shù)據(jù)剔除,而后再進(jìn)行一次人居指數(shù)與人口統(tǒng)計(jì)數(shù)據(jù)之間的線性回歸關(guān)系分析。圖11給出了剔除了鼓樓區(qū)、臺(tái)江區(qū)、晉安區(qū)修正后的人居指數(shù)與人口統(tǒng)計(jì)數(shù)據(jù)之間的線性回歸關(guān)系,這次回歸方程的決定系數(shù)較之前的(R2=0.534 9)有了非常顯著的提高(R2=0.916 9),這表明了人居指數(shù)與人口統(tǒng)計(jì)數(shù)據(jù)之間較強(qiáng)的相關(guān)性。通過修正后的線性回歸方程計(jì)算剩余各市(區(qū))縣的模擬人口數(shù)量,其中長樂區(qū)799 862人、閩侯縣785 280人、福清市1 166 929人、連江縣464 204人、閩清縣243 916人、平潭縣257 044人、倉山區(qū)625 162人、永泰縣278 428人、馬尾區(qū)287 393人、羅源縣280 044人。人口模擬相對誤差百分比絕對值為長樂區(qū)17.17%、閩侯縣18.60%、福清市5.50%、連江縣17.33%、閩清縣2.64%、平潭縣28.15%、倉山區(qū)18.04%、永泰縣11.61%、馬尾區(qū)23.91%、羅源縣34.85%,平均相對誤差為17.78%,可以看出模擬結(jié)果的精確度有了顯著的提高,整體的人口模擬誤差也出現(xiàn)了顯著的降低,如圖12所示。

    圖11 福州市各市(區(qū))縣剔除鼓樓區(qū)、臺(tái)江區(qū)、晉安區(qū)責(zé)任編輯:修正后的人居指數(shù)累計(jì)值與統(tǒng)計(jì)人口數(shù)量的回歸關(guān)系

    審圖號:閩S(2018)37號圖12 福州市各市(區(qū))縣人口模擬誤差分布

    3.4 人口數(shù)據(jù)的空間網(wǎng)格展布

    剔除中心城區(qū)3個(gè)區(qū)域后,人居指數(shù)與人口統(tǒng)計(jì)數(shù)據(jù)之間的線性回歸擬合關(guān)系已經(jīng)達(dá)到了一個(gè)較好的水準(zhǔn),因此可以對人口數(shù)據(jù)進(jìn)行空間網(wǎng)格的展布以實(shí)現(xiàn)人口分布的空間化。根據(jù)DMSP/OLS夜間燈光遙感影像、NDVI植被覆蓋指數(shù)影像、DEM高程影像等多源遙感數(shù)據(jù)所構(gòu)建的人居指數(shù)模型與人口統(tǒng)計(jì)數(shù)據(jù)之間的線性回歸關(guān)系,將人口統(tǒng)計(jì)數(shù)據(jù)展布到細(xì)致的柵格空間中,柵格格網(wǎng)的空間分辨具體為1 km×1 km,但由于模擬人口與統(tǒng)計(jì)人口之間仍然存在一些誤差尤其是此前予以剔除鼓樓區(qū)、臺(tái)江區(qū)、晉安區(qū)這3個(gè)中心城區(qū),故在這里還需要對數(shù)據(jù)進(jìn)行一些修正,采用計(jì)算統(tǒng)計(jì)人口與模擬人口之間的比值的方法,并用該比值對各市(區(qū))縣的模擬人口進(jìn)行修正,最終生成的最終生成1 km×1 km分辨率的人口網(wǎng)格分布密度圖如圖13所示。從圖13中不難看出,人口密度較高的地區(qū)也正是DMSP/OLS夜間燈光數(shù)據(jù)值較高的地區(qū)、同時(shí)也是NDVI植被覆蓋指數(shù)處于中間波段的地區(qū),說明了人口數(shù)據(jù)的空間網(wǎng)格分布與這兩種遙感數(shù)據(jù)的強(qiáng)相關(guān)性。福州市各市(區(qū))縣人口密度的平均值為606人/km2,最大值為25 961人/km2。人口密集區(qū)主要集中在以鼓樓區(qū)、臺(tái)江區(qū)、倉山區(qū)等的中心城區(qū)。從海拔對人口分布的影響來看,有686萬的人口分布在海拔500 m以下的地區(qū),占全市總?cè)丝诘?6%,而海拔在500 m以上地區(qū)的人口數(shù)量只有2.5萬人,占全市總?cè)丝诘?%,所以進(jìn)一步驗(yàn)證了海拔高度對人口的影響是顯著的,添加DEM高程對數(shù)據(jù)進(jìn)行修正也是十分有必要的,通過空間網(wǎng)格展布后的人口分布數(shù)據(jù)圖較此前以行政區(qū)劃為單位的人口分布圖明顯有了質(zhì)的提升。

    審圖號:閩S(2018)37號圖13 福州市各市(區(qū))縣1 km×1 km分辨率的人口密度

    4 結(jié)論與展望

    4.1 結(jié)論

    本文利用多源遙感數(shù)據(jù)實(shí)現(xiàn)人口空間網(wǎng)格化分布,將DMSP/OLS夜間燈光數(shù)據(jù)、NDVI等多源遙感數(shù)據(jù)構(gòu)建人居指數(shù)并與統(tǒng)計(jì)人口數(shù)據(jù)建立線性回歸關(guān)系,利用DEM高程數(shù)據(jù)對數(shù)據(jù)進(jìn)行修正、果斷剔除燈光過飽和的中心城區(qū)數(shù)據(jù),最終在1 km×1 km網(wǎng)格單元上實(shí)現(xiàn)2013年福州市各市(區(qū))縣的人口空間展布模擬。

    探究了空間遙感數(shù)據(jù)與統(tǒng)計(jì)人口數(shù)據(jù)的關(guān)系。采用融合DMSP/OLS夜間燈光數(shù)據(jù)與NDVI植被指數(shù)來構(gòu)建人居指數(shù)模型,并用該人居指數(shù)與統(tǒng)計(jì)人口數(shù)據(jù)進(jìn)行回歸分析建立線性回歸方程。該方法可以較為直觀地顯示遙感數(shù)據(jù)與統(tǒng)計(jì)人口數(shù)據(jù)的線性擬合關(guān)系,也有助于較快發(fā)現(xiàn)數(shù)據(jù)間存在的問題并進(jìn)行及時(shí)的探究和修正。

    探討了如何解決夜間燈光數(shù)據(jù)像元過飽和及溢出的問題。采用空間分辨率更高、數(shù)據(jù)更加精細(xì)的NDVI來對DMSP/OLS夜間燈光數(shù)據(jù)進(jìn)行修正,經(jīng)過修正后的夜間燈光數(shù)據(jù)在像元單元格網(wǎng)上表現(xiàn)得更加精準(zhǔn)和細(xì)膩,這對于人口數(shù)據(jù)的精細(xì)化空間展布起到了至關(guān)重要的作用。

    由于福州市位于東南丘陵地區(qū),地勢起伏明顯,較高海拔地區(qū)并不適宜人類居住。本文采用DEM高程數(shù)據(jù)對人居指數(shù)模型進(jìn)行修正,將不適宜人類居住的較高海拔地區(qū)予以數(shù)據(jù)矯正。該方法可降低只基于遙感數(shù)據(jù)所帶來的模擬誤差,可作為其他丘陵地區(qū)或地勢起伏較大地區(qū)進(jìn)行人口空間化研究時(shí)的參考。

    如何解決城市中心城區(qū)夜間燈光值過飽和問題,是本研究的一個(gè)障礙,同時(shí)也是一個(gè)創(chuàng)新點(diǎn)。本文通過人居指數(shù)與統(tǒng)計(jì)人口數(shù)據(jù)之間的線性回歸關(guān)系發(fā)現(xiàn)中心城區(qū)的數(shù)據(jù)對于模擬結(jié)果的誤差產(chǎn)生較大的影響。所以,本文選擇暫時(shí)擱置中心城區(qū)過飽和的夜間燈光數(shù)據(jù),通過探究表明剔除中心城區(qū)夜間燈光數(shù)據(jù)后,人居指數(shù)與統(tǒng)計(jì)人口數(shù)據(jù)之間建立了很強(qiáng)的線性擬合關(guān)系。因此,基于多源遙感數(shù)據(jù)實(shí)現(xiàn)人口空間網(wǎng)格的精細(xì)化展布是可行的。

    本文基于多源遙感數(shù)據(jù)實(shí)現(xiàn)人口空間網(wǎng)格化分布,這種方法較之前基于行政區(qū)劃來統(tǒng)計(jì)人口。從直觀性上,可更加直觀地看出人口的具體分布位置、密度;從經(jīng)濟(jì)性上,無須再費(fèi)時(shí)費(fèi)力的進(jìn)行以行政區(qū)為單位的人口調(diào)查,只需要利用到逐年的夜間燈光數(shù)據(jù)就能對人口進(jìn)行模擬及預(yù)測,更加省時(shí)省力節(jié)省成本;從實(shí)用性上,以像元單元網(wǎng)格為基礎(chǔ)的人口空間格網(wǎng)化,十分有利于城市建設(shè)總體規(guī)劃、防災(zāi)減災(zāi)的風(fēng)險(xiǎn)評估、社會(huì)可持續(xù)化發(fā)展等。

    4.2 不足之處

    1)并未找到很好的解決方案及模型,對中心城區(qū)過飽和的夜間燈光數(shù)據(jù)進(jìn)行修正使之建立與統(tǒng)計(jì)人口之間較強(qiáng)的線性回歸關(guān)系,該問題有待進(jìn)一步探究解決。

    2)由于受到時(shí)間與經(jīng)費(fèi)的限制,本研究并沒有對空間化后的人口網(wǎng)格數(shù)據(jù)進(jìn)行實(shí)地抽樣調(diào)查、檢驗(yàn)、評估,只是采用了傳統(tǒng)的誤差檢驗(yàn)方法。

    3)由于福州市整體海岸線較為曲折,導(dǎo)致沿海地區(qū)及周邊小島有數(shù)據(jù)丟失的情況,增加了數(shù)據(jù)的不穩(wěn)定性。

    4.3 研究展望

    1)未來可嘗試引入更加貼近于實(shí)際的分析模型用于解決城市中心城區(qū)夜間燈光過飽和及溢出等造成的無法與遙感數(shù)據(jù)建立較強(qiáng)相關(guān)性的問題。

    2)在時(shí)間、經(jīng)費(fèi)許可的情況下,可將傳統(tǒng)誤差檢驗(yàn)法與實(shí)地抽樣調(diào)查法相結(jié)合,使得所得到的模擬數(shù)據(jù)更加精確、合理,更加符合實(shí)際情況。

    3)在未來的研究中可將更多精度更高的遙感數(shù)據(jù)與夜間燈光數(shù)據(jù)相結(jié)合,以進(jìn)一步提高人口空間網(wǎng)格化的精細(xì)化程度。

    久久婷婷人人爽人人干人人爱| 三级毛片av免费| 亚洲美女搞黄在线观看 | 亚洲成人免费电影在线观看| 搡女人真爽免费视频火全软件 | 欧美成人一区二区免费高清观看| 亚洲av综合色区一区| 最近最新中文字幕大全电影3| 国产在线免费精品| 精品久久久久久久久亚洲| 国产亚洲最大av| 亚洲综合精品二区| 国产深夜福利视频在线观看| 伦理电影大哥的女人| 男人舔奶头视频| 1000部很黄的大片| 免费av不卡在线播放| 男女边摸边吃奶| 国产精品国产三级专区第一集| 精品人妻偷拍中文字幕| 制服丝袜香蕉在线| 国产高潮美女av| av播播在线观看一区| 日韩不卡一区二区三区视频在线| 日日撸夜夜添| 男女国产视频网站| 黄色怎么调成土黄色| videossex国产| 在线免费十八禁| 免费看日本二区| 精品亚洲成a人片在线观看 | 亚洲一级一片aⅴ在线观看| 国产精品国产三级国产av玫瑰| 欧美高清成人免费视频www| 欧美日韩亚洲高清精品| 日韩强制内射视频| 高清黄色对白视频在线免费看 | 欧美人与善性xxx| 黄色配什么色好看| 日本黄色日本黄色录像| 亚洲第一区二区三区不卡| 身体一侧抽搐| av国产免费在线观看| 精品久久久久久久久亚洲| 亚洲怡红院男人天堂| av免费在线看不卡| 激情 狠狠 欧美| 青青草视频在线视频观看| 人人妻人人添人人爽欧美一区卜 | 亚洲精品国产av蜜桃| 在线亚洲精品国产二区图片欧美 | 91久久精品国产一区二区三区| 成人影院久久| 好男人视频免费观看在线| 亚洲精品国产av成人精品| 一本一本综合久久| 国产69精品久久久久777片| 久久久久久人妻| 日韩,欧美,国产一区二区三区| 自拍偷自拍亚洲精品老妇| 亚洲精品国产成人久久av| 熟女人妻精品中文字幕| 久久毛片免费看一区二区三区| 91久久精品电影网| 狂野欧美激情性xxxx在线观看| 简卡轻食公司| 国产精品不卡视频一区二区| av播播在线观看一区| 噜噜噜噜噜久久久久久91| 最近2019中文字幕mv第一页| 80岁老熟妇乱子伦牲交| 精品人妻视频免费看| 国产精品女同一区二区软件| 免费看av在线观看网站| 国产黄片美女视频| 一本—道久久a久久精品蜜桃钙片| 久久午夜福利片| 亚洲国产精品国产精品| 永久网站在线| 成人特级av手机在线观看| 国产成人aa在线观看| 午夜免费观看性视频| 男的添女的下面高潮视频| 一级片'在线观看视频| 国产淫片久久久久久久久| 美女中出高潮动态图| 永久免费av网站大全| 大片电影免费在线观看免费| 久久99热6这里只有精品| 亚洲无线观看免费| 精品国产乱码久久久久久小说| 久久精品久久精品一区二区三区| 久久久久久久久久久丰满| 久久女婷五月综合色啪小说| 色吧在线观看| 一本一本综合久久| 色5月婷婷丁香| 91aial.com中文字幕在线观看| 视频中文字幕在线观看| 97在线人人人人妻| 国内少妇人妻偷人精品xxx网站| 麻豆国产97在线/欧美| 欧美xxxx性猛交bbbb| 国产午夜精品一二区理论片| 免费观看a级毛片全部| 亚洲欧美日韩卡通动漫| 中文在线观看免费www的网站| 全区人妻精品视频| 一级毛片电影观看| 国产精品蜜桃在线观看| 午夜福利在线观看免费完整高清在| 一本色道久久久久久精品综合| 国产片特级美女逼逼视频| 秋霞伦理黄片| 3wmmmm亚洲av在线观看| 日韩免费高清中文字幕av| 国产精品麻豆人妻色哟哟久久| 亚洲国产高清在线一区二区三| 免费不卡的大黄色大毛片视频在线观看| 久久韩国三级中文字幕| 亚洲第一av免费看| 久久久午夜欧美精品| 激情 狠狠 欧美| 三级国产精品片| 中国美白少妇内射xxxbb| 亚洲美女黄色视频免费看| 人人妻人人添人人爽欧美一区卜 | 成年美女黄网站色视频大全免费 | 成人毛片a级毛片在线播放| 一级a做视频免费观看| 国产成人aa在线观看| 蜜桃在线观看..| 街头女战士在线观看网站| 欧美日韩视频高清一区二区三区二| 性色av一级| 在线精品无人区一区二区三 | 国产在线一区二区三区精| 三级经典国产精品| 交换朋友夫妻互换小说| 久久久久久伊人网av| 亚洲国产欧美人成| 亚洲va在线va天堂va国产| 少妇 在线观看| 亚洲美女搞黄在线观看| 日韩在线高清观看一区二区三区| 国产免费一级a男人的天堂| 国产精品偷伦视频观看了| 久久久久久伊人网av| 搡老乐熟女国产| 王馨瑶露胸无遮挡在线观看| 国产精品伦人一区二区| 亚洲av成人精品一区久久| 欧美激情国产日韩精品一区| 联通29元200g的流量卡| 免费观看在线日韩| 一级毛片 在线播放| 色视频在线一区二区三区| 欧美成人午夜免费资源| 亚洲av欧美aⅴ国产| 我要看日韩黄色一级片| 最近最新中文字幕大全电影3| 中文天堂在线官网| 日日摸夜夜添夜夜爱| 老女人水多毛片| 国产永久视频网站| 国产成人a∨麻豆精品| 国产日韩欧美在线精品| 欧美bdsm另类| 国产黄色视频一区二区在线观看| h日本视频在线播放| 国产成人精品婷婷| 男人狂女人下面高潮的视频| 日韩欧美 国产精品| 欧美性感艳星| 亚洲第一av免费看| 成人国产麻豆网| 麻豆成人午夜福利视频| 免费av不卡在线播放| 有码 亚洲区| 嘟嘟电影网在线观看| 你懂的网址亚洲精品在线观看| 亚洲欧洲日产国产| 精品人妻一区二区三区麻豆| 国产精品福利在线免费观看| 日韩欧美精品免费久久| 麻豆成人av视频| 99久久精品热视频| 日韩 亚洲 欧美在线| 黄色配什么色好看| 亚洲精品国产色婷婷电影| 高清毛片免费看| 国内少妇人妻偷人精品xxx网站| 久久青草综合色| 国产成人a区在线观看| av在线观看视频网站免费| 日韩欧美精品免费久久| 日本黄色日本黄色录像| 2018国产大陆天天弄谢| av黄色大香蕉| 欧美日韩精品成人综合77777| 亚洲欧美清纯卡通| 欧美日本视频| 大码成人一级视频| 久久久久久人妻| 国产爱豆传媒在线观看| 一级毛片我不卡| 国产精品国产av在线观看| 久久精品久久久久久噜噜老黄| 久久久久网色| 人人妻人人澡人人爽人人夜夜| 免费av中文字幕在线| 中文字幕制服av| 人妻夜夜爽99麻豆av| av在线播放精品| 亚洲三级黄色毛片| 国产成人免费观看mmmm| 美女高潮的动态| 肉色欧美久久久久久久蜜桃| 91精品国产国语对白视频| 久久久a久久爽久久v久久| 中文欧美无线码| 国产精品久久久久成人av| 欧美激情国产日韩精品一区| 亚洲精品乱久久久久久| 欧美成人精品欧美一级黄| 亚洲电影在线观看av| 国产精品精品国产色婷婷| 日日啪夜夜爽| 国产精品无大码| 日韩免费高清中文字幕av| 久久亚洲国产成人精品v| 色婷婷av一区二区三区视频| 美女高潮的动态| 亚洲欧美一区二区三区黑人 | 午夜福利视频精品| 国产又色又爽无遮挡免| 久久久久视频综合| 亚洲不卡免费看| 男人舔奶头视频| 嫩草影院入口| 高清视频免费观看一区二区| 2018国产大陆天天弄谢| videos熟女内射| 国产精品国产三级专区第一集| 婷婷色综合www| 日日摸夜夜添夜夜添av毛片| 高清av免费在线| 欧美极品一区二区三区四区| 亚洲一区二区三区欧美精品| 热99国产精品久久久久久7| 十八禁网站网址无遮挡 | 毛片一级片免费看久久久久| 日日摸夜夜添夜夜添av毛片| 一级a做视频免费观看| 直男gayav资源| 日产精品乱码卡一卡2卡三| 能在线免费看毛片的网站| 国产又色又爽无遮挡免| 啦啦啦在线观看免费高清www| 久久久久久久大尺度免费视频| 久久精品人妻少妇| 女性被躁到高潮视频| 九草在线视频观看| 欧美一区二区亚洲| 亚洲丝袜综合中文字幕| 一本色道久久久久久精品综合| 女的被弄到高潮叫床怎么办| 亚洲精品,欧美精品| 我要看黄色一级片免费的| 亚洲国产日韩一区二区| 欧美性感艳星| 国产精品不卡视频一区二区| 成人特级av手机在线观看| tube8黄色片| 观看av在线不卡| 亚洲精品成人av观看孕妇| 天堂中文最新版在线下载| av在线播放精品| 免费看不卡的av| 18禁裸乳无遮挡动漫免费视频| 国产精品伦人一区二区| 国产午夜精品久久久久久一区二区三区| 久久久色成人| 欧美xxⅹ黑人| 一本久久精品| 国产av精品麻豆| 久久精品久久精品一区二区三区| 欧美人与善性xxx| 97在线视频观看| 国产高清三级在线| 国产精品秋霞免费鲁丝片| 天堂俺去俺来也www色官网| 中文天堂在线官网| 久久久久精品性色| 纵有疾风起免费观看全集完整版| 免费观看的影片在线观看| www.av在线官网国产| 国产精品精品国产色婷婷| 在线观看免费视频网站a站| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久精品古装| 国内少妇人妻偷人精品xxx网站| 国产免费一区二区三区四区乱码| 国产欧美日韩精品一区二区| 亚洲欧美成人综合另类久久久| 欧美丝袜亚洲另类| 亚洲av二区三区四区| 日韩av免费高清视频| 97在线人人人人妻| 国产成人免费观看mmmm| 国产日韩欧美亚洲二区| 亚洲av男天堂| av在线app专区| 亚洲熟女精品中文字幕| 热re99久久精品国产66热6| 亚洲精品亚洲一区二区| 国产中年淑女户外野战色| 久久av网站| 久久久精品94久久精品| 国产伦精品一区二区三区四那| 亚洲一级一片aⅴ在线观看| 欧美高清性xxxxhd video| 91aial.com中文字幕在线观看| 99热6这里只有精品| 男的添女的下面高潮视频| 日本欧美国产在线视频| 精品国产露脸久久av麻豆| 久久 成人 亚洲| 成人高潮视频无遮挡免费网站| 搡老乐熟女国产| av视频免费观看在线观看| 国产成人精品婷婷| 精品一品国产午夜福利视频| 一级a做视频免费观看| 成人综合一区亚洲| 欧美精品国产亚洲| 国产淫语在线视频| 日本猛色少妇xxxxx猛交久久| 少妇精品久久久久久久| 少妇人妻久久综合中文| 又爽又黄a免费视频| 国产精品av视频在线免费观看| 少妇熟女欧美另类| 成人毛片a级毛片在线播放| 日本欧美国产在线视频| 成人毛片60女人毛片免费| 一本久久精品| 国产成人精品久久久久久| 最近手机中文字幕大全| 你懂的网址亚洲精品在线观看| 一级毛片我不卡| 成人亚洲欧美一区二区av| 91久久精品国产一区二区成人| 在线精品无人区一区二区三 | 99精国产麻豆久久婷婷| 久久午夜福利片| 成人亚洲精品一区在线观看 | 伦精品一区二区三区| 99九九线精品视频在线观看视频| 国产成人一区二区在线| 国产黄片视频在线免费观看| 久久鲁丝午夜福利片| 最近中文字幕高清免费大全6| 欧美激情国产日韩精品一区| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧美一区二区三区国产| 内射极品少妇av片p| 少妇人妻 视频| 99久国产av精品国产电影| 精品久久久噜噜| 亚洲精品视频女| 久久97久久精品| 男人添女人高潮全过程视频| 亚洲欧美精品专区久久| 婷婷色综合www| 内射极品少妇av片p| 久久久国产一区二区| 欧美高清性xxxxhd video| 国产精品福利在线免费观看| 亚洲四区av| 岛国毛片在线播放| 黑人猛操日本美女一级片| 99久久人妻综合| 2022亚洲国产成人精品| 中文字幕精品免费在线观看视频 | 在线观看av片永久免费下载| 在线观看免费高清a一片| kizo精华| 久久久成人免费电影| 特大巨黑吊av在线直播| 日本vs欧美在线观看视频 | 亚洲第一区二区三区不卡| 大香蕉97超碰在线| 天天躁夜夜躁狠狠久久av| 亚洲精品第二区| 婷婷色麻豆天堂久久| 国产精品一二三区在线看| 黄片无遮挡物在线观看| 久久久久久久国产电影| 色吧在线观看| 秋霞伦理黄片| 久久ye,这里只有精品| 人妻一区二区av| 黄色怎么调成土黄色| www.av在线官网国产| 天天躁日日操中文字幕| 欧美三级亚洲精品| 人妻少妇偷人精品九色| 午夜福利在线在线| 亚洲av不卡在线观看| 亚洲美女视频黄频| 成人美女网站在线观看视频| 91久久精品国产一区二区三区| 久久av网站| 国精品久久久久久国模美| .国产精品久久| 最近手机中文字幕大全| 精品人妻偷拍中文字幕| 妹子高潮喷水视频| 男女下面进入的视频免费午夜| 直男gayav资源| 爱豆传媒免费全集在线观看| 秋霞伦理黄片| 观看美女的网站| 亚洲国产日韩一区二区| 一本久久精品| 亚洲色图av天堂| 中国国产av一级| 毛片女人毛片| 亚洲国产最新在线播放| 日本欧美视频一区| 国产乱人视频| 久久99热这里只频精品6学生| 欧美一级a爱片免费观看看| 少妇被粗大猛烈的视频| 国产黄色视频一区二区在线观看| 在线天堂最新版资源| 最后的刺客免费高清国语| 亚洲色图av天堂| 成人特级av手机在线观看| 亚洲av综合色区一区| 两个人的视频大全免费| 亚洲精品久久午夜乱码| 国产淫语在线视频| 美女国产视频在线观看| 婷婷色av中文字幕| 男女边摸边吃奶| 丝袜脚勾引网站| 黑丝袜美女国产一区| 最近中文字幕高清免费大全6| 久久久欧美国产精品| 精品午夜福利在线看| 高清在线视频一区二区三区| av国产免费在线观看| 欧美丝袜亚洲另类| 在线看a的网站| 久久午夜福利片| 少妇高潮的动态图| 18禁在线播放成人免费| 小蜜桃在线观看免费完整版高清| 亚洲精品色激情综合| 国产一级毛片在线| 国产一区二区在线观看日韩| 亚洲av在线观看美女高潮| 中文精品一卡2卡3卡4更新| 欧美变态另类bdsm刘玥| 久久精品久久精品一区二区三区| 高清午夜精品一区二区三区| 亚洲精品日韩在线中文字幕| 国产黄色视频一区二区在线观看| 纵有疾风起免费观看全集完整版| 男人爽女人下面视频在线观看| 一区二区三区精品91| 内地一区二区视频在线| 少妇的逼好多水| 国产精品99久久久久久久久| 一区二区三区乱码不卡18| 国产爱豆传媒在线观看| 国产精品蜜桃在线观看| 国产一区亚洲一区在线观看| 亚洲性久久影院| 王馨瑶露胸无遮挡在线观看| 一区二区三区乱码不卡18| 人体艺术视频欧美日本| 久久久久国产精品人妻一区二区| 26uuu在线亚洲综合色| 亚洲国产精品国产精品| 丰满乱子伦码专区| 性高湖久久久久久久久免费观看| 精品一区二区三区视频在线| 亚洲天堂av无毛| 国产日韩欧美亚洲二区| 国产亚洲最大av| 精品酒店卫生间| 国产精品嫩草影院av在线观看| 亚洲精品,欧美精品| 精品亚洲乱码少妇综合久久| 精品久久国产蜜桃| 岛国毛片在线播放| 亚洲国产色片| 国产一区二区在线观看日韩| 大又大粗又爽又黄少妇毛片口| 少妇裸体淫交视频免费看高清| 精品国产露脸久久av麻豆| 少妇 在线观看| 精品国产三级普通话版| 免费人成在线观看视频色| 色视频在线一区二区三区| 全区人妻精品视频| 久久精品久久久久久久性| 日日撸夜夜添| 欧美老熟妇乱子伦牲交| 性色avwww在线观看| 日日啪夜夜撸| 精品亚洲乱码少妇综合久久| 丰满乱子伦码专区| av福利片在线观看| 一级黄片播放器| 日本wwww免费看| 久久婷婷青草| 一区二区三区四区激情视频| 国产精品一区二区性色av| 亚洲精品456在线播放app| av又黄又爽大尺度在线免费看| 伦精品一区二区三区| 我要看黄色一级片免费的| 街头女战士在线观看网站| 丝袜喷水一区| 日韩成人av中文字幕在线观看| 亚洲电影在线观看av| 久久精品久久久久久久性| 十分钟在线观看高清视频www | a级一级毛片免费在线观看| 啦啦啦视频在线资源免费观看| 欧美xxⅹ黑人| 欧美 日韩 精品 国产| 亚洲精品色激情综合| 内射极品少妇av片p| 亚洲av免费高清在线观看| 黄片无遮挡物在线观看| 国产成人精品一,二区| 人人妻人人添人人爽欧美一区卜 | 好男人视频免费观看在线| 亚洲精品视频女| 久久精品国产鲁丝片午夜精品| 狂野欧美激情性xxxx在线观看| 国产av精品麻豆| 天堂俺去俺来也www色官网| 一本—道久久a久久精品蜜桃钙片| 人妻一区二区av| 大香蕉97超碰在线| 观看av在线不卡| 免费黄频网站在线观看国产| 一边亲一边摸免费视频| 另类亚洲欧美激情| 久久ye,这里只有精品| 内射极品少妇av片p| 免费观看性生交大片5| 伦精品一区二区三区| 日日摸夜夜添夜夜爱| 你懂的网址亚洲精品在线观看| 夜夜骑夜夜射夜夜干| 午夜福利网站1000一区二区三区| 国产精品久久久久成人av| 亚洲精品aⅴ在线观看| 欧美xxxx性猛交bbbb| 国产精品国产三级专区第一集| 91精品国产九色| 成人国产麻豆网| 在线观看人妻少妇| 日韩大片免费观看网站| 美女中出高潮动态图| 精品一区在线观看国产| 麻豆成人av视频| 久久 成人 亚洲| 亚洲电影在线观看av| 一级黄片播放器| 最黄视频免费看| 精品亚洲成国产av| 欧美精品国产亚洲| 91精品一卡2卡3卡4卡| 亚洲,一卡二卡三卡| 精品国产一区二区三区久久久樱花 | 又粗又硬又长又爽又黄的视频| 国产午夜精品一二区理论片| 人人妻人人添人人爽欧美一区卜 | 久久久久视频综合| 性高湖久久久久久久久免费观看| 久久6这里有精品| 亚洲电影在线观看av| 亚洲欧美一区二区三区国产| 精品久久久久久久久亚洲| 亚洲av在线观看美女高潮| 噜噜噜噜噜久久久久久91| 亚洲第一区二区三区不卡| 午夜精品国产一区二区电影| 97在线视频观看| 在线天堂最新版资源| 女的被弄到高潮叫床怎么办| 精品一区在线观看国产| 免费观看av网站的网址| 中文字幕免费在线视频6| 中国国产av一级| 高清黄色对白视频在线免费看 | 国产在线一区二区三区精| 日本猛色少妇xxxxx猛交久久| 免费观看性生交大片5| 黄色配什么色好看| 日韩欧美精品免费久久| 1000部很黄的大片| 国产亚洲一区二区精品| 亚洲国产av新网站| 日本欧美视频一区| 国产亚洲午夜精品一区二区久久| 欧美激情国产日韩精品一区|