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

    基于實驗室高光譜數據的大圍山森林土壤氧化鐵全量反演①

    2021-11-01 07:05:38周衛(wèi)軍崔浩杰
    土壤 2021年4期
    關鍵詞:大圍山包絡線全量

    譚 潔,陳 嚴,周衛(wèi)軍,崔浩杰,劉 沛

    基于實驗室高光譜數據的大圍山森林土壤氧化鐵全量反演①

    譚 潔1,陳 嚴1,周衛(wèi)軍2*,崔浩杰2,劉 沛1

    (1湖南農業(yè)大學風景園林與藝術設計學院,長沙 410128;2湖南農業(yè)大學資源環(huán)境學院,長沙 410128)

    氧化鐵是土壤中含鐵礦物的主體,是土壤發(fā)育和土壤分類最明顯和最有用的指標之一。本文以湖南省大圍山森林土壤為研究對象,通過實驗室化學成分測定和光譜采集,在光譜預處理及組合變換基礎上,采用相關性分析篩選土壤氧化鐵全量的敏感波段,并分別建立多元逐步回歸和偏最小二乘回歸反演模型。結果表明:不同土壤光譜曲線趨勢基本一致,均形似陡坎,且在420 ~ 580 nm波段,土壤氧化鐵全量與光譜反射率呈負相關關系;不同的光譜數據變換方式可以提高光譜與氧化鐵全量的相關性,Savitzky-Golay(S-G)平滑和去包絡線相結合優(yōu)于其他預處理方法;土壤氧化鐵全量的特征波段主要為392、427、529、523、549、559、565、570、994和1 040 nm,偏最小二乘回歸模型比多元逐步回歸模型具有更好的穩(wěn)定性,適合于快速估算紅黃壤區(qū)森林土壤氧化鐵全量。

    土壤光譜;氧化鐵;多元逐步回歸;偏最小二乘回歸

    鐵是巖石圈的主要元素之一,是土壤中重要的微量元素,不僅對植物的生長發(fā)育和農產品品質具有極其重要的作用,還對人體健康有著重要的意義。土壤中的氧化鐵是在地表特定水熱條件下,由硅酸鹽類礦物經徹底風化作用而形成的重要土壤物質。由于土壤氧化鐵的聚積和遷移活動受到化學風化、生物循環(huán)等因素影響,反映了土壤的淋溶過程、風化發(fā)育程度以及土壤的地帶性分布特征,常被作為描述土壤發(fā)育和土壤分類的最明顯和最有用的指標之一[1-3]。傳統(tǒng)土壤化學分析方法存在著數據獲取難、分析成本高、時效性差等問題,隨著土壤紅外光譜[4]、漫反射光譜[5]、高光譜[6]等光譜技術的快速發(fā)展,為快速、準確、定量獲取土壤氧化鐵全量信息提供了新的科學途徑。

    土壤是由不同理化性質的物質混合組成的系統(tǒng),其高光譜特性受到有機質、水分、黏土礦物、氧化鐵、母質及土壤質地等因素共同作用。氧化鐵在土壤中常處于還原淋溶和氧化淀積的交替過程,氧化鐵的形態(tài)、顆粒粗細以及是否以膜狀態(tài)出現,對土壤顏色產生了很大的影響,從而也顯著影響了土壤光譜特征[7]。Stone等人[8]的研究認為氧化鐵的吸收帶主要集中在反射光譜曲線的700 nm和 900 nm處;Hunt和Ashley[9]研究發(fā)現Fe2+在1 000 ~ 1 100 nm 附近可產生一個最常見的強而寬的譜帶,Fe3+形成的吸收帶主要位于 900 nm 附近;Galv?o等人[10]研究巴西熱帶土壤時發(fā)現土壤中氧化鐵含量與 1 710 nm 處的反射率的相關性顯著。何挺等[11]的研究發(fā)現土壤氧化鐵含量與1 550 ~ 2 320 nm 波段反射率的相關性很強,且含量越高,在400 ~ 1 200 nm 波段的反射率越低,利用紅波段和藍波段平均反射率的對數構成的土壤氧化鐵指數,在高光譜圖像解譯中具有應用潛力;彭杰等[12]認為隨游離氧化鐵含量的增加,350 ~ 570 nm 波段土壤反射率降低,570 ~ 2 500 nm 波段的反射率增加?;趯ν寥姥趸F的高光譜特征分析,學者們還采用多元逐步回歸模型[13]、神經網絡模型(BPNN)[14]、偏最小二乘回歸模型(PLSR)[15]、多層感知器神經網絡模型(MLP)[16]、隨機森林回歸算法(RFR)[17]、組合模型[18]和模糊識別[19]等方法,分別估算了土壤氧化鐵全量、游離鐵、針鐵礦和赤鐵礦含量。由此可見,土壤氧化鐵全量對土壤光譜存在明顯影響,可借助高光譜技術對土壤氧化鐵全量進行估測。然而,目前針對森林土壤氧化鐵全量的光譜反演研究相對較少,且高光譜測量易受到各種環(huán)境因素的影響,在一定程度上增加了微量元素反演的難度。因此,實驗室的可見光到近紅外區(qū)域高光譜反射率用于預測土壤氧化鐵全量的潛力還需要有效深入地挖掘。

    本研究以湖南省大圍山為研究區(qū),利用室內高光譜數據,通過運用多種平滑和數學變換形式對原始光譜進行預處理,分析土壤氧化鐵全量的高光譜特征和特征波段,厘清多元逐步回歸模型和偏最小二乘模型估算土壤氧化鐵全量的精度差異,明確模型適用性以及最佳預處理方法,為快速定量獲取紅黃壤區(qū)森林土壤氧化鐵全量提供高光譜模型依據。

    1 材料與方法

    1.1 研究區(qū)概況

    大圍山(114°2′ ~ 114°12′ E,28°21′ ~ 28°26′ N)位于湖南省瀏陽市東北部大圍山鎮(zhèn)與張坊鎮(zhèn)交界處,與江西省銅鼓縣接壤,屬中亞熱帶季風濕潤氣候,年平均氣溫11 ℃ ~ 17 ℃,年平均降雨量1 200 ~ 2 000 mm,年相對濕度高于83%。研究區(qū)屬中山地貌,海拔范圍為230 ~ 1 607.9 m,相對高差1 377.9 m,地質單元主要為燕山期花崗巖(γ)和中元古代冷家溪群板頁巖(Ptln)[20],山地土壤和植被具有明顯垂直地帶性。

    1.2 樣品采集與測試

    本研究采集的樣品為大圍山花崗巖發(fā)育的紅壤、黃紅壤、黃壤、黃棕壤。在海拔177、474、545、720、860、1 078、1 210、1 382 m選擇了8個具有代表性的土壤剖面設置樣地,按土壤發(fā)生學劃分原采集各土壤剖面樣品共計 26 個,樣品經自然風干,檢出肉眼可見的沙礫及植物殘體,研磨過100目篩。采用四分法將每個樣品分為2份,1份采用三酸消解–原子吸收分光光度法測定土壤氧化鐵全量[21],另一份用于土壤光譜測量。

    在人工暗室中利用SOC-710(surface optics cor-poration,USA)成像高光譜儀對土壤樣品進行了高光譜圖像采集。儀器光譜范圍為400 ~ 1 050 nm,光譜分辨率為1.3 nm,128個波段。將土壤樣品平鋪至黑色塑料盒(直徑10 cm,深1 cm),水平放置在成像高光譜儀的視野中,灰板也放置于同視場角內,寬光譜石英鹵素燈垂直照向土壤樣品,與土樣表面距離約100 cm。高光譜儀鏡頭光圈設置為5.6 mm,調整焦距直至顯示出清晰土壤樣品圖像,曝光時間為自動,采集高光譜數據。通過SRAnal710e軟件將土壤光譜數據進行輻射定標并將能量值轉化成真實反射率,借助ENVI軟件在不同的鋁盒中分別選取靠近鋁盒中心的50個像素點作為ROI區(qū)域建立波譜庫。

    1.3 高光譜數據處理與特征波段提取

    將 26個土壤樣品分別依照土壤氧化鐵全量進行聚類分析,并計算出每一類土壤氧化鐵全量的光譜均值,獲得土壤氧化鐵全量的光譜響應范圍。將土壤氧化鐵全量分別與原始光譜反射率、數學變換光譜反射率和去包絡線光譜反射率進行逐波段的相關分析,根據相關系數最值對應波長選擇特征波段。相關系數計算公式如下:

    按照相關系數的顯著程度可以分成以下3類:|y|<0.4為低度線性相關,0.4<|y|<0.7為顯著性相關,0.7<|y|<1為高度線性相關[23]。

    1.4 光譜預測模型與驗證

    為保證預測樣本和驗證樣本的不同土壤類型的土壤氧化鐵全量變化范圍一致,將不同土壤類型的土壤,根據氧化鐵全量的不同,劃分為若干個小范圍,統(tǒng)計每個范圍內的樣本數后,隨機選取20個建模樣本,6 個驗證樣本?;谌グj線全波段(380 ~ 1 050 nm)高光譜數據建立多元線性回歸估測模型(MLR);按照建模組與驗證組比例分別為0.5:0.5和0.75:0.25,建立了兩種土壤偏最小二乘回歸估測模型(PLSR)。模型精度評價為決定系數2、校正均方根誤差RMSEC、預測均方根RMSEP,2越大,RMSEC越小,RMSEP與RMSEC越接近,表明模型的校正及預測精度越高。根據Rossel和Webster[24]對殘余預測偏差的分析可知,RPD<1.4說明模型的預測能力很差;1.4≤RPD<2.0,表明模型初步具備定量預測可能;RPD≥2.0表明模型具備較好的定量預測能力。本研究PLSR 建模和驗證利用Unscrambler?X10.4完成,多元線性回歸建模、數據分析在 SPSS 22.0 完成,繪圖通過GraphPad Prism 8完成。

    2 結果與分析

    2.1 土壤氧化鐵全量變化特征

    由于海拔高低相差大,氣候條件差異較大,大圍山土壤呈明顯垂直地帶分布特征,低于海拔500 m多為山地紅壤,海拔500 ~ 800 m為山地黃紅壤,海拔800 ~ 1 200 m為山地黃壤,山地黃棕壤則主要分布在海拔1 200 ~ 1 600 m[25]。由圖1和圖2可知,8個土壤剖面的氧化鐵全量均較高,均值為46.33 ~ 79.61 g/kg,表現為紅壤>黃紅壤>黃壤>黃棕壤,海拔與氧化鐵全量之間呈極顯著負相關關系(= –0.861 5,<0.001,= 26)。這與隨著海拔升高,溫度降低,相對濕度升高,土壤脫硅富鐵鋁化程度降低,氧化鐵全量隨之不斷減少有關。剖面土壤氧化鐵全量以剖面紅壤1最高,為73.01 ~ 88.56 g/kg,剖面黃棕壤2最低,僅為46.48 ~ 50.41 g/kg。所有剖面土壤氧化鐵全量均表現出AB層、B層和BC層較A層高的特征。這反映出研究區(qū)雨水豐富,淋溶作用強烈,氧化鐵具有從表土層向下聚集遷移的趨勢。隨著土壤剖面深度的增加,外界環(huán)境和淋溶作用減弱,土壤的氧化鐵全量逐步下降并趨于穩(wěn)定狀態(tài),這與前人對太白山、廬山土壤氧化鐵全量的研究結論類似[26-27]。

    圖1 各土壤剖面氧化鐵全量變化特征

    圖2 土壤氧化鐵全量與海拔的相關性

    2.2 土壤光譜曲線特征

    根據樣本的土壤氧化鐵全量,將樣本分為 >70、60 ~ 70、50 ~ 60和 <50 g/kg 4個區(qū)間,分別計算樣本在 4 個區(qū)間內的反射率平均值得到平均光譜曲線(圖3)。雖然樣品的氧化鐵全量不同,但高光譜平均反射率曲線均呈上凸的拋物線型。在430 ~ 560 nm,隨著氧化鐵全量的增加,反射率呈現下降的趨勢,氧化鐵和反射率為負相關關系;在500 ~ 620 nm波段,反射率曲線斜率增加最快,反射率以遞增的速率增加,形似陡坎;在620 ~ 780 nm 波段反射率曲線斜率略微變小,反射率以遞減的速率增加;在波長780 nm處,反射率出現最大值,隨后波段反射率開始下降;在波長900 nm附近,高光譜曲線存在一個微弱的吸收峰,被認為是Fe2+、Fe3+的微弱吸收峰[28]。這與土壤氧化鐵全量增加時,可見光與近紅外部分吸收增強,使土壤光譜反射率朝著藍波段方向下降,導致土壤光譜反射率降低有關[7]。此結果與何挺等[11]、彭杰等[12]、郭穎等[29]認為氧化鐵全量與反射光譜率為負相關一致。

    圖3 不同氧化鐵全量的土壤反射率均值光譜

    2.3 土壤氧化鐵全量的光譜特征波段

    2.4 土壤氧化鐵全量多元逐步回歸模型

    基于去包絡線全波段高光譜數據,建立多元線性回歸估測模型及誤差分析見表1。391、565、994和1040為對應391、565、994及1 040 nm波段反射率經去包絡線后的數值。通過比較2及 RPD值可以看出,模型 D 的決定系數(2)及相對分析誤差(RPD)達到最大,所以模型D建立的多元逐步回歸模型精度相對較高。

    2.5 土壤氧化鐵全量的偏最小二乘回歸模型

    將26個土壤樣本分成建模組與驗證組的比例分別為0.5:0.5(樣本形式A)和0.75:0.25(樣本形式B),基于去包絡線高光譜數據建立了2種偏最小二乘反演模型。由表2可知,2種樣本形式的 RPD 值均大于1.4,表明樣本形式A、B模型都能較好地估測土壤氧化鐵全量。雖然樣本形式A的建模組的決定系數R和校正均方根誤差 RMSEC與樣本形式B相接近,分別為0.930 6和3.643 4,但樣本形式A的驗證組決定系數2僅為0.630 0,驗證組預測均方根誤差RMSEP為7.740 0,表明該模型雖然建模效果好,但預測能力較差。樣本形式 B 的建模組的決定系數2達到0.935 8,校正均方根誤差RMSEC為3.419 7,驗證組的2達到0.981 8,預測均方根誤差RMSEP為3.155 0,表明該模型不僅建模效果好,而且預測能力也很強??梢?,樣本形式 B 的建模效果優(yōu)于A,且形式 B 的預測效果也更佳,能夠更好地估測土壤氧化鐵全量。

    2.6 模型的驗證與優(yōu)選

    從多元逐步回歸模型和偏最小二乘回歸模型建模結果對比可見(表3),偏最小二乘回歸模型樣本形式B精度更高,其驗證組的決定系數(2)為0.981 8,預測均方根誤差RMSEP為3.155 0,預測相對分析誤差RPD為4.235 2。由兩種模型預測值和實測值擬合圖(圖5)也可看出,偏最小二乘回歸模型B的數據組相比多元逐步回歸模型,建模組和驗證組的預測效果更好,該結論和表2的分析結果一致。

    圖4 高光譜變換的相關性曲線

    表 1 基于去包絡線的多元逐步回歸模型

    表 2 基于去包絡線的偏最小二乘模型精度

    表 3 建模結果對比

    (A.多元逐步回歸模型預測值和實測值擬合;B.偏最小二乘回歸模型預測值和實測值擬合)

    3 結論

    1)不同森林土壤氧化鐵全量和不同土壤類型光譜曲線在可見光–近紅外波波段(400 ~ 1 000 nm)內變化趨勢基本一致,均形似陡坎。在420 ~ 580 nm波段,氧化鐵全量與光譜反射率呈負相關關系。

    2)土壤原始反射率經過S-G平滑,并去包絡線和微分處理后,光譜曲線特征差異明顯,其中去包絡線的處理反演效果最好。土壤氧化鐵高光譜的特征波段主要為392、427、529、523、549、559、565、570、994和1 040 nm。

    3)以去包絡線的土壤光譜反射率為自變量建立的土壤氧化鐵全量多元逐步回歸模型=–115.015×565–275.853×994+231.346×1040–76.7665×391+227.603精度較好,建模組與驗證組的比例為0.75:0.25優(yōu)于比例為0.5:0.5的偏最小二乘回歸模型。其中,建模組與驗證組的比例為0.75:0.25的偏最小二乘反演模型的2最高,為0.935 8;RMSEC最小,為3.419 7;RPD最高,為4.235 2,是較為適用于紅黃壤區(qū)森林土壤氧化鐵全量估測的高光譜模型。

    [1] Schwertmann U, Cornell R M.Iron oxides in the laboratory: Preparation and characterization[M].Wiley‐VCH Verlag GmbH: 2000.

    [2] 章明奎.土壤氧化鐵的形態(tài)在浙江紅壤和黃壤分類中的意義[J].浙江農業(yè)大學學報, 1990, 16(1): 42–50.

    [3] Terra F S, Demattê J A M, Viscarra Rossel R A.Proximal spectral sensing in pedological assessments: Vis-NIR spectra for soil classification based on weathering and pedogenesis[J].Geoderma, 2018, 318: 123–136.

    [4] Todorova M, Atanassova S, Sitaula B, et al.Application of NIRS as a rapid and alternative method for prediction of heavy metals content in soil [J].Agricultural Science and Technology, 2012, 4(4): 440-444.

    [5] Gupta A, Vasava H B, Das B S, et al.Local modeling approaches for estimating soil properties in selected Indian soils using diffuse reflectance data over visible to near-infrared region[J].Geoderma, 2018, 325: 59–71.

    [6] 史舟, 王乾龍, 彭杰, 等.中國主要土壤高光譜反射特性分類與有機質光譜預測模型[J].中國科學: 地球科學, 2014, 44(5): 978–988.

    [7] Fontes M P F, Carvalho I A Jr.Color attributes and mineralogical characteristics, evaluated by radiometry, of highly weathered tropical soils[J].Soil Science Society of America Journal, 2005, 69(4): 1162–1172.

    [8] Stoner E R, Baumgardner M F.Characteristic variations in reflectance of surface soils [J].Soil Science Society of America Journal, 1981, 45(6): 1161–1165.

    [9] Hunt G R, Ashley R P.Spectra of altered rocks in the visible and near infrared [J].Economic Geology, 1979, 74(7): 1613–1629.

    [10] Galv?o L S, Pizarro M A, Epiphanio J C N.Variations in reflectance of tropical soils: Spectral-chemical composition relationships from AVIRIS data[J].Remote Sensing of Environment, 2001, 75(2): 245–255.

    [11] 何挺, 王靜, 程燁, 等.土壤氧化鐵光譜特征研究[J].地理與地理信息科學, 2006, 22(2): 30–34.

    [12] 彭杰, 向紅英, 周清, 等.土壤氧化鐵的高光譜響應研究[J].光譜學與光譜分析, 2013, 33(2): 502–506.

    [13] de Aldana B R V, Criado B G, Ciudad A G, et al.Estimation of mineral content in natural grasslands by near infrared reflectance spectroscopy[J].Communications in Soil Science and Plant Analysis, 1995, 26(9/10): 1383–1396.

    [14] 魏昌龍, 趙玉國, 鄔登巍, 等.基于光譜分析的土壤游離鐵預測研究[J].土壤, 2014, 46(4): 678–683.

    [15] Liu Y L, Chen Y Y.Estimation of total iron content in floodplain soils using VNIR spectroscopy - a case study in the Le'an River floodplain, China[J].International Journal of Remote Sensing, 2012, 33(18): 5954–5972.

    [16] Magendran T, Sanjeevi S.Statistical and neural network analysis of hyperspectral radiometric data to characterise hematite of Singbhum iron ore belt, India[C]//SPIE Asia-Pacific Remote Sensing.Proc SPIE 9263, Multispectral, Hyperspectral, and Ultraspectral Remote Sensing Technology, Techniques and Applications V, Beijing, China.2014, 9263: 92630M.

    [17] 秦倩如, 齊雁冰, 吳娟, 等.基于高光譜的土壤游離鐵隨機森林模型估算研究[J].土壤通報, 2018, 49(6): 1286–1293.

    [18] 謝文, 趙小敏, 郭熙, 等.基于組合模型的廬山森林土壤有效鐵光譜間接反演研究[J].土壤學報, 2017, 54(3): 601–612.

    [19] 高開強, 徐泮林, 趙曉旭.基于模糊模式識別的土壤含鐵量高光譜估測[J].測繪與空間地理信息, 2020, 43(1): 127–129,132.

    [20] 趙振華, 劉偉, 趙亞輝, 等.大圍山第四紀冰川地質遺跡特征及成因探討[J].國土資源導刊, 2008, 5(6): 42–45.

    [21] 魯如坤.土壤農業(yè)化學分析方法[M].北京: 中國農業(yè)科技出版社, 1999.

    [22] 王濤, 喻彩麗, 張楠楠, 等.基于去包絡線和連續(xù)投影算法的棗園土壤電導率光譜檢測研究[J].干旱地區(qū)農業(yè)研究, 2019, 37(5): 193–199,217.

    [23] 魏建榕.經濟物理學中的金融數據分析:統(tǒng)計與建模[D].上海: 復旦大學, 2013.

    [24] Rossel R A V, Webster R.Predicting soil properties from the Australian soil visible-near infrared spectroscopic database [J].European Journal of Soil Science, 2012, 63(6): 848–860.

    [25] 羅蘭花, 王翠紅, 謝紅霞, 等.大圍山花崗巖風化物發(fā)育土壤抗蝕性垂直分異[J].水土保持研究, 2018, 25(1): 62–65, 71.

    [26] 雷梅, 常慶瑞, 馮立孝, 等.太白山土壤特性及氧化鐵發(fā)生學特征[J].地理研究, 2001, 20(1): 83–90.

    [27] 曹慶.廬山北坡土壤發(fā)生特性與系統(tǒng)分類[D].南京: 南京師范大學, 2012.

    [28] 龔紹琦, 王鑫, 沈潤平, 等.濱海鹽土重金屬含量高光譜遙感研究[J].遙感技術與應用, 2010, 25(2): 169–177.

    [29] 郭穎, 郭治興, 劉佳, 等.亞熱帶典型區(qū)域水稻土氧化鐵高光譜反演——以珠江三角洲為例[J].應用生態(tài)學報, 2017, 28(11): 3675–3683.

    [30] Baumgardner M F, Silva L F, Biehl L L, et al.Reflectance properties of soils [J].Advances in Agronomy, 1986, 38: 1–44.

    Inversion of Iron Oxide Contents in Forest Soils of Dawei Mountains Using Laboratory Hyperspectral Data

    TAN Jie1, CHEN Yan1, ZHOU Weijun2*, CUI Haojie2, LIU Pei1

    (1 College of Landscape Architecture and Art Design, Hunan Agricultural University, Changsha 410128, China; 2College of Resources and Environment, Hunan Agricultural University, Changsha 410128, China)

    Iron oxide is the main body of iron-bearing minerals in the soil and is one of the most obvious and useful indicators of soil development and soil classification.In this paper, forest soils in Dawei Mountain of Hunan Province were collected, iron oxide contents in soils were determined respectively by conventional chemical method and by hyperspectral inversion with models of multiple stepwise regression and partial least squares regression inversion established by screening sensitive bands after spectral preprocessing and combinatorial transformation.The results showed that soil spectral curves with different iron oxide contents all were in steep-hill shape in the whole band, iron oxide content was negatively correlated with spectral reflectance within 420-580 nm, different spectral transformation could improve the correlation, and the combination of Savitzky-Golay (S-G) smoothing and de-embedding lines was superior to other pretreatment methods in inversion.The characteristic bands of iron oxides were 392, 427, 529, 523, 549, 559, 565, 570, 994 and 1 040 nm.Partial least squares regression model had better stability than multiple stepwise regression model, and is suitable for rapid estimation of iron oxide contents in forest red and yellow soils.

    Soil spectra; Iron oxide; Multiple stepwise regression; Partial least squares regression

    S155.1

    A

    10.13758/j.cnki.tr.2021.04.025

    譚潔, 陳嚴, 周衛(wèi)軍, 等.基于實驗室高光譜數據的大圍山森林土壤氧化鐵全量反演.土壤, 2021, 53(4): 858–864.

    國家自然科學基金項目(41771272)資助。

    (wjzh0108@163.com)

    譚潔(1979—),女,湖南湘潭人,博士,副教授,主要從事土地/土壤環(huán)境過程及模擬研究。E-mail: 37080689@qq.com

    猜你喜歡
    大圍山包絡線全量
    湖南:瀏陽市大圍山鎮(zhèn)水果年產值超億元
    平面曲線運動中的包絡線問題探討
    成都市溫江區(qū)全力推進醫(yī)保全量數據采集試點工作
    拋體的包絡線方程的推導
    大圍山之戀(外一首)
    民族音樂(2019年6期)2019-01-09 17:52:15
    風景秀麗的大圍山
    創(chuàng)作(2018年1期)2018-11-12 21:10:04
    一種用于故障隔離的參數區(qū)間包絡線計算方法
    現代機械(2018年1期)2018-04-17 07:29:35
    土壤重金屬污染調查與評估的誤區(qū)及其改進方法
    麥秸全量還田下氮肥運籌對水稻產量及其產量構成的影響
    瀏陽大圍山客家山歌研究
    戲劇之家(2015年24期)2016-01-12 00:36:26
    老司机午夜十八禁免费视频| a 毛片基地| 亚洲国产av影院在线观看| 一级,二级,三级黄色视频| 男女下面插进去视频免费观看| 考比视频在线观看| 国产不卡av网站在线观看| 99re6热这里在线精品视频| 老司机影院毛片| 久久天躁狠狠躁夜夜2o2o| 日本撒尿小便嘘嘘汇集6| 亚洲综合色网址| 香蕉国产在线看| 欧美 日韩 精品 国产| 亚洲精品av麻豆狂野| 国产片内射在线| 日韩制服骚丝袜av| 久久久水蜜桃国产精品网| 搡老岳熟女国产| 人人妻人人爽人人添夜夜欢视频| 国产伦理片在线播放av一区| 亚洲成人国产一区在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 精品久久蜜臀av无| 少妇的丰满在线观看| 亚洲国产中文字幕在线视频| 精品一区二区三区av网在线观看 | 午夜91福利影院| 国产高清视频在线播放一区 | 精品久久久精品久久久| 无遮挡黄片免费观看| 美女脱内裤让男人舔精品视频| 一进一出抽搐动态| 亚洲av成人一区二区三| 亚洲精品成人av观看孕妇| 免费人妻精品一区二区三区视频| 国产精品一区二区在线观看99| 视频区欧美日本亚洲| 久久女婷五月综合色啪小说| 天堂中文最新版在线下载| 亚洲精品久久久久久婷婷小说| 美女主播在线视频| 麻豆国产av国片精品| 亚洲人成电影观看| 在线观看免费日韩欧美大片| 18在线观看网站| 国产极品粉嫩免费观看在线| 亚洲精品乱久久久久久| 成人影院久久| 91精品国产国语对白视频| 老汉色∧v一级毛片| 男女床上黄色一级片免费看| 99精品久久久久人妻精品| 我的亚洲天堂| 18禁观看日本| 国产成人精品久久二区二区免费| 老汉色∧v一级毛片| 老汉色av国产亚洲站长工具| 一区二区av电影网| 老司机深夜福利视频在线观看 | 免费人妻精品一区二区三区视频| 亚洲色图综合在线观看| 十八禁高潮呻吟视频| 老司机影院毛片| 极品人妻少妇av视频| 精品少妇内射三级| 亚洲 国产 在线| 成人亚洲精品一区在线观看| 美女主播在线视频| 天天躁日日躁夜夜躁夜夜| 人人妻,人人澡人人爽秒播| 日韩视频一区二区在线观看| 国产成人欧美在线观看 | 国产1区2区3区精品| 少妇粗大呻吟视频| 久久久国产精品麻豆| 日本一区二区免费在线视频| 免费在线观看视频国产中文字幕亚洲 | 国产黄频视频在线观看| 黄色毛片三级朝国网站| 免费黄频网站在线观看国产| 男女之事视频高清在线观看| 国产麻豆69| 各种免费的搞黄视频| 狂野欧美激情性xxxx| 国产激情久久老熟女| 久久精品国产综合久久久| 精品亚洲成国产av| 欧美激情极品国产一区二区三区| 1024香蕉在线观看| 久久av网站| 国产av精品麻豆| 人成视频在线观看免费观看| www.精华液| 精品少妇一区二区三区视频日本电影| 涩涩av久久男人的天堂| 亚洲 国产 在线| 宅男免费午夜| 中文欧美无线码| 国产日韩一区二区三区精品不卡| 一区二区三区乱码不卡18| 亚洲一区中文字幕在线| 亚洲精品成人av观看孕妇| 老司机靠b影院| 在线永久观看黄色视频| 性色av乱码一区二区三区2| 99re6热这里在线精品视频| av电影中文网址| 久久精品亚洲熟妇少妇任你| 在线观看舔阴道视频| 激情视频va一区二区三区| 欧美精品高潮呻吟av久久| 中文字幕人妻丝袜一区二区| 伊人亚洲综合成人网| 高清av免费在线| 亚洲专区国产一区二区| 十八禁网站网址无遮挡| bbb黄色大片| 搡老熟女国产l中国老女人| 啦啦啦免费观看视频1| 久久久久久久久久久久大奶| 久久青草综合色| 午夜视频精品福利| 性色av乱码一区二区三区2| √禁漫天堂资源中文www| 国产视频一区二区在线看| 人人妻人人爽人人添夜夜欢视频| 亚洲免费av在线视频| 亚洲国产精品成人久久小说| 亚洲精品美女久久久久99蜜臀| 后天国语完整版免费观看| 在线观看免费高清a一片| 精品亚洲乱码少妇综合久久| 中文字幕人妻丝袜一区二区| 国产高清视频在线播放一区 | 国产在线一区二区三区精| 国产亚洲av片在线观看秒播厂| 午夜福利,免费看| 国产伦人伦偷精品视频| 在线观看免费日韩欧美大片| 日本撒尿小便嘘嘘汇集6| 亚洲情色 制服丝袜| 欧美激情久久久久久爽电影 | 不卡一级毛片| 51午夜福利影视在线观看| 丰满人妻熟妇乱又伦精品不卡| 老汉色∧v一级毛片| 天天添夜夜摸| 日本猛色少妇xxxxx猛交久久| 亚洲精品在线美女| 黑丝袜美女国产一区| 色综合欧美亚洲国产小说| 亚洲中文字幕日韩| 午夜福利免费观看在线| 精品一区在线观看国产| 国产亚洲一区二区精品| 国产免费av片在线观看野外av| 国产精品av久久久久免费| tocl精华| 在线av久久热| 亚洲精品国产色婷婷电影| 亚洲专区字幕在线| 亚洲久久久国产精品| 美女主播在线视频| 九色亚洲精品在线播放| 国产精品自产拍在线观看55亚洲 | 麻豆国产av国片精品| 91九色精品人成在线观看| 久久亚洲精品不卡| 亚洲av男天堂| 亚洲精品国产区一区二| 亚洲欧美精品自产自拍| 精品少妇久久久久久888优播| 久久精品国产a三级三级三级| 男女下面插进去视频免费观看| 我的亚洲天堂| 久久亚洲精品不卡| e午夜精品久久久久久久| 韩国高清视频一区二区三区| 亚洲人成77777在线视频| 美女扒开内裤让男人捅视频| 成年av动漫网址| 国产免费福利视频在线观看| 99国产精品99久久久久| 精品视频人人做人人爽| 日本wwww免费看| a级毛片黄视频| 日本精品一区二区三区蜜桃| 熟女少妇亚洲综合色aaa.| 韩国高清视频一区二区三区| 国产av精品麻豆| 高清av免费在线| 12—13女人毛片做爰片一| 欧美xxⅹ黑人| 美女福利国产在线| 日韩欧美一区二区三区在线观看 | av有码第一页| 青春草亚洲视频在线观看| 嫩草影视91久久| 国产伦人伦偷精品视频| 国产成人欧美在线观看 | 90打野战视频偷拍视频| 日本av免费视频播放| 亚洲国产欧美网| 80岁老熟妇乱子伦牲交| 交换朋友夫妻互换小说| 丝袜人妻中文字幕| 国产亚洲午夜精品一区二区久久| 亚洲精品中文字幕在线视频| 久久精品国产亚洲av高清一级| 亚洲欧美日韩高清在线视频 | 男女之事视频高清在线观看| 亚洲欧洲日产国产| 国产av国产精品国产| 少妇 在线观看| 亚洲精品国产av蜜桃| 国产欧美亚洲国产| 别揉我奶头~嗯~啊~动态视频 | 人妻一区二区av| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲一区中文字幕在线| 国产三级黄色录像| 淫妇啪啪啪对白视频 | 一区二区av电影网| 啦啦啦啦在线视频资源| 日本a在线网址| 亚洲 欧美一区二区三区| 亚洲第一av免费看| 国产熟女午夜一区二区三区| 99热全是精品| 最近中文字幕2019免费版| 国产精品一区二区在线不卡| 日韩精品免费视频一区二区三区| 日本撒尿小便嘘嘘汇集6| av不卡在线播放| a 毛片基地| 老司机亚洲免费影院| 免费黄频网站在线观看国产| 亚洲性夜色夜夜综合| 少妇人妻久久综合中文| 一本久久精品| 国产又爽黄色视频| 亚洲精品中文字幕一二三四区 | 777米奇影视久久| 国产成人欧美| 丝袜人妻中文字幕| 水蜜桃什么品种好| 国产免费福利视频在线观看| 国产亚洲精品一区二区www | 精品乱码久久久久久99久播| 亚洲第一欧美日韩一区二区三区 | 日韩视频在线欧美| 精品国产乱码久久久久久小说| 欧美日韩精品网址| 国产精品自产拍在线观看55亚洲 | av天堂久久9| 亚洲国产毛片av蜜桃av| 久久这里只有精品19| 国产视频一区二区在线看| 免费观看av网站的网址| 99国产精品99久久久久| 国产一区二区三区综合在线观看| 亚洲精品国产色婷婷电影| 一个人免费看片子| 欧美日韩一级在线毛片| 999久久久精品免费观看国产| 交换朋友夫妻互换小说| 91九色精品人成在线观看| 一二三四社区在线视频社区8| 国产成人啪精品午夜网站| 精品人妻1区二区| 亚洲熟女精品中文字幕| 啦啦啦啦在线视频资源| 欧美亚洲 丝袜 人妻 在线| 欧美性长视频在线观看| 91av网站免费观看| 亚洲欧美一区二区三区久久| 五月开心婷婷网| 久久综合国产亚洲精品| 精品人妻一区二区三区麻豆| 国产在线视频一区二区| 少妇被粗大的猛进出69影院| 欧美老熟妇乱子伦牲交| 亚洲精品自拍成人| 男女国产视频网站| 国产1区2区3区精品| 久久人人97超碰香蕉20202| 婷婷成人精品国产| av不卡在线播放| 天天影视国产精品| 亚洲av欧美aⅴ国产| 国产精品一区二区在线不卡| 免费在线观看影片大全网站| 女性生殖器流出的白浆| 女人精品久久久久毛片| 欧美日韩av久久| 国产三级黄色录像| 啦啦啦在线免费观看视频4| 久久天堂一区二区三区四区| svipshipincom国产片| 欧美精品高潮呻吟av久久| 久久久精品国产亚洲av高清涩受| 久久久久久人人人人人| 如日韩欧美国产精品一区二区三区| 国产精品久久久久久精品电影小说| 精品国产乱码久久久久久男人| 午夜免费鲁丝| 99国产精品免费福利视频| 久久久精品94久久精品| 另类亚洲欧美激情| 国产男女内射视频| 成年动漫av网址| 18禁裸乳无遮挡动漫免费视频| 黑人巨大精品欧美一区二区蜜桃| 一区二区av电影网| 人人妻,人人澡人人爽秒播| 青草久久国产| 亚洲自偷自拍图片 自拍| 国产高清视频在线播放一区 | 精品久久久精品久久久| 欧美变态另类bdsm刘玥| 99热国产这里只有精品6| 五月天丁香电影| 亚洲精品国产精品久久久不卡| 亚洲欧洲日产国产| 精品亚洲成a人片在线观看| 久久天躁狠狠躁夜夜2o2o| 国产成人精品久久二区二区免费| 久久国产精品男人的天堂亚洲| 最近最新免费中文字幕在线| 精品少妇内射三级| 高清视频免费观看一区二区| av网站在线播放免费| 久久人人97超碰香蕉20202| 久久精品成人免费网站| 国产成+人综合+亚洲专区| 日韩三级视频一区二区三区| 成年美女黄网站色视频大全免费| 免费在线观看日本一区| 午夜免费观看性视频| 亚洲欧美成人综合另类久久久| 一级黄色大片毛片| 大陆偷拍与自拍| 精品第一国产精品| 日本欧美视频一区| 欧美日韩一级在线毛片| 韩国高清视频一区二区三区| 色94色欧美一区二区| 欧美精品av麻豆av| 欧美精品高潮呻吟av久久| 久久香蕉激情| 侵犯人妻中文字幕一二三四区| 国产一区二区三区在线臀色熟女 | 国产1区2区3区精品| 国产深夜福利视频在线观看| 最近最新免费中文字幕在线| 美女福利国产在线| 在线观看免费视频网站a站| 美女主播在线视频| 亚洲午夜精品一区,二区,三区| www.av在线官网国产| 亚洲欧洲日产国产| 亚洲视频免费观看视频| 91老司机精品| 亚洲国产日韩一区二区| 国产成人免费无遮挡视频| 老汉色av国产亚洲站长工具| 精品少妇久久久久久888优播| www.999成人在线观看| 久久久久久久精品精品| 日韩熟女老妇一区二区性免费视频| 91精品国产国语对白视频| 国产亚洲欧美在线一区二区| 男男h啪啪无遮挡| 日韩熟女老妇一区二区性免费视频| 久久国产亚洲av麻豆专区| 性高湖久久久久久久久免费观看| 亚洲av片天天在线观看| 黄频高清免费视频| 丝袜脚勾引网站| 最近最新免费中文字幕在线| 一区二区三区激情视频| 男人操女人黄网站| 成年av动漫网址| 日本wwww免费看| 97在线人人人人妻| 精品久久蜜臀av无| 精品亚洲乱码少妇综合久久| 国产黄频视频在线观看| 欧美国产精品va在线观看不卡| 韩国高清视频一区二区三区| 国产亚洲av高清不卡| 女性生殖器流出的白浆| 欧美另类亚洲清纯唯美| 欧美 日韩 精品 国产| 交换朋友夫妻互换小说| 91精品伊人久久大香线蕉| 国产不卡av网站在线观看| 亚洲免费av在线视频| 婷婷色av中文字幕| 一级片免费观看大全| a 毛片基地| 中国美女看黄片| 久久久久国产精品人妻一区二区| 国产一区二区三区av在线| 亚洲精品国产区一区二| 亚洲国产av新网站| 国产xxxxx性猛交| 多毛熟女@视频| 亚洲成人国产一区在线观看| 日日摸夜夜添夜夜添小说| 90打野战视频偷拍视频| 久久精品亚洲熟妇少妇任你| 亚洲激情五月婷婷啪啪| 婷婷色av中文字幕| 亚洲伊人久久精品综合| 日韩 欧美 亚洲 中文字幕| 欧美中文综合在线视频| 亚洲av美国av| 少妇猛男粗大的猛烈进出视频| 90打野战视频偷拍视频| 精品久久久精品久久久| av网站免费在线观看视频| 97人妻天天添夜夜摸| 男女边摸边吃奶| 激情视频va一区二区三区| 99香蕉大伊视频| 嫩草影视91久久| 一本色道久久久久久精品综合| 日韩电影二区| 久久久久国产一级毛片高清牌| 啦啦啦免费观看视频1| 欧美少妇被猛烈插入视频| 欧美激情 高清一区二区三区| 精品一区二区三区av网在线观看 | 在线精品无人区一区二区三| 中文字幕高清在线视频| 国产福利在线免费观看视频| 99热全是精品| 日本一区二区免费在线视频| 在线观看一区二区三区激情| 久久精品亚洲av国产电影网| 亚洲欧美精品综合一区二区三区| tocl精华| 精品亚洲乱码少妇综合久久| 国产男女内射视频| 99re6热这里在线精品视频| 亚洲中文日韩欧美视频| 日本a在线网址| 侵犯人妻中文字幕一二三四区| 国产av又大| 自线自在国产av| 女人精品久久久久毛片| 久久精品国产亚洲av香蕉五月 | 久久这里只有精品19| 亚洲欧美色中文字幕在线| 法律面前人人平等表现在哪些方面 | 99热国产这里只有精品6| 精品亚洲成国产av| 丰满人妻熟妇乱又伦精品不卡| 十八禁网站免费在线| 黑人操中国人逼视频| 国产免费福利视频在线观看| 日韩一区二区三区影片| 久久香蕉激情| videos熟女内射| √禁漫天堂资源中文www| 亚洲成av片中文字幕在线观看| 国产成人免费无遮挡视频| 国产伦人伦偷精品视频| 亚洲熟女毛片儿| 一进一出抽搐动态| 日韩视频在线欧美| 黄片播放在线免费| 久久中文字幕一级| 欧美日韩精品网址| 十八禁网站网址无遮挡| 一区二区三区乱码不卡18| 美女大奶头黄色视频| 又紧又爽又黄一区二区| 日本一区二区免费在线视频| 精品国产一区二区三区四区第35| 黄色a级毛片大全视频| 大陆偷拍与自拍| 男人爽女人下面视频在线观看| 成人手机av| 97在线人人人人妻| 欧美精品av麻豆av| 女性生殖器流出的白浆| 美女国产高潮福利片在线看| 青草久久国产| 午夜福利视频在线观看免费| 伦理电影免费视频| 一个人免费在线观看的高清视频 | 免费观看人在逋| 99热全是精品| 日韩一卡2卡3卡4卡2021年| h视频一区二区三区| 亚洲欧美日韩高清在线视频 | 女人精品久久久久毛片| 国产91精品成人一区二区三区 | 伊人久久大香线蕉亚洲五| 大片电影免费在线观看免费| 精品少妇黑人巨大在线播放| 欧美人与性动交α欧美软件| 久久天躁狠狠躁夜夜2o2o| 每晚都被弄得嗷嗷叫到高潮| 精品久久蜜臀av无| 十八禁网站免费在线| 国产在视频线精品| 日韩大码丰满熟妇| 久久人人97超碰香蕉20202| 成人三级做爰电影| 老汉色av国产亚洲站长工具| 一本大道久久a久久精品| 亚洲五月婷婷丁香| 在线观看人妻少妇| 色综合欧美亚洲国产小说| 欧美97在线视频| 日韩一卡2卡3卡4卡2021年| 丝袜喷水一区| 老司机午夜福利在线观看视频 | 老汉色∧v一级毛片| 一区在线观看完整版| 99热全是精品| 亚洲美女黄色视频免费看| 婷婷成人精品国产| 中文字幕色久视频| 亚洲久久久国产精品| 又黄又粗又硬又大视频| tocl精华| 美女中出高潮动态图| 淫妇啪啪啪对白视频 | 精品卡一卡二卡四卡免费| 一级,二级,三级黄色视频| 黑人巨大精品欧美一区二区蜜桃| 国产片内射在线| 国产精品自产拍在线观看55亚洲 | 黄色视频在线播放观看不卡| 最黄视频免费看| av欧美777| 99国产精品99久久久久| 中国美女看黄片| 乱人伦中国视频| 国产片内射在线| 国产精品一区二区在线观看99| 亚洲精品第二区| 日韩制服骚丝袜av| 欧美变态另类bdsm刘玥| 咕卡用的链子| 青草久久国产| 精品国产一区二区三区四区第35| 啦啦啦视频在线资源免费观看| 亚洲国产欧美一区二区综合| 国产在线一区二区三区精| 久久影院123| 丝袜美腿诱惑在线| 蜜桃国产av成人99| 国产男女内射视频| 午夜福利免费观看在线| 国产精品国产三级国产专区5o| 一本色道久久久久久精品综合| 十八禁网站免费在线| 亚洲国产欧美网| 亚洲精品中文字幕一二三四区 | 又紧又爽又黄一区二区| 亚洲五月婷婷丁香| 国产高清视频在线播放一区 | 亚洲国产欧美日韩在线播放| av福利片在线| 久久精品国产a三级三级三级| 久久久久视频综合| 人妻 亚洲 视频| 国产日韩一区二区三区精品不卡| 捣出白浆h1v1| 中文字幕高清在线视频| 欧美精品一区二区大全| 亚洲精品一区蜜桃| 天堂8中文在线网| 天天操日日干夜夜撸| 国产精品影院久久| 精品一区二区三区av网在线观看 | 久久香蕉激情| 午夜福利一区二区在线看| 国产免费一区二区三区四区乱码| 亚洲五月婷婷丁香| 久久狼人影院| 免费在线观看日本一区| 久久天躁狠狠躁夜夜2o2o| 午夜福利视频在线观看免费| 免费在线观看日本一区| 男女边摸边吃奶| 欧美av亚洲av综合av国产av| 中文字幕av电影在线播放| 亚洲av国产av综合av卡| 天天躁狠狠躁夜夜躁狠狠躁| xxxhd国产人妻xxx| 亚洲激情五月婷婷啪啪| 国产精品麻豆人妻色哟哟久久| 免费在线观看日本一区| 1024视频免费在线观看| 精品亚洲成国产av| 日日摸夜夜添夜夜添小说| 十八禁人妻一区二区| 久久国产精品男人的天堂亚洲| 丝袜美足系列| 久久九九热精品免费| 男女无遮挡免费网站观看| 久久精品国产综合久久久| 久9热在线精品视频| 欧美变态另类bdsm刘玥| 超碰成人久久| 久久这里只有精品19| 精品少妇黑人巨大在线播放| 亚洲专区国产一区二区|