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

    基于高分遙感數(shù)據(jù)的昌吉市棉花面積識別研究

    2017-11-04 07:48:51張佳琪武紅旗
    天津農(nóng)業(yè)科學 2017年10期
    關(guān)鍵詞:昌吉市棉花

    張佳琪++武紅旗

    摘 要:本文探討了如何利用高分遙感數(shù)據(jù)大范圍快速提取棉花種植面積。本研究以昌吉市為研究區(qū),基于高分一號(GF-1)遙感影像,利用實地調(diào)查的2015年昌吉市作物種植信息,選取不同的監(jiān)督分類方法,加入耕地掩膜,分析得出昌吉市棉花識別的最佳識別時相以及最佳識別方法。結(jié)果表明,棉花與其他種植作物分離程度最好月份為7月,棉花種植面積最佳識別月份為7月,支持向量機分類方法總體精度最高,總體精度為95.24%,Kappa系數(shù)0.935。棉花種植面積提取以小于6個像元為最小圖斑時結(jié)果最佳。高分一號遙感影像7月以支持向量機分類法解譯出的2015年昌吉市棉花面積為15 974 hm2,棉花面積提取精度為95.79%。

    關(guān)鍵詞:昌吉市;高分遙感數(shù)據(jù);棉花;面積識別

    中圖分類號:S127 文獻標識碼:A DOI 編碼:10.3969/j.issn.1006-6500.2017.10.012

    Research on Cotton Area Identification of Changji City based on High-revolution Remote Sensing Data

    ZHANG Jiaqi, WU Hongqi

    (College of Pratacultural and Environmental Sciences, Xinjiang Agricultural University, Urumqi, Xinjiang 830052, China)

    Abstract: This paper aimed to investigate how to use the high-revolution remote sensing data to extract cotton acreage on a large scale. This research selected Changji as the experimental area, and was based on the remote sensing images collected by GF-1. This research used the 2015 field survey of crop planting information, selected different supervised classification methods, and involved mask of arable land, to analyze and obtain the best distinguishing time phase and top identification method. The results showed that the best degree of separation of cotton from other crops was July, cotton acreage of best recognition month was July. Support vector machine method has been the supervise classification method with the highest interpretation accuracy so far, the overall accuracy was 95.24%, Kappa coefficient 0.935. Cotton acreage extraction with results was nicely when the minimum map spot was less than six pixel. GF-1 remote sensing image in July to support vector machine classification method of interpretation in 2015 Changji City, survey area was 15 974 hectares, the extraction precision of cotton area was 95.79%.

    Key words: Changji city; high-revolution remote sensing data; cotton; area identification

    新疆是中國最大的棉花生產(chǎn)基地,無論在國內(nèi)還是國際都占據(jù)重要地位,每年棉花種植面積及總產(chǎn)量受到全世界的關(guān)注[1-3]。傳統(tǒng)調(diào)查方法需要大量的人力和物力,且交通困難地區(qū)不方便到達、調(diào)查時間長,造成國家補貼分多次兌付且到位較晚,影響了棉農(nóng)還貸和準備次年生產(chǎn)。如何快速、準確地獲取棉花種植面積是解決農(nóng)業(yè)種植計劃的關(guān)鍵。隨著研究的不斷發(fā)展,國外Mathur[4]等人指出監(jiān)督分類的準確性在很大程度上取決于所使用的訓練樣本數(shù)據(jù)。同時在數(shù)據(jù)源方面,不同遙感數(shù)據(jù)源在農(nóng)作物遙感識別方面的研究應用廣泛,Landsat TM[5]、MODIS[6]、QuickBird[7]等影像數(shù)據(jù)在農(nóng)作物空間分布提取中發(fā)揮了重要作用。近年來,我國在作物遙感識別方面也開展了大量研究工作,包括單一影像的作物遙感識別、基于時間序列影像的作物遙感識別以及遙感數(shù)據(jù)與非遙感數(shù)據(jù)結(jié)合的作物遙感識別等方面[8-11]。作物遙感識別和種植面積提取的影像數(shù)據(jù)大部分來源于Landsat系列衛(wèi)星、ERS-1衛(wèi)星和SPOT系列衛(wèi)星,基于我國研發(fā)的高分一號衛(wèi)星的遙感識別研究還有待深入探索研究。本研究從面向?qū)ο笊a(chǎn)實踐角度出發(fā),結(jié)合實地調(diào)查的2015年作物種植信息,以昌吉市為研究區(qū),加入耕地掩膜,探討基于高分遙感數(shù)據(jù)的最佳識別月份、最佳識別方法,以期能大范圍快速提取棉花種植面積。

    1 研究區(qū)概況endprint

    昌吉市位于天山北麓、亞歐大陸腹地、準噶爾盆地南緣,地處東經(jīng)86°24′~87°37′,北緯43°06′~45°20′。昌吉市地貌主要由山地、平原、沙漠三部分組成。由南向北,地勢從高到低呈階梯狀,南北高差4 000多m。由于地形條件的影響,由南向北氣候差異較大,南部山區(qū)氣候特征明顯。積溫高,降雨少,夏季降水充沛,年際變化大,并且光熱資源豐富。昌吉市農(nóng)產(chǎn)品資源品質(zhì)較好主要源于獨特的光熱水土條件,非常有利于糧食作物的生長,平原地區(qū)年平均氣溫6.10 ℃,作物主要生長季節(jié)4—6月,日較差14.60 ℃左右,有利于作物干物質(zhì)合成的積累,有效積溫為3 400~3 584 ℃,平均無霜期160 d,平原地區(qū)年降水量180.1 mm,是最適宜種植番茄、葡萄、小麥的地區(qū)之一。

    2 研究方法

    本文利用提取棉花生長周期的2015年5月至9月高分一號遙感數(shù)據(jù),所獲得的數(shù)據(jù)通過ENVI 5.3進行輻射定標、大氣校正、幾何校正,并對遙感影像進行鑲嵌、裁剪,得到研究區(qū)經(jīng)過處理的影像;通過野外實際調(diào)查獲得棉花種植信息,并以此類信息建立ROI;選擇不同分類器對影像進行分類,對分類結(jié)果進行分析,從而實現(xiàn)對昌吉市棉花的面積的識別研究。

    2.1 研究區(qū)數(shù)據(jù)來源

    昌吉市作為研究區(qū),影像數(shù)據(jù)包括高分一號、高分二號的高分衛(wèi)星遙感數(shù)據(jù)來自新疆衛(wèi)星遙感應用工程中心;用于提取耕地掩膜的土地利用現(xiàn)狀圖數(shù)據(jù)來源于2010年全國第二次土地調(diào)查得到的比例尺為1∶10 000的土地利用現(xiàn)狀圖;積溫數(shù)據(jù)來自國家氣象中心頒布的1961—2010年全國積溫柵格圖。

    2.1.1 柵格數(shù)據(jù)介紹 本研究兩種柵格數(shù)據(jù)的數(shù)據(jù)源分別是高分一號(GF-1)和高分二號(GF-2)。高分一號衛(wèi)星搭載了2臺分辨率為2 m全色、8 m多光譜的高分辨率相機(PMS)和4臺分辨率為16 m的多光譜分辨率寬幅相機(WFV),現(xiàn)已下傳有效數(shù)據(jù)上萬景[12]。高分二號衛(wèi)星搭載有兩臺高分辨率1 m全色、4 m多光譜相機,觀測幅寬達到45 km。高分二號衛(wèi)星的發(fā)射標志著中國遙感衛(wèi)星進入亞米級“高分時代”[13]。本次研究所使用的高分一號衛(wèi)星遙感影像數(shù)據(jù)是2015年5月至9月云量為0的影像,范圍是新疆維吾爾自治區(qū)昌吉市,如表1所示。

    2.1.2 矢量數(shù)據(jù)介紹 通過實地調(diào)查了解2015年昌吉市的作物種植分布,并結(jié)合同一時期空間分辨率2 m的高分二號數(shù)據(jù),以2010年第二次全國土地調(diào)查中投影為UTM_WGS_1984的1∶10 000土地利用現(xiàn)狀圖為基礎(chǔ),更新該區(qū)的農(nóng)田分布圖。并根據(jù)研究區(qū)實地調(diào)查的2015年作物地塊的實際種植情況,獲取研究區(qū)2015年的土地利用現(xiàn)狀圖。根據(jù)昌吉市有效積溫圖,提取有效積溫大于3 000 ℃適宜棉花生長的耕地掩膜參與監(jiān)督分類計算,排除非耕地對監(jiān)督分類精度的影響。在后期監(jiān)督分類處理時加入掩膜,非耕地區(qū)域不參與監(jiān)督分類計算。

    2.1.3 實地調(diào)查數(shù)據(jù) (1)作物物候數(shù)據(jù)。昌吉市種植的作物主要有棉花、小麥、玉米和其他作物,棉花播種時期在4月中下旬,蕾期在6月下旬至7月下旬,花鈴期在7月下旬至8月下旬,吐絮期在9月上旬至10月中旬。小麥從拔節(jié)到成熟是從4月下旬至7月下旬。玉米從播種到成熟是從4月下旬至10月上旬[14]。

    (2)野外作物調(diào)查點數(shù)據(jù)。本研究分別于2015年8月、11月兩次實地調(diào)查采集昌吉市作物樣本經(jīng)緯度、調(diào)查時間、種植類型等信息在內(nèi)的共計1 130個野外調(diào)查點。綜合實地調(diào)查2015年昌吉市作物分布樣點信息與對遙感影像的目視解譯,由于不同月份的作物光譜特征不同,在遙感影像上表現(xiàn)出紋理不同,可以建立作物感興趣區(qū),分別有棉花、小麥、玉米和其他作物四類。選取作物樣本時一部分樣本選取作為訓練樣本,另一部分樣本選取作為驗證樣本,用來檢驗監(jiān)督分類時樣本的精度。選擇的樣本要具有代表性,能夠代表作物的顯著特征。以7月份的影像為參照,每種作物分類選取了100個訓練樣本和50個驗證樣本,共計600個作物樣本。見表2。

    2.2 數(shù)據(jù)預處理方法

    遙感影像的預處理過程就是糾正遙感成像的過程,也是對影像數(shù)據(jù)糾正與重建的過程[15-16]。受傳感器姿態(tài)、高度、速度等外在因素的影響,使遙感影像在一定程度上產(chǎn)生了幾何變形,很大程度會降低遙感影像的精度,需要對遙感圖像進行預處理工作[17-18]。本文對GF-1影像進行預處理,主要進行幾何校正、輻射定標、大氣校正和圖像裁剪等預處理。

    (1)利用ENVI 5.3提供的Radiometric Correction工具箱中的Radiometric Calibration工具完成影像的輻射定標。

    (2)以ENVI 5.3中FLASSH Atmospheric Correction對輻射定標后的影像進行大氣校正。

    (3)將得到的經(jīng)過大氣校正的影像進行幾何校正,本文使用空間分辨率2 m的GF-2影像為基準圖像,對GF-1影像進行校正,采用二次多項式、雙線性內(nèi)插法重采樣,地圖投影采用UTM_Zone_45_North,地理坐標系為WGS-84。

    (4)利用昌吉市行政區(qū)劃矢量數(shù)據(jù)對GF-1影像進行裁剪,裁剪時保持裁剪邊界與影像投影一致,利用ENVI 5.3的裁剪模塊對影像完成裁剪。

    3 結(jié)果與分析

    3.1 最佳識別月份的確定

    對不同月份的影像進行監(jiān)督分類之前需驗證樣本可分離性,根據(jù)昌吉市2015年5—9月高分一號遙感影像數(shù)據(jù)的棉花與其他三類作物可分離性報表可以看出,棉花與其他三類作物可分離性精度合格(分離度>1.8),可以進行監(jiān)督分類可分離性精度(表3)。

    根據(jù)研究區(qū)主要的四種作物分類類型分析,選擇2015年5—9月GF-1共5期遙感影像作為監(jiān)督分類的基礎(chǔ)數(shù)據(jù),基于提取的適宜棉花生長的耕地掩膜運用支持向量機的分類方法(SVM)進行監(jiān)督分類,結(jié)果如表4所示。endprint

    選取支持向量機分類方法對GF-1不同時期遙感影像進行監(jiān)督分類精度對比,其中支持向量機分類方法對GF1影像監(jiān)督分類結(jié)果進行精度分類驗證,分類精度為95.24%,Kappa系數(shù)為0.935分類精度最高的是7月影像。

    3.2 最佳識別方法的確定

    根據(jù)棉花與其他作物分類的可分離性分析(表3)可知,2015年7月GF-1遙感影像各待分地類之間的可分離性最高。研究中選擇2015年7月GF-1遙感影像進行監(jiān)督分類,以提取棉花面積種植信息。如表5所示。

    將可分離性最高7月GF-1遙感影像,分別用最大似然法、平行六面體、最小距離法、馬氏距離法、支持向量機這五種分類方法進行分類,從表5可以看出支持向量機分類法總體精度、Kappa系數(shù)均高于其他四類分類方法。

    3.3 提取棉花種植面積精度計算及驗證

    3.3.1 最小圖斑處理 在計算機上,通常面積大于6個像元的面狀地類、地物在遙感影像圖上可以解譯識別[19-20]。本研究在確定最佳識別月份和最佳識別方法后,對監(jiān)督分類后的結(jié)果先進行碎斑處理,主要目的是去除分類圖中比較均勻的圖斑上散亂分布的孤立像素。本文分別選取將小于3、6、9個像元點合并到其他圖斑,旨在研究合并最小圖斑像元點個數(shù)差異對于面積解譯精度的影響。

    3.3.2 野外實地調(diào)查點驗證解譯精度 根據(jù)監(jiān)督分類解譯的矢量數(shù)據(jù)結(jié)果圖與實地野外作物調(diào)查點數(shù)據(jù)疊加,以實地調(diào)查的野外調(diào)查點檢驗分析監(jiān)督分類精度。在ARCGIS 10.2中,通過Intersect工具對監(jiān)督分類結(jié)果矢量圖與1 130個野外調(diào)查點疊加,獲得野外調(diào)查點相對應地塊的解譯類型,根據(jù)屬性表篩選出監(jiān)督分類解譯與野外調(diào)查點信息一致的地塊數(shù)目。以野外調(diào)查點驗證監(jiān)督分類解譯精度(表6)。

    由表6可知,根據(jù)野外調(diào)查點進行精度驗證得到作物總體分類精度為75.93%,野外調(diào)查點中共計154個棉花野外調(diào)查點,其中解譯棉花部分的精度評價顯示,解譯結(jié)果與野外調(diào)查點一致的棉花野外調(diào)查點共136個,其中有18個漏分的棉花野外調(diào)查點,由于棉花與玉米生育期重疊時間比較長,其中12個漏分的棉花野外調(diào)查點被劃分為玉米野外調(diào)查點,棉花野外調(diào)查點漏分誤差11.69%。在驗證解譯精度的總體野外調(diào)查點中共有57個包括玉米、小麥、其他作物的野外調(diào)查點被錯分為棉花野外調(diào)查點,棉花錯分誤差29.53%,棉花解譯精度為88.31%。

    3.3.3 作物樣本驗證解譯精度 根據(jù)野外調(diào)查樣本數(shù)據(jù)建立的600個地表真實數(shù)據(jù)感興趣區(qū)(ROI)與監(jiān)督分類解譯的矢量數(shù)據(jù)結(jié)果圖疊加,可以獲得作物樣本相對應地塊的解譯類型。以作物樣本檢驗分析監(jiān)督分類精度,根據(jù)屬性表篩選出監(jiān)督分類解譯與作物樣本一致的地塊數(shù)目。以野外調(diào)查點為基礎(chǔ)建立的ROI驗證監(jiān)督分類解譯精度(表7)。

    由表7可知,根據(jù)野外調(diào)查點建立的作物樣本進行精度驗證得到作物總體分類精度為94.50%,其中解譯棉花部分的精度評價顯示,作物樣本中共計150個棉花作物樣本,其中解譯結(jié)果與作物樣本一致的棉花樣本共146個,由于棉花與玉米生育期重疊時間比較長,其中有4個漏分的棉花作物樣本,被劃分為玉米作物樣本,棉花作物漏分誤差2.67%。在驗證解譯精度的總體作物樣本中共有5個包括玉米、小麥、其他作物在內(nèi)的作物樣本被錯分為棉花作物樣本,棉花錯分誤差3.31%,棉花解譯精度為97.33%。

    本研究選取的作物樣本很合理,無論是對于作物總體分類精度還是棉花分類精度,選取的包括訓練樣本和驗證樣本的作物樣本對棉花種植分布進行精度檢驗顯示作物解譯精度和棉花解譯精度很好。

    3.3.4 提取棉花種植面積精度計算 在ENVI 5.3中得到的最佳監(jiān)督分類方法解譯出的最佳識別月份的遙感影像數(shù)據(jù),其監(jiān)督分類結(jié)果數(shù)據(jù)類型為柵格數(shù)據(jù),然后通過ARCGIS 10.2軟件中Raster to Polygon工具,將處理后的不同月份監(jiān)督分類結(jié)果圖矢量化。通過以上處理分析后,ARCGIS 10.2中,根據(jù)屬性表計算棉花分布區(qū)域的面積,并根據(jù)昌吉市統(tǒng)計局公布的棉花面積與解譯面積計算精度,提取精度=1-[ABS(調(diào)查面積-解譯面積)/調(diào)查面積],根據(jù)合并不同大小最小圖斑處理結(jié)果,選取棉花最佳識別面積方案。

    由表8可以看出,選擇不同像元點個數(shù)最小圖斑合并對解譯面積影響存在較為明顯影響,隨著合并像元圖斑個數(shù)的增加,合并0個像元圖斑對應的解譯面積相比合并3個像元圖斑對應的解譯面積差異較大,相差21 hm2。選擇合并9個像元圖斑對應的解譯面積與合并6個像元圖斑對應的解譯面積差異不顯著,相差1 hm2。選擇合并6個像元圖斑的棉花種植面積提取結(jié)果最準確。

    4 結(jié)論與討論

    4.1 結(jié) 論

    本研究選取2015年5—9月GF-1共5景遙感衛(wèi)星影像,對研究區(qū)昌吉市進行棉花種植面積信息提取研究,得出以下結(jié)論。

    (1)作物樣本可分離性最高影像為7月,GF-1支持向量機分類法確定的最佳時相也是7月。7月是識別棉花面積的最佳時相,可分離性高,分類精度高,提取棉花種植面積最為準確。

    (2)利用監(jiān)督分類方法中的馬氏距離、最小距離、平行六面體、基于支持向量機、最大似然法分別對GF-1遙感影像數(shù)據(jù)進行監(jiān)督分類,支持向量機分類法為解譯精度最高的監(jiān)督分類方法。

    (3)棉花面積解譯識別的最佳方案是,基于GF-1遙感數(shù)據(jù)7月為最佳識別月份的基于支持向量機為最佳識別方法,選取合并小于6個像元的最小圖斑處理為最佳,實現(xiàn)精度最高的棉花面積識別。2015年昌吉市棉花提取面積為15 974 hm2,提取精度為95.79%。

    4.2 討 論

    通過野外實地調(diào)查并結(jié)合2015年昌吉市作物物候歷可知,7月是棉花生長最旺盛的時期,棉花生物量達到飽和,而玉米處于營養(yǎng)生長期,小麥已收獲,通過影像上的紋理最容易區(qū)分棉花與其他農(nóng)作物。9月玉米還在成熟期,其他作物都已收獲,而棉花還在吐絮期,在此生長期棉花與玉米并不容易區(qū)分。面積提取結(jié)果表明,選取合并小于6個像元的最小圖斑處理為最佳,7月的識別精度最高最易于識別,提取精度高達95.79%,支持向量機是精度最高的監(jiān)督分類方法。endprint

    此外,本次研究也存在一些不足之處。主要有以下幾點。

    (1)本研究監(jiān)督分類選取空間分辨率為16 m的GF-1遙感影像數(shù)據(jù),未能與其他數(shù)據(jù)源進行對比研究,不同分辨率的遙感影像數(shù)據(jù)可能會影響監(jiān)督分類精度。

    (2)GF-1影像數(shù)據(jù)受時間分辨率影響,可選擇使用的影像時間范圍較小,所獲得的影像時間與作物生長期重合度有待提高。

    (3)本研究采用監(jiān)督分類的方法識別棉花種植面積。為更加快速解譯棉花分布范圍并獲得棉花種植面積,可以通過對數(shù)據(jù)的挖掘,獲得基于決策樹的作物分類方法,以提高地面調(diào)查工作量,提高棉花種植面積的提取周期。

    參考文獻:

    [1]蒙繼華, 杜鑫, 張淼, 等. 物候信息在大范圍作物長勢遙感監(jiān)測中的應用[J]. 遙感技術(shù)與應用, 2014, 29(2): 278-285.

    [2]李郁竹. 農(nóng)作物氣象衛(wèi)星遙感監(jiān)測和估產(chǎn)研究進展及前景探討[J]. 氣象科技, 1997, (3): 30-36.

    [3]江東, 王乃斌, 楊小喚. 我國糧食作物衛(wèi)星遙感估產(chǎn)的研究[J]. 自然雜志, 1999, 21(6): 351-355.

    [4]MATHUR A, FOODY G M.Crop classification by support vector machine with intelligently selected training data for an operational application[J]. International journal of remote sensing, 2008,29(8): 2227-2240.

    [5]ZHONG L, GONG P, BIGING G S. Efficient corn and soybean mapping with temporal extendability: A multi-year experiment using Landsat imagery[J]. Remote sensing of environment, 2014, 140(1): 1-13.

    [6]PAN Y, LI L, ZHANG J, et al. Crop area estimation based on MODIS-EVI time series according to distinct characteristics of key phenology phases: a case study of winter wheat area estimation in small-scale area[J]. Journal of remote sensing, 2011, 15(3): 578-594.

    [7]TURKER M, OZDARICI A. Field-based crop classification using SPOT4, SPOT5, IKONOS and QuickBird imagery for agricultural areas: a comparison study[J]. International journal of remote sensing, 2011, 32(24): 9735-9768.

    [8]王瓊, 王克如, 李少昆, 等. HJ衛(wèi)星數(shù)據(jù)在棉花種植面積提取中的應用研究[J]. 棉花學報, 2012, 23(6): 503-510.

    [9]熊勤學, 黃敬峰. 利用NDVI指數(shù)時序特征監(jiān)測秋收作物種植面積[J]. 農(nóng)業(yè)工程學報,2009,1(24): 144-148.

    [10]彭代亮. 基于統(tǒng)計與MODIS數(shù)據(jù)的水稻遙感估產(chǎn)方法研究[D]. 杭州: 浙江大學, 2009.

    [11]劉佳岐. 基于Landsat8遙感影像的扶風縣蘋果園地信息提取研究[D]. 楊凌: 西北農(nóng)林科技大學, 2015.

    [12]陳文倩, 丁建麗, 王嬌, 等. 基于高分一號影像的土地覆被分類方法初探[J]. 干旱區(qū)地理, 2016, 39(1): 182-189.

    [13]潘騰, 關(guān)暉, 賀瑋. “高分二號”衛(wèi)星遙感技術(shù)[J]. 航天返回與遙感, 2015, 36(4): 16-24.

    [14]趙良斌, 曹衛(wèi)彬, 唐春華, 等. 新疆棉花遙感識別最佳時相的選擇[J]. 新疆農(nóng)業(yè)科學, 2008, 45(4): 618-622.

    [15]劉排英. 基于光譜面積和IHS變換的水體提取的研究[D]. 長沙: 中南大學, 2010.

    [16]何同弟. 高光譜圖像的分類技術(shù)研究[D]. 重慶: 重慶大學, 2014.

    [17]YANG G, STWART C, SOFKA M, et al. Registration of challenging image pairs:initialization, estimation, and decision[J]. IEEE transactions on pattern analysis and machine intellingence, 2007, 29(11): 1973-1989.

    [18]王媛, 葉思菁, 岳彥利, 等. 面向高分一號遙感影像的自動幾何配準算法對比[J]. 農(nóng)業(yè)機械學報, 2015, 46(12): 260-266.

    [19]韋燕燕. 高分辨率遙感數(shù)據(jù)最小上圖圖斑面積的研究[D]. 焦作: 河南理工大學, 2012.

    [20]陳粲, 張忠啟, 史學正, 等. 不同分類級別土壤矢量圖與最小可分柵格的關(guān)系研究[J]. 土壤學報, 2011, 48(3): 470-477.endprint

    猜你喜歡
    昌吉市棉花
    新疆昌吉市水資源承載能力分析研究
    地下水(2022年6期)2022-11-29 10:28:52
    棉花是花嗎?
    昌吉市地下水水位動態(tài)變化特征研究
    地下水(2021年5期)2021-11-10 12:44:58
    棉花
    小讀者(2020年4期)2020-06-16 03:33:54
    快樂城的棉花糖
    小布老虎(2016年8期)2016-12-01 05:46:26
    昌吉市土地利用變化對生態(tài)環(huán)境影響分析
    心中的“棉花糖”
    第三講 棉花肥害診斷及其防治
    棉花-新陸早49號
    欧美一区二区国产精品久久精品| 我的老师免费观看完整版| 69人妻影院| 日韩高清综合在线| 级片在线观看| 国语自产精品视频在线第100页| 大型黄色视频在线免费观看| 精品久久久久久久久av| 人妻少妇偷人精品九色| 亚洲丝袜综合中文字幕| 亚洲国产高清在线一区二区三| 欧美日韩在线观看h| 国产精品人妻久久久影院| 亚洲婷婷狠狠爱综合网| 天堂影院成人在线观看| 亚洲专区国产一区二区| 日本与韩国留学比较| АⅤ资源中文在线天堂| 可以在线观看的亚洲视频| 亚洲在线自拍视频| 国产午夜福利久久久久久| 99热网站在线观看| 国产精华一区二区三区| 亚洲欧美日韩高清在线视频| 麻豆久久精品国产亚洲av| 成人国产麻豆网| 校园人妻丝袜中文字幕| 日本黄色片子视频| 国产精品99久久久久久久久| www日本黄色视频网| 超碰av人人做人人爽久久| 国产中年淑女户外野战色| 国产精品永久免费网站| 在线天堂最新版资源| 99久久精品热视频| 毛片女人毛片| 人妻丰满熟妇av一区二区三区| 免费观看在线日韩| 亚州av有码| 高清毛片免费观看视频网站| 天堂网av新在线| 你懂的网址亚洲精品在线观看 | 一本一本综合久久| 亚洲欧美日韩高清在线视频| 又粗又爽又猛毛片免费看| 国产黄a三级三级三级人| 亚洲av免费高清在线观看| 久久久久久久亚洲中文字幕| 亚洲成人久久性| 国产精品一区二区三区四区久久| 久久久国产成人免费| 亚洲在线自拍视频| 日韩精品青青久久久久久| 国国产精品蜜臀av免费| 麻豆乱淫一区二区| 日韩大尺度精品在线看网址| 看免费成人av毛片| 日日摸夜夜添夜夜添av毛片| 亚洲欧美日韩高清专用| 午夜爱爱视频在线播放| 老司机福利观看| 直男gayav资源| 精品久久久噜噜| 人人妻人人看人人澡| 国产精品乱码一区二三区的特点| 亚洲中文字幕日韩| 又粗又爽又猛毛片免费看| 国产精品av视频在线免费观看| 99久久中文字幕三级久久日本| 成年女人永久免费观看视频| 舔av片在线| av卡一久久| 午夜爱爱视频在线播放| 欧美成人精品欧美一级黄| 一本精品99久久精品77| 欧美日本视频| 大香蕉久久网| 日日撸夜夜添| 久久午夜亚洲精品久久| 97超碰精品成人国产| 婷婷精品国产亚洲av在线| 激情 狠狠 欧美| 干丝袜人妻中文字幕| 日韩国内少妇激情av| 97人妻精品一区二区三区麻豆| 午夜亚洲福利在线播放| 99热这里只有精品一区| 亚洲欧美精品自产自拍| 97在线视频观看| 国产美女午夜福利| 成人永久免费在线观看视频| 深夜a级毛片| 女人十人毛片免费观看3o分钟| 午夜福利成人在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 午夜免费男女啪啪视频观看 | 99热这里只有精品一区| 精品不卡国产一区二区三区| 有码 亚洲区| 一进一出抽搐动态| 91狼人影院| 日本色播在线视频| 搡老妇女老女人老熟妇| 99久久成人亚洲精品观看| 插逼视频在线观看| 国产精品爽爽va在线观看网站| 一级黄片播放器| 国产精品亚洲美女久久久| 美女免费视频网站| 国产成人a∨麻豆精品| 精品人妻一区二区三区麻豆 | 欧美激情在线99| 国产久久久一区二区三区| 看十八女毛片水多多多| 国产一区二区激情短视频| 最近的中文字幕免费完整| 91在线观看av| 精品一区二区免费观看| 成人av在线播放网站| 亚洲成av人片在线播放无| 99国产极品粉嫩在线观看| 亚洲av一区综合| 九九爱精品视频在线观看| 一级黄色大片毛片| 国产一区二区激情短视频| 成人鲁丝片一二三区免费| 午夜福利视频1000在线观看| 国内久久婷婷六月综合欲色啪| 女同久久另类99精品国产91| 久久久精品94久久精品| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人久久性| 国产精品一区二区三区四区久久| 91精品国产九色| 精品不卡国产一区二区三区| 亚洲性夜色夜夜综合| 99久国产av精品| av福利片在线观看| 九九热线精品视视频播放| 干丝袜人妻中文字幕| 亚洲欧美日韩东京热| 亚洲欧美成人综合另类久久久 | 岛国在线免费视频观看| 久久人妻av系列| 国内少妇人妻偷人精品xxx网站| 美女xxoo啪啪120秒动态图| 99国产极品粉嫩在线观看| 99热6这里只有精品| www.色视频.com| 无遮挡黄片免费观看| 长腿黑丝高跟| a级毛色黄片| 欧美最黄视频在线播放免费| 精品一区二区三区人妻视频| 在线观看美女被高潮喷水网站| 亚洲av免费高清在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲欧美精品自产自拍| 嫩草影院精品99| 亚洲人与动物交配视频| 免费看光身美女| 久久婷婷人人爽人人干人人爱| 亚洲精品国产av成人精品 | 51国产日韩欧美| 97超碰精品成人国产| 精品人妻视频免费看| 日日摸夜夜添夜夜爱| 国产 一区 欧美 日韩| 色哟哟哟哟哟哟| 欧美性猛交黑人性爽| 国产日本99.免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产高潮美女av| 亚洲av二区三区四区| 精品人妻一区二区三区麻豆 | 亚洲久久久久久中文字幕| 成人午夜高清在线视频| 久久精品综合一区二区三区| 亚洲综合色惰| 国产精品,欧美在线| 亚洲欧美日韩无卡精品| 日本与韩国留学比较| 搡老岳熟女国产| 亚洲美女黄片视频| 国产精品福利在线免费观看| 欧美日韩在线观看h| 久久久久久国产a免费观看| 午夜激情福利司机影院| 国产精品1区2区在线观看.| 长腿黑丝高跟| 日本三级黄在线观看| 国产成人精品久久久久久| 久久久久久国产a免费观看| 如何舔出高潮| 大型黄色视频在线免费观看| 夜夜爽天天搞| 在线国产一区二区在线| 免费av毛片视频| 亚洲欧美中文字幕日韩二区| 国产白丝娇喘喷水9色精品| 日韩三级伦理在线观看| 极品教师在线视频| 国产女主播在线喷水免费视频网站 | 搡老妇女老女人老熟妇| 青春草视频在线免费观看| 蜜臀久久99精品久久宅男| 丝袜美腿在线中文| 小蜜桃在线观看免费完整版高清| 草草在线视频免费看| 麻豆国产av国片精品| 一个人看视频在线观看www免费| 一进一出好大好爽视频| 国产麻豆成人av免费视频| 国产精品乱码一区二三区的特点| 亚洲熟妇熟女久久| 欧洲精品卡2卡3卡4卡5卡区| 日日啪夜夜撸| 亚洲最大成人手机在线| 久久久久国产网址| 我要看日韩黄色一级片| 国产综合懂色| 美女大奶头视频| 一级毛片我不卡| 中文字幕av在线有码专区| 哪里可以看免费的av片| 别揉我奶头 嗯啊视频| 两个人的视频大全免费| 夜夜爽天天搞| 欧美激情国产日韩精品一区| 变态另类丝袜制服| 国内久久婷婷六月综合欲色啪| 精品一区二区三区人妻视频| 嫩草影院入口| 国产精品不卡视频一区二区| 日韩欧美三级三区| 欧美成人免费av一区二区三区| 亚洲第一区二区三区不卡| 成人漫画全彩无遮挡| 精品久久久久久久久久久久久| 日韩av不卡免费在线播放| 欧美在线一区亚洲| 欧美国产日韩亚洲一区| 色在线成人网| 少妇人妻精品综合一区二区 | 寂寞人妻少妇视频99o| 精品免费久久久久久久清纯| 看片在线看免费视频| 久久精品久久久久久噜噜老黄 | 成人无遮挡网站| 人人妻人人看人人澡| 色综合亚洲欧美另类图片| 99热这里只有是精品在线观看| 日本成人三级电影网站| 嫩草影院新地址| 精品一区二区三区av网在线观看| 婷婷六月久久综合丁香| 美女高潮的动态| 国产午夜福利久久久久久| 一级毛片久久久久久久久女| 国产激情偷乱视频一区二区| 国产精品国产三级国产av玫瑰| 亚洲成a人片在线一区二区| 精品一区二区三区视频在线观看免费| 亚洲高清免费不卡视频| 国内精品一区二区在线观看| 国产精品女同一区二区软件| 联通29元200g的流量卡| 一进一出好大好爽视频| 国产精品久久电影中文字幕| 精品人妻视频免费看| 成人av在线播放网站| .国产精品久久| 免费大片18禁| 国产精品不卡视频一区二区| 97在线视频观看| 老女人水多毛片| 国产午夜精品久久久久久一区二区三区 | 国产久久久一区二区三区| 嫩草影视91久久| 天天躁日日操中文字幕| 日本熟妇午夜| 18禁在线播放成人免费| 成人精品一区二区免费| 国产亚洲91精品色在线| 亚洲天堂国产精品一区在线| 一区福利在线观看| 精品日产1卡2卡| 91午夜精品亚洲一区二区三区| 在线观看午夜福利视频| 99久久九九国产精品国产免费| 狂野欧美激情性xxxx在线观看| 精品人妻一区二区三区麻豆 | 国产真实伦视频高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 全区人妻精品视频| 亚洲乱码一区二区免费版| 搡女人真爽免费视频火全软件 | 别揉我奶头 嗯啊视频| 精品福利观看| 又粗又爽又猛毛片免费看| 亚洲精品国产成人久久av| 干丝袜人妻中文字幕| 国产探花在线观看一区二区| 尾随美女入室| 搡女人真爽免费视频火全软件 | 国产aⅴ精品一区二区三区波| 亚洲色图av天堂| 搡女人真爽免费视频火全软件 | 国内精品一区二区在线观看| 国产精品久久久久久久电影| 亚洲熟妇熟女久久| 成人鲁丝片一二三区免费| 少妇高潮的动态图| 国产私拍福利视频在线观看| 国产欧美日韩一区二区精品| 中文字幕av在线有码专区| 少妇人妻一区二区三区视频| 日日摸夜夜添夜夜爱| 真实男女啪啪啪动态图| 人妻夜夜爽99麻豆av| 亚洲无线在线观看| 美女高潮的动态| 91狼人影院| 成人三级黄色视频| 一个人免费在线观看电影| 麻豆乱淫一区二区| 成人永久免费在线观看视频| 深夜a级毛片| 卡戴珊不雅视频在线播放| 国产乱人偷精品视频| 国产成人精品久久久久久| 亚洲精品粉嫩美女一区| 久久人人爽人人片av| 波多野结衣巨乳人妻| 亚洲国产高清在线一区二区三| 久久精品国产自在天天线| 国产成人福利小说| 午夜爱爱视频在线播放| 最近最新中文字幕大全电影3| 亚洲经典国产精华液单| 国产 一区精品| 天美传媒精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 波野结衣二区三区在线| 久久精品国产亚洲av香蕉五月| 麻豆国产av国片精品| 国产在线精品亚洲第一网站| 中文亚洲av片在线观看爽| 长腿黑丝高跟| 99热全是精品| 亚洲精华国产精华液的使用体验 | 九九在线视频观看精品| 美女黄网站色视频| 插阴视频在线观看视频| 欧美潮喷喷水| 色在线成人网| 69av精品久久久久久| 色尼玛亚洲综合影院| 久久久久国产网址| 免费av毛片视频| 久久久久性生活片| 精品一区二区三区视频在线观看免费| 深夜精品福利| 国产亚洲精品久久久com| 成人美女网站在线观看视频| 看免费成人av毛片| 久久精品久久久久久噜噜老黄 | 色av中文字幕| 看免费成人av毛片| 黄片wwwwww| 亚洲精品影视一区二区三区av| 日本三级黄在线观看| 你懂的网址亚洲精品在线观看 | 日日摸夜夜添夜夜爱| 国产精品一二三区在线看| 免费av观看视频| 大型黄色视频在线免费观看| 国产男人的电影天堂91| 久久天躁狠狠躁夜夜2o2o| 午夜激情福利司机影院| 一级毛片aaaaaa免费看小| 婷婷精品国产亚洲av| 一区福利在线观看| 岛国在线免费视频观看| 国产亚洲欧美98| 精品久久久噜噜| 午夜视频国产福利| 在线观看免费视频日本深夜| 国产精品免费一区二区三区在线| 69人妻影院| 日本三级黄在线观看| 国产精品亚洲一级av第二区| 黑人高潮一二区| a级毛色黄片| 精品福利观看| 长腿黑丝高跟| 国产麻豆成人av免费视频| 特级一级黄色大片| 久久婷婷人人爽人人干人人爱| 亚洲国产精品成人综合色| 成熟少妇高潮喷水视频| 午夜视频国产福利| 日本与韩国留学比较| 国产av不卡久久| 国产一区二区激情短视频| 亚洲最大成人手机在线| 成人亚洲精品av一区二区| a级一级毛片免费在线观看| 国产亚洲91精品色在线| 九九爱精品视频在线观看| 国产精品99久久久久久久久| 久久婷婷人人爽人人干人人爱| 国产片特级美女逼逼视频| 精品国内亚洲2022精品成人| 性欧美人与动物交配| 亚洲欧美日韩高清专用| 综合色丁香网| 非洲黑人性xxxx精品又粗又长| 91在线观看av| 国产精品久久久久久精品电影| 1024手机看黄色片| 日韩大尺度精品在线看网址| 亚洲精品成人久久久久久| 韩国av在线不卡| 成人毛片a级毛片在线播放| 18禁在线无遮挡免费观看视频 | 日本精品一区二区三区蜜桃| 一进一出抽搐动态| 人人妻人人看人人澡| 国产三级在线视频| 日本免费a在线| 亚洲精品456在线播放app| 永久网站在线| 亚洲国产欧洲综合997久久,| 91精品国产九色| 精品人妻偷拍中文字幕| 女人十人毛片免费观看3o分钟| 久久久久久久久久久丰满| 2021天堂中文幕一二区在线观| 国产亚洲91精品色在线| 极品教师在线视频| 自拍偷自拍亚洲精品老妇| 美女内射精品一级片tv| 极品教师在线视频| 亚洲电影在线观看av| 校园人妻丝袜中文字幕| 成人永久免费在线观看视频| 变态另类丝袜制服| 亚洲久久久久久中文字幕| 久久久欧美国产精品| 美女黄网站色视频| 秋霞在线观看毛片| 白带黄色成豆腐渣| 国产精品,欧美在线| 男人狂女人下面高潮的视频| 欧美区成人在线视频| 免费看美女性在线毛片视频| 国产综合懂色| 久久久午夜欧美精品| 俺也久久电影网| 欧美一级a爱片免费观看看| 村上凉子中文字幕在线| 国产精品久久视频播放| 国产亚洲精品综合一区在线观看| 国产亚洲91精品色在线| 欧美日韩在线观看h| 中文字幕av在线有码专区| 蜜臀久久99精品久久宅男| 国产一区二区三区在线臀色熟女| 亚洲精品456在线播放app| 国内精品美女久久久久久| 精品免费久久久久久久清纯| 一级毛片aaaaaa免费看小| 国产精品人妻久久久影院| 欧美国产日韩亚洲一区| 国产三级在线视频| 夜夜夜夜夜久久久久| 国产精品久久久久久亚洲av鲁大| 久久午夜亚洲精品久久| 熟女电影av网| 亚洲性夜色夜夜综合| 天堂网av新在线| 深夜a级毛片| 韩国av在线不卡| 国产视频一区二区在线看| 日日摸夜夜添夜夜爱| 欧美激情国产日韩精品一区| 国国产精品蜜臀av免费| 十八禁国产超污无遮挡网站| 午夜a级毛片| 男人狂女人下面高潮的视频| 久久久午夜欧美精品| 能在线免费观看的黄片| 免费在线观看影片大全网站| 99热网站在线观看| 久久精品国产99精品国产亚洲性色| 午夜福利在线在线| 热99re8久久精品国产| 精品人妻偷拍中文字幕| 午夜免费激情av| 国产在线男女| 亚洲内射少妇av| 成人综合一区亚洲| 天堂动漫精品| 国产又黄又爽又无遮挡在线| 热99在线观看视频| 欧美成人免费av一区二区三区| 黄片wwwwww| 亚洲久久久久久中文字幕| 日韩欧美一区二区三区在线观看| 亚洲图色成人| 国产三级在线视频| 蜜桃久久精品国产亚洲av| 亚州av有码| 国产伦在线观看视频一区| 久久天躁狠狠躁夜夜2o2o| 天天躁夜夜躁狠狠久久av| 一卡2卡三卡四卡精品乱码亚洲| 免费高清视频大片| 看免费成人av毛片| 欧美成人一区二区免费高清观看| 女同久久另类99精品国产91| 亚洲七黄色美女视频| 一级黄片播放器| 人妻制服诱惑在线中文字幕| 听说在线观看完整版免费高清| 国产伦在线观看视频一区| 狠狠狠狠99中文字幕| 国产探花极品一区二区| 精品免费久久久久久久清纯| .国产精品久久| 久久午夜亚洲精品久久| 你懂的网址亚洲精品在线观看 | av在线亚洲专区| 午夜亚洲福利在线播放| 一个人看视频在线观看www免费| 嫩草影院新地址| 国产精品人妻久久久久久| 中国国产av一级| 成人永久免费在线观看视频| 嫩草影院新地址| 国产精品国产高清国产av| 卡戴珊不雅视频在线播放| 波野结衣二区三区在线| 国产成人影院久久av| 亚洲一区二区三区色噜噜| 麻豆成人午夜福利视频| 在线观看av片永久免费下载| 麻豆乱淫一区二区| 少妇人妻一区二区三区视频| 色av中文字幕| 99久久精品热视频| 国产色爽女视频免费观看| 日韩欧美三级三区| 欧美潮喷喷水| 亚洲性夜色夜夜综合| 成人高潮视频无遮挡免费网站| av卡一久久| 插逼视频在线观看| 麻豆国产av国片精品| 日本五十路高清| 男人狂女人下面高潮的视频| 久久久成人免费电影| 少妇的逼水好多| 少妇裸体淫交视频免费看高清| 精品熟女少妇av免费看| 女的被弄到高潮叫床怎么办| 国产av不卡久久| 亚洲色图av天堂| 亚洲精品国产av成人精品 | 国产淫片久久久久久久久| 高清午夜精品一区二区三区 | 男插女下体视频免费在线播放| 白带黄色成豆腐渣| 性色avwww在线观看| 亚洲中文字幕日韩| av女优亚洲男人天堂| 午夜福利在线观看免费完整高清在 | 亚洲精品在线观看二区| 99热6这里只有精品| 成人三级黄色视频| 九九爱精品视频在线观看| 国产私拍福利视频在线观看| 欧美国产日韩亚洲一区| 久久久久久久久久黄片| 国产av麻豆久久久久久久| 免费电影在线观看免费观看| 内地一区二区视频在线| 亚洲精品色激情综合| 精品国内亚洲2022精品成人| 最近视频中文字幕2019在线8| 18+在线观看网站| 丝袜美腿在线中文| 国产老妇女一区| 成人性生交大片免费视频hd| 亚洲色图av天堂| 深夜a级毛片| 亚洲av电影不卡..在线观看| 国产精品久久久久久精品电影| 国产精品一区二区免费欧美| 男女边吃奶边做爰视频| 中文字幕av在线有码专区| 亚洲av一区综合| 一级黄色大片毛片| 久久午夜亚洲精品久久| 日日啪夜夜撸| 欧美另类亚洲清纯唯美| 日韩三级伦理在线观看| 人人妻人人看人人澡| 免费人成视频x8x8入口观看| 精品一区二区三区视频在线观看免费| 三级毛片av免费| 美女免费视频网站| 久久久a久久爽久久v久久| 欧美性猛交╳xxx乱大交人|