[go: up one dir, main page]

CN114724029A - Remote sensing extraction method and system of mining goaf's influence on the surface - Google Patents

Remote sensing extraction method and system of mining goaf's influence on the surface Download PDF

Info

Publication number
CN114724029A
CN114724029A CN202210351884.XA CN202210351884A CN114724029A CN 114724029 A CN114724029 A CN 114724029A CN 202210351884 A CN202210351884 A CN 202210351884A CN 114724029 A CN114724029 A CN 114724029A
Authority
CN
China
Prior art keywords
area
index
goaf
remote sensing
crop
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210351884.XA
Other languages
Chinese (zh)
Inventor
杨邦会
王晓华
王春红
殷健
温莹莹
胡乔利
李冬阳
王琦文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongke Haihui Tianjin Technology Co ltd
Original Assignee
Zhongke Haihui Tianjin Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhongke Haihui Tianjin Technology Co ltd filed Critical Zhongke Haihui Tianjin Technology Co ltd
Priority to CN202210351884.XA priority Critical patent/CN114724029A/en
Publication of CN114724029A publication Critical patent/CN114724029A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a remote sensing extraction method and a system for the influence of a goaf of a mining area on the earth's surface, wherein the method comprises the following steps: determining a region to be monitored of a mining area, and acquiring remote sensing data of the region; acquiring a normalized vegetation index mean value of a time sequence of a region to be monitored and a crop identification index based on remote sensing data, and confirming that the region to be monitored is a farmland region through a first preset condition based on the crop identification index and the normalized vegetation index mean value of the time sequence; acquiring a temperature vegetation drought index and crop nitrogen content of the farmland region, preliminarily determining the farmland region as a region influenced by the goaf through a second preset condition based on the temperature vegetation drought index and the crop nitrogen content, and further finally determining the farmland region as the region influenced by the goaf. According to the method, the remote sensing extraction of the earth surface can be performed aiming at the goaf of the mining area through the remote sensing image, so that whether the area corresponding to the remote sensing data is the area influenced by the goaf or not is judged, and the purpose of accurately identifying the crop extraction under the influence of the goaf is achieved.

Description

矿区采空区对地表影响的遥感提取方法及系统Remote sensing extraction method and system of mining goaf's influence on the surface

技术领域technical field

本发明涉及农业和遥感技术领域,尤其涉及一种矿区采空区对地表影响的遥感提取方法及系统。The invention relates to the technical fields of agriculture and remote sensing, in particular to a remote sensing extraction method and system for the influence of a goaf in a mining area on the surface.

背景技术Background technique

采空区是由人为挖掘或者天然地质运动在地表下面产生的“空洞”,采空区的存在使得矿山的安全生产面临很大的安全问题,人员与机械设备都可能掉入采空区内部受到伤害。煤矿采空区:是指在煤矿作业过程中,将地下煤炭或煤矸石等开采完成后留下的空洞或空腔。The goaf is a "hole" created by artificial excavation or natural geological movement under the surface. The existence of the goaf makes the safety production of the mine face great safety problems. Personnel and mechanical equipment may fall into the goaf and be affected. harm. Coal mine goaf: refers to the cavity or cavity left after the mining of underground coal or coal gangue is completed during the operation of the coal mine.

矿区采空区如果治理不当,会导致地表沉降,直接影响是采空区上的房屋倾斜、采空区上的水体下泄、道路因为沉降也会导致塌陷、农田因为地下的镂空作用,导致地表不能保持水分导致干旱,并最终导致产量下降等问题,会发生系列地质灾害。因此,一般按照国家要求,需要绿色开采,一般要求边采矿、边治理。然而,采空区对农田作物的影响具有明显的隐蔽性,不容易被发现,成为目前采空区治理的难点。If the goaf in the mining area is not properly managed, it will lead to surface subsidence. The direct impact is that the houses on the gob are inclined, the water body on the gob is leaking, the road will also collapse due to subsidence, and the farmland will not be able to work on the surface due to the hollowing effect of the underground. A series of geological disasters can occur due to problems such as maintaining moisture, leading to droughts and, ultimately, reduced yields. Therefore, in accordance with national requirements, green mining is generally required, and mining and governance are generally required. However, the impact of goaf on farmland crops is obviously concealed, and it is not easy to be discovered, which has become a difficult point in the governance of goaf at present.

因此,如何通过遥感影像针对矿区采空区对地表的遥感提取,判断遥感数据对应的区域是否为采空区影响的区域,是当前课题亟需解决的技术问题。Therefore, how to use remote sensing images to extract the ground surface from the goaf in the mining area and determine whether the area corresponding to the remote sensing data is the area affected by the goaf is a technical problem that needs to be solved urgently.

发明内容SUMMARY OF THE INVENTION

本发明提供一种矿区采空区对地表影响的遥感提取方法及系统,用以解决现有技术中难以对采空区进行治理的缺陷,实现通过遥感影像针对矿区采空区对地表的遥感提取,判断遥感数据对应的区域是否为采空区影响的区域。The invention provides a remote sensing extraction method and system for the influence of a goaf in a mining area on the surface, so as to solve the defect that it is difficult to control the goaf in the prior art, and realize the remote sensing extraction of the surface of the goaf in the mining area through remote sensing images. , to judge whether the area corresponding to the remote sensing data is the area affected by the goaf.

本发明提供一种矿区采空区对地表影响的遥感提取方法,包括:The invention provides a remote sensing extraction method for the influence of a goaf in a mining area on the surface, comprising:

确定矿区的待监测区域,并获取待监测区域的遥感数据;Determine the area to be monitored in the mining area, and obtain remote sensing data of the area to be monitored;

在确认所述待监测区域为作物对降水敏感区的情况下,基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域;In the case of confirming that the to-be-monitored area is a crop sensitive area to precipitation, obtain the time-series normalized vegetation index mean value and crop identification index of the to-be-monitored area based on the remote sensing data, and based on the crop identification index and the mean value of the normalized vegetation index of the time series through the first preset condition to confirm that the area to be monitored is a farmland area;

获取所述农田区域的温度植被干旱指数以及作物氮素含量,并基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域。Obtain the temperature vegetation drought index and crop nitrogen content of the farmland area, and based on the temperature vegetation drought index and the crop nitrogen content, the farmland area is preliminarily confirmed as an area affected by a goaf through a second preset condition.

根据本发明提供的一种矿区采空区对地表影响的遥感提取方法,所述基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域,包括:According to a remote sensing extraction method for the influence of a goaf in a mining area on the surface provided by the present invention, the average value of the normalized vegetation index based on the crop identification index and the time series confirms that the to-be-monitored area is determined by a first preset condition. Farmland areas, including:

在所述作物识别指数不小于第一预设阈值,且所述时间序列的归一化植被指数均值不小于第二预设阈值的情况下,确认所述待监测区域为农田区域。When the crop identification index is not less than the first preset threshold, and the mean value of the normalized vegetation index of the time series is not less than the second preset threshold, it is confirmed that the area to be monitored is a farmland area.

根据本发明提供的一种矿区采空区对地表影响的遥感提取方法,所述基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域,包括:According to a remote sensing extraction method for the influence of a goaf in a mining area on the surface provided by the present invention, the farmland area is preliminarily confirmed as a goaf by a second preset condition based on the temperature, vegetation drought index and crop nitrogen content Affected areas include:

在所述温度植被干旱指数不小于第三预设阈值,且所述作物氮素含量不大于第四预设阈值的情况下,将所述农田区域初步确认为采空区影响的区域。Under the condition that the temperature vegetation drought index is not less than the third preset threshold, and the nitrogen content of the crop is not greater than the fourth preset threshold, the farmland area is preliminarily confirmed as the area affected by the goaf.

根据本发明提供的一种矿区采空区对地表影响的遥感提取方法,在所述温度植被干旱指数不小于第三预设阈值,且所述作物氮素含量不大于第四预设阈值的情况下,将所述农田区域初步确认为采空区影响的区域之后,所述方法还包括:According to a remote sensing extraction method for the influence of a goaf in a mining area on the surface provided by the present invention, when the temperature vegetation drought index is not less than a third preset threshold, and the nitrogen content of the crop is not greater than a fourth preset threshold Next, after the farmland area is preliminarily confirmed as the area affected by the goaf, the method further includes:

获取所述初步确认为采空区影响的区域内的时间序列的归一化植被指数均值以及作物生长周期,并获取矿区农田的归一化植被指数均值以及生长周期均值;Obtaining the mean value of the normalized vegetation index and the crop growth period of the time series in the area that is initially confirmed to be affected by the goaf, and obtaining the mean value of the normalized vegetation index and the mean value of the growth period of the farmland in the mining area;

在所述初步确认为采空区影响的区域内的时间序列的归一化植被指数均值与矿区农田的归一化植被指数均值的差值不小于第五预设阈值,且所述初步确认为采空区影响的区域内的作物生长周期不等于矿区农田的生长周期均值的情况下,将所述农田区域最终确认为采空区影响的区域。The difference between the mean value of the normalized vegetation index of the time series and the mean value of the normalized vegetation index of the farmland in the mining area in the area that is preliminarily confirmed to be affected by the goaf is not less than the fifth preset threshold, and the preliminary confirmation is If the crop growth cycle in the area affected by the goaf is not equal to the average growth cycle of the farmland in the mining area, the farmland area is finally confirmed as the area affected by the goaf.

根据本发明提供的一种矿区采空区对地表影响的遥感提取方法,所述基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,包括:According to a remote sensing extraction method for the impact of mined-out areas in mining areas on the surface provided by the present invention, the acquisition of the time-series normalized vegetation index mean value and crop identification index of the to-be-monitored area based on the remote sensing data includes:

在所述遥感数据的景数据中选择目标景数据,获取所述目标景数据的归一化植被指数以及时间序列的归一化植被指数均值;Selecting target scene data from the scene data of the remote sensing data, and obtaining the normalized vegetation index of the target scene data and the mean value of the normalized vegetation index of the time series;

基于所述目标景数据的归一化植被指数、时间序列的归一化植被指数均值以及遥感影像的景数据获取所述待监测区域的作物识别指数。The crop identification index of the to-be-monitored area is obtained based on the normalized vegetation index of the target scene data, the mean value of the normalized vegetation index of the time series, and the scene data of the remote sensing image.

根据本发明提供的一种矿区采空区对地表影响的遥感提取方法,确定矿区的待监测区域,并获取待监测区域的遥感数据之后,所述方法还包括:According to a remote sensing extraction method for the influence of a goaf in a mining area on the surface provided by the present invention, after determining the area to be monitored in the mining area and obtaining the remote sensing data of the area to be monitored, the method further includes:

获取所述待监测区域的年平均降水量以及作物种植季节降水量;Obtain the annual average precipitation of the area to be monitored and the precipitation in the crop planting season;

在所述年平均降水量不大于第六预设阈值,且所述作物种植季节降水量不大于第七预设阈值的情况下,将所述待监测区域判定为作物对降水敏感区。When the average annual precipitation is not greater than the sixth preset threshold, and the precipitation during the crop planting season is not greater than the seventh preset threshold, the to-be-monitored area is determined as a crop sensitive area to precipitation.

本发明还提供一种矿区采空区对地表影响的遥感提取系统,包括:The present invention also provides a remote sensing extraction system for the influence of the goaf in the mining area on the surface, including:

获取模块,用于确定矿区的待监测区域,并获取待监测区域的遥感数据;The acquisition module is used to determine the area to be monitored in the mining area and obtain the remote sensing data of the area to be monitored;

农田区域确认模块,用于在确认所述待监测区域为作物对降水敏感区的情况下,基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域;The farmland area confirmation module is used to obtain, based on the remote sensing data, the normalized vegetation index mean value and the crop identification index of the time series of the to-be-monitored area under the condition that the to-be-monitored area is confirmed as a crop sensitive area to precipitation, and confirming that the area to be monitored is a farmland area through a first preset condition based on the crop identification index and the mean value of the normalized vegetation index of the time series;

初步确认模块,用于获取所述农田区域的温度植被干旱指数以及作物氮素含量,并基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域。The preliminary confirmation module is used to obtain the temperature vegetation drought index and the crop nitrogen content of the farmland area, and based on the temperature vegetation drought index and the crop nitrogen content, the farmland area is preliminarily confirmed as harvesting through the second preset condition. The area affected by the empty zone.

本发明还提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述任一种所述矿区采空区对地表影响的遥感提取方法。The present invention also provides an electronic device, comprising a memory, a processor and a computer program stored in the memory and running on the processor, when the processor executes the program, the goaf of any one of the above-mentioned mining areas is realized Remote sensing extraction methods of regional impacts on the surface.

本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现如上述任一种所述矿区采空区对地表影响的遥感提取方法。The present invention also provides a non-transitory computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, implements the remote sensing extraction method for the influence of a goaf in a mining area on the surface as described above.

本发明还提供一种计算机程序产品,包括计算机程序,所述计算机程序被处理器执行时实现如上述任一种所述矿区采空区对地表影响的遥感提取方法。The present invention also provides a computer program product, including a computer program, which, when executed by a processor, implements the remote sensing extraction method for the impact of a goaf in a mining area on the surface as described above.

本发明提供的矿区采空区对地表影响的遥感提取方法及系统,通过确定矿区的待监测区域,并获取待监测区域的遥感数据,基于遥感数据获取待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认待监测区域为农田区域,获取农田区域的温度植被干旱指数以及作物氮素含量,并基于温度植被干旱指数和作物氮素含量通过第二预设条件将农田区域初步确认为采空区影响的区域。本发明能够通过遥感影像针对矿区采空区对地表的遥感提取,从而判断遥感数据对应的区域是否为采空区影响的区域,达到精准识别采空区影响下的作物提取的目的。The method and system for remote sensing extraction of the influence of mined-out areas in mining areas on the surface provided by the present invention, by determining the area to be monitored in the mining area, and obtaining the remote sensing data of the area to be monitored, based on the remote sensing data to obtain the normalized vegetation of the time series of the area to be monitored The average value of the index and the crop identification index, and based on the crop identification index and the mean value of the normalized vegetation index of the time series, the area to be monitored is confirmed as a farmland area through the first preset condition, and the temperature, vegetation drought index and crop nitrogen content of the farmland area are obtained, And based on the temperature vegetation drought index and crop nitrogen content, the farmland area is preliminarily confirmed as the area affected by the goaf through the second preset condition. The invention can use remote sensing images to extract the ground surface from the goaf in the mining area, so as to judge whether the area corresponding to the remote sensing data is the area affected by the goaf, and achieve the purpose of accurately identifying the crop extraction under the influence of the goaf.

附图说明Description of drawings

为了更清楚地说明本发明或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to explain the present invention or the technical solutions in the prior art more clearly, the following will briefly introduce the accompanying drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are the For some embodiments of the invention, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without any creative effort.

图1是本发明提供的矿区采空区对地表影响的遥感提取方法的流程示意图之一;Fig. 1 is one of the schematic flow charts of the remote sensing extraction method of the influence of mining goaf on the surface provided by the present invention;

图2是本发明提供的矿区采空区对地表影响的遥感提取方法的流程示意图之二;Fig. 2 is the second schematic flow chart of the remote sensing extraction method of the influence of mining goaf on the surface provided by the present invention;

图3是本发明提供的矿区采空区对地表影响的遥感提取系统的结构示意图;Fig. 3 is the structure schematic diagram of the remote sensing extraction system of the influence of the mining area goaf on the surface provided by the present invention;

图4是本发明提供的电子设备的结构示意图。FIG. 4 is a schematic structural diagram of an electronic device provided by the present invention.

具体实施方式Detailed ways

为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions in the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are part of the embodiments of the present invention. , not all examples. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

下面结合图1-图4描述本发明的矿区采空区对地表影响的遥感提取方法及系统。The method and system for remote sensing extraction of the influence of the mined goaf on the surface of the present invention will be described below with reference to FIGS. 1 to 4 .

参照图1,本发明提供一种矿区采空区对地表影响的遥感提取方法,包括以下步骤:Referring to Fig. 1, the present invention provides a remote sensing extraction method for the impact of a mining goaf on the surface, comprising the following steps:

步骤110、确定矿区的待监测区域,并获取待监测区域的遥感数据。Step 110: Determine the area to be monitored in the mining area, and acquire remote sensing data of the area to be monitored.

具体地,本实施例通过确定矿区的中心点坐标,通过专业地理信息软件做一定半径范围的缓冲区,并以此缓冲区为待监测区域。其中,覆盖矿区的区域为矿区主井、副井、矿区的煤仓生产区、生活区、公共设施以及运输管理等区域。确认待监测区域之后,获取待监测区域的遥感影像,并对遥感影像进行分析得到遥感数据,从而基于遥感数据对覆盖矿区的的待监测区域进行监测。Specifically, in this embodiment, the coordinates of the center point of the mining area are determined, and professional geographic information software is used to create a buffer area with a certain radius, and this buffer area is used as the area to be monitored. Among them, the areas covering the mining area are the main shaft, auxiliary shaft, coal bunker production area, living area, public facilities and transportation management areas of the mining area. After confirming the area to be monitored, the remote sensing image of the area to be monitored is obtained, and the remote sensing image is analyzed to obtain the remote sensing data, so as to monitor the area to be monitored covering the mining area based on the remote sensing data.

步骤120、在确认所述待监测区域为作物对降水敏感区的情况下,基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域。Step 120, in the case of confirming that the area to be monitored is a crop sensitive area to precipitation, obtain the mean value of the normalized vegetation index and the crop identification index of the time series of the area to be monitored based on the remote sensing data, and based on the The crop identification index and the mean value of the normalized vegetation index of the time series confirm that the to-be-monitored area is a farmland area through the first preset condition.

具体地,在覆盖矿区的待监测区域为作物对降水敏感区的基础上,根据遥感数据获取待监测区域的时间序列的归一化植被指数均值

Figure BDA0003580911680000061
和作物识别指数(CRI)作为判断指标,将待监测区域确认为农田区域,即该区域满足作物提取的首要目标。Specifically, on the basis that the to-be-monitored area covering the mining area is a crop sensitive to precipitation, the mean value of the normalized vegetation index of the time series of the to-be-monitored area is obtained according to the remote sensing data
Figure BDA0003580911680000061
The crop identification index (CRI) is used as the judgment index, and the area to be monitored is confirmed as the farmland area, that is, the area meets the primary goal of crop extraction.

步骤130、获取所述农田区域的温度植被干旱指数以及作物氮素含量,并基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域。Step 130: Obtain the temperature vegetation drought index and crop nitrogen content of the farmland area, and based on the temperature vegetation drought index and crop nitrogen content, preliminarily confirm that the farmland area is affected by the goaf through a second preset condition. Area.

具体地,本实施例在确认待监测区域为农田区域之后,获取该区域的温度植被干旱指数(TVDI)和作物氮素含量(TCARI/OSAVI)。当温度植被干旱指数和作物氮素含量满足预设条件的情况下,判定该区域为采空区影响的农田区域。Specifically, in this embodiment, after confirming that the area to be monitored is a farmland area, the temperature vegetation drought index (TVDI) and crop nitrogen content (TCARI/OSAVI) of the area are obtained. When the temperature vegetation drought index and crop nitrogen content meet the preset conditions, it is determined that the area is the farmland area affected by the goaf.

本发明提供的矿区采空区对地表影响的遥感提取方法,通过确定矿区的待监测区域,并获取待监测区域的遥感数据,基于遥感数据获取待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认待监测区域为农田区域,获取农田区域的温度植被干旱指数以及作物氮素含量,并基于温度植被干旱指数和作物氮素含量通过第二预设条件将农田区域初步确认为采空区影响的区域。本发明能够通过遥感影像针对矿区采空区对地表的遥感提取,从而判断遥感数据对应的区域是否为采空区影响的区域,达到精准识别采空区影响下的作物提取的目的。The present invention provides a remote sensing extraction method for the influence of a mined-out area of a mining area on the surface, by determining a to-be-monitored area of a mining area, and acquiring remote-sensing data of the to-be-monitored area, based on the remote-sensing data to obtain the mean value of the normalized vegetation index of the time series of the to-be-monitored area and the crop identification index, and based on the crop identification index and the mean value of the normalized vegetation index of the time series, the area to be monitored is confirmed as a farmland area through the first preset condition, and the temperature, vegetation drought index and crop nitrogen content of the farmland area are obtained, and based on The temperature vegetation drought index and the crop nitrogen content preliminarily confirm the farmland area as the area affected by the goaf through the second preset condition. The invention can use remote sensing images to extract the ground surface from the goaf in the mining area, so as to judge whether the area corresponding to the remote sensing data is the area affected by the goaf, and achieve the purpose of accurately identifying the crop extraction under the influence of the goaf.

基于以上实施例,所述基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域,包括:Based on the above embodiment, it is confirmed that the to-be-monitored area is a farmland area based on the crop identification index and the mean value of the normalized vegetation index of the time series through the first preset condition, including:

在所述作物识别指数不小于第一预设阈值,且所述时间序列的归一化植被指数均值不小于第二预设阈值的情况下,确认所述待监测区域为农田区域。When the crop identification index is not less than the first preset threshold, and the mean value of the normalized vegetation index of the time series is not less than the second preset threshold, it is confirmed that the area to be monitored is a farmland area.

具体地,本实施例是通过将作物识别指数(CRI)和时间序列的归一化植被指数

Figure BDA0003580911680000071
做一个范围限定,从而确定待监测区域为农田区域。Specifically, in this embodiment, the crop identification index (CRI) and the time series normalized vegetation index are
Figure BDA0003580911680000071
Make a scope limit to determine the area to be monitored as a farmland area.

即:CRI≥a1且

Figure BDA0003580911680000072
时,判定待监测区域为农田区域。That is: CRI≥a1 and
Figure BDA0003580911680000072
, it is determined that the area to be monitored is a farmland area.

需要说明的是,本实施例中的a1和a2根据实际情况可进行调整,在此不作具体限定。It should be noted that, a1 and a2 in this embodiment can be adjusted according to actual conditions, which are not specifically limited here.

参照图2,基于以上实施例,在所述温度植被干旱指数不小于第三预设阈值,且所述作物氮素含量不大于第四预设阈值的情况下,将所述农田区域初步确认为采空区影响的区域之后,所述方法还包括:Referring to FIG. 2, based on the above embodiment, under the condition that the temperature vegetation drought index is not less than the third preset threshold, and the nitrogen content of the crop is not greater than the fourth preset threshold, the farmland area is preliminarily confirmed as After the area affected by the gob, the method further includes:

步骤210、获取所述初步确认为采空区影响的区域内的时间序列的归一化植被指数均值以及作物生长周期,并获取矿区农田的归一化植被指数均值以及生长周期均值;Step 210, obtaining the mean value of the normalized vegetation index and the crop growth period of the time series in the area that is initially confirmed to be affected by the goaf, and obtaining the mean value of the normalized vegetation index and the mean value of the growth period of the farmland in the mining area;

步骤220、在所述初步确认为采空区影响的区域内的时间序列的归一化植被指数均值与矿区农田的归一化植被指数均值的差值不小于第五预设阈值,且所述初步确认为采空区影响的区域内的作物生长周期不等于矿区农田的生长周期均值的情况下,将所述农田区域最终确认为采空区影响的区域。Step 220: The difference between the mean value of the normalized vegetation index of the time series and the mean value of the normalized vegetation index of the farmland in the mining area is not less than the fifth preset threshold in the area that is initially confirmed to be affected by the goaf, and the If the crop growth cycle in the area preliminarily confirmed to be affected by the goaf is not equal to the average growth cycle of the farmland in the mining area, the farmland area is finally confirmed as the area affected by the goaf.

具体地,本实施在将待监测区域初步确认为采空区影响的区域之后,需要进一步确认,从而最终确认待监测区域为采空区影响的区域。Specifically, in this implementation, after the area to be monitored is preliminarily confirmed as the area affected by the goaf, further confirmation is required, so as to finally confirm that the area to be monitored is the area affected by the goaf.

首先获得初步确认为采空区影响的农田区域内的归一化植被指数

Figure BDA0003580911680000081
即时间序列的归一化植被指数均值,矿区农田的归一化植被指数均值
Figure BDA0003580911680000082
以及采空区影响的农田区域内的作物生长周期
Figure BDA0003580911680000083
和矿区农田的生长周期均值
Figure BDA0003580911680000084
First, obtain the normalized vegetation index in the farmland area that is preliminarily confirmed to be affected by the goaf
Figure BDA0003580911680000081
That is, the mean value of the normalized vegetation index of the time series, the mean value of the normalized vegetation index of the farmland in the mining area
Figure BDA0003580911680000082
and the crop growth cycle in the farmland area affected by the goaf
Figure BDA0003580911680000083
and the average growth cycle of farmland in the mining area
Figure BDA0003580911680000084

Figure BDA0003580911680000085
Figure BDA0003580911680000086
时,则可进一步判定该农田区域为最终确认的采空区影响的区域。when
Figure BDA0003580911680000085
and
Figure BDA0003580911680000086
, it can be further determined that the farmland area is the area affected by the final confirmed goaf.

本实施例通过对初步确认为采空区影响的农田区域内的归一化植被指数均值以及作物生长周期进行判断,从而将农田区域最终确认为采空区影响的农田区域,从而能够更加有效地、精准地进行识别采空区影响下的作物提取。In this embodiment, by judging the average value of the normalized vegetation index and the crop growth cycle in the farmland area that is initially confirmed to be affected by the goaf, the farmland area is finally confirmed as the farmland area affected by the goaf, so that it can be more effectively , Accurately identify crop extraction under the influence of goaf.

基于以上实施例,所述基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,包括:Based on the above embodiment, the acquisition of the time-series normalized vegetation index mean and crop identification index of the to-be-monitored area based on the remote sensing data includes:

在所述遥感数据的景数据中选择目标景数据,获取所述目标景数据的归一化植被指数以及时间序列的归一化植被指数均值;Selecting target scene data from the scene data of the remote sensing data, and obtaining the normalized vegetation index of the target scene data and the mean value of the normalized vegetation index of the time series;

基于所述目标景数据的归一化植被指数、时间序列的归一化植被指数均值以及遥感影像的景数据获取所述待监测区域的作物识别指数。The crop identification index of the to-be-monitored area is obtained based on the normalized vegetation index of the target scene data, the mean value of the normalized vegetation index of the time series, and the scene data of the remote sensing image.

具体地,本实施例从遥感数据的所有景数据n1中选择目标景数据作为计算景数据,记为i。获得目标景数据i的归一化植被指数,记为NDVIi。获取目标景数据i的时间序列的归一化植被指数均值,记为

Figure BDA0003580911680000087
作物识别指数通过以下公式实现:Specifically, in this embodiment, the target scene data is selected from all scene data n 1 of the remote sensing data as the calculation scene data, which is denoted as i. Obtain the normalized vegetation index of the target scene data i, denoted as NDVI i . Obtain the mean value of the normalized vegetation index of the time series of the target scene data i, denoted as
Figure BDA0003580911680000087
The crop identification index is achieved by the following formula:

Figure BDA0003580911680000091
Figure BDA0003580911680000091

其中,CRI为作物识别指数,n1为遥感数据的景数,i为用于计算的目标景数据,NDVIi为目标景数据i的归一化植被指数,

Figure BDA0003580911680000092
为目标景数据i的时间序列的归一化植被指数均值。Among them, CRI is the crop identification index, n 1 is the scene number of remote sensing data, i is the target scene data used for calculation, NDVI i is the normalized vegetation index of the target scene data i,
Figure BDA0003580911680000092
is the mean value of the normalized vegetation index of the time series of the target scene data i.

归一化植被指数均值,则由以下公式确定:The mean value of the normalized vegetation index is determined by the following formula:

Figure BDA0003580911680000093
Figure BDA0003580911680000093

其中,Rnir为近红外波段反射率,Rr为红波段的反射率。Among them, R nir is the reflectance in the near-infrared band, and R r is the reflectance in the red band.

本实施例通过对遥感数据的目标景数据进行计算,从而得到待监测区域内的时间序列的归一化植被指数均值、作物识别指数等参数,从而便于确认待监测区域是否为农田区域。In this embodiment, by calculating the target scene data of the remote sensing data, parameters such as the mean value of the normalized vegetation index and the crop identification index of the time series in the to-be-monitored area are obtained, thereby facilitating the confirmation of whether the to-be-monitored area is a farmland area.

基于以上实施例,所述基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域,包括:Based on the above embodiment, the said farmland area is preliminarily confirmed as the area affected by the goaf through the second preset condition based on the temperature vegetation drought index and the crop nitrogen content, including:

在所述温度植被干旱指数不小于第三预设阈值,且所述作物氮素含量不大于第四预设阈值的情况下,将所述农田区域初步确认为采空区影响的区域。Under the condition that the temperature vegetation drought index is not less than the third preset threshold, and the nitrogen content of the crop is not greater than the fourth preset threshold, the farmland area is preliminarily confirmed as the area affected by the goaf.

具体地,本实施例是通过将温度植被干旱指数(TVDI)和作物氮素含量(TCARI/OSAVI)做一个范围限定,从而将农田区域初步确认为采空区影响的区域。Specifically, in this embodiment, the temperature vegetation drought index (TVDI) and crop nitrogen content (TCARI/OSAVI) are defined within a range, thereby preliminarily confirming the farmland area as the area affected by the goaf.

即:TVDI≥a3且TCARI/OSAVI≤a4时,将农田区域初步确认为采空区影响的区域。That is, when TVDI≥a3 and TCARI/OSAVI≤a4, the farmland area is preliminarily confirmed as the area affected by the goaf.

需要说明的是,本实施例中的a3和a4根据实际情况可进行调整,在此不作具体限定。It should be noted that, a3 and a4 in this embodiment can be adjusted according to actual conditions, which are not specifically limited here.

一方面,本实施例中的TVDI是根据归一化植被指数NDVI和任意像元的地表温度Ts构成的二维矢量空间为三角形的基础上进行获得的,适用于研究特定年内某一时期某一区域的干旱程度。如以下公式所示:On the one hand, the TVDI in this embodiment is obtained on the basis that the two-dimensional vector space formed by the normalized vegetation index NDVI and the surface temperature Ts of any pixel is a triangle, which is suitable for studying a certain period of time in a specific year. the degree of drought in the area. As shown in the following formula:

Figure BDA0003580911680000101
Figure BDA0003580911680000101

其中,TVDI表示温度植被干旱指数,Tsmin表示最小地表最低温度,对应的是湿边,Ts是任意像元的地表温度,a和b为参数,即干边(dry edge)拟合方程的系数。Among them, TVDI represents the temperature vegetation drought index, Ts min represents the minimum surface temperature, corresponding to the wet edge, Ts is the surface temperature of any pixel, a and b are parameters, that is, the dry edge (dry edge) The coefficient of the fitting equation .

需要进一步说明的是,以上参数要求研究区域的范围足够大,地表覆盖从裸土变化到植被完全覆盖土壤表层含水量从萎蔫含水量变化到田间持水量。It should be further explained that the above parameters require that the scope of the study area is large enough, and the surface coverage changes from bare soil to fully covered with vegetation. The soil surface water content changes from wilting water content to field capacity.

另一方面,当作物受胁迫时,会导致叶片氮素会收到影响,因此需要结合氮素高光谱遥感诊断模型计算获得作物氮素含量(TCARI/OSAVI),如以下公式所示:On the other hand, when crops are under stress, leaf nitrogen will be affected. Therefore, it is necessary to combine the nitrogen hyperspectral remote sensing diagnostic model to calculate the crop nitrogen content (TCARI/OSAVI), as shown in the following formula:

Figure BDA0003580911680000102
Figure BDA0003580911680000102

其中,

Figure BDA0003580911680000103
表示作物氮素含量,Rx表示在x波段下对应的中心波长的反射率,Rx1、Rx2、Rx3和Rx4取值分别为700、670、550和800nm。in,
Figure BDA0003580911680000103
represents the nitrogen content of crops, Rx represents the reflectance of the corresponding central wavelength in the x-band, and the values of Rx 1 , Rx 2 , Rx 3 and Rx 4 are 700, 670, 550 and 800 nm, respectively.

基于以上实施例,确定矿区的待监测区域,并获取待监测区域的遥感数据之后,所述方法还包括:Based on the above embodiment, after determining the to-be-monitored area of the mining area and obtaining the remote sensing data of the to-be-monitored area, the method further includes:

获取所述待监测区域的年平均降水量以及作物种植季节降水量;Obtain the annual average precipitation of the area to be monitored and the precipitation in the crop planting season;

在所述年平均降水量不大于第六预设阈值,且所述作物种植季节降水量不大于第七预设阈值的情况下,将所述待监测区域判定为作物对降水敏感区。When the average annual precipitation is not greater than the sixth preset threshold, and the precipitation during the crop planting season is not greater than the seventh preset threshold, the to-be-monitored area is determined as a crop sensitive area to precipitation.

具体地,本实施例中通过获取待监测区域内的多年的年平均降水量

Figure BDA0003580911680000111
和作物种植季节降水量
Figure BDA0003580911680000112
从而判断待监测区域是否为作物对降水敏感区。Specifically, in this embodiment, the annual average precipitation in the area to be monitored for many years is obtained by obtaining
Figure BDA0003580911680000111
and precipitation during the planting season
Figure BDA0003580911680000112
Thereby, it can be judged whether the area to be monitored is a sensitive area of crops to precipitation.

即:当

Figure BDA0003580911680000113
Figure BDA0003580911680000114
时,将判定待监测区域为作物对降水敏感区。That is: when
Figure BDA0003580911680000113
and
Figure BDA0003580911680000114
, the to-be-monitored area will be determined as the crop sensitive area to precipitation.

需要说明的是,本实施例中的b1和b2根据实际情况可进行调整,在此不作具体限定。It should be noted that, b1 and b2 in this embodiment can be adjusted according to actual conditions, which are not specifically limited here.

通过对待监测区域进行降水敏感判断,从而便于获取待监测区域作物的作物识别指数以及归一化植被参数等数据,对判断待监测区域是否为采空区影响区域提供了基础。By judging the precipitation sensitivity of the area to be monitored, it is convenient to obtain the crop identification index and normalized vegetation parameters of the crops in the area to be monitored, which provides a basis for judging whether the area to be monitored is the affected area of the goaf.

下面对本发明提供的矿区采空区对地表影响的遥感提取系统进行描述,下文描述的矿区采空区对地表影响的遥感提取系统与上文描述的矿区采空区对地表影响的遥感提取方法可相互对应参照。The following describes the remote sensing extraction system for the influence of the mining area goaf on the surface provided by the present invention. refer to each other.

参照图3,本发明还提供一种矿区采空区对地表影响的遥感提取系统,包括以下模块:Referring to Fig. 3, the present invention also provides a remote sensing extraction system for the impact of a mining goaf on the surface, including the following modules:

获取模块310,用于确定矿区的待监测区域,并获取待监测区域的遥感数据;The acquisition module 310 is used to determine the area to be monitored in the mining area, and acquire remote sensing data of the area to be monitored;

农田区域确认模块320,用于在确认所述待监测区域为作物对降水敏感区的情况下,基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域;The farmland area confirmation module 320 is used to obtain the mean value of the normalized vegetation index and the crop identification index of the time series of the to-be-monitored area based on the remote sensing data when confirming that the area to be monitored is a crop sensitive to precipitation , and based on the crop identification index and the mean value of the normalized vegetation index of the time series, it is confirmed that the to-be-monitored area is a farmland area through a first preset condition;

初步确认模块330,用于获取所述农田区域的温度植被干旱指数以及作物氮素含量,并基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域。The preliminary confirmation module 330 is used to obtain the temperature vegetation drought index and the crop nitrogen content of the farmland area, and based on the temperature vegetation drought index and the crop nitrogen content, the farmland area is preliminarily confirmed as a second preset condition. The area affected by the gob.

基于以上实施例,所述农田区域确认模块具体用于:Based on the above embodiment, the farmland area confirmation module is specifically used for:

在所述作物识别指数不小于第一预设阈值,且所述时间序列的归一化植被指数均值不小于第二预设阈值的情况下,确认所述待监测区域为农田区域。When the crop identification index is not less than the first preset threshold, and the mean value of the normalized vegetation index of the time series is not less than the second preset threshold, it is confirmed that the area to be monitored is a farmland area.

基于以上实施例,所述初步确认模块具体用于:Based on the above embodiment, the preliminary confirmation module is specifically used for:

在所述温度植被干旱指数不小于第三预设阈值,且所述作物氮素含量不大于第四预设阈值的情况下,将所述农田区域初步确认为采空区影响的区域。Under the condition that the temperature vegetation drought index is not less than the third preset threshold, and the nitrogen content of the crop is not greater than the fourth preset threshold, the farmland area is preliminarily confirmed as the area affected by the goaf.

基于以上实施例,所述装置还包括:Based on the above embodiments, the device further includes:

初步区域参数获取模块,用于获取所述初步确认为采空区影响的区域内的时间序列的归一化植被指数均值以及作物生长周期,并获取矿区农田的归一化植被指数均值以及生长周期均值;The preliminary area parameter acquisition module is used to acquire the mean value of the normalized vegetation index and the crop growth period of the time series in the area that is initially confirmed to be affected by the goaf, and obtain the mean value of the normalized vegetation index and the growth period of the farmland in the mining area mean;

最终确认模块,用于在所述初步确认为采空区影响的区域内的时间序列的归一化植被指数均值与矿区农田的归一化植被指数均值的差值不小于第五预设阈值,且所述初步确认为采空区影响的区域内的作物生长周期不等于矿区农田的生长周期均值的情况下,将所述农田区域最终确认为采空区影响的区域。The final confirmation module is used for the difference between the mean value of the normalized vegetation index of the time series and the mean value of the normalized vegetation index of the farmland in the mining area in the area that is initially confirmed to be affected by the goaf is not less than the fifth preset threshold, And if the crop growth cycle in the area preliminarily confirmed to be affected by the goaf is not equal to the average growth cycle of the farmland in the mining area, the farmland area is finally confirmed as the area affected by the goaf.

基于以上实施例,所述农田区域确认模块具体用于:Based on the above embodiment, the farmland area confirmation module is specifically used for:

在所述遥感数据的景数据中选择目标景数据,获取所述目标景数据的归一化植被指数以及时间序列的归一化植被指数均值;Selecting target scene data from the scene data of the remote sensing data, and obtaining the normalized vegetation index of the target scene data and the mean value of the normalized vegetation index of the time series;

基于所述目标景数据的归一化植被指数、时间序列的归一化植被指数均值以及遥感影像的景数据获取所述待监测区域的作物识别指数。The crop identification index of the to-be-monitored area is obtained based on the normalized vegetation index of the target scene data, the mean value of the normalized vegetation index of the time series, and the scene data of the remote sensing image.

基于以上实施例,所述装置还包括:Based on the above embodiments, the device further includes:

降水量获取模块,用于获取所述待监测区域的年平均降水量以及作物种植季节降水量;a precipitation acquisition module, used to acquire the annual average precipitation of the to-be-monitored area and the precipitation in the crop planting season;

判定模块,用于在所述年平均降水量不大于第六预设阈值,且所述作物种植季节降水量不大于第七预设阈值的情况下,将所述待监测区域判定为作物对降水敏感区。A determination module, configured to determine the to-be-monitored area as a crop-to-rain precipitation under the condition that the annual average precipitation is not greater than the sixth preset threshold, and the precipitation during the crop planting season is not greater than the seventh preset threshold sensitive area.

图4示例了一种电子设备的实体结构示意图,如图4所示,该电子设备可以包括:处理器(processor)410、通信接口(Communications Interface)420、存储器(memory)430和通信总线440,其中,处理器410,通信接口420,存储器430通过通信总线440完成相互间的通信。处理器410可以调用存储器430中的逻辑指令,以执行矿区采空区对地表影响的遥感提取方法,该方法包括:FIG. 4 illustrates a schematic diagram of the physical structure of an electronic device. As shown in FIG. 4 , the electronic device may include: a processor (processor) 410, a communication interface (Communications Interface) 420, a memory (memory) 430, and a communication bus 440, The processor 410 , the communication interface 420 , and the memory 430 communicate with each other through the communication bus 440 . The processor 410 can invoke the logic instructions in the memory 430 to execute a remote sensing extraction method for the impact of a goaf in a mining area on the surface, the method comprising:

确定矿区的待监测区域,并获取待监测区域的遥感数据;Determine the area to be monitored in the mining area, and obtain remote sensing data of the area to be monitored;

在确认所述待监测区域为作物对降水敏感区的情况下,基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域;In the case of confirming that the to-be-monitored area is a crop sensitive area to precipitation, obtain the time-series normalized vegetation index mean value and crop identification index of the to-be-monitored area based on the remote sensing data, and based on the crop identification index and the mean value of the normalized vegetation index of the time series through the first preset condition to confirm that the area to be monitored is a farmland area;

获取所述农田区域的温度植被干旱指数以及作物氮素含量,并基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域。Obtain the temperature vegetation drought index and crop nitrogen content of the farmland area, and based on the temperature vegetation drought index and the crop nitrogen content, the farmland area is preliminarily confirmed as an area affected by a goaf through a second preset condition.

此外,上述的存储器430中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。In addition, the above-mentioned logic instructions in the memory 430 can be implemented in the form of software functional units and can be stored in a computer-readable storage medium when sold or used as an independent product. Based on such understanding, the technical solution of the present invention can be embodied in the form of a software product in essence, or the part that contributes to the prior art or the part of the technical solution. The computer software product is stored in a storage medium, including Several instructions are used to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, Read-Only Memory (ROM, Read-Only Memory), Random Access Memory (RAM, Random Access Memory), magnetic disk or optical disk and other media that can store program codes .

另一方面,本发明还提供一种计算机程序产品,所述计算机程序产品包括计算机程序,计算机程序可存储在非暂态计算机可读存储介质上,所述计算机程序被处理器执行时,计算机能够执行上述各方法所提供的矿区采空区对地表影响的遥感提取方法,该方法包括:In another aspect, the present invention also provides a computer program product, the computer program product includes a computer program, the computer program can be stored on a non-transitory computer-readable storage medium, and when the computer program is executed by a processor, the computer can Execute the remote sensing extraction method for the impact of mined-out areas in mining areas on the surface provided by the above methods, the method includes:

确定矿区的待监测区域,并获取待监测区域的遥感数据;Determine the area to be monitored in the mining area, and obtain remote sensing data of the area to be monitored;

在确认所述待监测区域为作物对降水敏感区的情况下,基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域;In the case of confirming that the area to be monitored is a crop sensitive area to precipitation, obtain the mean value of the normalized vegetation index and the crop identification index of the time series of the area to be monitored based on the remote sensing data, and based on the crop identification index and the mean value of the normalized vegetation index of the time series through the first preset condition to confirm that the area to be monitored is a farmland area;

获取所述农田区域的温度植被干旱指数以及作物氮素含量,并基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域。Obtain the temperature vegetation drought index and crop nitrogen content of the farmland area, and based on the temperature vegetation drought index and the crop nitrogen content, the farmland area is preliminarily confirmed as an area affected by a goaf through a second preset condition.

又一方面,本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各方法提供的矿区采空区对地表影响的遥感提取方法,该方法包括:In yet another aspect, the present invention also provides a non-transitory computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, the computer program is implemented to perform the effects of the mine goaf on the surface provided by the above methods. Remote sensing extraction method, the method includes:

确定矿区的待监测区域,并获取待监测区域的遥感数据;Determine the area to be monitored in the mining area, and obtain remote sensing data of the area to be monitored;

在确认所述待监测区域为作物对降水敏感区的情况下,基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域;In the case of confirming that the to-be-monitored area is a crop sensitive area to precipitation, obtain the time-series normalized vegetation index mean value and crop identification index of the to-be-monitored area based on the remote sensing data, and based on the crop identification index and the mean value of the normalized vegetation index of the time series through the first preset condition to confirm that the area to be monitored is a farmland area;

获取所述农田区域的温度植被干旱指数以及作物氮素含量,并基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域。Obtain the temperature vegetation drought index and crop nitrogen content of the farmland area, and based on the temperature vegetation drought index and the crop nitrogen content, the farmland area is preliminarily confirmed as an area affected by a goaf through a second preset condition.

以上所描述的系统实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。The system embodiments described above are only illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in One place, or it can be distributed over multiple network elements. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution in this embodiment. Those of ordinary skill in the art can understand and implement it without creative effort.

通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。From the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and certainly can also be implemented by hardware. Based on this understanding, the above-mentioned technical solutions can be embodied in the form of software products in essence or the parts that make contributions to the prior art, and the computer software products can be stored in computer-readable storage media, such as ROM/RAM, magnetic A disc, an optical disc, etc., includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform the methods described in various embodiments or some parts of the embodiments.

最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, but not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that it can still be Modifications are made to the technical solutions described in the foregoing embodiments, or some technical features thereof are equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions depart from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (10)

1.一种矿区采空区对地表影响的遥感提取方法,其特征在于,包括:1. a kind of remote sensing extraction method of mining goaf impact on surface, is characterized in that, comprises: 确定矿区的待监测区域,并获取待监测区域的遥感数据;Determine the area to be monitored in the mining area, and obtain remote sensing data of the area to be monitored; 在确认所述待监测区域为作物对降水敏感区的情况下,基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域;In the case of confirming that the to-be-monitored area is a crop sensitive area to precipitation, obtain the time-series normalized vegetation index mean value and crop identification index of the to-be-monitored area based on the remote sensing data, and based on the crop identification index and the mean value of the normalized vegetation index of the time series through the first preset condition to confirm that the area to be monitored is a farmland area; 获取所述农田区域的温度植被干旱指数以及作物氮素含量,并基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域。Acquire the temperature vegetation drought index and crop nitrogen content of the farmland area, and based on the temperature vegetation drought index and the crop nitrogen content, the farmland area is preliminarily confirmed as an area affected by a goaf through a second preset condition. 2.根据权利要求1所述的矿区采空区对地表影响的遥感提取方法,其特征在于,所述基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域,包括:2. The remote sensing extraction method of mining goaf according to claim 1, characterized in that, the normalized vegetation index mean value based on the crop identification index and time series is confirmed by the first preset condition The to-be-monitored area is a farmland area, including: 在所述作物识别指数不小于第一预设阈值,且所述时间序列的归一化植被指数均值不小于第二预设阈值的情况下,确认所述待监测区域为农田区域。When the crop identification index is not less than the first preset threshold, and the mean value of the normalized vegetation index of the time series is not less than the second preset threshold, it is confirmed that the area to be monitored is a farmland area. 3.根据权利要求1所述的矿区采空区对地表影响的遥感提取方法,其特征在于,所述基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域,包括:3. The remote sensing extraction method of mining goaf impact on surface according to claim 1, characterized in that, described farmland area based on the temperature vegetation drought index and crop nitrogen content by the second preset condition. Areas initially identified as goaf affected include: 在所述温度植被干旱指数不小于第三预设阈值,且所述作物氮素含量不大于第四预设阈值的情况下,将所述农田区域初步确认为采空区影响的区域。Under the condition that the temperature vegetation drought index is not less than the third preset threshold, and the nitrogen content of the crop is not greater than the fourth preset threshold, the farmland area is preliminarily confirmed as the area affected by the goaf. 4.根据权利要求3所述的矿区采空区对地表影响的遥感提取方法,其特征在于,在所述温度植被干旱指数不小于第三预设阈值,且所述作物氮素含量不大于第四预设阈值的情况下,将所述农田区域初步确认为采空区影响的区域之后,所述方法还包括:4. The remote sensing extraction method of mining goaf impact on the surface according to claim 3, characterized in that, at the temperature, the vegetation drought index is not less than the third preset threshold, and the nitrogen content of the crops is not greater than the third preset threshold. In the case of four preset thresholds, after the farmland area is preliminarily confirmed as the area affected by the goaf, the method further includes: 获取所述初步确认为采空区影响的区域内的时间序列的归一化植被指数均值以及作物生长周期,并获取矿区农田的归一化植被指数均值以及生长周期均值;Obtaining the mean value of the normalized vegetation index and the crop growth period of the time series in the area that is initially confirmed to be affected by the goaf, and obtaining the mean value of the normalized vegetation index and the mean value of the growth period of the farmland in the mining area; 在所述初步确认为采空区影响的区域内的时间序列的归一化植被指数均值与矿区农田的归一化植被指数均值的差值不小于第五预设阈值,且所述初步确认为采空区影响的区域内的作物生长周期不等于矿区农田的生长周期均值的情况下,将所述农田区域最终确认为采空区影响的区域。The difference between the mean value of the normalized vegetation index of the time series and the mean value of the normalized vegetation index of the farmland in the mining area in the area that is preliminarily confirmed to be affected by the goaf is not less than the fifth preset threshold, and the preliminary confirmation is If the crop growth cycle in the area affected by the goaf is not equal to the average growth cycle of the farmland in the mining area, the farmland area is finally confirmed as the area affected by the goaf. 5.根据权利要求1所述的矿区采空区对地表影响的遥感提取方法,其特征在于,所述基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,包括:5. the remote sensing extraction method of mining goaf according to claim 1, is characterized in that, described based on described remote sensing data obtains the normalized vegetation index mean value and crop of the time series of described area to be monitored Identify indices, including: 在所述遥感数据的景数据中选择目标景数据,获取所述目标景数据的归一化植被指数以及时间序列的归一化植被指数均值;Selecting target scene data from the scene data of the remote sensing data, and obtaining the normalized vegetation index of the target scene data and the mean value of the normalized vegetation index of the time series; 基于所述目标景数据的归一化植被指数、时间序列的归一化植被指数均值以及遥感影像的景数据获取所述待监测区域的作物识别指数。The crop identification index of the to-be-monitored area is obtained based on the normalized vegetation index of the target scene data, the mean value of the normalized vegetation index of the time series, and the scene data of the remote sensing image. 6.根据权利要求1-5任一所述的矿区采空区对地表影响的遥感提取方法,其特征在于,确定矿区的待监测区域,并获取待监测区域的遥感数据之后,所述方法还包括:6. The remote sensing extraction method for the influence of mined-out areas in mining areas on the surface according to any one of claims 1-5, characterized in that, after determining the area to be monitored in the mining area, and after obtaining the remote sensing data of the area to be monitored, the method further includes: include: 获取所述待监测区域的年平均降水量以及作物种植季节降水量;Obtain the annual average precipitation of the area to be monitored and the precipitation in the crop planting season; 在所述年平均降水量不大于第六预设阈值,且所述作物种植季节降水量不大于第七预设阈值的情况下,将所述待监测区域判定为作物对降水敏感区。When the average annual precipitation is not greater than the sixth preset threshold, and the precipitation during the crop planting season is not greater than the seventh preset threshold, the to-be-monitored area is determined as a crop sensitive area to precipitation. 7.一种矿区采空区对地表影响的遥感提取系统,其特征在于,包括:7. A remote sensing extraction system for the influence of a goaf in a mining area on the surface, comprising: 获取模块,用于确定矿区的待监测区域,并获取待监测区域的遥感数据;The acquisition module is used to determine the area to be monitored in the mining area, and obtain the remote sensing data of the area to be monitored; 农田区域确认模块,用于在确认所述待监测区域为作物对降水敏感区的情况下,基于所述遥感数据获取所述待监测区域的时间序列的归一化植被指数均值和作物识别指数,并基于所述作物识别指数和时间序列的归一化植被指数均值通过第一预设条件确认所述待监测区域为农田区域;a farmland area confirmation module, configured to obtain the mean value of the normalized vegetation index and the crop identification index of the time series of the to-be-monitored area based on the remote sensing data under the condition that the to-be-monitored area is confirmed as a crop sensitive area to precipitation, and confirming that the to-be-monitored area is a farmland area through the first preset condition based on the crop identification index and the mean value of the normalized vegetation index of the time series; 初步确认模块,用于获取所述农田区域的温度植被干旱指数以及作物氮素含量,并基于所述温度植被干旱指数和作物氮素含量通过第二预设条件将所述农田区域初步确认为采空区影响的区域。The preliminary confirmation module is used to obtain the temperature vegetation drought index and the crop nitrogen content of the farmland area, and based on the temperature vegetation drought index and the crop nitrogen content, the farmland area is preliminarily confirmed as harvesting through the second preset condition. The area affected by the empty zone. 8.一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至6任一项所述矿区采空区对地表影响的遥感提取方法。8. An electronic device comprising a memory, a processor and a computer program stored on the memory and running on the processor, wherein the processor implements the program as claimed in claim 1 when executing the program A remote sensing extraction method for the impact of goafs in mining areas on the surface of any one of to 6. 9.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至6任一项所述矿区采空区对地表影响的遥感提取方法。9. A non-transitory computer-readable storage medium having a computer program stored thereon, is characterized in that, when the computer program is executed by the processor, it is characterized in that, when the computer program is executed by the processor, the mining area goaf as described in any one of claims 1 to 6 is implemented. Remote sensing extraction methods for surface impacts. 10.一种计算机程序产品,包括计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至6任一项所述矿区采空区对地表影响的遥感提取方法。10 . A computer program product, comprising a computer program, characterized in that, when the computer program is executed by a processor, the remote sensing extraction method for the influence of a goaf in a mining area on the ground surface according to any one of claims 1 to 6 is implemented.
CN202210351884.XA 2022-04-02 2022-04-02 Remote sensing extraction method and system of mining goaf's influence on the surface Pending CN114724029A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210351884.XA CN114724029A (en) 2022-04-02 2022-04-02 Remote sensing extraction method and system of mining goaf's influence on the surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210351884.XA CN114724029A (en) 2022-04-02 2022-04-02 Remote sensing extraction method and system of mining goaf's influence on the surface

Publications (1)

Publication Number Publication Date
CN114724029A true CN114724029A (en) 2022-07-08

Family

ID=82242733

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210351884.XA Pending CN114724029A (en) 2022-04-02 2022-04-02 Remote sensing extraction method and system of mining goaf's influence on the surface

Country Status (1)

Country Link
CN (1) CN114724029A (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102706328A (en) * 2012-04-28 2012-10-03 中国神华能源股份有限公司 Method and device for detecting subsidence area of ground and data processing equipment
WO2018173622A1 (en) * 2017-03-23 2018-09-27 日本電気株式会社 Vegetation index calculation device, vegetation index calculation method, and computer readable recording medium
CN110472357A (en) * 2019-08-21 2019-11-19 华北理工大学 Assess construction method and the application of the remote sensing comprehensive ecological model RSIEI of the different effect of Mining Development compact district earth's surface thermal environment point
CN110567892A (en) * 2019-09-17 2019-12-13 中国水利水电科学研究院 A Hyperspectral Prediction Method of Nitrogen in Summer Maize Based on Critical Nitrogen Concentration
CN110598553A (en) * 2019-08-09 2019-12-20 中国科学院南京地理与湖泊研究所 Original true landform mining damaged area detection method based on remote sensing image and topographic data
CN113984212A (en) * 2021-10-27 2022-01-28 中国气象科学研究院 Agricultural irrigated area extraction method and system
CN114021973A (en) * 2021-11-04 2022-02-08 国能神东煤炭集团有限责任公司 Evaluation method and system of ecosystem stability in coal mining subsidence area

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102706328A (en) * 2012-04-28 2012-10-03 中国神华能源股份有限公司 Method and device for detecting subsidence area of ground and data processing equipment
WO2018173622A1 (en) * 2017-03-23 2018-09-27 日本電気株式会社 Vegetation index calculation device, vegetation index calculation method, and computer readable recording medium
CN110598553A (en) * 2019-08-09 2019-12-20 中国科学院南京地理与湖泊研究所 Original true landform mining damaged area detection method based on remote sensing image and topographic data
CN110472357A (en) * 2019-08-21 2019-11-19 华北理工大学 Assess construction method and the application of the remote sensing comprehensive ecological model RSIEI of the different effect of Mining Development compact district earth's surface thermal environment point
CN110567892A (en) * 2019-09-17 2019-12-13 中国水利水电科学研究院 A Hyperspectral Prediction Method of Nitrogen in Summer Maize Based on Critical Nitrogen Concentration
CN113984212A (en) * 2021-10-27 2022-01-28 中国气象科学研究院 Agricultural irrigated area extraction method and system
CN114021973A (en) * 2021-11-04 2022-02-08 国能神东煤炭集团有限责任公司 Evaluation method and system of ecosystem stability in coal mining subsidence area

Similar Documents

Publication Publication Date Title
Muchoney et al. Change detection for monitoring forest defoliation
Johansen et al. Using Unmanned Aerial Vehicles to assess the rehabilitation performance of open cut coal mines
Mutlu et al. Mapping surface fuel models using lidar and multispectral data fusion for fire behavior
McMorrow Linear regression modelling for the estimation of oil palm age from Landsat TM
Eckert et al. Identification and classification of structural soil conservation measures based on very high resolution stereo satellite data
CN112712038A (en) Method and system for monitoring wheat lodging condition based on multispectral satellite image
CN117233123B (en) Large-scale remote sensing monitoring method and device for rice bacterial blight based on Sentinel-2
Nguyen et al. Mapping Land use/land cover using a combination of Radar Sentinel-1A and Sentinel-2A optical images
CN108416508A (en) The extracting method and system on a kind of earth's surface ecology influence boundary
Monnet et al. Tree top detection using local maxima filtering: a parameter sensitivity analysis
CN115690580A (en) Corn lodging remote sensing monitoring method and system, electronic device and medium
Kennaway et al. Mapping land cover and estimating forest structure using satellite imagery and coarse resolution lidar in the Virgin Islands
Feldpausch et al. Biomass, harvestable area, and forest structure estimated from commercial timber inventories and remotely sensed imagery in southern Amazonia
CN114724029A (en) Remote sensing extraction method and system of mining goaf's influence on the surface
Benseghir et al. Shortwave infrared vegetation index-based modelling for aboveground vegetation biomass assessment in the arid steppes of Algeria
Rinnamang et al. Estimation of aboveground biomass using aerial photogrammetry from unmanned aerial vehicle in teak (Tectona grandis) plantation in Thailand
Zhang et al. A novel index for mapping crop residue covered cropland using remote sensing data
CN115775354B (en) Grouting period rice harvest-stop extraction method based on fusion of remote sensing indexes
Xu et al. Advanced exploratory data analysis for mapping regional canopy cover
CN115935111A (en) Tobacco field waterlogging early identification method based on hyperspectrum of unmanned aerial vehicle
Noonan et al. A method for mapping the distribution of willow at a catchment scale using bi‐seasonal SPOT5 imagery
Peştemalci et al. Acreage estimation of wheat and barley fields in the province of Adana, Turkey
Vafidis et al. Mapping scrub vegetation cover from photogrammetric point clouds: A useful tool in reserve management
Seidl IDENTIFICATION OF RETENTION AREAS USING AIRBORNE LIDAR DATA. A CASE STUDY FROM CENTRAL SWEDEN.
Gessner et al. Towards an automated estimation of vegetation cover fractions on multiple scales: Examples of Eastern and Southern Africa

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination