CN-121725359-B - Nickel ore target area prediction method and system based on hyperspectral remote sensing data analysis

Abstract

The invention provides a nickel mine target area prediction method and a system based on hyperspectral remote sensing data analysis, which relate to the technical field of mineral resource exploration, firstly acquire hyperspectral remote sensing data bodies of areas to be predicted, then execute alteration characteristic space clustering construction processing on the hyperspectral remote sensing data bodies to form alteration characteristic space cluster bodies, identify the space distribution of various geological alteration types and bind the spectrum characteristics, and then carrying out ore formation alteration combined sequential spatial evolution treatment to obtain an ore formation alteration sequential evolution body, then carrying out target area indication feature anchoring shaping construction treatment on the ore formation alteration sequential evolution body to form a target area indication spatial shaping body, positioning nickel ore target area key indication features, and finally carrying out nickel ore target area space normalization circle setting treatment to generate a nickel ore target area prediction result. According to the method, the spectrum and space information of the hyperspectral remote sensing data are fully utilized, and the accuracy and reliability of nickel ore target area prediction are effectively improved.

Inventors

• QIN SHUKAI
• Wen Zhenqin
• LING YAJUN
• CHENG WENQI
• WU DEQIANG
• KONG DECAI
• JIAN CHU
• TANG JUNHUA
• PENG SONGLIN
• LONG BO

Assignees

• 四川省第四地质大队

Dates

Publication Date

20260505

Application Date

20260212

Claims (10)

1. 1. The method for predicting the target area of the nickel ore based on hyperspectral remote sensing data analysis is characterized by comprising the following steps of: Acquiring a hyperspectral remote sensing data body of a region to be predicted, wherein the hyperspectral remote sensing data body comprises full-band spectral response characteristics of each geographic detection point in the region to be predicted and geographic space spread characteristics of each geographic detection point, the full-band spectral response characteristics are continuous response characteristics of spectral reflection and spectral absorption formed by the geographic detection points in a hyperspectral remote sensing detection band range, and the geographic space spread characteristics are geographic coordinates of each geographic detection point and spatial adjacent and spatial extension relations among the points; Performing clustered construction treatment on the hyperspectral remote sensing data body to form a clustered body of the changed characteristic space, wherein the clustered body of the changed characteristic space is a continuous clustered distribution structure formed by various geological change types depending on geographical space spread characteristics, and the continuous clustered distribution structure binds full-band spectral response characteristics of the corresponding geological change types; Performing ore-forming alteration combination sequential spatial evolution processing on the alteration feature space cluster body to form an ore-forming alteration sequential evolution body, wherein the ore-forming alteration sequential evolution body is a spatial sequential extension structure formed by geological alteration combination matched with nickel ore-forming features according to the occurrence sequence of the ore-forming alterations, and the spatial sequential extension structure binds full-band spectral response feature evolution content of the corresponding alteration combination; Performing target region indication feature anchoring forming construction treatment on the ore forming alteration sequential evolution body to form a target region indication space forming body, wherein the target region indication space forming body is a formed space distribution structure constructed by taking nickel ore forming core alteration features as anchoring centers, and the formed space distribution structure binds multi-stage target region indication alteration features; And carrying out nickel ore target space normalization circle setting treatment on the target space indication space excipient to form a nickel ore target prediction result.
2. 2. The method for predicting a target area of a nickel ore based on hyperspectral remote sensing data analysis as claimed in claim 1, wherein the performing the processing of the clustered construction of the altered feature space on the hyperspectral remote sensing data body to form the clustered body of the altered feature space comprises: Extracting full-band spectral response characteristics in the hyperspectral remote sensing data body, and carrying out point-by-point binding on the full-band spectral response characteristics of each geographic detection point and geographic coordinates of the corresponding geographic detection point to form a spectral characteristic geographic point binding set, wherein each geographic coordinate in the spectral characteristic geographic point binding set uniquely corresponds to one group of full-band spectral response characteristics; a standard spectral response library of the geological alteration type related to nickel ore mineralization is prepared, standard full-wave band spectral response characteristics corresponding to various geological alteration types in the standard spectral response library are extracted, and all the standard full-wave band spectral response characteristics are subjected to wave band synchronous normalization according to the wave band division rules of hyperspectral remote sensing detection, so that a normalized standard spectral characteristic set is formed; The single standard full-band spectral response characteristic in the normalized standard spectral characteristic set and the single group of detection full-band spectral response characteristic in the spectral characteristic geographic point position binding set are subjected to band-by-band characteristic comparison, spectral fragments matched with the standard full-band spectral response characteristic are screened from the detection full-band spectral response characteristic, and a single-point position spectral characteristic matching set is formed; Sequentially completing band-by-band feature comparison and spectrum fragment screening of detection full-band spectrum response features of all geographic detection points in the spectrum feature geographic point binding set, collecting all geographic detection points forming a single-point spectrum feature matching set, and classifying the corresponding geographic detection points according to the matched geological alteration type to form an alteration type geographic point classification set; Extracting geographic space spread features of all geographic detection points corresponding to a single type of geological alteration type in the alteration type geographic point classification set, analyzing the spatial adjacent relation among the geographic detection points, and identifying geographic detection point groups with continuous spatial adjacent relation to form a single type of alteration continuous point group set; extracting geographic coordinates of all geographic detection points in a single geographic detection point group in the single continuous point group set, continuously linking the geographic coordinates of all the points according to the spatial extension direction of the geographic detection points to form single-type geological changed spatially continuous distribution lines, and constructing a single-type changed spatially line set; Extracting full-band spectral response characteristics of all geographic detection points on a single spatially continuous distribution line in the single type of alteration space line set, continuously arranging the full-band spectral response characteristics of each point according to the spatially extending direction to form spectral characteristic continuous arrangement content of the single line, and binding the spectral characteristic continuous arrangement content to the corresponding spatially continuous distribution line to form an alteration line spectral binding set; Integrating all the spatial continuous distribution lines and the continuous distribution content of the bound spectral features under the same geological alteration type in the spectral binding set of the alteration lines, and constructing an integral cluster distribution structure of single geological alteration by combining the spatial position relation among the lines to form a single alteration cluster structure set; Arranging all the cluster distribution structures of the geological alteration type in the single type alteration cluster structure set according to the actual geographical space distribution relation in the region to be predicted, and reserving the actual spatial position correlation among various cluster distribution structures to form a global alteration cluster arrangement set; And performing global fusion on all cluster distribution structures in the global alteration cluster arrangement set, full-band spectral response characteristics and geospatial spread characteristics of corresponding geological alteration types to form alteration characteristic space cluster bodies, wherein the alteration characteristic space cluster bodies comprise all cluster distribution structures of nickel ore mineralization related geological alteration types in a region to be predicted and corresponding bound characteristic contents.
3. 3. The method for predicting a target region of nickel ore based on hyperspectral remote sensing data analysis as claimed in claim 1, wherein the performing the ore-forming alteration combined sequential spatial evolution processing on the altered feature space cluster body to form an ore-forming alteration sequential evolution body comprises: Extracting various geology-changed cluster-shaped distribution structures and corresponding geospatial distribution ranges in the changed characteristic space cluster-shaped body, and integrally distributing all cluster-shaped distribution structures according to the geospatial trend of the area to be predicted to form a changed cluster-shaped space trend distribution set; Analyzing the position connection relation of different geological alteration cluster distribution structures in the geographical space in the alteration cluster space trend arrangement set, and identifying geological alteration cluster distribution structure combinations with adjacent or direct connection on the spatial distribution to form an alteration cluster space connection combination set; Taking a sequential occurrence rule of geological changes in the nickel ore mineralization process, extracting the occurrence sequence of various geological changes in the sequential occurrence rule and the evolution association relation among the changes, and extracting core geological change types and combination forms corresponding to each stage of nickel ore mineralization to form an ore-formation-change sequential rule set; Comparing the types of the changes contained in the single changes in the cluster-shaped space connection combination and the arrangement sequence of the single changes in the space with the main sequence of the changes defined in the mine-forming changes in the sequence rule set, if the types of the changes in the cluster-shaped space connection combination and the arrangement sequence are one continuous subsequence of the main sequence, judging that the changes are matched, and returning all the matched changes to the mine-related changes in the cluster-shaped space connection combination; Extracting various geological alteration cluster-shaped distribution structures in a single ore-forming related alteration combination in the initial selection of the ore-forming related alteration combination, and sequentially arranging the cluster-shaped distribution structures according to the alteration occurrence sequence of nickel ore forming to form a single-combination alteration cluster-shaped sequential arrangement structure; Extracting geospatial coordinates of cluster-shaped distribution structures subjected to geological alteration at each stage in the single-combination alteration cluster-shaped sequential distribution structure, continuously connecting the cluster-shaped distribution structures subjected to the previous alteration and the subsequent alteration according to a space adjacent relation to form a space sequential extension structure subjected to ore formation alteration combination, and constructing a single-combination space sequential extension set; Extracting full-band spectral response characteristics of geological alteration of each stage in the single-combination spatial sequential extension centralized mineral alteration combined spatial sequential extension structure, and continuously arranging the full-band spectral response characteristics of each stage according to an alteration sequential evolution direction to form spectral characteristic sequential evolution arrangement content, wherein the spectral characteristic sequential evolution arrangement content is bound to a corresponding spatial sequential extension structure; extracting the alteration cluster distribution structure and the bound spectral feature sequential evolution arrangement content corresponding to the nickel ore forming core stage in the single-combination spatial sequential extension set, carrying out key extraction on the feature content of the core stage to form a core ore forming stage alteration feature set, and binding to the spatial sequential extension structure corresponding to the ore forming alteration combination; the spatial sequential extension structures of all the ore-forming related alteration combinations in the initial selection of the ore-forming related alteration combinations and the sequential evolution arrangement content of the bound spectrum features are integrally arranged according to the actual geographic spatial distribution of the region to be predicted, so that a global ore-forming alteration sequential arrangement set is formed; And performing global fusion on the spatial sequential extension structure of all the ore-forming alteration combinations in the global ore-forming alteration sequential arrangement and various bound characteristic contents to form an ore-forming alteration sequential evolution body, wherein the ore-forming alteration sequential evolution body comprises the dual characteristics of sequential evolution and spatial distribution of all the ore-forming related alteration combinations in a region to be predicted.
4. 4. The method of claim 1, wherein performing a target indication feature anchor shape construction process on the mine-forming alteration sequential evolution body to form a target indication space shape comprises: Extracting a spatial sequential extension structure of each ore-forming alteration combination in the ore-forming alteration sequential evolution body, identifying alteration cluster-shaped distribution structures corresponding to a nickel ore-forming core stage in each ore-forming alteration combination, determining a geographic space center coordinate of the alteration cluster-shaped distribution structures, and taking the geographic space center coordinate as a core anchoring point position indicated by a target area to form a core anchoring point position set; extracting spatial sequential extension structures of single diagenetic changes in the diagenetic changes sequential evolution body, identifying all previous alteration cluster-shaped distribution structures positioned in front of a core alteration cluster-shaped distribution structure corresponding to a core anchoring point, and all subsequent alteration cluster-shaped distribution structures positioned behind the previous alteration cluster-shaped distribution structures, and extracting a geographic space range covered by an external polygon of the alteration cluster-shaped distribution structure to form an initial alteration combination spatial association range; fitting according to a preset geometric template corresponding to the typical nickel ore forming space morphology based on the external polygon of the initial alteration combination space association range by taking the core anchoring point as the center to generate a single combination forming space distribution structure; extracting geological alteration cluster distribution structures and full-band spectral response characteristics corresponding to geometric parts in the single combined shaping space distribution structure, and binding various characteristic contents with specific geometric parts of the shaping space point by point to form a shaping space characteristic part binding set; Binding the formed space feature parts, marking the changed features corresponding to the core anchoring points as first-stage anchoring features, dividing a plurality of concentric annular bands according to a preset distance radius by taking the core anchoring points as circle centers, sequentially marking the changed features in different annular bands as multi-stage anchoring features, and forming a multi-stage target area indication anchoring feature set; extracting geospatial coordinates of each level of anchoring features in the multi-level target region indication anchoring feature set, constructing spatial association links between each level of anchoring features and each level of anchoring features, and binding all the spatial association links with corresponding shaped spatial distribution structures to form a shaped spatial feature association link set; Binding the forming space distribution structure corresponding to all the forming change combinations in the forming change sequential evolution body and the forming space feature part binding set, and integrally arranging the multi-level target region indication anchoring feature set and the forming space feature associated link set according to the actual geographic space position of the region to be predicted to form a global target region indication forming arrangement set; extracting boundary characteristics of each forming space distribution structure in the global target region indication forming arrangement set, and carrying out characterization marking on the boundary of the forming space by combining the corresponding multi-level target region indication anchoring characteristic set so as to ensure that the distribution range of boundary marking content and anchoring characteristics is kept consistent to form a forming space characterization boundary set; Extracting all the forming space distribution structures and various bound characteristic contents in the global target region indication forming arrangement set, removing characteristic contents and space links which are not directly related to nickel ore forming, reserving target region indication characteristics and space association relations of cores, and forming a simplified target region indication forming set; And performing global fusion on all the formed spatial distribution structures in the reduced target region indication forming set and the bound multi-level target region indication anchoring feature set, the formed spatial feature association link set and the formed spatial feature boundary set to form a target region indication spatial forming body, wherein the target region indication spatial forming body comprises anchoring contents of all target region indication features in a region to be predicted and the formed spatial distribution structures.
5. 5. The method for predicting a nickel mine target area based on hyperspectral remote sensing data analysis as claimed in claim 1, wherein the performing the nickel mine target area space normalization processing on the target area indication space excipient to form a nickel mine target area prediction result comprises: Extracting each formed space distribution structure in the target region indication space forming body, extracting a characteristic boundary of each formed space distribution structure, and obtaining geographic coordinate information corresponding to the characteristic boundary to form a formed space boundary coordinate extraction set; extracting boundary geographic coordinate points of a single forming space distribution structure in the forming space boundary coordinate extraction set, and continuously linking discrete boundary coordinate points according to geographic space trend to form a single forming space geographic range with continuous boundaries, so as to construct a single forming space continuous range set; Analyzing the geographic space range of the adjacent forming space distribution structure in the single forming space continuous range, extracting continuous boundary coordinates of two forming spaces with space adjacent or space overlapping, and carrying out fusion treatment on the boundary coordinates of adjacent or overlapping parts to form a fused continuous geographic boundary; Extracting all geographic space ranges covered by the fused continuous geographic boundaries, taking the geographic space ranges as an integral mineral potential space range, and extracting a multi-level target area indication anchoring feature set and a forming space feature association link set bound in the mineral potential space range to form a fused mineral potential space feature set; Extracting the geographic space ranges of the independent forming spatial distribution structures which are not formed in the single forming spatial continuous range set and are adjacent or overlapped in space, and reserving the original continuous geographic boundaries and the indication characteristic contents of various bound target areas to form an independent mineralization potential spatial characteristic set; The fusion type ore-forming potential space feature set and the independent ore-forming potential space feature set are subjected to global integration, all ore-forming potential space ranges are distributed according to geographical space partitions of a region to be predicted, and a global ore-forming potential space primary integral set is formed; Extracting complete geographic boundary coordinate sequences of all mineral formation potential space ranges in the whole domain mineral formation potential space initial integration, thinning the boundary coordinate sequences by using a Targes-Prak algorithm, and reserving characteristic coordinate points within a set tolerance range to form a regularized mineral formation potential space boundary set; Accurately matching the geographic boundary of each mineral formation potential space range in the regularized mineral formation potential space boundary set with the corresponding multi-level target region indication anchoring feature set and the mineral formation alteration combination characterization content to form a one-to-one correspondence of the space range and the feature characterization, and constructing a spatial feature accurate matching set; Extracting geographic coordinate ranges, core anchoring features, multi-stage anchoring features and ore-forming alteration combination evolution relations of all ore-forming potential spatial ranges in the spatial feature accurate matching set, and carrying out region-by-region complete recording on all contents to form a global ore-forming potential spatial feature record set; And performing global summarization on the global mineral formation potential space feature record set to form a nickel mine target area prediction result, wherein the nickel mine target area prediction result comprises all geographic space ranges with nickel mine mineral formation potential in the area to be predicted and complete mineral formation feature characterization contents bound by all the geographic space ranges.
6. 6. The method for predicting a nickel ore target area based on hyperspectral remote sensing data analysis according to claim 2, wherein the step of comparing the single standard full-band spectral response characteristic in the normalized standard spectral characteristic set with the single group of detected full-band spectral response characteristics in the spectral characteristic geographic point binding set, and screening out spectral fragments matched with the standard full-band spectral response characteristics from the detected full-band spectral response characteristics to form a single-point spectral characteristic matching set comprises the steps of: extracting standard full-band spectral response characteristics of single-class geological alteration in the normalized standard spectral characteristics, splitting the standard full-band spectral response characteristics into a plurality of continuous band units according to a band division rule of hyperspectral remote sensing detection, and reserving complete spectral reflection and spectral absorption response characteristics of each band unit to form a single-class standard spectral band unit set; Extracting the detection full-band spectral response characteristics of a single geographic detection point in the binding set of the geographic point positions of the spectral characteristics, splitting the detection full-band spectral response characteristics into a plurality of band units corresponding to the standard spectrum according to the same band division rule to form a single-point detection spectrum band unit set, and synchronously splitting the detection spectrum and the band units of the standard spectrum; Extracting single standard band units in the single standard spectrum band unit set and single detection band units corresponding to the single point detection spectrum band unit set, analyzing the change trend of spectrum reflection and spectrum absorption in the two band units, and comparing the feature form matching of the two band units to form a feature comparison result of the single band unit; Aiming at the characteristic comparison result of the single-band unit with consistent spectrum variation trend, calculating the spectrum angle of the detection band unit and the corresponding standard band unit in the same band range, and if the spectrum angle is smaller than a set threshold value, judging that characteristic details are matched to form details confirmation content of the matched band unit; sequentially completing unit-by-unit comparison of all standard band units in the single-class standard spectrum band unit set and detection band units corresponding to the single-point detection spectrum band unit set, and collecting all band units which form characteristic form matching and detail matching to form a single-point matching band unit set; Extracting the band numbers of all the matched band units in the single-point position matched band unit set, and continuously connecting the matched band units according to the sequence of the band numbers to form continuous spectrum fragments corresponding to single geographic detection points and single type standard spectrums to form single-point position single type spectrum matched fragments; Counting the number of band units contained in the single-point single-class spectrum matching segment, and confirming that the spectrum matching segment is in a continuous band distribution form according to the continuous distribution state of the band units to form a single-point single-class spectrum matching segment validity confirmation content; Detecting full-band spectral response characteristics of single geographic detection points, and sequentially completing the band-by-band comparison and the spectrum matching fragment extraction with standard full-band spectral response characteristics of all geological alteration types in the normalized standard spectral feature set to form a single-point multi-class spectrum matching fragment set; Extracting all spectrum matching fragments in continuous wave band distribution form in the single-point multi-class spectrum matching fragment set, binding each spectrum matching fragment with a corresponding geological alteration type, and forming a single-point spectrum characteristic matching initial set; And sorting the single-point spectral feature matching initial set, determining the geological alteration type and the corresponding continuous spectral matching segment corresponding to the single geographic detection point, and removing the content of the spectral segments which are discontinuously distributed to form the single-point spectral feature matching set.
7. 7. The method for predicting a target area of a nickel ore based on hyperspectral remote sensing data analysis as claimed in claim 3, wherein the feature matching of the single alteration combination in the alteration cluster-shaped space engagement combination set and the ore-forming alteration sequential rule set is performed, and the alteration combination matched with the ore-forming alteration sequence of the nickel ore is screened out to form an ore-forming related alteration combination primary set, which comprises the following steps: Extracting all geological alteration types contained in a single alteration combination in the alteration cluster space connection combination, and sequencing the alteration types along a preset space reference direction according to the geometric center point coordinates of various alteration cluster distribution structures to form a single combination alteration space distribution sequential structure; Extracting a main sequence of geological changes corresponding to a nickel ore forming core stage in the ore forming change sequence rule set, and determining the occurrence sequence and type combination requirements of various geological changes in the main sequence of geological changes to form an ore forming core change main sequence set; carrying out overall feature matching on the single-combination alteration space distribution sequential structure and the ore-forming core alteration main sequential set, and analyzing overall matching of combination and arrangement sequence of geological alteration types in the single-combination alteration space distribution sequential structure and the ore-forming core alteration main sequential set to form a single-combination overall sequential matching result; Extracting the space extending direction of each geological alteration type cluster distribution structure in the single-combination alteration space distribution sequence aiming at the single-combination whole sequential matching result forming the whole matching property, analyzing the matching property of the space extending direction and the alteration evolution direction of the ore-forming core alteration main sequential set, and forming a single-combination space evolution direction matching result; Extracting full-band spectral response characteristics of each geological alteration type in the single-combination alteration spatial distribution sequential structure, analyzing evolution association relations among the characteristics, and comparing the evolution association relations with the matching of spectral evolution characteristics of the concentrated alteration of the ore-forming alteration sequential rules to form a single-combination spectral evolution characteristic matching result; For a single-combination integral sequential matching result without integral matching, extracting geological alteration type fragments which are locally matched with the ore-forming core alteration main sequential set in the single-combination alteration spatial distribution sequence, and analyzing whether the geological alteration type fragments are alteration combinations of the nickel ore-forming core stage or not to form a single-combination local sequential matching result; Extracting the rest geological alteration types in the alteration combination according to the single combination local sequential matching result forming the local matching, analyzing whether the corresponding alteration type is the associated alteration type of the nickel ore mineralization alteration or not, and whether the corresponding alteration type is in spatial distribution and evolution association with the core alteration type or not, so as to form a single combination associated alteration matching result; the method comprises the steps of collecting the alteration combinations forming overall matching, the alteration combinations forming local matching and accompanying alteration matching to form a candidate set of the ore-forming related alteration combinations; extracting cluster distribution structures corresponding to each alteration combination in the ore-forming related alteration combination candidate set, calculating the ratio of the total length of continuous distribution lines in each cluster distribution structure to the area surrounded by the lines, and judging that the alteration cluster distribution structures are combinations with poor space extensibility and rejecting if the ratio is lower than a set threshold value, so as to form an ore-forming related alteration combination screening set; And taking the ore-forming related alteration combination screening set as an ore-forming related alteration combination primary screening set, and determining all geological alteration combinations matched with the ore-forming characteristics of the nickel ore in the region to be predicted.
8. 8. The method of claim 4, wherein the extracting spatially sequential extension structures of single mineforming changes in the mineforming changes sequential evolution body, identifying all preceding changed cluster-shaped distribution structures located before and all following changed cluster-shaped distribution structures located after a core changed cluster-shaped distribution structure corresponding to a core anchoring point, extracting a geospatial range covered by circumscribed polygons of the changed cluster-shaped distribution structures, and forming an initial changed combination spatial association range, comprises: Extracting a spatial sequential extension structure of a single ore forming alteration combination in the ore forming alteration sequential evolution body, extracting a core alteration cluster-shaped distribution structure corresponding to a core anchoring point, and obtaining the complete geospatial distribution coordinates of the cluster-shaped distribution structure to form a core alteration cluster-shaped space feature set; According to the sequence of the alteration sequences defined in the space sequence extending structure, extracting a geological alteration cluster distribution structure positioned in front of the core alteration cluster distribution structure in the sequence as a first-order precursor alteration cluster distribution structure, and obtaining geographic space distribution coordinates of the geological alteration cluster distribution structure to form a first-order precursor alteration space feature set; According to the sequence of the alteration sequences defined in the space sequence extending structure, extracting a geological alteration cluster distribution structure positioned behind the core alteration cluster distribution structure in the sequence as a first-stage subsequent alteration cluster distribution structure, and acquiring geographic space distribution coordinates of the geological alteration cluster distribution structure to form a first-stage subsequent alteration space feature set; Continuously extracting a multi-level precursor geological alteration cluster distribution structure corresponding to the first-level precursor alteration space feature set along the precursor evolution direction until an alteration type which is not directly related to nickel ore formation is extracted, and recording geographic space distribution coordinates of precursor alteration of each level and space distance of core alteration to form a multi-level precursor alteration space feature set; Continuously extracting multi-level post-processing geological alteration cluster distribution structures corresponding to the first-level post-processing alteration space feature sets along the post-processing evolution direction until alteration types which are not directly related to nickel ore formation are extracted, and recording geographic space distribution coordinates of post-processing alteration of each level and space distances of post-processing alteration of each level and core alteration to form the multi-level post-processing alteration space feature sets; Integrating the geospatial distribution coordinates of all the changed cluster-shaped distribution structures in the core changed cluster-shaped space feature set, the multi-stage preamble changed space feature set and the multi-stage follow-up changed space feature set, and drawing the whole space distribution profile of all the changed cluster-shaped distribution structures to form a changed combined global space profile map; taking the geographic space center of the core alteration cluster-shaped space feature set as an origin, and defining the minimum geographic space range which can contain all the ore-forming related alteration cluster-shaped distribution structures by combining the spatial distances between the precursor and the subsequent alteration cluster-shaped distribution structures of each level and the origin, so as to form an alteration combined space association range primary outline; extracting the geographic coordinate boundary of the primary contour of the etching combination space association range, and expanding the boundary of the primary contour according to the space extension edges of all the mining-related etching cluster-shaped distribution structures to ensure that the whole of all the mining-related etching cluster-shaped distribution structures is in an expanded range, thereby forming an etching combination space association range expansion contour; Extracting all the alteration cluster distribution structures in the initial outline of the alteration combination space association range, if the non-relevant alteration cluster distribution structures which are not consistent with the definition type of the ore-forming alteration sequential rule set are identified, removing the vertex closest to the geometric center of the non-relevant alteration cluster distribution structure from the vertexes of the circumscribed polygon, and reconstructing the convex polygon boundary to form an alteration combination space association range adjustment outline; And extracting the geospatial distribution range of the modification combination space association range adjustment profile and the corresponding geographic coordinate boundary to form an initial modification combination space association range, and determining all the mining-related modification cluster distribution structures contained in the geospatial distribution range.
9. 9. The method for predicting a nickel mine target area based on hyperspectral remote sensing data analysis according to claim 5, wherein the analyzing the geospatial range of the single-shaped space continuous range centralizing adjacent shaped space distribution structure, extracting continuous boundary coordinates of two shaped spaces with space adjacent or space overlapping, and performing fusion processing on the boundary coordinates of adjacent or overlapping parts to form a fused continuous geographic boundary comprises: extracting geographic space ranges of two adjacent forming space distribution structures in the single forming space continuous range set, obtaining complete continuous boundary geographic coordinates of the two adjacent forming space distribution structures, determining boundary coordinate points and space trend of the two space ranges, and forming a double forming space boundary coordinate set; Comparing geographic coordinate boundaries of two shaping spaces in the coordinate set of the boundary of the double shaping space, identifying coordinate segments with adjacent contact or coordinate areas with space overlapping of the boundary, marking specific geographic coordinates of adjacent or overlapping parts, and forming a double shaping space adjacent overlapping mark set; Extracting a pair of coordinate points closest to adjacent boundaries of two space ranges as connection points aiming at coordinate segments in adjacent contact with the adjacent boundaries of the double-shaped space adjacent overlapping mark set, and connecting the pair of coordinate points by using straight line segments to form adjacent boundary joint coordinate segments; Extracting all geographic coordinate points of two space ranges in the overlapping region aiming at the coordinate regions of the adjacent overlapping marks of the double forming space, calculating the density clustering center of the geographic coordinate points, and reserving coordinate points belonging to main density clusters to form a core coordinate set of the overlapping region; Integrating the core coordinate set of the overlapping region with the boundary coordinate points of the non-overlapping region in two space ranges, continuously connecting all coordinate points according to the natural trend of the geographic space, constructing a continuous boundary coordinate segment after the overlapping region is fused, and forming an overlapping boundary fusion coordinate segment; Integrating the adjacent boundary joint coordinate segments or the overlapping boundary fusion coordinate segments with boundary coordinates of non-adjacent and non-overlapping areas of two space ranges to form a preliminary continuous geographic boundary after two shaping space fusion to form a double shaping space fusion primary boundary set; Extracting all boundary coordinate points in the double-shaped space fusion primary boundary set, sequencing the boundary coordinate points in a clockwise or anticlockwise direction, interpolating or resampling the sequence of the sequenced coordinate points to enable the interval distance between adjacent coordinate points to be smaller than a set threshold value, and forming a double-shaped space fusion regular boundary set; Comparing the fused geographic boundary corresponding to the double-shaped space fusion regular boundary set with geographic space ranges of other surrounding shaped space distribution structures again, if adjacent or overlapped geographic space ranges still exist, continuing to extract boundary coordinates and performing fusion processing to form a multi-shaped space fusion boundary set; And (3) extracting all boundary coordinate points in the multi-shaped space fusion boundary set, sequencing the boundary coordinate points in a clockwise or anticlockwise direction, calculating an included angle formed by continuous three points, deleting intermediate points if the included angle is within a set threshold value range, or calculating the angle change of each vertex of the polygon, reserving the vertex with curvature change larger than the set threshold value to form a fusion space core boundary coordinate set, continuously linking the fusion space core boundary coordinate set according to the geographic space trend to form a fused continuous geographic boundary, and determining the complete coordinate trend and coverage of the geographic boundary.
10. 10. Nickel mine target area prediction system based on hyperspectral remote sensing data analysis, which is characterized by comprising: A processor; a machine-readable storage medium storing machine-executable instructions for the processor; Wherein the processor is configured to perform the method of nickel mine target zone prediction based on hyperspectral remote sensing data analysis of any one of claims 1 to 9 via execution of the machine executable instructions.

Description

Nickel ore target area prediction method and system based on hyperspectral remote sensing data analysis Technical Field The invention relates to the technical field of mineral resource exploration, in particular to a nickel mine target area prediction method and system based on hyperspectral remote sensing data analysis. Background In the field of mineral resource exploration, nickel ore is used as an important metal mineral, and efficient and accurate exploration has important significance for national economic development and resource safety guarantee. Traditional nickel ore exploration methods mainly rely on geological surveys, geophysical surveys, geochemical surveys and the like. Geological investigation is to infer the possible existence area of nickel ore by analyzing geological information such as rock, stratum and the like through field investigation, but the method is limited by investment of manpower and material resources, has limited investigation range and is not deep enough to know the deep geological condition. Geophysical exploration utilizes physical methods, such as gravity, magnetic methods, electric methods and the like, to detect physical property differences of underground geologic bodies to search for nickel ores, however, similarity of physical properties of different geologic bodies may exist, so that interpretation results have multiple solutions, and accuracy is difficult to guarantee. Geochemical exploration defines a nickel mine target area by analyzing the content of chemical elements in the media such as soil, rock, water system sediments and the like, but the method is easily influenced by surface environmental factors such as weathering, carrying and the like, so that element distribution is changed, and the reliability of exploration results is influenced. In recent years, hyperspectral remote sensing technology is gradually applied to the field of mineral resource exploration. The hyperspectral remote sensing data has rich spectral information, can reflect the fine spectral characteristic difference of the ground object, and provides a new means for mineral exploration. However, most of the existing mineral exploration methods based on hyperspectral remote sensing data simply use spectrum information for classification or identification, do not fully consider geospatial information contained in hyperspectral data and complex relations between geological alteration and mineralization, are difficult to accurately and comprehensively predict a nickel ore target area, and cannot meet the requirements of modern nickel ore exploration. Disclosure of Invention In view of the above-mentioned problems, in combination with the first aspect of the present invention, an embodiment of the present invention provides a method for predicting a target area of nickel ore based on hyperspectral remote sensing data analysis, the method comprising: Acquiring a hyperspectral remote sensing data body of a region to be predicted, wherein the hyperspectral remote sensing data body comprises full-band spectral response characteristics of each geographic detection point in the region to be predicted and geographic space spread characteristics of each geographic detection point, the full-band spectral response characteristics are continuous response characteristics of spectral reflection and spectral absorption formed by the geographic detection points in a hyperspectral remote sensing detection band range, and the geographic space spread characteristics are geographic coordinates of each geographic detection point and spatial adjacent and spatial extension relations among the points; Performing clustered construction treatment on the hyperspectral remote sensing data body to form a clustered body of the changed characteristic space, wherein the clustered body of the changed characteristic space is a continuous clustered distribution structure formed by various geological change types depending on geographical space spread characteristics, and the continuous clustered distribution structure binds full-band spectral response characteristics of the corresponding geological change types; Performing ore-forming alteration combination sequential spatial evolution processing on the alteration feature space cluster body to form an ore-forming alteration sequential evolution body, wherein the ore-forming alteration sequential evolution body is a spatial sequential extension structure formed by geological alteration combination matched with nickel ore-forming features according to the occurrence sequence of the ore-forming alterations, and the spatial sequential extension structure binds full-band spectral response feature evolution content of the corresponding alteration combination; Performing target region indication feature anchoring forming construction treatment on the ore forming alteration sequential evolution body to form a target region indication space forming body, wherein the target region indic