李恒凯, 雷军, 吴娇. 基于多源时序NDVI的稀土矿区土地毁损与恢复过程分析[J]. 农业工程学报, 2018, 34(1): 232-240. DOI: 10.11975/j.issn.1002-6819.2018.01.032
    引用本文: 李恒凯, 雷军, 吴娇. 基于多源时序NDVI的稀土矿区土地毁损与恢复过程分析[J]. 农业工程学报, 2018, 34(1): 232-240. DOI: 10.11975/j.issn.1002-6819.2018.01.032
    Li Hengkai, Lei Jun, Wu Jiao. Analysis of land damage and recovery process in rare earth mining area based on multi-source sequential NDVI[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(1): 232-240. DOI: 10.11975/j.issn.1002-6819.2018.01.032
    Citation: Li Hengkai, Lei Jun, Wu Jiao. Analysis of land damage and recovery process in rare earth mining area based on multi-source sequential NDVI[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(1): 232-240. DOI: 10.11975/j.issn.1002-6819.2018.01.032

    基于多源时序NDVI的稀土矿区土地毁损与恢复过程分析

    Analysis of land damage and recovery process in rare earth mining area based on multi-source sequential NDVI

    • 摘要: 为分析稀土开采的时空分布及矿区土地毁损与恢复过程,该文以岭北稀土矿区为例,以1990-2016年的HJ-1B CCD、Landsat 5和Landsat 8遥感数据为数据源,结合回归分析法、遥感时序NDVI分析方法,对岭北稀土矿区的稀土开采状况及土地毁损与恢复情况进行分析。为减少不同数据由于传感器自身原因而带来的NDVI误差,采用回归分析法构建HJ-1B CCD、Landsat5/8数据的NDVI转换方程并利用均方根误差对转换方程的精度进行检验。结果表明:HJ-1B CCD、Landsat5/8数据的NDVI 影像对应转换方程的R2值均超过0.9以及模拟影像与真实影像之间的均方根误差值均小于0.05,说明HJ-1B CCD、Landsat5/8的NDVI数据之间存在较为显著的线性正相关,所求的转换方程具有较高精度,可以消除数据来源差异对NDVI的影响;稀土开采时空分布分析表明岭北稀土矿区在2001-2004年、2006年的开采面积较大,均超过1 km2,其中2006的开采面积最大,达到2.546 1 km2;稀土开采在空间分布上较为分散,不仅增加了监管的困难,造成了资源的浪费,还一定程度上增加了治理困难;土地毁损与恢复分析表明受到人为扰动的区域占矿区总面积的54.59%,其中受森林砍伐扰动影响的面积最大,为86.5341 km2;稀土开采扰动的土地面积为11.354 4 km2,其植被恢复的平均加权时长为11年。其中复垦恢复的面积为5.004 9 km2,仍有6.349 5 km2的区域植被未恢复,急需治理。该研究能较好的反映稀土干扰下的土地毁损与恢复过程,为岭北稀土矿区的生态环境治理及矿区的可持续发展提供借鉴。

       

      Abstract: Abstract: The problem of land and ecological environment damage caused by ion mining has attracted more and more attention. Taking the Lingbei rare earth mining area as an example, combined with remote sensing normalized difference vegetation index (NDVI) time sequence analysis method, we used HJ-1B CCD (charge-coupled device), Landsat 5 and Landsat 8 remote sensing data during 1990-2016, and analyzed land distribution and rare earth mining area damage and recovery process. In order to reduce the NDVI error caused by the sensor itself, NDVI conversion equation of HJ-1B CCD and Landsat 5/8 data was constructed by using the regression analysis method, NDVI image of HJ-1B CCD was transformed into the NDVI image corresponding to the years of Landsat 5/8 (the analog image) by using conversion equation, and the root mean square error (RMSE) between the simulated image and the real Landsat 5/8 NDVI image was calculated and the accuracy of the transformation equation was verified. The results showed that the R2 values of NDVI equations from HJ-1B CCD and Landsat 5, HJ-1B CCD and Landsat 8 data were 0.9319 and 0.9552, both more than 0.9, indicating that there is obvious linear positive correlation between the NDVI of HJ-1B CCD and Landsat 5/8. The RMSE values of NDVI imaging between simulation images and real Landsat 5, and simulation images and real Landsat 8 were 0.0293 and 0.040, which were less than 0.05, showing that the obtained conversion equation has high accuracy and can eliminate the influence of different data sources on the NDVI. Combining NDVI threshold of mining and non-mining, and threshold of non-rare earth mining interference and rare earth mining interference, Lingbei rare earth mining area had annual rare earth mining activity from 1990 to 2016, rare earth mining area of which was smaller in 1990, 1991,1995 and 2000, mining area was more than 1 km2 in 2001-2006, and the mining area reached 2.546 1 km2 in 2006; after 2006, to maintain the rare earth mining scale at the lower level, the country began to implement control of rare earth mining, and the spatial distribution of rare earth mining was more dispersed, which not only caused difficulties in supervision and waste of resources, but also increased governance difficulties to a certain extent. Testing the extraction accuracy of land destruction and restoration type, extraction accuracy of various types was found above 85%, most of which reached more than 90%, satisfying the requirement of accuracy. Analysis of land destruction and restoration type showed that the area that was not disturbed was 97.082 1 km2, accounting for 45.41% of the total mining area; anthropogenic disturbance area was 116.709 3 km2, accounting for 54.59% of the total mining area; the area affected by deforestation disturbance was the largest, which was 86.534 1 km2, and the average time length of vegetation restoration was 3 years; the arable land that was not abandoned occupied 5.404 5 km2; cultivated land abandoned for many years and that for returning to forest occupied 9.390 6 km2; orchard development area was 3.148 2 km2, and the average time length of vegetation restoration for 7 years; rare earth mining occupied 11.354 4 km2, whose average recovery time was 11 years, and the recovery of the reclamation area was 5.004 9 km2, so there was still 6.349 5 km2 area with the regional vegetation that did not be recovered, and relevant departments should pay attention to the situation.

       

    /

    返回文章
    返回