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

    武漢市1991~2013年城市擴(kuò)張與地表覆蓋變化研究

    2016-12-28 08:41:49孫安昌張敘葭牛瑞卿
    地理空間信息 2016年12期
    關(guān)鍵詞:用地光譜水體

    孫安昌,陳 濤,張敘葭,牛瑞卿

    (1.中國(guó)地質(zhì)大學(xué)(武漢)地球物理與空間信息學(xué)院,430074)

    武漢市1991~2013年城市擴(kuò)張與地表覆蓋變化研究

    孫安昌1,陳 濤1,張敘葭1,牛瑞卿1

    (1.中國(guó)地質(zhì)大學(xué)(武漢)地球物理與空間信息學(xué)院,430074)

    利用1991~2013年間6個(gè)時(shí)相的武漢市Landsat遙感數(shù)據(jù),使用多端元光譜混合分析方法進(jìn)行混合像元分解,得到各時(shí)相的地表覆蓋豐度圖。根據(jù)豐度信息,計(jì)算土地變化動(dòng)態(tài)度與土地覆蓋變化強(qiáng)度,得到地表覆蓋變化情況,分析武漢地區(qū)城市擴(kuò)張與土地轉(zhuǎn)移特征。結(jié)果表明,武漢市主城區(qū)建設(shè)用地源自對(duì)城市外圍的綠地與內(nèi)部湖泊地區(qū)的侵占,湖泊與綠地面積總體呈下降趨勢(shì)。

    城市擴(kuò)張;地表覆蓋變化;多端元光譜混合分析;武漢市;遙感

    城市化過程將自然景觀轉(zhuǎn)化為建設(shè)用地,導(dǎo)致城市地表覆蓋物理特性變化迅速,造成城市內(nèi)部環(huán)境發(fā)生改變[1]。例如,建設(shè)用地的快速增長(zhǎng)會(huì)導(dǎo)致城市熱島效應(yīng)加劇[2-3]。衛(wèi)星遙感技術(shù)的發(fā)展為周期性大范圍城市變化監(jiān)測(cè)提供了支持[4]。對(duì)于中、低空間分辨率的衛(wèi)星影像,因城市地表復(fù)雜,場(chǎng)景中的部分地物小于影像像元的尺寸,不能被直接檢測(cè)[5],傳感器獲取的不是單一地物的地表反射率,而是不同地物按照一定比例加權(quán)的總和,形成混合像元[6]。對(duì)于混合像元,使用傳統(tǒng)的圖像硬分類,直接將像元判為某一地物,必然會(huì)影響地物識(shí)別與圖像分類精度。

    光譜混合分析(spectral mixture analysis,SMA)可將像元分解為不同的端元組分,求得其所占像元的百分含量,解決了傳統(tǒng)分類法對(duì)混合像元分析精度下降以及信息丟失的問題,得到更加準(zhǔn)確的結(jié)果[4],可用于解決混合像元在土地分類中遇到的問題。SMALL采用線性光譜混合模型研究城市植被空間與時(shí)間變化[7-8];LU等人運(yùn)用SMA技術(shù)提取城市不透水面,分類精度達(dá)到83.78%[9]。ROBERTS等人根據(jù)SMA提出了多端元光譜分解(multiple endmember spectral mixture analysis, MESMA),該方法通過調(diào)整端元數(shù)量,評(píng)估和選擇最優(yōu)端元組合,以動(dòng)態(tài)方式處理混合像元,被證明是一種更合理、精度更高的方法[6,10-13]。ROBERTS等人運(yùn)用MESMA分別對(duì)美國(guó)Los Angeles、巴西Manaus的TM圖像進(jìn)行城市土地覆蓋像元分解并測(cè)定城市形態(tài)物質(zhì)組成,結(jié)果表明MESMA在亞像元水平定量獲取地表信息是有效的[10,13]。

    MESMA在分析城市地表成分上應(yīng)用廣泛,但少有使用該方法對(duì)基于時(shí)間序列的遙感影像進(jìn)行定性定量的地表覆蓋變化研究。本文以武漢市為例,使用6個(gè)時(shí)相的Landsat TM遙感影像,運(yùn)用MESMA進(jìn)行混合像元分解,通過分解結(jié)果獲取地表信息,以此驗(yàn)證該方法在多時(shí)相城市地表定量變化監(jiān)測(cè)中的應(yīng)用;同時(shí)分析總結(jié)武漢市近20 a來城市擴(kuò)張的特點(diǎn),為城市發(fā)展與規(guī)劃提供參考。

    1 研究區(qū)概況

    武漢市位于江漢平原東部,長(zhǎng)江和漢水交匯處(113°41'E~115°05'E,29°50'N~31°22'N),市區(qū)總面積為8 494.41 km2,是我國(guó)中部地區(qū)第一大城市,也是重要的金融、商業(yè)、貿(mào)易、文化中心。武漢屬于亞熱帶季風(fēng)氣候,雨量充沛、四季分明;在地形上北高南低,屬于殘丘性河湖沖積平原,北部為山地丘陵,其余均屬江漢平原,地勢(shì)平坦、河道縱橫、湖泊星羅棋布。

    2 數(shù)據(jù)獲取與預(yù)處理

    2.1 數(shù)據(jù)準(zhǔn)備

    本文使用的遙感數(shù)據(jù)為L(zhǎng)andsat影像,共6景(表1);另外還有全國(guó)行政省、市(縣)邊界圖及其他相關(guān)資料,所用圖像處理軟件為Envi5.0。

    表1 遙感數(shù)據(jù)介紹

    2.2 數(shù)據(jù)預(yù)處理

    本文使用的Landsat數(shù)據(jù)為L(zhǎng)1T級(jí),已經(jīng)過輻射校正與幾何校正。在6景影像內(nèi)均勻選取同名像點(diǎn)進(jìn)行配準(zhǔn),誤差控制在0.5個(gè)像素內(nèi),并用矢量數(shù)據(jù)進(jìn)行裁剪。考慮到數(shù)據(jù)限制以及研究需求,僅選取武漢的主城區(qū)作為研究區(qū)域,即江岸區(qū)、江漢區(qū)、硚口區(qū)、漢陽(yáng)區(qū)、武昌區(qū)、青山區(qū)、洪山區(qū)。

    圖1 研究區(qū)遙感影像縮略圖(Landsat8 R、Band4 G、Band3 B、Band2)

    3 算法及模型

    3.1 線性光譜混合模型

    光譜混合分析是將混合像元表現(xiàn)為地表純凈端元光譜組合形式[4]。其常用線性光譜模型定義為:像元在某一光譜波段的反射率是由構(gòu)成像元基本組分的反射率及其所占像元面積比例為權(quán)重系數(shù)的線性組合??捎靡韵鹿奖磉_(dá):

    其中,Riλ為第λ波段第i像元的反射率(已知);N為端元數(shù)目;fi為對(duì)應(yīng)第i像元某一端元的豐度;ρiλ為各端元的反射率;ελ為殘余誤差值。

    評(píng)價(jià)模型用殘差ελ或均方根誤差RMSE表示:

    其中,M是傳感器波段數(shù)。

    線性模型從混合像元中分離和提取各組分的平均光譜響應(yīng),通過求解線性方程來反解端元在像元中所占的比例(豐度),從而將所有像元分解成這些基本組分的分量。模型計(jì)算的結(jié)果表現(xiàn)為各端元的分量值和以均方根誤差表示的殘余誤差圖像。

    3.2 端元提取

    端元必須準(zhǔn)確全面地代表地面材料(地表覆蓋類型)[15]。建立應(yīng)用于混合像元分解的端元光譜庫(kù)的挑戰(zhàn)是雙面的。首先,每類地物的光譜庫(kù)應(yīng)該包含足夠數(shù)量的光譜,這樣才能充分代表地面材料的光譜變化。然而,隨著端元總數(shù)的增加,又會(huì)導(dǎo)致計(jì)算效率下降[14]。為達(dá)到以上目的,更好地選擇出能夠代表土地覆蓋類別的端元,本文參考RIDD提出的概念[17],結(jié)合研究區(qū)域的實(shí)際情況,定義端元類別為水體、植被、建設(shè)用地、裸地。具體說明如下:

    1)水體包括江河湖泊、水田、漁場(chǎng)。

    2)植被包括林地、草場(chǎng)、農(nóng)田。

    3)建設(shè)用地包括所有建成的人工建筑、道路等。4)裸地包括土壤、未建成的建設(shè)用地。

    3.3 最優(yōu)端元光譜篩選

    將端元集制作成光譜庫(kù),并通過參數(shù)EAR(endmember average RMSE)從庫(kù)中篩選出最優(yōu)端元光譜。EAR由DENNISON等在2003年提出[15],該參數(shù)能夠評(píng)判出地表覆蓋類型中最具有代表性的端元。計(jì)算公式如下:

    式中,A代表端元類別;n為類別A的光譜數(shù)量;Ai為具體的光譜。

    EAR是用來評(píng)價(jià)每個(gè)端元模擬分解其地表類型能力的參數(shù),用于衡量端元在該類別分解過程中的表現(xiàn)能力,不能用來衡量端元的純凈程度。EAR值較低的端元比EAR值較高的端元在分解過程中有更好的表現(xiàn)能力,即EAR值較低的端元是最能代表這一地物類別的端元。因此,篩選的最優(yōu)端元應(yīng)具有較低的EAR值??紤]到數(shù)據(jù)空間分辨率和光譜分辨率的限制,以及研究中對(duì)研究區(qū)域影像的判讀和了解,不同時(shí)期數(shù)據(jù)中同一地表覆蓋類別所篩選的光譜數(shù)不相同。各端元光譜數(shù)量見表2。

    表2 不同時(shí)相地表覆蓋模型數(shù)與像元成功分解比例

    3.4 混合像元分解

    本文根據(jù)線性光譜混合模型,采用雙端元模型,將像元分解為2個(gè)組分。增加一個(gè)端元個(gè)數(shù)雖然可以減少RMSE,但對(duì)結(jié)果的提升不大[15]。過多端元會(huì)造成端元之間光譜的混淆,隨著端元個(gè)數(shù)增加,模型數(shù)量也隨之增加,計(jì)算效率就會(huì)下降。雙端元即代表實(shí)際地物的端元光譜與考慮到傳感器照明的變化而加上的陰影(shade,SHD)端元。本次研究中,設(shè)置分解的RMSE小于 0.025,若RMSE大于該閾值,則劃為未識(shí)別區(qū)域。

    3.5 地表信息參數(shù)

    引入兩個(gè)地表信息參數(shù):土地變化動(dòng)態(tài)度和土地覆蓋變化強(qiáng)度。

    1)土地變化動(dòng)態(tài)度。土地變化動(dòng)態(tài)度公式如下:

    式中,AV表示某土地利用類型的變化率;Aε和As分別表示前后兩個(gè)時(shí)相的面積,該面積為像元的豐度乘以該像元所代表的面積,T表示兩個(gè)時(shí)相的時(shí)間差(單位:年)。該系數(shù)能定量描述土地利用變化狀況,土地變化動(dòng)態(tài)度越大說明該地表覆蓋類別在相應(yīng)時(shí)間段內(nèi)的變化越劇烈。

    2)土地覆蓋變化強(qiáng)度。參考MICHISHITA等人[16]的做法,引入土地覆蓋變化強(qiáng)度(LCCI)。土地覆蓋變化以像素為單位,反映了兩景影像時(shí)間內(nèi)每個(gè)土地類別每天的變化情況。公式如下:

    式中,LCFt1和LCFt2為兩個(gè)時(shí)相的LCF圖;?t為兩個(gè)時(shí)相間隔的天數(shù);A為像元所代表的面積,即900 m2,LCCI的單位為m2/day。

    4 結(jié)果與討論

    4.1 精度評(píng)價(jià)

    分解成功像元所占影像的比例如表2所示。1991~2009年期間5個(gè)時(shí)相的數(shù)據(jù)分解成功的像元所占比例大于99%,2013年成功率為94.14%。圖2為1991年分解結(jié)果。

    由于缺乏地表覆蓋豐度數(shù)據(jù),無法從豐度上給出精度評(píng)價(jià)。因此,本文將每個(gè)地表覆蓋類別在影像上隨機(jī)選擇500個(gè)像元,計(jì)算得到1991~2013年混淆矩陣。結(jié)果表明,總體精度最高達(dá)到90%,最低為75%。這4個(gè)地表覆蓋類型中,裸地的混淆情況比較嚴(yán)重,在該光譜分辨率下裸地與某些建設(shè)用地的光譜相近。但總體而言,結(jié)果滿足進(jìn)一步研究需求。

    圖2 分解結(jié)果示例(1991)

    4.2 地表覆蓋變化分析

    根據(jù)土地變化動(dòng)態(tài)度的計(jì)算結(jié)果(表3),從1991~2013年,研究區(qū)域內(nèi)建設(shè)用地面積增長(zhǎng)迅速,年均增長(zhǎng)率為4.36%,植被、水體以及裸地呈負(fù)增長(zhǎng),變化率分別為-2.33%、-1.1%、-3.45%。

    建設(shè)用地在1996~2000年變化最快,達(dá)到6.82%,隨后增長(zhǎng)速度放緩,2005~2009年在研究區(qū)域出現(xiàn)了負(fù)增長(zhǎng),減少了1.69%。建設(shè)用地的整體變化率是所有土地類型中最大的,其快速增長(zhǎng)也反映了城市20多年來的迅猛發(fā)展。

    在綠地面積減少的大趨勢(shì)下,2000~2005年變化率減小,僅為-0.74%;在2005~2009年,則有小幅度的上升。對(duì)比綠地與建設(shè)用地,變化相反且強(qiáng)度相對(duì)一致。

    水體面積也呈下降趨勢(shì),但在1996~2000年有明顯的恢復(fù)。經(jīng)查閱資料,武漢市為了經(jīng)濟(jì)發(fā)展,采取了“退耕還漁”的相關(guān)策略,養(yǎng)殖水面增加,導(dǎo)致湖泊也有較大增長(zhǎng),這與曹麗琴[2]的研究結(jié)果一致。至2009年,水體面積與1996年基本持平,2009~2013年則降到更低。

    裸地變化雜亂,其整體呈減少趨勢(shì),但在1991~1996年和2005~2009年間又有明顯增長(zhǎng)。

    表3 1991~2013年武漢主城區(qū)土地變化動(dòng)態(tài)度/%

    4.3 城區(qū)擴(kuò)張與土地轉(zhuǎn)移分析

    為了探究城市擴(kuò)張以及土地轉(zhuǎn)移的情況,本文提取研究區(qū)域2013年建設(shè)用地部分,并對(duì)其他影像進(jìn)行掩膜處理,制成時(shí)間序列的土地覆蓋豐度變化圖(圖3)。其中,紅色表示建設(shè)用地,綠色、藍(lán)色和灰色分別表示綠地、水體和裸地在2013年前轉(zhuǎn)換為建設(shè)用地的部分;地表覆蓋的豐度越高,對(duì)應(yīng)的顏色越深。

    圖3 武漢主城區(qū)土地覆蓋豐度變化圖

    該時(shí)間序列圖主要闡述了武漢地區(qū)建設(shè)用地的增長(zhǎng)模式以及土地來源。從空間上看,建設(shè)用地的擴(kuò)張是以武漢的地理中心逐層向外,以圈層式的結(jié)構(gòu)發(fā)展。漢口、漢陽(yáng)和武昌的擴(kuò)張?zhí)攸c(diǎn)各不相同(圖3a)。漢口地區(qū)依托老城區(qū),呈弧線向西北地區(qū)逐層發(fā)展;漢陽(yáng)地區(qū)的變化緩慢,主要沿著長(zhǎng)江與漢江方向擴(kuò)張;武昌地區(qū)變化較大,在2009年以前,一方面向東擴(kuò)張,另一方面青山工業(yè)地區(qū)向南和武昌建成的老城區(qū)相連,2009~2013年向外圍郊區(qū)擴(kuò)張明顯。

    建設(shè)用地?cái)U(kuò)張侵占了植被以及水體面積。對(duì)于圖中藍(lán)色被侵占的水體,主要為城市內(nèi)部湖泊,且存在于2005年以前,尤其是2000~2005年(圖3c、3d),大量藍(lán)色區(qū)域消失,填湖情況嚴(yán)重;但在2005年以后,建設(shè)用地對(duì)水體的侵占不明顯。

    4.4 城市擴(kuò)張模式分析

    根據(jù)土地覆蓋變化強(qiáng)度,將建設(shè)用地的發(fā)展分為兩種形式:一種是新建設(shè)用地的擴(kuò)張;另一種是城區(qū)內(nèi)的改造或再發(fā)展。新建設(shè)用地的開發(fā)是其他土地變化為建設(shè)用地的覆蓋類型,而再發(fā)展是將已有的建設(shè)用地進(jìn)行改造,造成其豐度變化。

    由表4可以看出,建設(shè)用地變化最快的時(shí)間段為1996~2000年,植被、水體、裸地快速減少,其中植被和水體的變化強(qiáng)度在歷年中最大;1991~1996年和2009~2013年建設(shè)用地增長(zhǎng)速度適中,相應(yīng)的其他土地覆蓋變化強(qiáng)度為負(fù),這3個(gè)時(shí)間段主要為新城區(qū)的擴(kuò)張。2000~2005年建設(shè)用地的變化強(qiáng)度僅次于1996~2000年,但植被和水體變化強(qiáng)度弱,即該段時(shí)間武漢城市的變化主要為城區(qū)內(nèi)的改造和再發(fā)展。在2005~2009年,建設(shè)用地的面積出現(xiàn)了負(fù)增長(zhǎng),本文認(rèn)為這也與城市內(nèi)部的改造與再發(fā)展有關(guān)系。

    表4 土地覆蓋變化強(qiáng)度表/m2/day

    5 結(jié) 語

    本文使用6個(gè)時(shí)相的武漢Landsat遙感影像為數(shù)據(jù)源,使用多端元光譜混合分析方法進(jìn)行雙端元模型的混合像元分解,得到結(jié)論如下:

    1)通過目視解譯直接選取端元,使用EAR參數(shù)進(jìn)行端元光譜篩選,并采用雙端元模型進(jìn)行多端元光譜混合分析,成功提取了地表覆蓋類別與豐度信息,精度能夠滿足研究需求。根據(jù)結(jié)果計(jì)算土地利用動(dòng)態(tài)度與LCCI,繪制時(shí)間序列LCF圖與LCCI圖,將土地覆蓋類別、土地覆蓋豐度與空間信息相結(jié)合,較好地監(jiān)測(cè)并分析了1991~2013年武漢城市擴(kuò)張情況與特點(diǎn)。

    2)1991~2013年這23 a中,武漢市主城區(qū)地表變化劇烈,建設(shè)用地城市區(qū)域擴(kuò)張明顯,植被和水體總體減少,裸地則隨年份不同變化較大。其中,建設(shè)用地和植被區(qū)域在2005~2009年分別減少和增加,水體在1996~2000年間也有增長(zhǎng),這3個(gè)變化與整體趨勢(shì)相背。

    3)武漢的城市擴(kuò)張以長(zhǎng)江與漢江的交匯處為中心,逐層向外以圈層式的結(jié)構(gòu)發(fā)展。在1991~2000年和2009~2013年主要為城市向外擴(kuò)張,在2000~2009年更多表現(xiàn)為內(nèi)部重建。漢口、漢陽(yáng)和武昌三鎮(zhèn)擴(kuò)張?zhí)攸c(diǎn)不同,漢口依托建成的老城區(qū),呈弧線向西北地區(qū)逐層發(fā)展;漢陽(yáng)地區(qū)的變化緩慢,主要沿著長(zhǎng)江與漢江擴(kuò)張;武昌地區(qū)和青山地區(qū)逐漸發(fā)展相連,隨后沿東軸呈跳躍式增長(zhǎng)。參考文獻(xiàn)

    [1] HUANG Wei, ZENG Yongnian,LI Songnian. An Ananlysis of Urban Expansion and Its Associated Thermal Characteristics Using Landsat Imagery [J]. Geocarto International, 2015(30):93-103

    [2] 曹麗琴.城市熱環(huán)境及其影像因子的定量遙感研究[D].武漢:武漢大學(xué),2009

    [3] 陳濤,孫安昌,王鑫,等.基于遙感的武漢地區(qū)LUCC對(duì)熱島效應(yīng)的影響研究[J].地理空間信息,2015,13(1):62-65

    [4] 趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京:科學(xué)出版社,2003

    [5] ALAN H S, CURTIS E W, JAMES A S. On the Nature of Models in Remote Sensing [J].Remote Sensing of Environment, 1986, 20(2): 121-139

    [6] TAREK R, JOHN R W, DAR R, et al. Measuring the Physical Composition of Urban. Morphology Using Multiple Endmember Spectral Mixture Models[J].Photogrammetric Engineering and Remote Sensing, 2003(69):1 011-1 020

    [7] SMALL C. Multitemporal Analysis of Urban Reflectance [J].Remote Sensing of Environment, 2002(81):427-442

    [8] SMALL C. Estimation of Urban Vegetation Abundance by Spectral Mixture Analysis [J]. International Journal of Remote Sensing, 2001, 22(7):1 305-1 334

    [9] LU Dengsheng, WENG Qihao. Urban Classification Using Full Spectral Information of Landsat ETM+Imagery in Marion County, Indiana [J]. Photogrammetric Engineering & Remote Sensing, 2005,71(11):1 275-1 284

    [10] POWDLL R L,ROBERTS D A,DENNISON P E, et al. Subpixel Mapping of Urban Land Cover Using Multiple Endmember Spectral Mixture Analysis: Manaus, Brazil [J]. Remote Sensing of Environment, 2007(106):253-267

    [11] WU C. Normalized Spectral Mixture Analysis for Monitoring Urban Composition Using ETM+Imagery [J] .Remote Sensing of Environment, 2004(93): 480-492

    [12] FRANKE J, ROBERTS D A, HALLIGAN K, et al.Hierarchical Multiple Endmembers Spectral Mixture Analysis (MESMA) of Hyperspectral Imagery for Urban Environments [J].Remote Sensing of Environment,2009(113):1 712-1 723

    [13] ROBERTS D A, GARDNER M, CHURCH R, et al.Mapping Chaparral in the Santa Monica Mountains Using Multiple Endmember Spectral Mixture Models [J]. Remote Sensing of Environment, 1998, 65(3): 267-279

    [14] ROBERTS D A, NUMATA I, HOLMES K, et al. Large Area Mapping of Land-cover Change in Rondonia Using Multitemporal Spectral Mixture Analysis and Decision Tree Classifiers [J]. Journal of Geophysical Research Atmospheres, 2002, 107(D20): LBA 40-1-LBA 40-18

    [15] DENNISON P E, ROBERTS D A. Endmember Selection for Multiple Endmember Spectral Mixture Analysis Using Endmember Average RMSE [J]. Remote Sensing of Environment, 2003(87):123-135

    [16] RYO M, JIANG Zhiben, XU Bing. Monitoring Two Decades of Urbanization in the Poyang Lake Area, China Through Spectral Unmixing [J].Remote Sensing of Environment, 2012(117):3-18

    [17] 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].International Journal of Remote Sensing, 1995(16): 2 165-2 185

    P273

    B

    1672-4623(2016)12-0033-04

    10.3969/j.issn.1672-4623.2016.12.012

    孫安昌,碩士研究生,研究方向?yàn)榄h(huán)境遙感、遙感地質(zhì)、高光譜遙感。

    2015-08-10。

    項(xiàng)目來源:國(guó)家高技術(shù)研究發(fā)展計(jì)劃資助項(xiàng)目(2012AA121303);數(shù)字制圖與國(guó)土信息應(yīng)用工程國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室開放研究基金資助項(xiàng)目(GCWD201202);湖北省自然科學(xué)基金資助項(xiàng)目(2012FFB06501);長(zhǎng)江科學(xué)院開放研究基金資助項(xiàng)目 (CKWV2013221/KY)。

    猜你喜歡
    用地光譜水體
    基于三維Saab變換的高光譜圖像壓縮方法
    農(nóng)村黑臭水體治理和污水處理淺探
    生態(tài)修復(fù)理念在河道水體治理中的應(yīng)用
    城鄉(xiāng)建設(shè)用地增減掛鉤研究進(jìn)展綜述
    廣元:治理黑臭水體 再還水清岸美
    城鄉(xiāng)建設(shè)用地增減掛鉤政策的演變
    城鄉(xiāng)建設(shè)用地增減掛鉤的實(shí)踐與認(rèn)識(shí)
    城鄉(xiāng)建設(shè)用地增減掛鉤的實(shí)踐與認(rèn)識(shí)
    星載近紅外高光譜CO2遙感進(jìn)展
    苦味酸與牛血清蛋白相互作用的光譜研究
    亚洲精品国产一区二区精华液| 十八禁高潮呻吟视频| 国产一卡二卡三卡精品| 母亲3免费完整高清在线观看| 无人区码免费观看不卡| 在线观看免费高清a一片| 精品国产亚洲在线| 久热这里只有精品99| 大陆偷拍与自拍| 亚洲人成电影免费在线| 多毛熟女@视频| av电影中文网址| 亚洲专区中文字幕在线| 久久香蕉国产精品| 热99re8久久精品国产| 一级毛片女人18水好多| 国产视频一区二区在线看| 成年动漫av网址| 免费女性裸体啪啪无遮挡网站| 午夜老司机福利片| av电影中文网址| 亚洲精品在线观看二区| 成人精品一区二区免费| 欧美成人午夜精品| 欧美老熟妇乱子伦牲交| 制服诱惑二区| 中文字幕av电影在线播放| 身体一侧抽搐| 大香蕉久久网| 亚洲欧美激情综合另类| 老司机深夜福利视频在线观看| 日韩欧美国产一区二区入口| 香蕉丝袜av| 精品久久久久久电影网| 国产免费男女视频| 在线观看免费视频网站a站| 90打野战视频偷拍视频| 女人精品久久久久毛片| 久久久久国产精品人妻aⅴ院 | 美女扒开内裤让男人捅视频| 国产精品1区2区在线观看. | 亚洲精品国产一区二区精华液| 久久精品亚洲精品国产色婷小说| 久久精品国产a三级三级三级| 黄色视频不卡| 午夜老司机福利片| 不卡av一区二区三区| 精品久久久久久,| 久久人妻福利社区极品人妻图片| 精品午夜福利视频在线观看一区| 水蜜桃什么品种好| 亚洲,欧美精品.| 一边摸一边抽搐一进一小说 | 中文字幕人妻丝袜一区二区| 免费观看a级毛片全部| 亚洲在线自拍视频| 久久人人97超碰香蕉20202| 在线天堂中文资源库| 在线观看日韩欧美| 91精品国产国语对白视频| 飞空精品影院首页| 国产av精品麻豆| 99久久精品国产亚洲精品| 国产成人精品久久二区二区免费| 午夜福利视频在线观看免费| 国产精品 国内视频| 午夜久久久在线观看| 曰老女人黄片| 国产蜜桃级精品一区二区三区 | 亚洲专区中文字幕在线| 久久中文字幕人妻熟女| 日韩欧美在线二视频 | 久99久视频精品免费| 亚洲 欧美一区二区三区| 黄色 视频免费看| 19禁男女啪啪无遮挡网站| 后天国语完整版免费观看| 国产99久久九九免费精品| 悠悠久久av| 大香蕉久久网| 久久久国产一区二区| 成年版毛片免费区| 高清视频免费观看一区二区| 黄色片一级片一级黄色片| 成人精品一区二区免费| 久久国产精品男人的天堂亚洲| 高清黄色对白视频在线免费看| 校园春色视频在线观看| 欧美成人午夜精品| 一个人免费在线观看的高清视频| 操美女的视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产又爽黄色视频| 又大又爽又粗| xxxhd国产人妻xxx| videos熟女内射| 欧美激情高清一区二区三区| 成人18禁在线播放| 男人舔女人的私密视频| 中文字幕另类日韩欧美亚洲嫩草| 在线播放国产精品三级| 美女视频免费永久观看网站| 亚洲九九香蕉| 成熟少妇高潮喷水视频| 热re99久久精品国产66热6| 亚洲精品美女久久久久99蜜臀| 国产欧美日韩一区二区精品| 无限看片的www在线观看| 国产真人三级小视频在线观看| 免费人成视频x8x8入口观看| 人人妻人人添人人爽欧美一区卜| 久久精品国产a三级三级三级| 日韩欧美免费精品| 伦理电影免费视频| 999精品在线视频| 国产精品九九99| 免费在线观看亚洲国产| a级毛片黄视频| 91麻豆精品激情在线观看国产 | 人妻丰满熟妇av一区二区三区 | a在线观看视频网站| 久久久久久久久免费视频了| 亚洲成人免费av在线播放| 亚洲精品中文字幕一二三四区| 国产男女内射视频| 999精品在线视频| 精品久久久久久久毛片微露脸| 亚洲成国产人片在线观看| 亚洲午夜精品一区,二区,三区| 丰满的人妻完整版| 精品熟女少妇八av免费久了| 好看av亚洲va欧美ⅴa在| 精品国产一区二区三区四区第35| 99riav亚洲国产免费| av国产精品久久久久影院| www.自偷自拍.com| 一夜夜www| 欧美中文综合在线视频| 悠悠久久av| 99精国产麻豆久久婷婷| www.精华液| 亚洲熟女毛片儿| 岛国毛片在线播放| 成熟少妇高潮喷水视频| 不卡一级毛片| 国产av精品麻豆| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文字幕日韩| 日韩欧美免费精品| 50天的宝宝边吃奶边哭怎么回事| 色在线成人网| 大香蕉久久成人网| 久久精品亚洲精品国产色婷小说| 天堂动漫精品| 老司机午夜福利在线观看视频| 国产精品一区二区精品视频观看| 一二三四社区在线视频社区8| 国产成人av教育| 成人永久免费在线观看视频| 亚洲成人国产一区在线观看| 国产97色在线日韩免费| 极品教师在线免费播放| 制服人妻中文乱码| 午夜影院日韩av| av超薄肉色丝袜交足视频| 午夜免费鲁丝| 国产一区有黄有色的免费视频| 免费人成视频x8x8入口观看| 美女福利国产在线| 国产精品久久久人人做人人爽| 在线观看日韩欧美| 国产精品免费视频内射| 俄罗斯特黄特色一大片| 久久人人爽av亚洲精品天堂| 亚洲国产毛片av蜜桃av| 国产欧美日韩一区二区精品| 国产成+人综合+亚洲专区| 两个人免费观看高清视频| 欧美日韩视频精品一区| 真人做人爱边吃奶动态| 一区福利在线观看| 免费在线观看亚洲国产| 侵犯人妻中文字幕一二三四区| 国产欧美日韩精品亚洲av| 欧美丝袜亚洲另类 | 国产亚洲精品一区二区www | 看黄色毛片网站| 精品一区二区三卡| 国产国语露脸激情在线看| www日本在线高清视频| 男女床上黄色一级片免费看| 久久久国产一区二区| 国产成人欧美| 亚洲一区二区三区欧美精品| 亚洲专区字幕在线| 人人妻人人爽人人添夜夜欢视频| 精品欧美一区二区三区在线| 自线自在国产av| 亚洲av熟女| 久久青草综合色| 天天影视国产精品| bbb黄色大片| 国产乱人伦免费视频| 黄色视频不卡| 最新美女视频免费是黄的| 十八禁高潮呻吟视频| 18在线观看网站| 精品卡一卡二卡四卡免费| 一区二区日韩欧美中文字幕| 久久久久久亚洲精品国产蜜桃av| 美女视频免费永久观看网站| 久久精品成人免费网站| 欧美性长视频在线观看| 亚洲av日韩精品久久久久久密| 亚洲avbb在线观看| 女人高潮潮喷娇喘18禁视频| 波多野结衣一区麻豆| 757午夜福利合集在线观看| 国产精品二区激情视频| 久久草成人影院| 美女视频免费永久观看网站| 在线观看午夜福利视频| 久久久久国产精品人妻aⅴ院 | 热99国产精品久久久久久7| 国产精品二区激情视频| 大片电影免费在线观看免费| av中文乱码字幕在线| 中国美女看黄片| 91在线观看av| 亚洲精品国产精品久久久不卡| 老司机福利观看| 亚洲国产欧美日韩在线播放| 99riav亚洲国产免费| 国产成+人综合+亚洲专区| 亚洲成av片中文字幕在线观看| 久久精品国产99精品国产亚洲性色 | 日韩熟女老妇一区二区性免费视频| 法律面前人人平等表现在哪些方面| 久久人人97超碰香蕉20202| 亚洲一区二区三区欧美精品| 成年人午夜在线观看视频| 国产精品九九99| 高清在线国产一区| 亚洲少妇的诱惑av| 国产欧美日韩综合在线一区二区| 下体分泌物呈黄色| 亚洲av电影在线进入| 制服人妻中文乱码| 人妻 亚洲 视频| 精品无人区乱码1区二区| 亚洲五月色婷婷综合| 少妇裸体淫交视频免费看高清 | 亚洲自偷自拍图片 自拍| 99国产精品免费福利视频| 777米奇影视久久| 日韩免费av在线播放| 最近最新中文字幕大全免费视频| netflix在线观看网站| 免费在线观看完整版高清| 欧美一级毛片孕妇| 色在线成人网| 日本vs欧美在线观看视频| 亚洲中文字幕日韩| 这个男人来自地球电影免费观看| 美女国产高潮福利片在线看| 欧美激情高清一区二区三区| 久久精品国产综合久久久| 午夜福利,免费看| 欧美一级毛片孕妇| 欧美av亚洲av综合av国产av| 精品一区二区三区四区五区乱码| 欧美日韩亚洲高清精品| 欧美日韩视频精品一区| 激情视频va一区二区三区| 久久精品国产亚洲av高清一级| 成人精品一区二区免费| 19禁男女啪啪无遮挡网站| 国产高清国产精品国产三级| 精品人妻在线不人妻| 老司机午夜十八禁免费视频| 嫩草影视91久久| 久久精品亚洲av国产电影网| 国产精品亚洲一级av第二区| 久久人妻熟女aⅴ| 啦啦啦视频在线资源免费观看| 一进一出好大好爽视频| 国产激情久久老熟女| 国产在线一区二区三区精| 欧美精品一区二区免费开放| 水蜜桃什么品种好| 老鸭窝网址在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产欧美亚洲国产| 午夜免费鲁丝| 天堂动漫精品| 国产成人av激情在线播放| 一区二区日韩欧美中文字幕| 美女扒开内裤让男人捅视频| 国产又色又爽无遮挡免费看| 91麻豆av在线| 搡老岳熟女国产| 精品久久久久久电影网| 999久久久国产精品视频| 91国产中文字幕| 亚洲专区国产一区二区| av超薄肉色丝袜交足视频| 免费在线观看视频国产中文字幕亚洲| 中文字幕另类日韩欧美亚洲嫩草| 在线观看舔阴道视频| 国产一区二区激情短视频| 国产在线精品亚洲第一网站| 国产主播在线观看一区二区| 欧洲精品卡2卡3卡4卡5卡区| 日本欧美视频一区| 欧美成人免费av一区二区三区 | 亚洲人成电影免费在线| 每晚都被弄得嗷嗷叫到高潮| 久久九九热精品免费| av免费在线观看网站| 狠狠婷婷综合久久久久久88av| 成人18禁在线播放| 国产蜜桃级精品一区二区三区 | 亚洲精品美女久久久久99蜜臀| 亚洲精品国产色婷婷电影| 免费黄频网站在线观看国产| 国产成人欧美在线观看 | 国产成人啪精品午夜网站| 午夜福利视频在线观看免费| 久久久国产成人精品二区 | 大陆偷拍与自拍| 中文字幕高清在线视频| 精品久久久久久久毛片微露脸| 国产激情欧美一区二区| 国产亚洲欧美在线一区二区| 真人做人爱边吃奶动态| 露出奶头的视频| 黄色毛片三级朝国网站| 国产亚洲欧美在线一区二区| 国产99久久九九免费精品| 欧美黑人欧美精品刺激| 怎么达到女性高潮| 久久天堂一区二区三区四区| 国产日韩一区二区三区精品不卡| 美女高潮喷水抽搐中文字幕| 久久人人爽av亚洲精品天堂| a级片在线免费高清观看视频| 国产精品 欧美亚洲| 一a级毛片在线观看| 成年女人毛片免费观看观看9 | 国产极品粉嫩免费观看在线| 久久久久久久国产电影| 99热国产这里只有精品6| 女人被躁到高潮嗷嗷叫费观| 99热国产这里只有精品6| 欧美 日韩 精品 国产| 制服诱惑二区| 国产欧美日韩一区二区三区在线| 三级毛片av免费| 天天躁夜夜躁狠狠躁躁| 亚洲欧美激情综合另类| 免费黄频网站在线观看国产| 91麻豆av在线| 人妻 亚洲 视频| √禁漫天堂资源中文www| 亚洲国产看品久久| 国产精品 欧美亚洲| 精品人妻在线不人妻| 欧美国产精品一级二级三级| 亚洲精品一二三| 国产高清videossex| 在线永久观看黄色视频| 欧美精品av麻豆av| 一本一本久久a久久精品综合妖精| 久久ye,这里只有精品| 国产精品国产av在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产高清激情床上av| √禁漫天堂资源中文www| 一区二区日韩欧美中文字幕| 国产精品永久免费网站| 久久中文看片网| 91大片在线观看| 久久天堂一区二区三区四区| 精品一区二区三区av网在线观看| 亚洲精品粉嫩美女一区| 99久久99久久久精品蜜桃| 国产精品自产拍在线观看55亚洲 | 母亲3免费完整高清在线观看| 欧美日韩精品网址| 桃红色精品国产亚洲av| 欧美日韩亚洲高清精品| 日韩一卡2卡3卡4卡2021年| 乱人伦中国视频| e午夜精品久久久久久久| 午夜影院日韩av| 久久久久视频综合| 久久精品熟女亚洲av麻豆精品| 国产亚洲欧美精品永久| 午夜福利一区二区在线看| 久久久久精品国产欧美久久久| 少妇被粗大的猛进出69影院| 久久人妻av系列| 91精品三级在线观看| av国产精品久久久久影院| 国产亚洲欧美在线一区二区| 香蕉久久夜色| 99国产精品一区二区蜜桃av | 亚洲一卡2卡3卡4卡5卡精品中文| 99久久人妻综合| 性少妇av在线| 精品国产一区二区三区四区第35| 亚洲欧洲精品一区二区精品久久久| 久久中文看片网| 高潮久久久久久久久久久不卡| 日本黄色日本黄色录像| 9191精品国产免费久久| 在线看a的网站| 下体分泌物呈黄色| 久久精品亚洲av国产电影网| 久久久久久久国产电影| 首页视频小说图片口味搜索| 亚洲国产看品久久| 最近最新中文字幕大全电影3 | 欧美中文综合在线视频| 日韩欧美三级三区| 国产成人免费无遮挡视频| 免费高清在线观看日韩| 色婷婷久久久亚洲欧美| 18禁裸乳无遮挡免费网站照片 | 亚洲精华国产精华精| 一级a爱视频在线免费观看| a在线观看视频网站| 国产成人欧美| e午夜精品久久久久久久| 日韩免费av在线播放| 国产国语露脸激情在线看| 亚洲少妇的诱惑av| 国产精品免费大片| 男女床上黄色一级片免费看| 精品卡一卡二卡四卡免费| 男女之事视频高清在线观看| 亚洲精品成人av观看孕妇| 亚洲色图 男人天堂 中文字幕| 男人舔女人的私密视频| 国产单亲对白刺激| 一进一出好大好爽视频| 日日爽夜夜爽网站| 国产亚洲精品久久久久5区| 成年版毛片免费区| 亚洲色图 男人天堂 中文字幕| 亚洲 国产 在线| 狂野欧美激情性xxxx| 久久精品国产亚洲av香蕉五月 | 国产一区二区三区在线臀色熟女 | av在线播放免费不卡| 久久久久精品国产欧美久久久| 91精品三级在线观看| 色播在线永久视频| 满18在线观看网站| 午夜两性在线视频| 国产精品秋霞免费鲁丝片| tube8黄色片| 国产成人啪精品午夜网站| 一本一本久久a久久精品综合妖精| 久久国产亚洲av麻豆专区| 欧美日韩乱码在线| 变态另类成人亚洲欧美熟女 | 在线观看午夜福利视频| 三上悠亚av全集在线观看| 两人在一起打扑克的视频| 99re在线观看精品视频| 淫妇啪啪啪对白视频| 亚洲av日韩精品久久久久久密| 欧美性长视频在线观看| 在线观看免费日韩欧美大片| 婷婷精品国产亚洲av在线 | 捣出白浆h1v1| 搡老熟女国产l中国老女人| 精品福利永久在线观看| 黄色 视频免费看| 国产成人免费观看mmmm| 精品国产一区二区三区四区第35| 淫妇啪啪啪对白视频| 国产精品自产拍在线观看55亚洲 | 一级作爱视频免费观看| 亚洲欧美色中文字幕在线| 高清在线国产一区| 黄色视频不卡| 亚洲国产精品合色在线| tocl精华| 亚洲成人免费av在线播放| av一本久久久久| 国产一区二区激情短视频| 亚洲va日本ⅴa欧美va伊人久久| 日韩大码丰满熟妇| 久久精品国产a三级三级三级| 校园春色视频在线观看| 久久久久精品人妻al黑| 不卡av一区二区三区| 亚洲五月婷婷丁香| 涩涩av久久男人的天堂| 成人影院久久| 国产成人精品在线电影| 亚洲欧美精品综合一区二区三区| 日本一区二区免费在线视频| 啦啦啦免费观看视频1| 久久热在线av| 香蕉丝袜av| 国产视频一区二区在线看| 国产成人av教育| 在线观看免费视频网站a站| 人妻丰满熟妇av一区二区三区 | 一级毛片高清免费大全| 91精品国产国语对白视频| 国产成人精品无人区| 女同久久另类99精品国产91| 精品一区二区三区四区五区乱码| 黄色女人牲交| 亚洲国产欧美网| 亚洲av熟女| 国产成人欧美| 成人三级做爰电影| 成年人免费黄色播放视频| 亚洲va日本ⅴa欧美va伊人久久| 久久精品亚洲av国产电影网| 村上凉子中文字幕在线| av天堂久久9| 国产精品自产拍在线观看55亚洲 | 视频区图区小说| av欧美777| 黄色毛片三级朝国网站| 夜夜爽天天搞| 丝瓜视频免费看黄片| 免费av中文字幕在线| 大香蕉久久网| 国产精品98久久久久久宅男小说| 操出白浆在线播放| 久久精品aⅴ一区二区三区四区| 高潮久久久久久久久久久不卡| 俄罗斯特黄特色一大片| 日本黄色视频三级网站网址 | 久久精品国产清高在天天线| 欧美日韩国产mv在线观看视频| 最新的欧美精品一区二区| 亚洲精品在线观看二区| 欧美大码av| 丝袜美腿诱惑在线| 亚洲中文字幕日韩| 国产成人欧美在线观看 | 欧美乱色亚洲激情| 在线观看免费日韩欧美大片| 在线视频色国产色| 免费观看精品视频网站| 亚洲黑人精品在线| 亚洲av成人不卡在线观看播放网| 欧美精品av麻豆av| 久久中文字幕一级| 日本黄色视频三级网站网址 | 日韩欧美三级三区| 老司机深夜福利视频在线观看| 亚洲色图 男人天堂 中文字幕| 巨乳人妻的诱惑在线观看| 一级毛片女人18水好多| 精品福利永久在线观看| 国产深夜福利视频在线观看| 国产高清videossex| 国产精品av久久久久免费| 久久香蕉精品热| 国产精品亚洲av一区麻豆| 天堂动漫精品| 曰老女人黄片| 亚洲人成伊人成综合网2020| 亚洲成人国产一区在线观看| 国产免费男女视频| 中文字幕人妻熟女乱码| 国产蜜桃级精品一区二区三区 | 这个男人来自地球电影免费观看| 午夜影院日韩av| 一级毛片女人18水好多| 老鸭窝网址在线观看| 婷婷精品国产亚洲av在线 | 久久中文字幕一级| 不卡一级毛片| 午夜久久久在线观看| 国产男女超爽视频在线观看| 一区二区三区激情视频| 自拍欧美九色日韩亚洲蝌蚪91| 成年人免费黄色播放视频| 91九色精品人成在线观看| 久久精品国产亚洲av高清一级| 激情视频va一区二区三区| 色综合欧美亚洲国产小说| 国产欧美亚洲国产| 高清欧美精品videossex| 亚洲一区二区三区欧美精品| 欧美+亚洲+日韩+国产| 亚洲欧洲精品一区二区精品久久久| 少妇裸体淫交视频免费看高清 | 一级a爱片免费观看的视频| 免费一级毛片在线播放高清视频 | 在线看a的网站| 欧美+亚洲+日韩+国产| 丰满迷人的少妇在线观看| 亚洲一码二码三码区别大吗| 伊人久久大香线蕉亚洲五| 国产高清国产精品国产三级| 亚洲欧美一区二区三区黑人| 男人舔女人的私密视频| 精品第一国产精品| 男女下面插进去视频免费观看| av线在线观看网站| 黑人欧美特级aaaaaa片| av网站免费在线观看视频| 新久久久久国产一级毛片| 亚洲 欧美一区二区三区|