CN105488488B - City black and odorous water remote sensing recognition method and device - Google Patents
City black and odorous water remote sensing recognition method and device Download PDFInfo
- Publication number
- CN105488488B CN105488488B CN201510917811.2A CN201510917811A CN105488488B CN 105488488 B CN105488488 B CN 105488488B CN 201510917811 A CN201510917811 A CN 201510917811A CN 105488488 B CN105488488 B CN 105488488B
- Authority
- CN
- China
- Prior art keywords
- identified
- water surface
- surface point
- point
- chromaticity coordinate
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
The embodiment of the invention discloses a kind of city black and odorous water remote sensing recognition method and devices, obtain the remotely-sensed data of water surface point to be identified；The spectral tristimulus value of water surface point to be identified is calculated based on acquired remotely-sensed data；The saturation degree of water surface point to be identified is calculated based on spectral tristimulus value；If saturation degree is less than preset threshold, determine that water surface point to be identified is black and odorous water.It can be seen that, black and odorous water remote sensing recognition method and device in city provided in an embodiment of the present invention, saturation degree is calculated based on remotely-sensed data, it whether is that black and odorous water identifies to water body based on saturation degree, inherently disclose the difference of black and odorous water and other water bodys in city, the blank of city black and odorous water remote sensing recognition method has been filled up in the identification for realizing city black and odorous water.
Description
Technical field
The present invention relates to remote sensing technology field, more specifically to a kind of city black and odorous water remote sensing recognition method and
Device.
Background technique
In the process of city prosperity, industry, agricultural and sanitary wastewater make the river water in completed region of the city no longer limpid,
It is even black, and stench is distributed, become city black and odorous water.Traditional research for black and odorous water be mainly environment,
The relevant research of ecology, lays particular emphasis on biochemical genetic analysis and differentiation, the research ten about remote sensing monitoring city black and odorous water
Divide rare.
In terms of using remote sensing image identification black and odorous water, has identify " lake is general " using remote sensing images at present, lake is general can be with
It is considered a kind of special black and odorous water, but is not belonging to typical city black and odorous water.
Inventor has found in the implementation of the present invention, identifies the poor universality of the general method in lake at present, can only identify
Lake pan class black and odorous water out, and cannot identify city black and odorous water.Therefore, how to identify that city black and odorous water becomes urgently to be resolved
The problem of.
Summary of the invention
The object of the present invention is to provide a kind of city black and odorous water remote sensing recognition method and devices, to city black and odorous water
It is identified.
To achieve the above object, the present invention provides the following technical scheme that
A kind of city black and odorous water remote sensing recognition method, comprising:
Obtain the remotely-sensed data of water surface point to be identified；
The spectral tristimulus value of water surface point to be identified is calculated based on the remotely-sensed data；
The saturation degree of water surface point to be identified is calculated based on the spectral tristimulus value；
If the saturation degree is less than preset threshold, determine that water surface point to be identified is black and odorous water.
The above method, it is preferred that the remotely-sensed data of the water surface point to be identified include: in water surface remote sensing images with the water
The tristimulus values of the corresponding pixel of face point to be identified；
The spectral tristimulus value for calculating water surface point to be identified based on the remotely-sensed data includes:
The tristimulus values are converted to spectral tristimulus value by the transformational relation according to tristimulus values and tristimulus values.
The above method, it is preferred that the remotely-sensed data of the water surface point to be identified includes: each of water surface point to be identified
The actual measurement Remote Sensing Reflectance of wave band；
The spectral tristimulus value for calculating water surface point to be identified based on the remotely-sensed data includes:
The actual measurement Remote Sensing Reflectance of each wave band is converted into satellite band equivalent reflectivity；
It is obtaining water surface point to be identified with trichromatic wave band pair according to the corresponding relationship of equivalent reflectivity and wave band
The tristimulus values answered；
The tristimulus values are converted to spectral tristimulus value by the transformational relation according to tristimulus values and tristimulus values.
The above method, it is preferred that the saturation degree that water surface point to be identified is calculated based on the spectral tristimulus value
Include:
The spectral tristimulus value is converted to the first chromaticity coordinate of water surface point to be identified；
Calculate the second chromaticity coordinate of the dominant wavelength of the characterized color of the first chromaticity coordinate；
Determine the first distance between the equal-energy white point in first chromaticity coordinate and chromatic diagram and second color
Spend the second distance between the equal-energy white point in coordinate and chromatic diagram；
The ratio of the first distance and the second distance is determined as to the saturation degree of water surface point to be identified.
The above method, it is preferred that the dominant wavelength corresponding second for calculating the characterized color of the first chromaticity coordinate
Chromaticity coordinate includes:
It calculates using the equal-energy white point in chromatic diagram as starting point, using first chromaticity coordinate as the vector of terminal and the color
Spend the angle of preset coordinate axis in figure；
According to the corresponding relationship of preset angle and dominant wavelength, dominant wavelength corresponding with the angle is determined；
According to the corresponding relationship of preset dominant wavelength and chromaticity coordinate, determination is corresponding with identified dominant wavelength
Chromaticity coordinate.
A kind of city black and odorous water remote sensing recognition device, comprising:
Module is obtained, for obtaining the remotely-sensed data of water surface point to be identified；
First computing module, for calculating the spectral tristimulus value of water surface point to be identified based on the remotely-sensed data；
Second computing module, for calculating the saturation degree of water surface point to be identified based on the spectral tristimulus value；
Determining module determines that water surface point to be identified is black and odorous water if being less than preset threshold for the saturation degree.
Above-mentioned apparatus, it is preferred that the remotely-sensed data of the water surface point to be identified include: in water surface remote sensing images with the water
The tristimulus values of the corresponding pixel of face point to be identified；
First computing module includes:
First converting unit converts the tristimulus values for the transformational relation according to tristimulus values and tristimulus values
For spectral tristimulus value.
Above-mentioned apparatus, it is preferred that the remotely-sensed data of the water surface point to be identified includes: each of water surface point to be identified
The actual measurement Remote Sensing Reflectance of wave band；
First computing module includes:
Second converting unit, for the actual measurement Remote Sensing Reflectance of each wave band to be converted to satellite band equivalent reflective
Rate；
Acquiring unit, for according to equivalent reflectivity and wave band corresponding relationship, obtain water surface point to be identified with
The corresponding tristimulus values of trichromatic wave band；
Third converting unit, for being converted to the tristimulus values according to the transformational relation of tristimulus values and tristimulus values
Spectral tristimulus value.
Above-mentioned apparatus, it is preferred that second computing module includes:
4th converting unit, the first coloration for the spectral tristimulus value to be converted to water surface point to be identified are sat
Mark；
Computing unit, the second chromaticity coordinate of the dominant wavelength for calculating the characterized color of the first chromaticity coordinate；
First determination unit, for determine between the equal-energy white point in first chromaticity coordinate and chromatic diagram first away from
From and second chromaticity coordinate and chromatic diagram in equal-energy white point between second distance；
Second determination unit waits knowing for the ratio of the first distance and the second distance to be determined as the water surface
The saturation degree of other point.
Above-mentioned apparatus, it is preferred that the computing unit includes:
Computation subunit is eventually with first chromaticity coordinate for calculating using the equal-energy white point in chromatic diagram as starting point
The angle of preset coordinate axis in the vector and the chromatic diagram of point；
First determines subelement, for the corresponding relationship according to preset angle and dominant wavelength, the determining and folder
The corresponding dominant wavelength in angle；
Second determines subelement, for the corresponding relationship according to preset dominant wavelength and chromaticity coordinate, determining and institute
The corresponding chromaticity coordinate of determining dominant wavelength.
By above scheme it is found that a kind of city black and odorous water remote sensing recognition method and device provided by the present application, obtains
The remotely-sensed data of water surface point to be identified；The spectral tristimulus value of water surface point to be identified is calculated based on acquired remotely-sensed data；Base
The saturation degree of water surface point to be identified is calculated in spectral tristimulus value；If saturation degree is less than preset threshold, water surface point to be identified is determined
For black and odorous water.As it can be seen that black and odorous water remote sensing recognition method and device in city provided in an embodiment of the present invention, is based on remotely-sensed data
Saturation degree is calculated, whether is that black and odorous water identifies to water body based on saturation degree, inherently discloses black smelly water in city
The difference of body and other water bodys realizes the identification of city black and odorous water, has filled up city black and odorous water remote sensing recognition method
Blank.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.
Fig. 1 is a kind of implementation flow chart of black and odorous water remote sensing recognition method in city provided in an embodiment of the present invention；
Fig. 2 is one kind of the saturation degree provided in an embodiment of the present invention that water surface point to be identified is calculated based on spectral tristimulus value
Implementation flow chart；
Fig. 3 is the schematic diagram for the CIE1931 chromatic diagram that International Commission on Illumination formulates；
Fig. 4 is characterized corresponding second color of dominant wavelength of color by first chromaticity coordinate of calculating provided in an embodiment of the present invention
Spend a kind of implementation flow chart of coordinate；
Fig. 5 is a kind of structural schematic diagram of black and odorous water remote sensing recognition device in city provided in an embodiment of the present invention；
Fig. 6 is a kind of structural schematic diagram of the first computing module provided in an embodiment of the present invention；
Fig. 7 is another structural schematic diagram of the first computing module provided in an embodiment of the present invention；
Fig. 8 is a kind of structural schematic diagram of the second computing module provided in an embodiment of the present invention；
Fig. 9 is a kind of structural schematic diagram of computing unit provided in an embodiment of the present invention.
Specification and claims and term " first " in above-mentioned attached drawing, " second ", " third " " the 4th " etc. (if
In the presence of) it is part for distinguishing similar, without being used to describe a particular order or precedence order.It should be understood that using in this way
Data be interchangeable under appropriate circumstances, so that embodiments herein described herein can be in addition to illustrating herein
Sequence in addition is implemented.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art are obtained every other under that premise of not paying creative labor
Embodiment shall fall within the protection scope of the present invention.
Referring to Fig. 1, Fig. 1 is a kind of realization stream of black and odorous water remote sensing recognition method in city provided in an embodiment of the present invention
Cheng Tu may include:
Step S11: the remotely-sensed data of water surface point to be identified is obtained；
Water surface point to be identified refers to any point in waters to be identified in completed region of the city, the remote sensing number of water surface point to be identified
According to three primary colors (the i.e. red, green and blue of pixel corresponding with the point to be identified in the remote sensing images that can be waters to be identified
Color) value (e.g., trichromatic reflectivity), it is also possible to the actual measurement Remote Sensing Reflectance of each wave band of the point to be identified.
Step S12: the spectral tristimulus value of water surface point to be identified is calculated based on acquired remotely-sensed data；
In order to carry out accurate digitized representations, International Commission on Illumination (The Commission to color
Internationale de L ' E^clairage, CIE) develop a set of standard color system CIE-XYZ color system.CIE-
Select three ideal primary colors [X], [Y], [Z] instead of three primary color reds [R], green in CIE-RGB system in XYZ system
[G], blue [B], so that the tristimulus values XYZ of spectrum is all positive value in colorimeter system, and two primary colors of X, Z only represents coloration, does not have
Have brightness, brightness only with tristimulus values Y at than.
Step S13: the saturation degree of water surface point to be identified is calculated based on spectral tristimulus value；
Step S14: if saturation degree is less than preset threshold, determine that water surface point to be identified is black and odorous water.
Wherein, preset threshold can be determined according to the experience of measured data.For example, by Shenyang black and odorous water reality
Ground measurement determines that preset threshold corresponding with the region is 0.1.Certainly, after having changed new survey region, this threshold value may
It can be varied.
Black and odorous water remote sensing recognition method in city provided in an embodiment of the present invention obtains the remote sensing number of water surface point to be identified
According to；The spectral tristimulus value of water surface point to be identified is calculated based on acquired remotely-sensed data；Water is calculated based on spectral tristimulus value
The saturation degree of face point to be identified；If saturation degree is less than preset threshold, determine that water surface point to be identified is black and odorous water.As it can be seen that this hair
The city black and odorous water remote sensing recognition method that bright embodiment provides calculates saturation degree based on remotely-sensed data, based on saturation degree to water
Whether body is that black and odorous water is identified, inherently discloses the difference of city black and odorous water Yu other water bodys, realizes city
The blank of city black and odorous water remote sensing recognition method has been filled up in the identification of city's black and odorous water.
Optionally, if the remotely-sensed data of water surface point to be identified is picture corresponding with water surface point to be identified in water surface remote sensing images
The tristimulus values of member, then the implementation that the spectral tristimulus value of water surface point to be identified is calculated based on remotely-sensed data can be with are as follows:
Tristimulus values are converted to spectral tristimulus value by the transformational relation according to tristimulus values and tristimulus values.
The relation formula of CIE-RGB and CIE-XYZ that CIE is defined are as follows:
In the embodiment of the present invention, tristimulus values can be converted into spectral tristimulus value according to formula (1).
Optionally, if the remotely-sensed data of water surface point to be identified is the actual measurement remote sensing reflection of each wave band of water surface point to be identified
Rate；The implementation that the spectral tristimulus value of water surface point to be identified is then calculated based on remotely-sensed data can be with are as follows:
The actual measurement Remote Sensing Reflectance of each wave band is converted into satellite band equivalent reflectivity；
Optionally, actual measurement Remote Sensing Reflectance can be converted by following formula:
Wherein, ReqFor satellite band equivalent reflectivity, Rrs(λ) is actual measurement Remote Sensing Reflectance, fSRF(λ) is satellite band light
Compose receptance function.
It is obtaining water surface point to be identified with trichromatic wave band pair according to the corresponding relationship of equivalent reflectivity and wave band
The tristimulus values answered；And the corresponding tristimulus values of trichromatic wave band are equivalent reflectivity corresponding with trichromatic wave band.
Different wave bands corresponds to different reflectivity.Therefore, according to trichromatic wavelength and pair of wavelength and reflectivity
It should be related to, can determine trichromatic reflectivity, in the embodiment of the present invention, using trichromatic reflectivity as tristimulus values.
Trichromatic grayvalue transition is spectral tristimulus value by transformational relation according to tristimulus values and tristimulus values.
In the embodiment of the present invention, tristimulus values can be converted into spectral tristimulus value according to formula (1).
Optionally, the one of the saturation degree provided in an embodiment of the present invention that water surface point to be identified is calculated based on spectral tristimulus value
Kind implementation flow chart is as shown in Fig. 2, may include:
Step S21: spectral tristimulus value is converted to the first chromaticity coordinate of water surface point to be identified；
Tristimulus values in CIE-XYZ system is particularly useful to color is defined, and is but not easy its visual result.For this purpose, CIE in
Define two-dimensional chromaticity figure within 1931, for depicting color on X-Y scheme, but color therein is unrelated with brightness, chromatic diagram
Upper two-dimensional coordinate x and y is to calculate the chromaticity coordinate got from tristimulus values XYZ, and calculation formula is as follows:
Due to x+y+z=1, with x, two values of y are assured that a color, therefore can use CIE-xy chromatic diagram (figure
3) all colours in visible-range are indicated, each color all corresponds to a chromaticity coordinate (x, y), all falls within horseshoe-shaped light
In the range of spectrum track surrounds.S is referred to as equal-energy white luminous point in figure, and chromaticity coordinate is (0.3333,0.3333), indicates three kinds of originals
Color mixed in equal amounts；C point is a color chromaticity coordinates point, and coordinate is (x, y), and the extended line of SC hands over spectrum locus (i.e. monochromatic rail
Mark) in D point, the wavelength of D point spectrum colour is exactly the dominant wavelength of C point color.Dominant wavelength is the important indicator of color quantizing, it
It will be seen that light color is indicated from 380nm to 700nm with 1 nanometer intervals, a kind of color can be indicated in the form of specific wavelength
Tone.
Step S22: the second chromaticity coordinate of the dominant wavelength of the characterized color of the first chromaticity coordinate is calculated；
If the first chromaticity coordinate is C point coordinate in Fig. 3, the second color of the dominant wavelength of the characterized color of the first chromaticity coordinate
Spending coordinate is D point coordinate.
Step S23: the first distance and the second color between the equal-energy white point in the first chromaticity coordinate and chromatic diagram are determined
Spend the second distance between the equal-energy white point in coordinate and chromatic diagram；
Step S24: the ratio of first distance and second distance is determined as to the saturation degree of water surface point to be identified.
Taking what is shown in fig. 3 as an example, the saturation degree of C point be C point to S point distance SC and D point to S point distance SD ratio
SC/SD。
Optionally, the dominant wavelength corresponding second provided in an embodiment of the present invention for calculating the characterized color of the first chromaticity coordinate
A kind of implementation flow chart of chromaticity coordinate is as shown in figure 4, may include:
Step S41: it calculates using the equal-energy white point in chromatic diagram as starting point, using the first chromaticity coordinate as the vector of terminal and color
Spend the angle of preset coordinate axis in figure；
For ease of description, as shown in figure 3, establishing a new rectangular coordinate system o'x'y' in chromatic diagram, wherein coordinate
It is that equal-energy white point S, x' axis that origin o' is located in chromatic diagram is parallel with chromatic diagram y-axis and direction is consistent, x in y' axis and chromatic diagram
Axis is parallel and direction is consistent, and x' axis is vertical with y' axis to be met at o' point chromatic diagram is divided into four quadrants.
The chromaticity coordinate (x, y) of each nano wave length of chromatic diagram object-line (i.e. spectrum locus) is converted into new coordinate
Coordinate (x', y') under system, conversion formula are as follows:
In o'x'y' coordinate system, it is assumed that vector (x', y') and x' positive axis angle are α, then angle α calculation formula are as follows: α=
ARCTAN2 (y', x')=ARCTAN2 (x-0.3333, y-0.3333) (4)
Wherein, ARCTAN2 function representation bivariate arctan function, the codomain of bivariate tangent function be (- 180 °,
180 °), i.e., vector (x', y') is rotated counterclockwise by x' positive axis from the negative axis of x' again to the process of the negative axis of x', arctan function from-
180 ° change to -180 °.
Since the x' axis in o'x'y' coordinate system is parallel with the y-axis in chromatic diagram, vector (x', y') and x' positive axis
Angle is equal to the angle of the y-axis in vector (x', y') and chromatic diagram.
After coordinate is converted, ARCTAN2 function can make the angle α with dominant wavelength be incremented by, i.e., from 380nm to
The corresponding angle α 700nm is gradually increased, in this way, the α angle that can establish a corresponding incremental order from 380nm to 700nm is looked into
Look for table, the i.e. mapping table of angle and dominant wavelength.
In the embodiment of the present invention, wave spectrum chromaticity coordinate each in the spectrum locus of CIE-xy chromatic diagram is also established one
From 380nm to 700nm, wavelength corresponds to the look-up table of spectral chromaticity coordinates, the i.e. corresponding relationship of dominant wavelength and chromaticity coordinate
Table.
Step S42: according to the corresponding relationship of preset angle and dominant wavelength, face corresponding with determined angle is determined
Color dominant wavelength；
Specifically, being selected if identified angle value is not present in the corresponding relationship of preset angle and dominant wavelength
It selects and the corresponding dominant wavelength of angle that identified angle is nearest is dominant wavelength corresponding with identified angle.Example
Such as, it is assumed that identified angle α is 6 degree, and in the corresponding relationship of preset angle and dominant wavelength, 558nm is corresponding to be
5.51 degree.559nm corresponding is 7.01 degree, and 6 degree closer from 5.51 degree, it is thus determined that color master corresponding with determined angle
Wavelength is 558nm.
Step S43: according to the corresponding relationship of preset dominant wavelength and chromaticity coordinate, determining and identified color master
The corresponding chromaticity coordinate of wavelength.
Corresponding with embodiment of the method, the embodiment of the present invention also provides a kind of city black and odorous water remote sensing recognition device, this
A kind of structural schematic diagram for the city black and odorous water remote sensing recognition device that inventive embodiments provide is as shown in figure 5, may include:
Obtain module 51, the first computing module 52, the second computing module 53 and determining module 54；Wherein,
Obtain the remotely-sensed data that module 51 is used to obtain water surface point to be identified；
Water surface point to be identified refers to any point in waters to be identified, and the remotely-sensed data of water surface point to be identified can be wait know
Three primary colors (i.e. red, green and blue) value of pixel corresponding with the point to be identified in the remote sensing images in other waters, can also be with
It is the actual measurement Remote Sensing Reflectance of each wave band of the point to be identified.
First computing module 52 is used to calculate the spectral tristimulus value of water surface point to be identified based on remotely-sensed data；
In order to carry out accurate digitized representations, International Commission on Illumination (The Commission to color
Internationale de L ' E^clairage, CIE) develop a set of standard color system CIE-XYZ color system.CIE-
Select three ideal primary colors [X], [Y], [Z] instead of three primary color reds [R], green in CIE-RGB system in XYZ system
[G], blue [B], so that the tristimulus values XYZ of spectrum is all positive value in colorimeter system, and two primary colors of X, Z only represents coloration, does not have
Have brightness, brightness only with tristimulus values Y at than.
Second computing module 53 is used to calculate the saturation degree of water surface point to be identified based on spectral tristimulus value；
If determining module 54 is less than preset threshold for saturation degree, determine that water surface point to be identified is black and odorous water.
Wherein, preset threshold can be determined according to the experience of measured data.For example, by Shenyang black and odorous water reality
Ground measurement determines that preset threshold corresponding with the region is 0.1.Certainly, after having changed new survey region, this threshold value may
It can be varied.
Black and odorous water remote sensing recognition device in city provided in an embodiment of the present invention obtains the remote sensing number of water surface point to be identified
According to；The spectral tristimulus value of water surface point to be identified is calculated based on acquired remotely-sensed data；Water is calculated based on spectral tristimulus value
The saturation degree of face point to be identified；If saturation degree is less than preset threshold, determine that water surface point to be identified is black and odorous water.As it can be seen that this hair
The city black and odorous water remote sensing recognition device that bright embodiment provides calculates saturation degree based on remotely-sensed data, based on saturation degree to water
Whether body is that black and odorous water is identified, inherently discloses the difference of black and odorous water and other water bodys in city, realizes
The blank of city black and odorous water remote sensing recognition method has been filled up in the identification of city black and odorous water.
Optionally, a kind of structural schematic diagram of the first computing module 52 provided in an embodiment of the present invention is as shown in fig. 6, can be with
Include:
First converting unit 61, if the remotely-sensed data for water surface point to be identified is to wait knowing with the water surface in water surface remote sensing images
The tristimulus values of the corresponding pixel of other point, according to the transformational relation of tristimulus values (specially gray value) and tristimulus values, by three
Primary color value is converted to spectral tristimulus value.
Tristimulus values can be converted into spectral tristimulus value according to formula (1).
Optionally, another structural schematic diagram of the first computing module 52 provided in an embodiment of the present invention is as shown in fig. 7, can
To include:
Second conversion module 71, if the remotely-sensed data for water surface point to be identified is each wave band of water surface point to be identified
Remote Sensing Reflectance is surveyed, the actual measurement Remote Sensing Reflectance of each wave band is converted into satellite band equivalent reflectivity；
The actual measurement Remote Sensing Reflectance of each wave band can be converted into satellite band equivalent reflectivity according to formula (2).
Module 72 is obtained, for the corresponding relationship according to equivalent reflectivity and wave band, obtains water surface point to be identified
Tristimulus values corresponding with trichromatic wave band；Tristimulus values corresponding with trichromatic wave band be and trichromatic wave band pair
The equivalent reflectivity answered.
Different wave bands corresponds to different reflectivity.Therefore, according to trichromatic wavelength and pair of wavelength and reflectivity
It should be related to, can determine trichromatic reflectivity, in the embodiment of the present invention, using trichromatic reflectivity as tristimulus values.
Third conversion module 73 turns trichromatic gray value for the transformational relation according to tristimulus values and tristimulus values
It is changed to spectral tristimulus value.
Tristimulus values can be converted into spectral tristimulus value according to formula (1).
Optionally, a kind of structural schematic diagram of the second computing module 53 provided in an embodiment of the present invention is as shown in figure 8, can be with
Include:
4th converting unit 81, computing unit 82, the first determination unit 83 and the second determination unit 84；Wherein,
4th converting unit 81 is used to be converted to spectral tristimulus value the first chromaticity coordinate of water surface point to be identified；
Tristimulus values in CIE-XYZ system is particularly useful to color is defined, and is but not easy its visual result.For this purpose, CIE in
Define two-dimensional chromaticity figure within 1931, for depicting color on X-Y scheme, but color therein is unrelated with brightness, chromatic diagram
Upper two-dimensional coordinate x and y is to calculate the chromaticity coordinate got from tristimulus values XYZ, and calculation formula is as follows:
Due to x+y+z=1, with x, two values of y are assured that a color, therefore can use CIE-xy chromatic diagram (figure
3) all colours in visible-range are indicated, each color all corresponds to a chromaticity coordinate (x, y), all falls within horseshoe-shaped light
In the range of spectrum track surrounds.S is referred to as equal-energy white luminous point in figure, and chromaticity coordinate is (0.3333,0.3333), indicates three kinds of originals
Color mixed in equal amounts；C point is a color chromaticity coordinates point, and coordinate is (x, y), and the extended line of SC hands over spectrum locus (i.e. monochromatic rail
Mark) in D point, the wavelength of D point spectrum colour is exactly the dominant wavelength of C point color.Dominant wavelength is the important indicator of color quantizing, it
It will be seen that light color is indicated from 380nm to 700nm with 1 nanometer intervals, a kind of color can be indicated in the form of specific wavelength
Tone.
Computing unit 82 is used to calculate the second chromaticity coordinate of the dominant wavelength of the characterized color of the first chromaticity coordinate；
If the first chromaticity coordinate is C point coordinate in Fig. 3, the second color of the dominant wavelength of the characterized color of the first chromaticity coordinate
Spending coordinate is D point coordinate.
First determination unit 83 is used to determine the first distance between equal-energy white point in the first chromaticity coordinate and chromatic diagram,
And the second distance between the second equal-energy white point in chromaticity coordinate and chromatic diagram；
Second determination unit 84 is used to for the ratio of first distance and second distance being determined as the saturation of water surface point to be identified
Degree.
Taking what is shown in fig. 3 as an example, the saturation degree of C point be C point to S point distance SC and D point to S point distance SD ratio
SC/SD。
Optionally, a kind of structural schematic diagram of computing unit 82 provided in an embodiment of the present invention is as shown in figure 9, can wrap
It includes:
Computation subunit 91, first determines that subelement 92 and second determines subelement 93；Wherein,
Computation subunit 91 is for calculating using the equal-energy white point in chromatic diagram as starting point, using the first chromaticity coordinate as terminal
The angle of preset coordinate axis in vector and chromatic diagram；
For ease of description, as shown in figure 3, establishing a new rectangular coordinate system o'x'y' in chromatic diagram, wherein coordinate
It is that equal-energy white point S, x' axis that origin o' is located in chromatic diagram is parallel with chromatic diagram y-axis and direction is consistent, x in y' axis and chromatic diagram
Axis is parallel and direction is consistent, and x' axis is vertical with y' axis to be met at o' point chromatic diagram is divided into four quadrants.
The chromaticity coordinate (x, y) of each nano wave length of chromatic diagram object-line (i.e. spectrum locus) is converted into new coordinate
Coordinate (x', y') under system, conversion formula are as follows:
In o'x'y' coordinate system, it is assumed that vector (x', y') and x' positive axis angle are α, then angle α calculation formula are as follows: α=
ARCTAN2 (y', x')=ARCTAN2 (x-0.3333, y-0.3333) (4)
Wherein, ARCTAN2 function representation bivariate arctan function, the codomain of bivariate tangent function be (- 180 °,
180 °), i.e., vector (x', y') is rotated counterclockwise by x' positive axis from the negative axis of x' again to the process of the negative axis of x', arctan function from-
180 ° change to -180 °.
Since the x' axis in o'x'y' coordinate system is parallel with the y-axis in chromatic diagram, vector (x', y') and x' positive axis
Angle is equal to the angle of the y-axis in vector (x', y') and chromatic diagram.
After coordinate is converted, ARCTAN2 function can make the angle α with dominant wavelength be incremented by, i.e., from 380nm to
The corresponding angle α 700nm is gradually increased, in this way, the α angle that can establish a corresponding incremental order from 380nm to 700nm is looked into
Look for table, the i.e. mapping table of angle and dominant wavelength.
In the embodiment of the present invention, wave spectrum chromaticity coordinate each in the spectrum locus of CIE-xy chromatic diagram is also established one
From 380nm to 700nm, wavelength corresponds to the look-up table of spectral chromaticity coordinates, the i.e. corresponding relationship of dominant wavelength and chromaticity coordinate
Table.
First determines subelement 92 for the corresponding relationship according to preset angle and dominant wavelength, determining and angle pair
The dominant wavelength answered；
Specifically, being selected if identified angle value is not present in the corresponding relationship of preset angle and dominant wavelength
It selects and the corresponding dominant wavelength of angle that identified angle is nearest is dominant wavelength corresponding with identified angle.Example
Such as, it is assumed that identified angle α is 6 degree, and in the corresponding relationship of preset angle and dominant wavelength, 558nm is corresponding to be
5.51 degree.559nm corresponding is 7.01 degree, and 6 degree closer from 5.51 degree, it is thus determined that color master corresponding with determined angle
Wavelength is 558nm.
Second determines subelement 93 for the corresponding relationship according to preset dominant wavelength and chromaticity coordinate, determining and institute
The corresponding chromaticity coordinate of determining dominant wavelength.
Those of ordinary skill in the art may be aware that list described in conjunction with the examples disclosed in the embodiments of the present disclosure
Member and algorithm steps can be realized with the combination of electronic hardware or computer software and electronic hardware.These functions are actually
It is implemented in hardware or software, the specific application and design constraint depending on technical solution.Professional technician
Each specific application can be used different methods to achieve the described function, but this realization is it is not considered that exceed
The scope of the present invention.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description
The specific work process of (if present), device and unit, can refer to corresponding processes in the foregoing method embodiment, herein no longer
It repeats.
In several embodiments provided herein, it should be understood that disclosed system (if present), device and side
Method may be implemented in other ways.For example, the apparatus embodiments described above are merely exemplary, for example, described
The division of unit, only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple units
Or component can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, institute
Display or the mutual coupling, direct-coupling or communication connection discussed can be through some interfaces, device or unit
Indirect coupling or communication connection can be electrical property, mechanical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme
's.
It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit
It is that each unit physically exists alone, can also be integrated in one unit with two or more units.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.
The foregoing description of the disclosed embodiments enables those skilled in the art to implement or use the present invention.
Various modifications to these embodiments will be readily apparent to those skilled in the art, as defined herein
General Principle can be realized in other embodiments without departing from the spirit or scope of the present invention.Therefore, of the invention
It is not intended to be limited to the embodiments shown herein, and is to fit to and the principles and novel features disclosed herein phase one
The widest scope of cause.
Claims (8)
1. a kind of city black and odorous water remote sensing recognition method characterized by comprising
Obtain the remotely-sensed data of water surface point to be identified；The remotely-sensed data of the water surface point to be identified includes: that the water surface is to be identified
The actual measurement Remote Sensing Reflectance of each wave band of point；
The spectral tristimulus value of water surface point to be identified is calculated based on the remotely-sensed data, comprising: will be described each by formula
The actual measurement Remote Sensing Reflectance of a wave band is converted to satellite band equivalent reflectivity；According to equivalent reflectivity pass corresponding with wave band
System, obtains the tristimulus values corresponding with trichromatic wave band of water surface point to be identified；According to tristimulus values and tristimulus values
Transformational relation the tristimulus values are converted into spectral tristimulus value；The formula are as follows:
Wherein, ReqFor satellite band equivalent reflectivity, Rrs(λ) is actual measurement Remote Sensing Reflectance, fSRF(λ) rings for satellite band spectrum
Answer function；
The saturation degree of water surface point to be identified is calculated based on the spectral tristimulus value；
If the saturation degree is less than preset threshold corresponding with identified water body affiliated area, determine that water surface point to be identified is
Black and odorous water.
2. the method according to claim 1, wherein the remotely-sensed data of water surface point to be identified includes: the water surface
The tristimulus values of pixel corresponding with water surface point to be identified in remote sensing images；
The spectral tristimulus value for calculating water surface point to be identified based on the remotely-sensed data includes:
The tristimulus values are converted to spectral tristimulus value by the transformational relation according to tristimulus values and tristimulus values.
3. the method according to claim 1, wherein described calculate the water surface based on the spectral tristimulus value
The saturation degree of point to be identified includes:
The spectral tristimulus value is converted to the first chromaticity coordinate of water surface point to be identified；
Calculate the second chromaticity coordinate of the dominant wavelength of the characterized color of the first chromaticity coordinate；
Determine that first distance and second coloration between the equal-energy white point in first chromaticity coordinate and chromatic diagram are sat
The second distance between equal-energy white point in mark and chromatic diagram；
The ratio of the first distance and the second distance is determined as to the saturation degree of water surface point to be identified.
4. according to the method described in claim 3, it is characterized in that, described calculate the characterized color of the first chromaticity coordinate
Corresponding second chromaticity coordinate of dominant wavelength includes:
It calculates using the equal-energy white point in chromatic diagram as starting point, using first chromaticity coordinate as the vector of terminal and the chromatic diagram
The angle of middle preset coordinate axis；
According to the corresponding relationship of preset angle and dominant wavelength, dominant wavelength corresponding with the angle is determined；
According to the corresponding relationship of preset dominant wavelength and chromaticity coordinate, color corresponding with identified dominant wavelength is determined
Spend coordinate.
5. a kind of city black and odorous water remote sensing recognition device characterized by comprising
Module is obtained, for obtaining the remotely-sensed data of water surface point to be identified；The remotely-sensed data of the water surface point to be identified includes: institute
State the actual measurement Remote Sensing Reflectance of each wave band of water surface point to be identified；
First computing module, for calculating the spectral tristimulus value of water surface point to be identified based on the remotely-sensed data；Include:
Second converting unit, for the actual measurement Remote Sensing Reflectance of each wave band to be converted to satellite band equivalent reflective by formula
Rate；Acquiring unit, for the corresponding relationship according to equivalent reflectivity and wave band, obtain water surface point to be identified and three primary colors
The corresponding tristimulus values of wave band；Third converting unit, will be described for the transformational relation according to tristimulus values and tristimulus values
Tristimulus values are converted to spectral tristimulus value；The formula are as follows:
Wherein, ReqFor satellite band equivalent reflectivity, Rrs(λ) is actual measurement Remote Sensing Reflectance, fSRF(λ) rings for satellite band spectrum
Answer function；
Second computing module, for calculating the saturation degree of water surface point to be identified based on the spectral tristimulus value；
Determining module, if for the saturation degree be less than preset threshold corresponding with identified water body affiliated area, determine described in
Water surface point to be identified is black and odorous water.
6. device according to claim 5, which is characterized in that the remotely-sensed data of the water surface point to be identified includes: the water surface
The tristimulus values of pixel corresponding with water surface point to be identified in remote sensing images；
First computing module includes:
The tristimulus values are converted to light for the transformational relation according to tristimulus values and tristimulus values by the first converting unit
Compose tristimulus values.
7. device according to claim 5, which is characterized in that second computing module includes:
4th converting unit, for the spectral tristimulus value to be converted to the first chromaticity coordinate of water surface point to be identified；
Computing unit, the second chromaticity coordinate of the dominant wavelength for calculating the characterized color of the first chromaticity coordinate；
First determination unit, for determining the first distance between the equal-energy white point in first chromaticity coordinate and chromatic diagram,
And the second distance between the equal-energy white point in second chromaticity coordinate and chromatic diagram；
Second determination unit, for the ratio of the first distance and the second distance to be determined as water surface point to be identified
Saturation degree.
8. device according to claim 7, which is characterized in that the computing unit includes:
Computation subunit, for calculating using the equal-energy white point in chromatic diagram as starting point, using first chromaticity coordinate as terminal
The angle of preset coordinate axis in vector and the chromatic diagram；
First determines subelement, for the corresponding relationship according to preset angle and dominant wavelength, the determining and angle pair
The dominant wavelength answered；
Second determines subelement, for the corresponding relationship according to preset dominant wavelength and chromaticity coordinate, determines and determines
The corresponding chromaticity coordinate of dominant wavelength.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510917811.2A CN105488488B (en) | 2015-12-10 | 2015-12-10 | City black and odorous water remote sensing recognition method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510917811.2A CN105488488B (en) | 2015-12-10 | 2015-12-10 | City black and odorous water remote sensing recognition method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105488488A CN105488488A (en) | 2016-04-13 |
| CN105488488B true CN105488488B (en) | 2018-12-21 |
Family
ID=55675460
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510917811.2A Active CN105488488B (en) | 2015-12-10 | 2015-12-10 | City black and odorous water remote sensing recognition method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105488488B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109142296A (en) * | 2018-08-16 | 2019-01-04 | 中国科学院合肥物质科学研究院 | The black smelly quick identification measuring method of urban water-body based on multi-source optical spectrum feature |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106874948A (en) * | 2017-02-08 | 2017-06-20 | 武汉海卓科科技有限公司 | A kind of black smelly water automatic identification and appraisal procedure |
| CN106950197B (en) * | 2017-03-03 | 2017-12-08 | 环境保护部卫星环境应用中心 | The Remotely sensed acquisition methods, devices and systems of sewage draining exit polluted-water |
| CN106918559B (en) * | 2017-05-04 | 2019-03-12 | 环境保护部卫星环境应用中心 | A kind of method and device detecting Water quality |
| CN107167431B (en) * | 2017-05-26 | 2019-07-05 | 中国科学院遥感与数字地球研究所 | A kind of black and odorous water recognition methods and system based on spectral index model |
| CN108333148B (en) * | 2018-01-11 | 2020-12-01 | 中国科学院遥感与数字地球研究所 | Method and device for identifying black and odorous water body |
| CN109269993A (en) * | 2018-10-15 | 2019-01-25 | 广州地理研究所 | A kind of recognition methods of dissolved oxygen, device, storage medium and equipment |
| CN110390277B (en) * | 2019-07-04 | 2021-07-06 | 中科卫星应用德清研究院 | Complex underlying surface water body identification method and black and odorous water body prediction method |
| CN111272662B (en) * | 2019-11-18 | 2022-07-26 | 深圳市深水水务咨询有限公司 | Urban black and odorous water body identification method based on remote sensing spectrum |
| CN110987825B (en) * | 2019-12-05 | 2022-09-02 | 南京师范大学 | Urban black and odorous water body grading method based on spectrum matching |
| CN113780071B (en) * | 2021-07-30 | 2024-01-30 | 中南安全环境技术研究院股份有限公司 | Black and odorous water body remote sensing identification system and method |
| CN117152602A (en) * | 2023-09-12 | 2023-12-01 | 广东省水利水电科学研究院 | Multisource satellite remote sensing image water body extraction method and system based on chromaticity index |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102542259A (en) * | 2011-12-29 | 2012-07-04 | 宁波工程学院 | Identification method for near-shore on-land water body |
| CN104361602A (en) * | 2014-11-26 | 2015-02-18 | 中国科学院遥感与数字地球研究所 | Water color detecting method and device based on MODIS image |
-
2015
- 2015-12-10 CN CN201510917811.2A patent/CN105488488B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102542259A (en) * | 2011-12-29 | 2012-07-04 | 宁波工程学院 | Identification method for near-shore on-land water body |
| CN104361602A (en) * | 2014-11-26 | 2015-02-18 | 中国科学院遥感与数字地球研究所 | Water color detecting method and device based on MODIS image |
Non-Patent Citations (2)
| Title |
|---|
| A Study on Retrieval Algorithm of Black Water Aggregation in Taihu Lake Based on HJ-1 Satellite Images;ZOU Lei et al;《IOP Conf. Series: Earth and Environmental Science》;20141231;第1-6页 * |
| 基于CIE1931的颜色识别方法研究;徐洋 等;《电视技术》;20121231;第36卷(第5期);第126-130页 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109142296A (en) * | 2018-08-16 | 2019-01-04 | 中国科学院合肥物质科学研究院 | The black smelly quick identification measuring method of urban water-body based on multi-source optical spectrum feature |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105488488A (en) | 2016-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105488488B (en) | City black and odorous water remote sensing recognition method and device | |
| Mokrzycki et al. | Colour difference∆ EA survey | |
| US8634640B2 (en) | Method and apparatus for selecting a color palette | |
| CN113766087A (en) | Color changing device | |
| CN105282396B (en) | Image processing apparatus and image processing method | |
| CN107702798B (en) | Method for measuring color accuracy of display and method for correcting color expression of display | |
| KR100441595B1 (en) | Color conversion system and method using the same | |
| US8861053B2 (en) | Method and apparatus for automated generation of a white ink separation out of CMYK data or RGB print data | |
| US20050219380A1 (en) | Digital camera for determining chromaticity coordinates and related color temperature of an object and method thereof | |
| US10728425B2 (en) | Color gamut mapping method and color gamut mapping device | |
| CN108020519A (en) | A kind of virtual multiple light courcess spectrum reconstruction method based on color constancy | |
| CN101478698A (en) | Image quality estimation device and method | |
| CN104361602B (en) | A kind of water body method for detecting color and device based on MODIS images | |
| CN103108469A (en) | Method and device for controlling output light of illuminating device and illuminating system | |
| CN103236251B (en) | Method and device for calibrating uniformity of splicing wall | |
| CN105070269B (en) | The color display and system of a kind of mobile terminal | |
| EP1058461A3 (en) | Colour image projector | |
| CN108401148A (en) | Method for executing automatic white balance to image | |
| WO2007078522A3 (en) | Perceptual color matching method between two different polychromatic displays | |
| Li et al. | Comparison of real colour gamuts using a new reflectance database | |
| CN106506902B (en) | Color management for light and color management method | |
| CN110277076A (en) | A kind of Color Mapping Approach and device, display equipment, medium | |
| CN203104900U (en) | Lighting control equipment and illumination system | |
| CN101147390B (en) | Color material mixing ratio calculating device, color material selection assisting device, and their methods and programs | |
| CN101585257A (en) | Complex frequency spectrum printing ink color matching method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |