CN111157459A - Hyperspectral image processing method for mineral information extraction - Google Patents
Hyperspectral image processing method for mineral information extraction Download PDFInfo
- 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
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 26
- 238000003672 processing method Methods 0.000 title claims abstract description 19
- 229910052500 inorganic mineral Inorganic materials 0.000 title claims description 15
- 239000011707 mineral Substances 0.000 title claims description 15
- 229910052891 actinolite Inorganic materials 0.000 claims abstract description 15
- 229910052934 alunite Inorganic materials 0.000 claims abstract description 15
- 239000010424 alunite Substances 0.000 claims abstract description 15
- 229910001649 dickite Inorganic materials 0.000 claims abstract description 15
- 229910052602 gypsum Inorganic materials 0.000 claims abstract description 15
- 239000010440 gypsum Substances 0.000 claims abstract description 15
- 229910052935 jarosite Inorganic materials 0.000 claims abstract description 15
- XMWCXZJXESXBBY-UHFFFAOYSA-L manganese(ii) carbonate Chemical compound [Mn+2].[O-]C([O-])=O XMWCXZJXESXBBY-UHFFFAOYSA-L 0.000 claims abstract description 15
- 229910052903 pyrophyllite Inorganic materials 0.000 claims abstract description 15
- KPZTWMNLAFDTGF-UHFFFAOYSA-D trialuminum;potassium;hexahydroxide;disulfate Chemical compound [OH-].[OH-].[OH-].[OH-].[OH-].[OH-].[Al+3].[Al+3].[Al+3].[K+].[O-]S([O-])(=O)=O.[O-]S([O-])(=O)=O KPZTWMNLAFDTGF-UHFFFAOYSA-D 0.000 claims abstract description 15
- 229910001679 gibbsite Inorganic materials 0.000 claims abstract description 12
- 229910052899 lizardite Inorganic materials 0.000 claims abstract description 12
- IBPRKWGSNXMCOI-UHFFFAOYSA-N trimagnesium;disilicate;hydrate Chemical compound O.[Mg+2].[Mg+2].[Mg+2].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-] IBPRKWGSNXMCOI-UHFFFAOYSA-N 0.000 claims abstract description 12
- 238000002310 reflectometry Methods 0.000 claims abstract description 7
- 238000012952 Resampling Methods 0.000 claims abstract description 4
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 6
- 241000722273 Saururus cernuus Species 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image 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
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.
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)
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)
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 |
-
2019
- 2019-12-20 CN CN201911327172.9A patent/CN111157459A/en active Pending
Patent Citations (10)
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)
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 |