CN111157459A - Hyperspectral image processing method for mineral information extraction - Google Patents

Hyperspectral image processing method for mineral information extraction Download PDF

Info

Publication number
CN111157459A
CN111157459A CN201911327172.9A CN201911327172A CN111157459A CN 111157459 A CN111157459 A CN 111157459A CN 201911327172 A CN201911327172 A CN 201911327172A CN 111157459 A CN111157459 A CN 111157459A
Authority
CN
China
Prior art keywords
image
band image
pixel value
wave band
condition
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
CN201911327172.9A
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.)
Beijing Research Institute of Uranium Geology
Original Assignee
Beijing Research Institute of Uranium Geology
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 Beijing Research Institute of Uranium Geology filed Critical Beijing Research Institute of Uranium Geology
Priority to CN201911327172.9A priority Critical patent/CN111157459A/en
Publication of CN111157459A publication Critical patent/CN111157459A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

The invention relates to a hyperspectral image processing method for information extraction of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite, in particular to a hyperspectral image processing method for information extraction of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite, which reduces the influence of other ground objects or noise and reduces the data processing amount. The invention comprises the following steps in sequence: step 1, acquiring a hyperspectral image; preprocessing the hyperspectral image, and performing atmospheric correction to obtain image data of ground reflectivity; step 2, resampling the image data of the ground reflectivity, and extracting images of different wave bands; step 3, comparing pixel values of images of different wave bands, and marking a comparison result; and 4, determining abundance values of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite according to the step 3.

Description

Hyperspectral image processing method for mineral information extraction
Technical Field
The invention relates to a hyperspectral image processing method for information extraction of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite, in particular to a hyperspectral image processing method for information extraction of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite, which reduces the influence of other ground objects or noise and reduces the data processing amount.
Background
The existing extraction method of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite information of hyperspectral remote sensing images mainly adopts full-band spectrum matching or partial continuous-band spectrum matching, specific algorithms include spectral angles, mixed demodulation filtering and the like, and because the material composition of the ground surface is rarely composed of single mineral, the methods are easily influenced by the spectrum or noise of other ground objects in the information extraction process, and the extraction accuracy is relatively low. Secondly, the existing spectrum extraction method has many manual operation steps, and the manual judgment error is increased. Thirdly, the hyperspectral data has more wave bands and large data volume, the existing method has long processing time, and the data processing speed and the application scale are reduced. Therefore, how to reduce the influence of other ground objects or noise, manual operation steps and data processing amount in the process of extracting information of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite becomes one of the leading edges of the current hyperspectral remote sensing image processing.
Disclosure of Invention
The invention aims to provide a hyperspectral image processing method for information extraction of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite, which can reduce the influence of other ground objects or noise and reduce the data volume to be processed. The hyperspectral image is resampled, specific wave bands are extracted, a series of judgments and calculations are carried out, and abundance values of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite in different areas in the image range are calculated.
The technical scheme adopted by the invention is as follows:
a hyperspectral image processing method for mineral information extraction sequentially comprises the following steps: step 1, acquiring a hyperspectral image; preprocessing the hyperspectral image, and performing atmospheric correction to obtain image data of ground reflectivity; step 2, resampling the image data of the ground reflectivity, and extracting images of different wave bands; step 3, comparing pixel values of images of different wave bands, and marking a comparison result; and 4, determining abundance values of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite according to the step 3.
In step 2, the wavelength range of the extracted image includes 980nm, 1040nm, 1145nm, 1175nm, 1325nm, 1340nm, 1355nm, 1415nm, 1430nm, 1445nm, 1460nm, 1490nm, 1520nm, 1535nm, 1550nm, 1595nm, 1625nm, 1700nm, 1715nm, 1745nm, 1760nm, 1775nm, 1790nm, 1805nm, 1820nm, 1940nm, 1940nm, 1955nm, 1985nm, 2015nm, 2045nm, 2060nm, 2075nm, 2090nm, 2075nm, 2120nm, 2135nm, 2150nm, 2155 nm, 2180nm, 2195nm, 2210nm, 2225nm, 2230nm, 2240nm, 2250 nm, 2270nm, 2285nm, 2300nm, 2315nm, 2330nm, 2425 nm, 2360nm, 2390nm, 235 nm, 230nm, 235 nm, 230nm and 235 nm.
The step 3 includes the steps of,
when each pixel value of the 1145nm band image is greater than the corresponding pixel value of the 1175nm band image, judging that the condition I is met; obtaining a result image H1 of subtracting the 1175nm wave band image from the 1145nm wave band image;
when each pixel value of the 1325nm band image is larger than the pixel value corresponding to the 1445nm band image, judging that a second condition is met; obtaining a result image H2 of subtracting the 1445nm wave band image from the 1325nm wave band image;
when each pixel value of the 1340nm waveband image is larger than the pixel value corresponding to the 1355nm waveband image, judging that the condition three is met; obtaining a result image H3 of subtracting 1455nm wave band image from 1340nm wave band image;
when each pixel value of the 1355nm band image is larger than the corresponding pixel value of the 980nm band image, judging that the condition four is met; obtaining a result image H4 of subtracting the 980nm wave band image from the 1355nm wave band image;
when each pixel value of the 1355nm band image is larger than the corresponding pixel value of the 1415nm band image, judging that the condition five is met; obtaining a result image H5 of subtracting the 1415nm wave band image from the 1355nm wave band image;
when each pixel value of the 1355nm band image is larger than the corresponding pixel value of the 1445nm band image, judging that a sixth condition is met; obtaining a result image H5 of subtracting the 1445nm wave band image from the 1355nm wave band image;
when each pixel value of the 1415nm band image is greater than the corresponding pixel value of the 1430nm band image, judging that the condition seven is met; obtaining a result image H7 obtained by subtracting the 1430nm wave band image from the 1415nm wave band image;
when each pixel value of the 1415nm band image is greater than the corresponding pixel value of the 1460nm band image, judging that the condition eight is met; obtaining a result image H8 of subtracting the 1460nm wave band image from the 1415nm wave band image;
when each pixel value of the 1445nm band image is larger than the corresponding pixel value of the 1415nm band image, judging that the condition nine is met; obtaining a result image H9 of subtracting the 1415nm wave band image from the 1445nm wave band image;
when each pixel value of the 1460nm band image is larger than the pixel value corresponding to the 1445nm band image, judging that the condition ten is met; obtaining a 1460nm wave band image minus a 1445nm wave band image as a result image H10;
when each pixel value of the 1490nm band image is larger than the pixel value corresponding to the 1445nm band image, judging that the condition eleven is met; obtaining a result image H11 of subtracting the 1445nm wave band image from the 1490nm wave band image;
when each pixel value of the 1490nm band image is larger than the pixel value corresponding to the 1520nm band image, judging that the condition twelve is met; obtaining a result image H12 of subtracting the 1520nm wave band image from the 1490nm wave band image;
when each pixel value of the 1520nm band image is larger than the pixel value corresponding to the 1490nm band image, judging that the condition thirteen is met; obtaining a result image H13 of subtracting the 1490nm wave band image from the 1520nm wave band image;
when each pixel value of the 1535nm wave band image is larger than the corresponding pixel value of the 1445nm wave band image, judging that the condition is fourteen; obtaining a result image H14 of subtracting the 1445nm wave band image from the 1535nm wave band image;
when each pixel value of the 1595nm band image is larger than the corresponding pixel value of the 1550nm band image, judging that the condition fifteen is met; obtaining a result image H15 obtained by subtracting the 1550nm waveband image from the 1595nm waveband image;
when each pixel value of the 1625nm waveband image is larger than the corresponding pixel value of the 1550nm waveband image, judging that the condition is sixteen; obtaining a result image H16 obtained by subtracting the 1550nm waveband image from the 1625nm waveband image;
when each pixel value of the 1700nm band image is larger than the corresponding pixel value of the 1040nm band image, judging that the condition is seventeen; obtaining a result image H17 of subtracting the 1040nm waveband image from the 1700nm waveband image;
when each pixel value of the 1700nm waveband image is larger than the corresponding pixel value of the 1745nm waveband image, judging that the condition is met by eighteen; obtaining a result image H18 of subtracting the 1745nm band image from the 1700nm band image;
when each pixel value of the 1715nm waveband image is larger than the pixel value corresponding to the 1760nm waveband image, judging that the condition is met; obtaining a result image H19 obtained by subtracting the 1760nm wave band image from the 1715nm wave band image;
when each pixel value of the 1775nm wave band image is larger than the pixel value corresponding to the 1940nm wave band image, judging that the condition twenty is met; obtaining a result image H20 of subtracting the 1940nm wave band image from the 1775nm wave band image;
when each pixel value of the 1790nm band image is larger than the pixel value corresponding to the 2015nm band image, judging that the condition twenty-one is met; obtaining a result image H21 of subtracting the 2015nm waveband image from the 1790nm waveband image;
when each pixel value of the 1805nm wave band image is larger than the corresponding pixel value of the 1775nm wave band image, judging that the condition twenty-two is met; obtaining a 1805nm wave band image minus a 1775nm wave band image H22;
when each pixel value of the 1820nm band image is larger than the corresponding pixel value of the 1745nm band image, judging that the condition twenty-three is met; obtaining a result image H23 of subtracting the 1745nm band image from the 1820nm band image;
when each pixel value of the 1820nm wave band image is larger than the pixel value corresponding to the 1940nm wave band image, judging that the condition twenty-four is met; obtaining a result image H24 of subtracting the 1940nm wave band image from the 1820nm wave band image;
when each pixel value of the 1820nm band image is larger than the pixel value corresponding to the 1955nm band image, judging that the condition is satisfied by twenty-five; obtaining a 1820nm wave band image minus the 1955nm wave band image as a result H25;
when each pixel value of the 2045nm band image is greater than the corresponding pixel value of the 1985nm band image, judging that the condition twenty-six is met; obtaining a 2045nm wave band image and subtracting a 1985nm wave band image to obtain a result image H26;
when each pixel value of the 2045nm band image is larger than the pixel value corresponding to the 2015nm band image, judging that the condition is met by twenty-seven; obtaining a result image H27 of subtracting the 2015nm waveband image from the 2045nm waveband image;
when each pixel value of the 2045nm band image is larger than the corresponding pixel value of the 2075nm band image, judging that the condition is satisfied twenty-eight; obtaining a result image H28 of subtracting the 2075nm wave band image from the 2045nm wave band image;
when each pixel value of the 2060nm wave band image is larger than the pixel value corresponding to the 1940nm wave band image, judging that the condition is satisfied twenty-nine; obtaining a result image H29 of subtracting the 1940nm wave band image from the 2060nm wave band image;
when each pixel value of the 2060nm wave band image is larger than the pixel value corresponding to the 2210nm wave band image, judging that the condition is thirty; obtaining a result image H30 of subtracting the 2210nm wave band image from the 2060nm wave band image;
when each pixel value of the 2075nm band image is larger than the pixel value corresponding to the 2120nm band image, judging that the condition is satisfied thirty-one; obtaining a result image H31 of subtracting the 2120nm wave band image from the 2075nm wave band image;
when each pixel value of the 2105nm wave band image is larger than the pixel value corresponding to the 2090nm wave band image, judging that the condition is satisfied thirty-two; obtaining a result image H32 obtained by subtracting the 2090nm wave band image from the 2105nm wave band image;
when each pixel value of the 2105nm wave band image is larger than the corresponding pixel value of the 2165nm wave band image, judging that the condition is satisfied; obtaining a result image H33 obtained by subtracting the 2165nm wave band image from the 2105nm wave band image;
when each pixel value of the 2120nm band image is larger than the corresponding pixel value of the 2165nm band image, judging that the thirty-four conditions are met; obtaining a result image H34 obtained by subtracting the 2165nm wave band image from the 2120nm wave band image;
when each pixel value of the 2135nm wave band image is larger than the corresponding pixel value of the 2165nm wave band image, judging that the condition is satisfied thirty-five; obtaining a result image H35 obtained by subtracting the 2165nm wave band image from the 2135nm wave band image;
when each pixel value of the 2135nm wave band image is larger than the pixel value corresponding to the 2180nm wave band image, judging that the condition is satisfied thirty-six; obtaining a result image H36 of subtracting the 2180nm waveband image from the 2135nm waveband image;
when each pixel value of the 2150nm band image is larger than the corresponding pixel value of the 2255nm band image, judging that the condition thirty-seven is met; obtaining a result image H37 obtained by subtracting the 2255nm wave band image from the 2150nm wave band image;
when each pixel value of the 2195nm wave band image is larger than the corresponding pixel value of the 2165nm wave band image, judging that the condition is satisfied thirty-eight; obtaining a result image H38 of subtracting the 2165nm wave band image from the 2195nm wave band image;
when each pixel value of the 2195nm band image is larger than the corresponding pixel value of the 2180nm band image, judging that the condition is satisfied by thirty-nine; obtaining a result image H39 of subtracting the 2180nm waveband image from the 2195nm waveband image;
when each pixel value of the 2195nm band image is larger than the pixel value corresponding to the 2210nm band image, judging that the condition is met by forty; obtaining a result image H40 of subtracting the 2210nm wave band image from the 2195nm wave band image;
when each pixel value of the 2195nm waveband image is larger than the corresponding pixel value of the 2225nm waveband image, judging that the condition is forty one; obtaining a result image H41 obtained by subtracting the 2225nm waveband image from the 2195nm waveband image;
when each pixel value of the 2225nm band image is greater than the corresponding pixel value of the 2120nm band image, judging that the condition is satisfied forty-two; obtaining a result image H42 of subtracting the 2120nm wave band image from the 2225nm wave band image;
when each pixel value of the 2225nm band image is larger than the corresponding pixel value of the 2255nm band image, judging that the condition is satisfied forty-three; obtaining a result image H43 obtained by subtracting the 2255nm wave band image from the 2225nm wave band image;
when each pixel value of the 2225nm wave band image is larger than the corresponding pixel value of the 2300nm wave band image, judging that the condition fourteen is met; obtaining a result image H44 of subtracting the 2300nm wave band image from the 2225nm wave band image;
when each pixel value of the 2225nm band image is larger than the corresponding pixel value of the 2360nm band image, judging that the condition forty-five is met; obtaining a result image H45 obtained by subtracting the 2360nm wave band image from the 2225nm wave band image;
when each pixel value of the 2240nm waveband image is larger than the pixel value corresponding to the 2210nm waveband image, judging that the condition is satisfied forty-six; obtaining a result image H46 of subtracting the 2210nm wave band image from the 2240nm wave band image;
when each pixel value of the 2240nm waveband image is larger than the corresponding pixel value of the 2255nm waveband image, judging that the condition is forty-seven; obtaining a result image H47 obtained by subtracting the 2255nm wave band image from the 2240nm wave band image;
when each pixel value of the 2240nm waveband image is larger than the corresponding pixel value of the 2315nm waveband image, judging that the condition is satisfied forty-eight; obtaining a result image H48 obtained by subtracting the 2315nm wave band image from the 2240nm wave band image;
when each pixel value of the 2255nm band image is greater than the pixel value corresponding to the 2210nm band image, judging that the condition forty-nine is met; obtaining a result image H49 of subtracting the 2210nm wave band image from the 2255nm wave band image;
when each pixel value of the 2270nm band image is larger than the corresponding pixel value of the 2300nm band image, judging that the condition is met; obtaining a result image H50 of subtracting the 2300nm wave band image from the 2270nm wave band image;
when each pixel value of the 2285nm band image is larger than the corresponding pixel value of the 2300nm band image, judging that the condition is met; obtaining a result image H51 of subtracting the 2300nm wave band image from the 2285nm wave band image;
when each pixel value of the 2285nm band image is larger than the corresponding pixel value of the 2315nm band image, judging that the condition is met by fifty-two; obtaining a result image H52 of subtracting the 2315nm wave band image from the 2285nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2405nm band image, judging that the condition is satisfied; obtaining a result image H53 obtained by subtracting the 2405nm wave band image from the 2330nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2255nm band image, judging that the condition fifty-four is met; obtaining a result image H54 obtained by subtracting the 2255nm wave band image from the 2330nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2315nm band image, judging that the condition fifty-five is met; obtaining a result image H55 obtained by subtracting the 2315nm wave band image from the 2330nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2390nm band image, judging that the condition is met fifty-six; obtaining a result image H56 obtained by subtracting the 2390nm wave band image from the 2330nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2405nm band image, judging that the condition is met; obtaining a result image H57 obtained by subtracting the 2405nm wave band image from the 2330nm wave band image;
when each pixel value of the 2345nm wave band image is larger than the corresponding pixel value of the 2315nm wave band image, judging that the condition is satisfied fifty-eight; obtaining a result image H58 of subtracting the 2315nm wave band image from the 2345nm wave band image;
when each pixel value of the 2345nm wave band image is larger than the pixel value corresponding to the 2330nm wave band image, judging that the condition is met by fifty-nine; obtaining a result image H59 of subtracting the 2330nm waveband image from the 2345nm waveband image;
when each pixel value of the 2360nm band image is larger than the corresponding pixel value of the 2375nm band image, judging that the condition sixty is met; obtaining a result image H60 obtained by subtracting the 2375nm wave band image from the 2360nm wave band image;
when each pixel value of the 2360nm band image is larger than the corresponding pixel value of the 2390nm band image, judging that the condition is sixty-one satisfied; obtaining a result image H61 obtained by subtracting the 2390nm wave band image from the 2360nm wave band image;
when each pixel value of the 2360nm band image is larger than the corresponding pixel value of the 2405nm band image, judging that the condition is satisfied sixty-two; obtaining a result image H62 obtained by subtracting the 2405nm wave band image from the 2360nm wave band image;
when each pixel value of the 2375nm band image is larger than the corresponding pixel value of the 2390nm band image, judging that the condition sixty-three is met; obtaining a result image H63 obtained by subtracting the 2390nm wave band image from the 2375nm wave band image;
when each pixel value of the 2405nm waveband image is larger than the corresponding pixel value of the 2360nm waveband image, judging that the condition sixty-four is met; obtaining a result image H64 obtained by subtracting the 2360nm wave band image from the 2405nm wave band image;
when each pixel value of the 2405nm band image is larger than the corresponding pixel value of the 2390nm band image, judging that the condition sixty-five is met; obtaining a result image H65 obtained by subtracting the 2390nm wave band image from the 2405nm wave band image;
when each pixel value of the 2420nm wave band image is larger than the corresponding pixel value of the 2390nm wave band image, judging that the condition is sixty-six met; obtaining a result image H66 obtained by subtracting the 2390nm wave band image from the 2420nm wave band image;
in the step 4, pixels satisfying conditions five, nine, thirty-six, thirty-nine, forty-six, forty-seven, fifty-one, fifty-six, and sixty-five are selected at the same time, and a sum H69 of H5, H9, H36, H39, H40, H46, H47, H51, H56, and H65 in the above pixel range is obtained, and a higher H69 value represents a higher dickite abundance value in the pixel.
In the step 4, pixels which simultaneously satisfy conditions four, eight, fifteen, twenty-one, thirty-one, forty-two, forty-three, fifty-four and fifty-seven are selected, a sum H70 of sums of H4, H8, H15, H21, H31, H42, H43, H54 and H57 in the pixel range is obtained, and the higher the value of H70 is, the higher the abundance value of jarosite in the pixels is.
In the step 4, pixels satisfying the conditions three, twenty-five, twenty-six, forty-one, forty-four, fifty-nine and sixty are selected at the same time, and the sum H71 of the sums of H3, H25, H26, H41, H44, H59 and H60 in the pixel range is obtained. The higher the value of H71, the higher the value of the richness of the lizardtail within the representative picture element.
In the step 4, selecting the pixels which simultaneously meet the conditions of fourteen, twenty-four, twenty-seven, thirty-five, forty-five and sixty-four, and obtaining the sum H72 of the sums of H14, H24, H27, H35, H45 and H64 in the pixel range, wherein the higher the H72 value is, the higher the abundance value of rhodochrosite in the pixels is.
In the step 4, pixels which simultaneously satisfy the conditions of the first condition, the second condition, the tenth condition, the eighteen condition, the twenty-third condition, the twenty-fourth condition, the twenty-ninth condition, the thirty condition and the fifty-third condition are selected, a sum H73 of the sums of H1, H2, H10, H18, H23, H24, H29, H30 and H53 in the pixel range is obtained, and the higher the value of H73 is, the higher the gypsum abundance value in the pixel is represented.
In the step 4, selecting the pixels which simultaneously satisfy the conditions six, eleven, twelve, sixteen, twenty and thirty-seven, and obtaining the sum H74 of the sum of H6, H11, H12, H16, H20 and H37 in the pixel range, wherein the higher the H74 value is, the higher the abundance value of the actinolite in the pixels is.
In the step 4, pixels which simultaneously satisfy the conditions of twenty-eight, thirty-two, thirty-three, thirty-eight, forty-eight, fifty-five and sixty-three are selected to obtain the sum H75 of the sums of H28, H32, H33, H38, H48, H55 and H63 in the pixel range, and the higher the value of H75, the higher the pyrophyllite abundance value in the pixel is.
In the step 4, pixels satisfying the conditions seven, thirteen, nineteen, twenty-two, twenty-four, thirty-nine, forty-nine, fifty-two, and sixty-two are selected at the same time, and the sum H76 of the sums of H7, H13, H19, H22, H24, H34, H39, H49, H52, and H62 in the pixel range is obtained. The higher the H76 value, the higher the alunite abundance value in the representative pixel.
The method of the invention is computer automatic extraction, reduces the operation steps of principal component transformation, end member spectrum selection and the like, and can greatly improve the operation speed. Because most of wave bands which are not relevant to information extraction are removed, the interference of other substances or noise on the spectrum of the wave bands is reduced, and the accuracy of information extraction is improved. The method has good effect and significance in quickly extracting information of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite in hyperspectral image data.
Detailed Description
The invention comprises the following steps in sequence:
firstly, acquiring a hyperspectral image; preprocessing the hyperspectral image, and performing atmospheric correction to obtain image data of ground reflectivity;
secondly, resampling the image data of ground reflectance to extract images with wavebands of 980nm, 1040nm, 1145nm, 1175nm, 1325nm, 1340nm, 1355nm, 1415nm, 1430nm, 1445nm, 1460nm, 1490nm, 1520nm, 1535nm, 1550nm, 1595nm, 1625nm, 1700nm, 1715nm, 1745nm, 1760nm, 1775nm, 1790nm, 1805nm, 1820nm, 1940nm, 1940nm, 1955nm, 1985nm, 2015nm, 2045nm, 2060nm, 2075nm, 2090nm, 2105nm, 2120nm, 2135nm, 2150nm, 2165nm, 2180nm, 2195nm, 2210nm, 2225nm, 2230nm, 2240nm, 2255nm, 225nm, 2270nm, 2285nm, 2300nm, 2315nm, 2420nm, 2425 nm, 2360nm, 2390nm, 235 nm, and 235 nm;
thirdly, when each pixel value of the 1145nm waveband image is larger than the corresponding pixel value of the 1175nm waveband image, judging that a first condition is met; obtaining a result image H1 of subtracting the 1175nm wave band image from the 1145nm wave band image;
fourthly, when each pixel value of the 1325nm wave band image is larger than the pixel value corresponding to the 1445nm wave band image, judging that a second condition is met; obtaining a result image H2 of subtracting the 1445nm wave band image from the 1325nm wave band image;
fifthly, when each pixel value of the 1340nm wave band image is larger than the pixel value corresponding to the 1355nm wave band image, judging that a third condition is met; obtaining a result image H3 of subtracting 1455nm wave band image from 1340nm wave band image;
sixthly, when each pixel value of the 1355 nm-band image is larger than the corresponding pixel value of the 980 nm-band image, judging that a fourth condition is met; obtaining a result image H4 of subtracting the 980nm wave band image from the 1355nm wave band image;
seventhly, when each pixel value of the 1355nm band image is larger than the corresponding pixel value of the 1415nm band image, judging that a fifth condition is met; obtaining a result image H5 of subtracting the 1415nm wave band image from the 1355nm wave band image;
eighthly, when each pixel value of the 1355 nm-band image is larger than the corresponding pixel value of the 1445 nm-band image, judging that a sixth condition is met; obtaining a result image H5 of subtracting the 1445nm wave band image from the 1355nm wave band image;
ninthly, when each pixel value of the 1415nm band image is larger than the corresponding pixel value of the 1430nm band image, judging that a seventh condition is met; obtaining a result image H7 obtained by subtracting the 1430nm wave band image from the 1415nm wave band image;
step ten, when each pixel value of the 1415nm band image is larger than the corresponding pixel value of the 1460nm band image, judging that the condition eight is met; obtaining a result image H8 of subtracting the 1460nm wave band image from the 1415nm wave band image;
eleventh, when each pixel value of the 1445nm band image is greater than the pixel value corresponding to the 1415nm band image, determining that the condition nine is satisfied; obtaining a result image H9 of subtracting the 1415nm wave band image from the 1445nm wave band image;
twelfth, when each pixel value of the 1460nm band image is larger than the pixel value corresponding to the 1445nm band image, judging that the condition ten is met; obtaining a 1460nm wave band image minus a 1445nm wave band image as a result image H10;
thirteenth, when each pixel value of the 1490nm band image is larger than the corresponding pixel value of the 1445nm band image, judging that the condition eleven is met; obtaining a result image H11 of subtracting the 1445nm wave band image from the 1490nm wave band image;
fourteenth, when each pixel value of the 1490nm band image is greater than the pixel value corresponding to the 1520nm band image, judging that the condition twelve is met; obtaining a result image H12 of subtracting the 1520nm wave band image from the 1490nm wave band image;
fifteenth step, when each pixel value of the 1520nm wave band image is larger than the pixel value corresponding to the 1490nm wave band image, judging that the condition thirteen is met; obtaining a result image H13 of subtracting the 1490nm wave band image from the 1520nm wave band image;
sixthly, when each pixel value of the 1535nm wave band image is larger than the pixel value corresponding to the 1445nm wave band image, judging that the condition fourteen is met; obtaining a result image H14 of subtracting the 1445nm wave band image from the 1535nm wave band image;
seventeenth, when each pixel value of the 1595nm band image is larger than the corresponding pixel value of the 1550nm band image, judging that a condition fifteen is met; obtaining a result image H15 obtained by subtracting the 1550nm waveband image from the 1595nm waveband image;
eighteen, when each pixel value of the 1625nm band image is larger than the corresponding pixel value of the 1550nm band image, judging that the condition is sixteen; obtaining a result image H16 obtained by subtracting the 1550nm waveband image from the 1625nm waveband image;
nineteenth step, when each pixel value of the 1700nm waveband image is larger than the corresponding pixel value of the 1040nm waveband image, judging that the condition is seventeen; obtaining a result image H17 of subtracting the 1040nm waveband image from the 1700nm waveband image;
twentieth, when each pixel value of the 1700nm band image is larger than the corresponding pixel value of the 1745nm band image, judging that the condition is met eighteen; obtaining a result image H18 of subtracting the 1745nm band image from the 1700nm band image;
twenty-first step, when each pixel value of the 1715nm wave band image is larger than the pixel value corresponding to the 1760nm wave band image, judging that the condition is met; obtaining a result image H19 obtained by subtracting the 1760nm wave band image from the 1715nm wave band image;
twenty-second step, when each pixel value of the 1775nm wave band image is larger than the pixel value corresponding to the 1940nm wave band image, judging that the condition is met twenty; obtaining a result image H20 of subtracting the 1940nm wave band image from the 1775nm wave band image;
a twenty-first step, when each pixel value of the 1790nm band image is larger than the pixel value corresponding to the 2015nm band image, judging that the condition is met; obtaining a result image H21 of subtracting the 2015nm waveband image from the 1790nm waveband image;
twenty-fourth, when each pixel value of the 1805nm wave band image is larger than the corresponding pixel value of the 1775nm wave band image, judging that the condition twenty-two is met; obtaining a 1805nm wave band image minus a 1775nm wave band image H22;
twenty-fifth step, when each pixel value of the 1820nm wave band image is larger than the corresponding pixel value of the 1745nm wave band image, judging that the condition twenty-third is met; obtaining a result image H23 of subtracting the 1745nm band image from the 1820nm band image;
twenty-sixth step, when each pixel value of the 1820nm wave band image is larger than the pixel value corresponding to the 1940nm wave band image, judging that the condition twenty-four is met; obtaining a result image H24 of subtracting the 1940nm wave band image from the 1820nm wave band image;
twenty-seventh step, when each pixel value of the 1820nm wave band image is larger than the pixel value corresponding to the 1955nm wave band image, judging that the condition is satisfied twenty-fifth; obtaining a 1820nm wave band image minus the 1955nm wave band image as a result H25;
twenty-eighth step, when each pixel value of the 2045nm wave band image is larger than the corresponding pixel value of the 1985nm wave band image, judging that the condition is satisfied with twenty-six; obtaining a 2045nm wave band image and subtracting a 1985nm wave band image to obtain a result image H26;
a twenty-ninth step of judging that the condition is satisfied when each pixel value of the 2045nm band image is greater than the pixel value corresponding to the 2015nm band image; obtaining a result image H27 of subtracting the 2015nm waveband image from the 2045nm waveband image;
thirty-eighth, when each pixel value of the 2045nm band image is larger than the corresponding pixel value of the 2075nm band image, judging that the condition is satisfied; obtaining a result image H28 of subtracting the 2075nm wave band image from the 2045nm wave band image;
thirty-ninth step, when each pixel value of the 2060nm wave band image is larger than the pixel value corresponding to the 1940nm wave band image, judging that the condition is satisfied; obtaining a result image H29 of subtracting the 1940nm wave band image from the 2060nm wave band image;
thirty-second step, when each pixel value of the 2060nm wave band image is greater than the pixel value corresponding to the 2210nm wave band image, judging that the condition is satisfied with thirty; obtaining a result image H30 of subtracting the 2210nm wave band image from the 2060nm wave band image;
thirty-first, when each pixel value of the 2075nm band image is larger than the pixel value corresponding to the 2120nm band image, judging that the condition is satisfied thirty-first; obtaining a result image H31 of subtracting the 2120nm wave band image from the 2075nm wave band image;
thirty-fourth step, when each pixel value of the 2105nm wave band image is larger than the pixel value corresponding to the 2090nm wave band image, judging that the condition is satisfied thirty-two; obtaining a result image H32 obtained by subtracting the 2090nm wave band image from the 2105nm wave band image;
thirty-fifth step, when each pixel value of the 2105nm wave band image is larger than the pixel value corresponding to the 2165nm wave band image, judging that the condition is satisfied; obtaining a result image H33 obtained by subtracting the 2165nm wave band image from the 2105nm wave band image;
sixthly, when each pixel value of the 2120nm band image is larger than the pixel value corresponding to the 2165nm band image, judging that the thirty-four conditions are met; obtaining a result image H34 obtained by subtracting the 2165nm wave band image from the 2120nm wave band image;
thirty-seventh step, when each pixel value of the 2135nm wave band image is larger than the corresponding pixel value of the 2165nm wave band image, judging that the condition is satisfied thirty-five; obtaining a result image H35 obtained by subtracting the 2165nm wave band image from the 2135nm wave band image;
thirty-eighth step, when each pixel value of the 2135nm wave band image is larger than the pixel value corresponding to the 2180nm wave band image, determining that the condition is satisfied thirty-six; obtaining a result image H36 of subtracting the 2180nm waveband image from the 2135nm waveband image;
thirty-ninth step, when each pixel value of the 2150nm band image is larger than the pixel value corresponding to the 2255nm band image, judging that the condition is satisfied thirty-seven; obtaining a result image H37 obtained by subtracting the 2255nm wave band image from the 2150nm wave band image;
forty step, when each pixel value of the 2195nm wave band image is larger than the pixel value corresponding to the 2165nm wave band image, judging that the condition is satisfied thirty-eight; obtaining a result image H38 of subtracting the 2165nm wave band image from the 2195nm wave band image;
forty-ninth step, when each pixel value of the 2195nm waveband image is larger than the corresponding pixel value of the 2180nm waveband image, judging that the condition is satisfied; obtaining a result image H39 of subtracting the 2180nm waveband image from the 2195nm waveband image;
step forty, when each pixel value of the 2195nm wave band image is larger than the pixel value corresponding to the 2210nm wave band image, judging that the condition forty is met; obtaining a result image H40 of subtracting the 2210nm wave band image from the 2195nm wave band image;
step forty, when each pixel value of the 2195nm waveband image is larger than the corresponding pixel value of the 2225nm waveband image, judging that the condition forty one is met; obtaining a result image H41 obtained by subtracting the 2225nm waveband image from the 2195nm waveband image;
forty-fourth step, when each pixel value of the 2225nm waveband image is larger than the corresponding pixel value of the 2120nm waveband image, judging that the condition is met by forty-two; obtaining a result image H42 of subtracting the 2120nm wave band image from the 2225nm wave band image;
forty-fifth step, when each pixel value of the 2225nm waveband image is larger than the corresponding pixel value of the 2255nm waveband image, judging that the condition is met by forty-three; obtaining a result image H43 obtained by subtracting the 2255nm wave band image from the 2225nm wave band image;
sixthly, when each pixel value of the 2225nm waveband image is larger than the corresponding pixel value of the 2300nm waveband image, judging that the condition is forty-four satisfied; obtaining a result image H44 of subtracting the 2300nm wave band image from the 2225nm wave band image;
forty-seventh step, when each pixel value of the 2225nm waveband image is larger than the corresponding pixel value of the 2360nm waveband image, judging that the condition is forty-five satisfied; obtaining a result image H45 obtained by subtracting the 2360nm wave band image from the 2225nm wave band image;
forty-eight step, when each pixel value of the 2240nm waveband image is larger than the pixel value corresponding to the 2210nm waveband image, judging that the condition is satisfied forty-six; obtaining a result image H46 of subtracting the 2210nm wave band image from the 2240nm wave band image;
forty-ninth step, when each pixel value of the 2240nm waveband image is larger than the corresponding pixel value of the 2255nm waveband image, judging that the condition is forty-seven; obtaining a result image H47 obtained by subtracting the 2255nm wave band image from the 2240nm wave band image;
fifthly, when each pixel value of the 2240nm waveband image is larger than the pixel value corresponding to the 2315nm waveband image, judging that the condition is satisfied forty-eight; obtaining a result image H48 obtained by subtracting the 2315nm wave band image from the 2240nm wave band image;
fifthly, when each pixel value of the 2255nm band image is greater than the pixel value corresponding to the 2210nm band image, judging that the condition forty-nine is met; obtaining a result image H49 of subtracting the 2210nm wave band image from the 2255nm wave band image;
fifthly, when each pixel value of the 2270nm band image is larger than the corresponding pixel value of the 2300nm band image, judging that the condition is fifty; obtaining a result image H50 of subtracting the 2300nm wave band image from the 2270nm wave band image;
fiftieth, when each pixel value of the 2285nm band image is larger than the corresponding pixel value of the 2300nm band image, judging that the condition is met; obtaining a result image H51 of subtracting the 2300nm wave band image from the 2285nm wave band image;
fiftieth, when each pixel value of the 2285nm band image is larger than the corresponding pixel value of the 2315nm band image, judging that the condition is met by fifty-two; obtaining a result image H52 of subtracting the 2315nm wave band image from the 2285nm wave band image;
fiftieth, when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2405nm band image, judging that the condition is met; obtaining a result image H53 obtained by subtracting the 2405nm wave band image from the 2330nm wave band image;
fiftieth, when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2255nm band image, judging that the condition fifty-four is met; obtaining a result image H54 obtained by subtracting the 2255nm wave band image from the 2330nm wave band image;
fiftieth and seventh steps, when each pixel value of the 2330nm wave band image is larger than the corresponding pixel value of the 2315nm wave band image, judging that the condition is met by fifty-five; obtaining a result image H55 obtained by subtracting the 2315nm wave band image from the 2330nm wave band image;
fiftieth and eighth step, when each pixel value of the 2330nm wave band image is larger than the corresponding pixel value of the 2390nm wave band image, judging that the condition is met by fifty-six; obtaining a result image H56 obtained by subtracting the 2390nm wave band image from the 2330nm wave band image;
fiftieth, when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2405nm band image, judging that the condition is met; obtaining a result image H57 obtained by subtracting the 2405nm wave band image from the 2330nm wave band image;
sixthly, when each pixel value of the 2345nm wave band image is larger than the corresponding pixel value of the 2315nm wave band image, judging that the condition is met, namely fifty-eight; obtaining a result image H58 of subtracting the 2315nm wave band image from the 2345nm wave band image;
sixthly, when each pixel value of the 2345nm wave band image is larger than the pixel value corresponding to the 2330nm wave band image, judging that the condition is met by fifty-nine; obtaining a result image H59 of subtracting the 2330nm waveband image from the 2345nm waveband image;
sixty, when each pixel value of the 2360nm wave band image is larger than the corresponding pixel value of the 2375nm wave band image, judging that the condition sixty is met; obtaining a result image H60 obtained by subtracting the 2375nm wave band image from the 2360nm wave band image;
sixty-one step, when each pixel value of the 2360nm band image is larger than the corresponding pixel value of the 2390nm band image, judging that the condition sixty-one is met; obtaining a result image H61 obtained by subtracting the 2390nm wave band image from the 2360nm wave band image;
sixty-fourth step, when each pixel value of the 2360nm wave band image is larger than the corresponding pixel value of the 2405nm wave band image, judging that the condition is met sixty-two; obtaining a result image H62 obtained by subtracting the 2405nm wave band image from the 2360nm wave band image;
sixty-fifth step, when each pixel value of the 2375nm wave band image is larger than the corresponding pixel value of the 2390nm wave band image, judging that the condition sixty-three is met; obtaining a result image H63 obtained by subtracting the 2390nm wave band image from the 2375nm wave band image;
sixty-sixth, when each pixel value of the 2405nm wave band image is larger than the corresponding pixel value of the 2360nm wave band image, judging that the condition sixty-four is met; obtaining a result image H64 obtained by subtracting the 2360nm wave band image from the 2405nm wave band image;
sixty-seventh step, when each pixel value of the 2405nm wave band image is larger than the corresponding pixel value of the 2390nm wave band image, judging that the condition sixty-five is met; obtaining a result image H65 obtained by subtracting the 2390nm wave band image from the 2405nm wave band image;
sixty-eight steps, when each pixel value of the 2420nm wave band image is larger than the corresponding pixel value of the 2390nm wave band image, judging that the condition is sixty-six met; obtaining a result image H66 obtained by subtracting the 2390nm wave band image from the 2420nm wave band image;
sixty-ninth step, selecting the pixels which simultaneously meet the conditions of five, nine, thirty-six, thirty-nine, forty-six, forty-seven, fifty-one, fifty-six and sixty-five to obtain the sum H69 of the summation of H5, H9, H36, H39, H40, H46, H47, H51, H56 and H65 in the pixel range. The higher the value of H69, the higher the value of dickite abundance in the representative pixel.
Seventy, selecting the pixels which simultaneously meet the conditions of four, eight, fifteen, twenty-one, thirty-one, forty-two, forty-three, fifty-four and fifty-seven to obtain the sum H70 of the sum of H4, H8, H15, H21, H31, H42, H43, H54 and H57 in the pixel range. The higher the value of H70, the higher the abundance value of jarosite in the representative pixel.
Seventy, selecting the pixels which simultaneously meet the conditions of three, twenty five, twenty six, forty-one, forty-four, fifty-nine and sixty to obtain the sum H71 of the sums of H3, H25, H26, H41, H44, H59 and H60 in the pixel range. The higher the value of H71, the higher the value of the richness of the lizardtail within the representative picture element.
Seventy-second step, selecting the pixels satisfying the conditions fourteen, twenty-four, twenty-seven, thirty-five, forty-five and sixty-four simultaneously to obtain the sum H72 of the sums of H14, H24, H27, H35, H45 and H64 in the pixel range. The higher the H72 value is, the higher the abundance value of rhodochrosite in the representative picture element is.
Seventy-third step, selecting the pixels which simultaneously meet the conditions of the first, the second, the tenth, the eighteen, the twenty-third, the twenty-fourth, the twenty-ninth, the thirty and the fifty-third, and obtaining the sum H73 of the summation of H1, H2, H10, H18, H23, H24, H29, H30 and H53 in the pixel range. The higher the value of H73, the higher the gypsum abundance value within the representative pixel.
Seventy-fourth step, selecting the pixels which simultaneously satisfy the conditions six, eleven, twelve, sixteen, twenty and thirty-seven, and obtaining the sum H74 of the sums of H6, H11, H12, H16, H20 and H37 in the pixel range. The higher the value of H74, the higher the value representing the abundance of actinolite within the pixel.
Seventy-fifth step, selecting the pixels which simultaneously satisfy the conditions of twenty-eight, thirty-two, thirty-three, thirty-eight, forty-eight, fifty-five and sixty-three to obtain the sum H75 of the sums of H28, H32, H33, H38, H48, H55 and H63 in the pixel range. The higher the value of H75, the higher the value of pyrophyllite abundance in the representative pixel.
Seventy-sixth step, selecting the pixels which simultaneously meet the conditions of seven, thirteen, nineteen, twenty-two, twenty-four, thirty-nine, forty-nine, fifty-two and sixty-two, and obtaining the sum H76 of the sum of H7, H13, H19, H22, H24, H34, H39, H49, H52 and H62 in the pixel range. The higher the H76 value, the higher the alunite abundance value in the representative pixel.

Claims (11)

1. A hyperspectral image processing method for mineral information extraction is characterized by comprising the following steps: the method sequentially comprises the following steps: step (1), acquiring a hyperspectral image; preprocessing the hyperspectral image, and performing atmospheric correction to obtain image data of ground reflectivity; resampling the image data of the ground reflectivity, and extracting images of different wave bands; step (3), comparing pixel values of images of different wave bands, and marking a comparison result; and (4) determining abundance values of dickite, jarosite, lizardite, rhodochrosite, alunite, gypsum, gibbsite, actinolite and pyrophyllite according to the step (3).
2. The hyperspectral image processing method for mineral information extraction according to claim 1, wherein: in the step (2), the wavelength range of the extracted image includes 980nm, 1040nm, 1145nm, 1175nm, 1325nm, 1340nm, 1355nm, 1415nm, 1430nm, 1445nm, 1460nm, 1490nm, 1520nm, 1535nm, 1550nm, 1595nm, 1625nm, 1700nm, 1715nm, 1745nm, 1760nm, 1775nm, 1790nm, 1805nm, 1820nm, 1940nm, 1940nm, 1955nm, 1985nm, 2015nm, 2045nm, 2060nm, 2075nm, 2090nm, 2105nm, 2120nm, 2135nm, 2150nm, 2165nm, 2180nm, 2195nm, 2210nm, 2225nm, 2230nm, 2240nm, 2250 nm, 2255nm, 2270nm, 2285nm, 2300nm, 2315nm, 2330nm, 2345nm, 2360nm, 24290 nm, 235 nm, 230nm, 2405nm, and 2400 nm.
3. The hyperspectral image processing method for mineral information extraction according to claim 1, wherein: the step (3) includes the steps of,
when each pixel value of the 1145nm band image is greater than the corresponding pixel value of the 1175nm band image, judging that the condition I is met; obtaining a result image H1 of subtracting the 1175nm wave band image from the 1145nm wave band image;
when each pixel value of the 1325nm band image is larger than the pixel value corresponding to the 1445nm band image, judging that a second condition is met; obtaining a result image H2 of subtracting the 1445nm wave band image from the 1325nm wave band image;
when each pixel value of the 1340nm waveband image is larger than the pixel value corresponding to the 1355nm waveband image, judging that the condition three is met; obtaining a result image H3 of subtracting 1455nm wave band image from 1340nm wave band image;
when each pixel value of the 1355nm band image is larger than the corresponding pixel value of the 980nm band image, judging that the condition four is met; obtaining a result image H4 of subtracting the 980nm wave band image from the 1355nm wave band image;
when each pixel value of the 1355nm band image is larger than the corresponding pixel value of the 1415nm band image, judging that the condition five is met; obtaining a result image H5 of subtracting the 1415nm wave band image from the 1355nm wave band image;
when each pixel value of the 1355nm band image is larger than the corresponding pixel value of the 1445nm band image, judging that a sixth condition is met; obtaining a result image H5 of subtracting the 1445nm wave band image from the 1355nm wave band image;
when each pixel value of the 1415nm band image is greater than the corresponding pixel value of the 1430nm band image, judging that the condition seven is met; obtaining a result image H7 obtained by subtracting the 1430nm wave band image from the 1415nm wave band image;
when each pixel value of the 1415nm band image is greater than the corresponding pixel value of the 1460nm band image, judging that the condition eight is met; obtaining a result image H8 of subtracting the 1460nm wave band image from the 1415nm wave band image;
when each pixel value of the 1445nm band image is larger than the corresponding pixel value of the 1415nm band image, judging that the condition nine is met; obtaining a result image H9 of subtracting the 1415nm wave band image from the 1445nm wave band image;
when each pixel value of the 1460nm band image is larger than the pixel value corresponding to the 1445nm band image, judging that the condition ten is met; obtaining a 1460nm wave band image minus a 1445nm wave band image as a result image H10;
when each pixel value of the 1490nm band image is larger than the pixel value corresponding to the 1445nm band image, judging that the condition eleven is met; obtaining a result image H11 of subtracting the 1445nm wave band image from the 1490nm wave band image;
when each pixel value of the 1490nm band image is larger than the pixel value corresponding to the 1520nm band image, judging that the condition twelve is met; obtaining a result image H12 of subtracting the 1520nm wave band image from the 1490nm wave band image;
when each pixel value of the 1520nm band image is larger than the pixel value corresponding to the 1490nm band image, judging that the condition thirteen is met; obtaining a result image H13 of subtracting the 1490nm wave band image from the 1520nm wave band image;
when each pixel value of the 1535nm wave band image is larger than the corresponding pixel value of the 1445nm wave band image, judging that the condition is fourteen; obtaining a result image H14 of subtracting the 1445nm wave band image from the 1535nm wave band image;
when each pixel value of the 1595nm band image is larger than the corresponding pixel value of the 1550nm band image, judging that the condition fifteen is met; obtaining a result image H15 obtained by subtracting the 1550nm waveband image from the 1595nm waveband image;
when each pixel value of the 1625nm waveband image is larger than the corresponding pixel value of the 1550nm waveband image, judging that the condition is sixteen; obtaining a result image H16 obtained by subtracting the 1550nm waveband image from the 1625nm waveband image;
when each pixel value of the 1700nm band image is larger than the corresponding pixel value of the 1040nm band image, judging that the condition is seventeen; obtaining a result image H17 of subtracting the 1040nm waveband image from the 1700nm waveband image;
when each pixel value of the 1700nm waveband image is larger than the corresponding pixel value of the 1745nm waveband image, judging that the condition is met by eighteen; obtaining a result image H18 of subtracting the 1745nm band image from the 1700nm band image;
when each pixel value of the 1715nm waveband image is larger than the pixel value corresponding to the 1760nm waveband image, judging that the condition is met; obtaining a result image H19 obtained by subtracting the 1760nm wave band image from the 1715nm wave band image;
when each pixel value of the 1775nm wave band image is larger than the pixel value corresponding to the 1940nm wave band image, judging that the condition twenty is met; obtaining a result image H20 of subtracting the 1940nm wave band image from the 1775nm wave band image;
when each pixel value of the 1790nm band image is larger than the pixel value corresponding to the 2015nm band image, judging that the condition twenty-one is met; obtaining a result image H21 of subtracting the 2015nm waveband image from the 1790nm waveband image;
when each pixel value of the 1805nm wave band image is larger than the corresponding pixel value of the 1775nm wave band image, judging that the condition twenty-two is met; obtaining a 1805nm wave band image minus a 1775nm wave band image H22;
when each pixel value of the 1820nm band image is larger than the corresponding pixel value of the 1745nm band image, judging that the condition twenty-three is met; obtaining a result image H23 of subtracting the 1745nm band image from the 1820nm band image;
when each pixel value of the 1820nm wave band image is larger than the pixel value corresponding to the 1940nm wave band image, judging that the condition twenty-four is met; obtaining a result image H24 of subtracting the 1940nm wave band image from the 1820nm wave band image;
when each pixel value of the 1820nm band image is larger than the pixel value corresponding to the 1955nm band image, judging that the condition is satisfied by twenty-five; obtaining a 1820nm wave band image minus the 1955nm wave band image as a result H25;
when each pixel value of the 2045nm band image is greater than the corresponding pixel value of the 1985nm band image, judging that the condition twenty-six is met; obtaining a 2045nm wave band image and subtracting a 1985nm wave band image to obtain a result image H26;
when each pixel value of the 2045nm band image is larger than the pixel value corresponding to the 2015nm band image, judging that the condition is met by twenty-seven; obtaining a result image H27 of subtracting the 2015nm waveband image from the 2045nm waveband image;
when each pixel value of the 2045nm band image is larger than the corresponding pixel value of the 2075nm band image, judging that the condition is satisfied twenty-eight; obtaining a result image H28 of subtracting the 2075nm wave band image from the 2045nm wave band image;
when each pixel value of the 2060nm wave band image is larger than the pixel value corresponding to the 1940nm wave band image, judging that the condition is satisfied twenty-nine; obtaining a result image H29 of subtracting the 1940nm wave band image from the 2060nm wave band image;
when each pixel value of the 2060nm wave band image is larger than the pixel value corresponding to the 2210nm wave band image, judging that the condition is thirty; obtaining a result image H30 of subtracting the 2210nm wave band image from the 2060nm wave band image;
when each pixel value of the 2075nm band image is larger than the pixel value corresponding to the 2120nm band image, judging that the condition is satisfied thirty-one; obtaining a result image H31 of subtracting the 2120nm wave band image from the 2075nm wave band image;
when each pixel value of the 2105nm wave band image is larger than the pixel value corresponding to the 2090nm wave band image, judging that the condition is satisfied thirty-two; obtaining a result image H32 obtained by subtracting the 2090nm wave band image from the 2105nm wave band image;
when each pixel value of the 2105nm wave band image is larger than the corresponding pixel value of the 2165nm wave band image, judging that the condition is satisfied; obtaining a result image H33 obtained by subtracting the 2165nm wave band image from the 2105nm wave band image;
when each pixel value of the 2120nm band image is larger than the corresponding pixel value of the 2165nm band image, judging that the thirty-four conditions are met; obtaining a result image H34 obtained by subtracting the 2165nm wave band image from the 2120nm wave band image;
when each pixel value of the 2135nm wave band image is larger than the corresponding pixel value of the 2165nm wave band image, judging that the condition is satisfied thirty-five; obtaining a result image H35 obtained by subtracting the 2165nm wave band image from the 2135nm wave band image;
when each pixel value of the 2135nm wave band image is larger than the pixel value corresponding to the 2180nm wave band image, judging that the condition is satisfied thirty-six; obtaining a result image H36 of subtracting the 2180nm waveband image from the 2135nm waveband image;
when each pixel value of the 2150nm band image is larger than the corresponding pixel value of the 2255nm band image, judging that the condition thirty-seven is met; obtaining a result image H37 obtained by subtracting the 2255nm wave band image from the 2150nm wave band image;
when each pixel value of the 2195nm wave band image is larger than the corresponding pixel value of the 2165nm wave band image, judging that the condition is satisfied thirty-eight; obtaining a result image H38 of subtracting the 2165nm wave band image from the 2195nm wave band image;
when each pixel value of the 2195nm band image is larger than the corresponding pixel value of the 2180nm band image, judging that the condition is satisfied by thirty-nine; obtaining a result image H39 of subtracting the 2180nm waveband image from the 2195nm waveband image;
when each pixel value of the 2195nm band image is larger than the pixel value corresponding to the 2210nm band image, judging that the condition is met by forty; obtaining a result image H40 of subtracting the 2210nm wave band image from the 2195nm wave band image;
when each pixel value of the 2195nm waveband image is larger than the corresponding pixel value of the 2225nm waveband image, judging that the condition is forty one; obtaining a result image H41 obtained by subtracting the 2225nm waveband image from the 2195nm waveband image;
when each pixel value of the 2225nm band image is greater than the corresponding pixel value of the 2120nm band image, judging that the condition is satisfied forty-two; obtaining a result image H42 of subtracting the 2120nm wave band image from the 2225nm wave band image;
when each pixel value of the 2225nm band image is larger than the corresponding pixel value of the 2255nm band image, judging that the condition is satisfied forty-three; obtaining a result image H43 obtained by subtracting the 2255nm wave band image from the 2225nm wave band image;
when each pixel value of the 2225nm wave band image is larger than the corresponding pixel value of the 2300nm wave band image, judging that the condition fourteen is met; obtaining a result image H44 of subtracting the 2300nm wave band image from the 2225nm wave band image;
when each pixel value of the 2225nm band image is larger than the corresponding pixel value of the 2360nm band image, judging that the condition forty-five is met; obtaining a result image H45 obtained by subtracting the 2360nm wave band image from the 2225nm wave band image;
when each pixel value of the 2240nm waveband image is larger than the pixel value corresponding to the 2210nm waveband image, judging that the condition is satisfied forty-six; obtaining a result image H46 of subtracting the 2210nm wave band image from the 2240nm wave band image;
when each pixel value of the 2240nm waveband image is larger than the corresponding pixel value of the 2255nm waveband image, judging that the condition is forty-seven; obtaining a result image H47 obtained by subtracting the 2255nm wave band image from the 2240nm wave band image;
when each pixel value of the 2240nm waveband image is larger than the corresponding pixel value of the 2315nm waveband image, judging that the condition is satisfied forty-eight; obtaining a result image H48 obtained by subtracting the 2315nm wave band image from the 2240nm wave band image;
when each pixel value of the 2255nm band image is greater than the pixel value corresponding to the 2210nm band image, judging that the condition forty-nine is met; obtaining a result image H49 of subtracting the 2210nm wave band image from the 2255nm wave band image;
when each pixel value of the 2270nm band image is larger than the corresponding pixel value of the 2300nm band image, judging that the condition is met; obtaining a result image H50 of subtracting the 2300nm wave band image from the 2270nm wave band image;
when each pixel value of the 2285nm band image is larger than the corresponding pixel value of the 2300nm band image, judging that the condition is met; obtaining a result image H51 of subtracting the 2300nm wave band image from the 2285nm wave band image;
when each pixel value of the 2285nm band image is larger than the corresponding pixel value of the 2315nm band image, judging that the condition is met by fifty-two; obtaining a result image H52 of subtracting the 2315nm wave band image from the 2285nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2405nm band image, judging that the condition is satisfied; obtaining a result image H53 obtained by subtracting the 2405nm wave band image from the 2330nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2255nm band image, judging that the condition fifty-four is met; obtaining a result image H54 obtained by subtracting the 2255nm wave band image from the 2330nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2315nm band image, judging that the condition fifty-five is met; obtaining a result image H55 obtained by subtracting the 2315nm wave band image from the 2330nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2390nm band image, judging that the condition is met fifty-six; obtaining a result image H56 obtained by subtracting the 2390nm wave band image from the 2330nm wave band image;
when each pixel value of the 2330nm band image is larger than the corresponding pixel value of the 2405nm band image, judging that the condition is met; obtaining a result image H57 obtained by subtracting the 2405nm wave band image from the 2330nm wave band image;
when each pixel value of the 2345nm wave band image is larger than the corresponding pixel value of the 2315nm wave band image, judging that the condition is satisfied fifty-eight; obtaining a result image H58 of subtracting the 2315nm wave band image from the 2345nm wave band image;
when each pixel value of the 2345nm wave band image is larger than the pixel value corresponding to the 2330nm wave band image, judging that the condition is met by fifty-nine; obtaining a result image H59 of subtracting the 2330nm waveband image from the 2345nm waveband image;
when each pixel value of the 2360nm band image is larger than the corresponding pixel value of the 2375nm band image, judging that the condition sixty is met; obtaining a result image H60 obtained by subtracting the 2375nm wave band image from the 2360nm wave band image;
when each pixel value of the 2360nm band image is larger than the corresponding pixel value of the 2390nm band image, judging that the condition is sixty-one satisfied; obtaining a result image H61 obtained by subtracting the 2390nm wave band image from the 2360nm wave band image;
when each pixel value of the 2360nm band image is larger than the corresponding pixel value of the 2405nm band image, judging that the condition is satisfied sixty-two; obtaining a result image H62 obtained by subtracting the 2405nm wave band image from the 2360nm wave band image;
when each pixel value of the 2375nm band image is larger than the corresponding pixel value of the 2390nm band image, judging that the condition sixty-three is met; obtaining a result image H63 obtained by subtracting the 2390nm wave band image from the 2375nm wave band image;
when each pixel value of the 2405nm waveband image is larger than the corresponding pixel value of the 2360nm waveband image, judging that the condition sixty-four is met; obtaining a result image H64 obtained by subtracting the 2360nm wave band image from the 2405nm wave band image;
when each pixel value of the 2405nm band image is larger than the corresponding pixel value of the 2390nm band image, judging that the condition sixty-five is met; obtaining a result image H65 obtained by subtracting the 2390nm wave band image from the 2405nm wave band image;
when each pixel value of the 2420nm wave band image is larger than the corresponding pixel value of the 2390nm wave band image, judging that the condition is sixty-six met; and obtaining a result image H66 of subtracting the 2390nm wave band image from the 2420nm wave band image.
4. The hyperspectral image processing method for mineral information extraction according to claim 3, wherein: in the step (4), pixels satisfying conditions five, nine, thirty-six, thirty-nine, forty-six, forty-seven, fifty-one, fifty-six, and sixty-five are selected at the same time, and a sum H69 of the sums of H5, H9, H36, H39, H40, H46, H47, H51, H56, and H65 in the above pixel range is obtained, and a higher value of H69 represents a higher value of dickite abundance in the pixel.
5. The hyperspectral image processing method for mineral information extraction according to claim 3, wherein: in the step (4), selecting the pixels which simultaneously satisfy the conditions of four, eight, fifteen, twenty-one, thirty-one, forty-two, forty-three, fifty-four and fifty-seven, obtaining the sum H70 of the sum of H4, H8, H15, H21, H31, H42, H43, H54 and H57 in the above pixel range, wherein the higher the value of H70 is, the higher the abundance value of jarosite in the pixels is.
6. The hyperspectral image processing method for mineral information extraction according to claim 3, wherein: in the step (4), pixels satisfying the conditions three, twenty-five, twenty-six, forty-one, forty-four, fifty-nine and sixty are selected to obtain a sum H71 of the sums of H3, H25, H26, H41, H44, H59 and H60 in the above pixel range. The higher the value of H71, the higher the value of the richness of the lizardtail within the representative picture element.
7. The hyperspectral image processing method for mineral information extraction according to claim 3, wherein: in the step (4), selecting the pixels which simultaneously meet the conditions of fourteen, twenty-four, twenty-seven, thirty-five, forty-five and sixty-four, and obtaining the sum H72 of the H14, H24, H27, H35, H45 and H64 in the pixel range, wherein the higher the H72 value is, the higher the abundance value of rhodochrosite in the pixels is.
8. The hyperspectral image processing method for mineral information extraction according to claim 3, wherein: in the step (4), pixels which simultaneously satisfy the conditions of the first condition, the second condition, the tenth condition, the eighteen condition, the twenty-third condition, the twenty-fourth condition, the twenty-ninth condition, the thirty condition and the fifty-third condition are selected, a sum H73 of the sums of H1, H2, H10, H18, H23, H24, H29, H30 and H53 in the pixel range is obtained, and the higher the value of H73 is, the higher the gypsum abundance value in the pixel is represented.
9. The hyperspectral image processing method for mineral information extraction according to claim 3, wherein: in the step (4), selecting the pixels which simultaneously satisfy the conditions six, eleven, twelve, sixteen, twenty and thirty-seven, and obtaining the sum H74 of the sum of H6, H11, H12, H16, H20 and H37 in the pixel range, wherein the higher the H74 value is, the higher the abundance value of the actinolite in the pixels is.
10. The hyperspectral image processing method for mineral information extraction according to claim 3, wherein: in the step (4), pixels which simultaneously satisfy the conditions of twenty-eight, thirty-two, thirty-three, thirty-eight, forty-eight, fifty-five and sixty-three are selected to obtain a sum H75 of the sums of H28, H32, H33, H38, H48, H55 and H63 in the pixel range, and the higher the value of H75, the higher the pyrophyllite abundance value in the pixel is.
11. The hyperspectral image processing method for mineral information extraction according to claim 3, wherein: in the step (4), pixels satisfying the conditions seven, thirteen, nineteen, twenty-two, twenty-four, thirty-nine, forty-nine, fifty-two, and sixty-two are selected at the same time, and the sum H76 of the sums of H7, H13, H19, H22, H24, H34, H39, H49, H52, and H62 in the above pixel range is obtained. The higher the H76 value, the higher the alunite abundance value in the representative pixel.
CN201911327172.9A 2019-12-20 2019-12-20 Hyperspectral image processing method for mineral information extraction Pending CN111157459A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911327172.9A CN111157459A (en) 2019-12-20 2019-12-20 Hyperspectral image processing method for mineral information extraction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911327172.9A CN111157459A (en) 2019-12-20 2019-12-20 Hyperspectral image processing method for mineral information extraction

Publications (1)

Publication Number Publication Date
CN111157459A true CN111157459A (en) 2020-05-15

Family

ID=70557505

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911327172.9A Pending CN111157459A (en) 2019-12-20 2019-12-20 Hyperspectral image processing method for mineral information extraction

Country Status (1)

Country Link
CN (1) CN111157459A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111562273A (en) * 2020-06-05 2020-08-21 大连工业大学 Hyperspectrum-based fish water jet descaling slight damage visualization method

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103900965A (en) * 2012-12-25 2014-07-02 核工业北京地质研究院 Hyperspectral image processing method used for extracting calcite information
CN103903225A (en) * 2012-12-25 2014-07-02 核工业北京地质研究院 Hyperspectral image processing method for dolomite information extraction
CN104574283A (en) * 2013-10-23 2015-04-29 核工业北京地质研究院 Hyperspectral image processing method for extracting pyrophyllite information
CN104732488A (en) * 2013-12-20 2015-06-24 核工业北京地址研究院 Hyperspectral image processing method for actinolite information extraction
CN105405102A (en) * 2014-08-28 2016-03-16 核工业北京地质研究院 High-spectral image processing method for gibbsite information extraction
CN105787915A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 Hyperspectral image processing method for extracting information of jarosite
CN105787483A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 Hyperspectral image processing method for extracting information of lizardite
CN105787916A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 Hyperspectral image processing method for extracting information of dickite
CN105784602A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 High spectral image processing method for extracting information of rhodochrosite
CN105787914A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 Hyperspectral image processing method for alunite information extraction

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103900965A (en) * 2012-12-25 2014-07-02 核工业北京地质研究院 Hyperspectral image processing method used for extracting calcite information
CN103903225A (en) * 2012-12-25 2014-07-02 核工业北京地质研究院 Hyperspectral image processing method for dolomite information extraction
CN104574283A (en) * 2013-10-23 2015-04-29 核工业北京地质研究院 Hyperspectral image processing method for extracting pyrophyllite information
CN104732488A (en) * 2013-12-20 2015-06-24 核工业北京地址研究院 Hyperspectral image processing method for actinolite information extraction
CN105405102A (en) * 2014-08-28 2016-03-16 核工业北京地质研究院 High-spectral image processing method for gibbsite information extraction
CN105787915A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 Hyperspectral image processing method for extracting information of jarosite
CN105787483A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 Hyperspectral image processing method for extracting information of lizardite
CN105787916A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 Hyperspectral image processing method for extracting information of dickite
CN105784602A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 High spectral image processing method for extracting information of rhodochrosite
CN105787914A (en) * 2014-12-22 2016-07-20 核工业北京地质研究院 Hyperspectral image processing method for alunite information extraction

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111562273A (en) * 2020-06-05 2020-08-21 大连工业大学 Hyperspectrum-based fish water jet descaling slight damage visualization method

Similar Documents

Publication Publication Date Title
CN111157459A (en) Hyperspectral image processing method for mineral information extraction
Kim et al. Ginsenoside profiles and related gene expression during foliation in Panax ginseng Meyer
EP4033881A4 (en) Treating plants using feature values and ground planes extracted from a single image
CA2206157A1 (en) Method of extraction of commercially valuable fractions of fenugreek
Fitriansyah et al. Correlation of total phenolic, flavonoid and carotenoid content of Sesbania sesban (L. Merr) leaves extract with DPPH scavenging activities
WO2007075580A3 (en) Recovery of residual plant components after distillation of essential oils
WO2010104603A3 (en) Radiation curable resin composition and rapid three-dimensional imaging process using the same
WO2012163851A2 (en) Easily removable nail varnish composition
EP2120094A3 (en) Black photosensitive resin composition, and color filter and method of producing the same
EP3598233A1 (en) Liquid photoinitiating compound and use of the same
Eftekhari et al. Chemodiversity evaluation of grape (Vitis vinifera) vegetative parts during summer and early fall
Curl Bell pepper carotenoids, carotenoids of green bell peppers
EP3747918A4 (en) Polymer including polymerizable triptycene derivative compound as component
MXPA04007027A (en) Method of making flavedo powder for enhancement of orange juice and products thereof.
Francis et al. Ultrastructural characterization of yellow shoulder disorder in a uniform ripening tomato genotype
EP3922654A4 (en) Method for producing polyphenol derivative, polyphenol derivative, and polyphenol derivative-containing resin composition material
Xu et al. Occurrence of long-chain n-alkenols, diols, keto-ols and sec-alkanols in a sediment core from a hypereutrophic, freshwater lake
EP3992713A4 (en) Production method for radiation-sensitive resin composition
DE102013101035A1 (en) Easily removable nail polish composition
JP2019066291A (en) Method for identifying plant species of licorice
Habib et al. Efficient oil palm total RNA extraction with a total RNA extraction kit
Warnasih et al. Phytochemical screening and antioxidant activity of date (Phoenix dactylifera) seed extracts
DE102010030995A1 (en) Use of beta-isophorone as a solvent
CN104109454A (en) Anti-solarization powdery paint
DE2607760A1 (en) Ski boot sole scraper - made from plastics or light metal including number of scraping edges to remove ice or snow

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200515

RJ01 Rejection of invention patent application after publication