CN105528760A - Cloud edge height matching method and system based on thermal infrared data - Google Patents
Cloud edge height matching method and system based on thermal infrared data Download PDFInfo
- Publication number
- CN105528760A CN105528760A CN201510885852.8A CN201510885852A CN105528760A CN 105528760 A CN105528760 A CN 105528760A CN 201510885852 A CN201510885852 A CN 201510885852A CN 105528760 A CN105528760 A CN 105528760A
- Authority
- CN
- China
- Prior art keywords
- data
- cloud
- impact point
- thermal infrared
- pixel
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000001514 detection method Methods 0.000 claims abstract description 53
- 238000010845 search algorithm Methods 0.000 claims abstract description 10
- 238000012952 Resampling Methods 0.000 claims description 13
- 238000013481 data capture Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 abstract description 6
- 230000004304 visual acuity Effects 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000000053 physical method Methods 0.000 description 2
- 230000036760 body temperature Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000001757 thermogravimetry curve Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
- G06T3/4061—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution by injecting details from different spectral ranges
-
- 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/10—Image acquisition modality
- G06T2207/10048—Infrared image
-
- 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
- G06T2207/30192—Weather; Meteorology
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Radiation Pyrometers (AREA)
Abstract
The present invention provides a cloud edge height matching method and system based on thermal infrared data. The method comprises a step of obtaining the visible light cloud detection data and the thermal infrared cloud height data in a same time area, a step of obtaining a target image according to the visible light cloud detection data, a step of obtaining a background image according to the thermal infrared cloud height data, a step of selecting one pixel in the target image as a target point, carrying out Euclidean distance transform on the pixel in the background image according to the target point, and obtaining Euclidean distance transform data, a step of obtaining the height value of the target point through a preset search algorithm according to the Euclidean distance transform data, and matching the target point. According to the method and the system, the matching of low-resolution thermal infrared cloud height data to the corresponding high-resolution visible light cloud detection data can be realized, a clout height distribution result with a high resolution, rich details and diversity is obtained, and the limit of obtaining the cloud height data by the thermal infrared technology in resolution is solved to some extent.
Description
Technical field
The present invention relates to image processing field, particularly relate to a kind of cloud brim height matching process based on Thermal Infrared Data and system.
Background technology
In traditional sensing techniques, the acquisition of the cloud level has multiple method, as utilized Laser-ceilometer, millimetre-wave radar and thermal infrared ceilometer etc.Above-mentioned measuring method comparatively wastes time and energy, and the acquisition cloud level that cannot be real-time in a big way, there is certain limitation.Along with the continuous maturation of satellite remote sensing technology, utilize the various parameters of remotely-sensed data inverting earth's surface, air and cloud, acquisition cloud level degrees of data that can be real-time on a large scale.
Remote sensing technology obtains, Retrieval of Cloud is high comprises physical method and method of geometry.Employing physical method obtains, namely calculate with vertical desuperheat rate according to Atmosphere and humidity profiles, the difference obtaining Black body temperature and reference temperature during clear sky as utilized remotely-sensed data than vertical desuperheat rate to estimate the cloud level, the method relies on the priori of cloud radiation characteristic and atmospheric outline etc., has certain error when inverting.Employing method of geometry obtains, and as utilized remotely-sensed data, take Fourier transform to find the corresponding point of cloud and cloud shade, and utilize geometric relationship to calculate the cloud level, the method is applicable to the calculating of Image of Flat Ground and computation process is comparatively complicated; Or utilize geometry method to calculate cloud shade, and set threshold value according to computational shadowgraph and the actual similarity detected between shade of remotely-sensed data, the cloud level that iteration is different obtains last cloud level degree.
Above based on the acquisition methods of remote sensing technology, the cloud level of acquisition is all more single cloud top or the height of cloud base, does not consider that cloud body is inner, the otherness of edge or different cloud height, and this otherness is more outstanding when high resolving power.Although thermal infrared technology conveniently can obtain the cloud height with otherness, but its resolution is lower, there is limitation compared to high-resolution visible data in tradition resampling result, compare geometric techniques be subject to a definite limitation when applying in the shape and detail characteristic of cloud.
Summary of the invention
For defect of the prior art, the invention provides a kind of cloud brim height matching process based on Thermal Infrared Data and system, for realizing low resolution thermal infrared cloud level Data Matching to corresponding High Resolution Visible Light cloud detection data.
First aspect, the invention provides a kind of cloud brim height matching process based on Thermal Infrared Data, described method comprises:
Obtain visible ray cloud detection data and the thermal infrared cloud level data of same time zone;
According to described visible ray cloud detection data, obtain the pixel set that target figure, described target figure form for visible ray cloud detection data edges pixel;
According to described thermal infrared cloud level data, background extraction figure, the pixel set that described Background forms for thermal infrared cloud level data edges pixel;
Choose a pixel on described target figure as impact point, and according to described impact point, Euclidean Distance Transform is carried out to the pixel in Background, obtain Euclidean Distance Transform data;
According to described Euclidean Distance Transform data, by preset search algorithm, obtain the height value of described impact point, described impact point is mated.
Preferably, before background extraction figure, described method also comprises:
Obtain the resolution of described visible ray cloud detection data;
Non-cloud region in described thermal infrared cloud level data is set as the territory, cloud sector of preset height value, obtains the thermal infrared cloud level data after setting;
Thermal infrared cloud level data after described setting are carried out resampling, and reduces non-cloud region after resampling, obtain and the thermal infrared cloud level data of described visible ray cloud detection data equal resolution.
Preferably, the height value of the described impact point of described acquisition, mates described impact point, comprising:
According to described Euclidean Distance Transform data, obtain Euclidean Distance Transform figure;
Choose region on described Euclidean Distance Transform figure as contour district, described contour district be with described impact point for the center of circle, there is the circle of preset search radius;
Calculate the Distance geometry of each candidate point and described impact point, described candidate point is the set of described contour district pixel;
According to the Distance geometry of described each candidate point and described impact point, calculate the distance weights of each candidate point and described impact point;
According to described each candidate point and the Distance geometry of described impact point and the distance weights of described each candidate point and described impact point, calculate the height value of described impact point, described impact point is mated.
Preferably, the Distance geometry of each candidate point of described calculating and described impact point, comprising:
Wherein, i is candidate point numbering, and n is candidate point number, d
ibe the distance of i-th candidate point and described impact point, Sd is the Distance geometry of each candidate point and described impact point.
Preferably, the distance weights of each candidate point of described calculating and described impact point, comprising:
w
i＝d
i/Sd，
Wherein, w
ibe the distance weights of i-th candidate point and described impact point.
Preferably, the height value of the described impact point of described calculating, comprising:
Wherein, h
ibe the height value of i-th candidate point, H is the height value of impact point.
Preferably, after the height value calculating described impact point, described method also comprises:
Choose one other pixel point on described target figure as another impact point, and according to another impact point described, Euclidean Distance Transform is carried out to the pixel in Background, obtain another Euclidean Distance Transform data;
According to another Euclidean Distance Transform data described, searched for by contour, obtain the height value of another impact point described, another impact point described is mated;
The step of the one other pixel point on described target figure as another impact point is chosen, until complete the pixel coupling on all described target figure described in execution.
Second aspect, the invention provides a kind of cloud brim height matching system based on Thermal Infrared Data, described system comprises:
Data capture unit, for obtaining visible ray cloud detection data and the thermal infrared cloud level data of same time zone;
Data extracting unit, for the visible ray cloud detection data that obtain according to described data capture unit and thermal infrared cloud level data, obtain target figure and Background, the pixel set that described target figure forms for visible ray cloud detection data edges pixel, the pixel set that described Background forms for thermal infrared cloud level data edges pixel;
Data matching unit, for the target figure that obtains according to data extracting unit and Background, chooses a pixel on described target figure as impact point, and obtains the height value of described impact point, mate described impact point.
As shown from the above technical solution, the invention provides a kind of cloud brim height matching process based on Thermal Infrared Data and system, obtain visible ray cloud detection data and the thermal infrared cloud level data of same time zone, contour search is carried out according to preset search algorithm, by low resolution thermal infrared cloud level Data Matching to corresponding High Resolution Visible Light cloud detection data by Euclidean Distance Transform.The present invention can realize low resolution thermal infrared cloud level Data Matching to corresponding High Resolution Visible Light cloud detection data, acquisition high resolving power, details are enriched and have the cloud level degree distribution results of otherness, solve the limitation of thermal infrared technical limit spacing cloud level degrees of data in resolution to a certain extent.
Accompanying drawing explanation
In order to be illustrated more clearly in disclosure embodiment or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only embodiments more of the present disclosure, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these figure.
The schematic flow sheet of a kind of cloud brim height matching process based on Thermal Infrared Data that Fig. 1 provides for one embodiment of the invention;
The schematic flow sheet of a kind of cloud brim height matching process based on Thermal Infrared Data that Fig. 2 provides for another embodiment of the present invention;
The structural representation of a kind of cloud brim height matching system based on Thermal Infrared Data that Fig. 3 provides for one embodiment of the invention;
Fig. 4 is preset search algorithm principle figure of the present invention.
Embodiment
Below in conjunction with drawings and Examples, the specific embodiment of the present invention is described in further detail.Following examples for illustration of the present invention, but are not used for limiting the scope of the invention.
Fig. 1 shows the schematic flow sheet of a kind of cloud brim height matching process based on Thermal Infrared Data that one embodiment of the invention provides, as shown in Figure 1, a kind of cloud brim height matching process based on Thermal Infrared Data of the present invention, described method comprises:
S11, the visible ray cloud detection data obtaining same time zone and thermal infrared cloud level data.
Be understandable that, thermal infrared cloud level data conveniently can obtain the cloud level degree with otherness, but its resolution is lower.And visible ray cloud detection data resolution is higher, but it cannot obtain that cloud body is inner, the otherness of edge or the different cloud level.Therefore, by low resolution thermal infrared cloud level Data Matching to corresponding High Resolution Visible Light cloud detection data, high resolving power can be obtained, details enriches and have the cloud level degree distribution results of otherness.
In the present embodiment, visible ray cloud detection data and thermal infrared cloud level data are all obtained by the form of view data, and its data message is embodied by the Pixel Information of view data.According to shape facility and geographic position, visible ray cloud detection data and thermal infrared cloud level data are carried out registration.
S12, according to described visible ray cloud detection data, obtain target figure.
Be understandable that, pixel is the elementary cell forming picture, is commonly used to the size representing photo resolution.Specifically, the pixel set that forms for visible ray cloud detection data edges pixel of target figure.
S13, according to described thermal infrared cloud level data, background extraction figure.
Specifically, the pixel set that forms for thermal infrared cloud level data edges pixel of Background.
S14, choose a pixel on described target figure as impact point, and according to described impact point, Euclidean Distance Transform is carried out to the pixel in Background, obtain Euclidean Distance Transform data.
Be understandable that, the pixel set that target figure forms for visible ray cloud detection data edges pixel, choose an impact point and be the pixel chosen in a target figure.
In prior art, range conversion is that also identifier space point is to the process of impact point distance in calculating, and it is finally transformed to gray level image bianry image, and wherein, the gray-scale value of each grid equals the distance that it arrives nearest impact point.Specifically, Euclidean Distance Transform (EDT) is high due to precision, conforms to actual range, applies more extensive.
S15, according to described Euclidean Distance Transform data, by preset search algorithm, obtain the height value of described impact point, described impact point is mated.
Specifically, as shown in Figure 4, preset search algorithm is contour search procedure, is namely the center of circle with impact point, and certain pixel is wide is search radius formation contour district, and relative target point from the close-by examples to those far off progressively expands hunting zone.For the pixel in contour district, according to its apart from distance of impact point by its height value weighting assignment to impact point.
The present embodiment can realize, according to shape and position feature, visible ray cloud detection data and thermal infrared cloud level data are carried out registration, by low resolution thermal infrared cloud level Data Matching to corresponding High Resolution Visible Light cloud detection data, acquisition high resolving power, details are enriched and have the cloud level degree distribution results of otherness, solve the limitation of thermal infrared technical limit spacing cloud level degrees of data in resolution to a certain extent.
Fig. 2 shows the schematic flow sheet of a kind of cloud brim height matching process based on Thermal Infrared Data that another embodiment of the present invention provides, as shown in Figure 2, a kind of cloud brim height matching process based on Thermal Infrared Data of the present invention, described method comprises:
S21, the visible ray cloud detection data obtaining same time zone and thermal infrared cloud level data.
Be understandable that, thermal infrared cloud level data are obtained by thermal infrared imager, namely the infrared energy distribution pattern utilizing infrared eye and optical imagery object lens to accept cloud layer is reflected on the light activated element of infrared eye, thus acquisition Infrared Thermogram, this thermal infrared picture figure is corresponding with the heat distribution field on cloud layer surface.
S22, resampling is carried out to thermal infrared cloud level data, obtain the thermal infrared cloud level data with described visible ray cloud detection data equal resolution.
Be understandable that, having equal resolution is the prerequisite of carrying out better matched data, therefore, needs by carrying out resampling to thermal infrared cloud level data and obtains the resolution identical with described visible ray cloud detection data.Resampling is image recovery process in essence, and it represents original image two dimension continuous function with the discrete digital image reconstruction of input, then samples by new image space-between and pixel position.
In addition, before resampling is carried out to thermal infrared cloud level data, in order to reduce the impact on cloud brim height in resampling process of non-cloud earth's surface, also comprise: the territory, the cloud sector non-cloud region in described thermal infrared cloud level data being set as preset height value, obtain the thermal infrared cloud level data after setting.Described non-cloud region is earth surface area, and preset height generally chooses this region cloud level minimum value.It should be noted that after resampling is carried out to thermal infrared cloud level data, also need to reduce non-cloud region.
S23, according to described visible ray cloud detection data, obtain target figure.
Specifically, the pixel set that forms for visible ray cloud detection data edges pixel of target figure.
S24, according to the thermal infrared cloud level data after described resampling, background extraction figure.
Specifically, the pixel set that forms for the thermal infrared cloud level data edges pixel after resampling of Background.
S25, choose a pixel on described target figure as impact point, and according to described impact point, Euclidean Distance Transform is carried out to the pixel in Background, obtain Euclidean Distance Transform data.
Be understandable that, the pixel set that target figure forms for visible ray cloud detection data edges pixel, choose an impact point and be the pixel chosen in a target figure.
S26, according to described Euclidean Distance Transform data, obtain Euclidean Distance Transform figure, and the region chosen on described Euclidean Distance Transform figure is as contour district.Be understandable that, as shown in Figure 4, contour district is be the center of circle with impact point, has the circle of preset search radius.The preset search radius deterministic process in described contour district is as follows: choose initial contour search radius r
1=1.5a, expands contour search radius according to default Changing Pattern gradually until the candidate point in described contour region reaches preset number.Described default Changing Pattern is for working as 1.5a≤r
iduring≤7.5a, r
i+1=r
i+ 0.5a; As 8a≤r
iduring≤12a, r
i+1=r
i+ a; Wherein, r
ibe the contour radius in i-th contour district, a is the width of a pixel.
Specifically, as shown in Figure 4, candidate point is the pixel set in contour district, and in the matching process, the number of candidate point needs to reach 2-8, and the contour search radius meeting this condition is suitable search radius.
S27, calculate the Distance geometry of each candidate point and described impact point.
Be understandable that, described candidate point is the set of described contour district pixel.Calculate the Distance geometry of each candidate point and described impact point, comprising:
Wherein, i is candidate point numbering, and n is candidate point number, d
ibe the distance of i-th candidate point and described impact point, Sd is the Distance geometry of each candidate point and described impact point.
S28, Distance geometry according to described each candidate point and described impact point, calculate the distance weights of each candidate point and described impact point.
Be understandable that, calculate the distance weights of each candidate point and described impact point, comprising:
w
i＝d
i/Sd，
Wherein, w
ibe the distance weights of i-th candidate point and described impact point.
S29, according to described each candidate point and the Distance geometry of described impact point and the distance weights of described each candidate point and described impact point, calculate the height value of described impact point, described impact point is mated.
Be understandable that, calculate the height value of described impact point, comprising:
Wherein, h
ibe the height value of i-th candidate point, H is the height value of impact point.
Meanwhile, the differing greatly in resolution due to visible ray cloud detection data and thermal infrared cloud level data, some cloud masses in small, broken bits in target figure there will be in corresponding Background does not have corresponding situation.Therefore, for some cloud masses in small, broken bits exceeding certain contour hunting zone, then the height value of assignment impact point is the mean value of the cloud level in Background.Wherein, exceed certain contour hunting zone and be set as 20a-25a, a is the width of a pixel.
Consider that detection range crosses conference under-represented, need detection range to limit.Thus, also need to mate again for exceeding the impact point limiting distance, described in specific as follows: the relevance of closing on pixel height is large, in order to obtain better matching effect, new pixel set is combined into the set of pixels of not mating in target figure, new background pixel set is combined into, repeated matching process with the set of pixels of having mated in target figure.
Be understandable that, realize low resolution thermal infrared cloud level Data Matching to need impact points all on target figure to mate to corresponding High Resolution Visible Light cloud detection data.And the coupling completing all impact points on target figure is a process choosing that different pixels point on target figure repeats step S25 to step S29 as impact point.
Specifically, choose one other pixel point on described target figure as another impact point, and according to another impact point described, Euclidean Distance Transform is carried out to the pixel in Background, obtain another Euclidean Distance Transform data; According to another Euclidean Distance Transform data described, searched for by contour, obtain the height value of another impact point described, another impact point described is mated; The step of the one other pixel point on described target figure as another impact point is chosen, until complete the pixel coupling on all described target figure described in execution.
The present embodiment can realize, according to shape facility and geographic position, visible ray cloud detection data and thermal infrared cloud level data are carried out registration, specifically, the present embodiment can realize, first according to geographical location information, secondly according to shape facility, visible ray cloud detection data and thermal infrared cloud level data being carried out registration.Thus, by low resolution thermal infrared cloud level Data Matching to corresponding High Resolution Visible Light cloud detection data, acquisition high resolving power, details are enriched and have the cloud level degree distribution results of otherness, solve the limitation of thermal infrared technical limit spacing cloud level degrees of data in resolution to a certain extent.
Fig. 3 shows the structural representation of a kind of cloud brim height matching system based on Thermal Infrared Data that one embodiment of the invention provides, and as shown in Figure 3, a kind of cloud brim height matching system 30 based on Thermal Infrared Data of the present invention, comprising:
Data capture unit 31, for obtaining visible ray cloud detection data and the thermal infrared cloud level data of same time zone.
Data extracting unit 32, for the visible ray cloud detection data that obtain according to data capture unit 31 and thermal infrared cloud level data, obtains target figure and Background.
Be understandable that, the pixel set that target figure forms for visible ray cloud detection data edges pixel, the pixel set that Background forms for thermal infrared cloud level data edges pixel.
Data matching unit 33, for the target figure that obtains according to data extracting unit 32 and Background, chooses a pixel on described target figure as impact point, obtains the height value of described impact point, mate described impact point.
Be understandable that, specifically, data matching unit 33 for choosing a pixel on target figure as impact point, and carries out Euclidean Distance Transform according to described impact point to the pixel in Background, obtains Euclidean Distance Transform data; According to described Euclidean Distance Transform data, by preset search algorithm, obtain the height value of described impact point, described impact point is mated.
Specifically, as shown in Figure 4, preset search algorithm is contour search procedure, is namely the center of circle with impact point, and certain pixel is wide is search radius formation contour district, and relative target point from the close-by examples to those far off progressively expands hunting zone.For the pixel in contour district, according to its apart from distance of impact point by its height value weighting assignment to impact point.
The present embodiment can realize low resolution thermal infrared cloud level Data Matching to corresponding High Resolution Visible Light cloud detection data, acquisition high resolving power, details are enriched and have the cloud level degree distribution results of otherness, solve the limitation of thermal infrared technical limit spacing cloud level degrees of data in resolution to a certain extent.
In sum, the invention provides a kind of cloud brim height matching process based on Thermal Infrared Data and system, obtain visible ray cloud detection data and the thermal infrared cloud level data of same time zone, contour search is carried out according to preset search algorithm, by low resolution thermal infrared cloud level Data Matching to corresponding High Resolution Visible Light cloud detection data by Euclidean Distance Transform.The present invention can realize low resolution thermal infrared cloud level Data Matching to corresponding High Resolution Visible Light cloud detection data, acquisition high resolving power, details are enriched and have the cloud level degree distribution results of otherness, solve the limitation of thermal infrared technical limit spacing cloud level degrees of data in resolution to a certain extent.
One of ordinary skill in the art will appreciate that: above each embodiment, only in order to technical scheme of the present invention to be described, is not intended to limit; Although with reference to foregoing embodiments to invention has been detailed description, those of ordinary skill in the art is to be understood that: it still can be modified to the technical scheme described in foregoing embodiments, or carries out equivalent replacement to wherein some or all of technical characteristic; And these amendments or replacement, do not make the essence of appropriate technical solution depart from the scope of the claims in the present invention.
Claims (8)
1., based on a cloud brim height matching process for Thermal Infrared Data, it is characterized in that, described method comprises:
Obtain visible ray cloud detection data and the thermal infrared cloud level data of same time zone;
According to described visible ray cloud detection data, obtain the pixel set that target figure, described target figure form for visible ray cloud detection data edges pixel;
According to described thermal infrared cloud level data, background extraction figure, the pixel set that described Background forms for thermal infrared cloud level data edges pixel;
Choose a pixel on described target figure as impact point, and according to described impact point, Euclidean Distance Transform is carried out to the pixel in Background, obtain Euclidean Distance Transform data;
According to described Euclidean Distance Transform data, by preset search algorithm, obtain the height value of described impact point, described impact point is mated.
2. the method for claim 1, is characterized in that, before background extraction figure, described method also comprises:
Obtain the resolution of described visible ray cloud detection data;
Non-cloud region in described thermal infrared cloud level data is set as the territory, cloud sector of preset height value, obtains the thermal infrared cloud level data after setting;
Thermal infrared cloud level data after described setting are carried out resampling, and reduces non-cloud region after resampling, obtain and the thermal infrared cloud level data of described visible ray cloud detection data equal resolution.
3. method as claimed in claim 2, is characterized in that the height value of the described impact point of described acquisition mates described impact point, comprising:
According to described Euclidean Distance Transform data, obtain Euclidean Distance Transform figure;
Choose region on described Euclidean Distance Transform figure as contour district, described contour district be with described impact point for the center of circle, there is the circle of preset search radius;
Calculate the Distance geometry of each candidate point and described impact point, described candidate point is the set of described contour district pixel;
According to the Distance geometry of described each candidate point and described impact point, calculate the distance weights of each candidate point and described impact point;
According to described each candidate point and the Distance geometry of described impact point and the distance weights of described each candidate point and described impact point, calculate the height value of described impact point, described impact point is mated.
4. method as claimed in claim 3, it is characterized in that, the Distance geometry of each candidate point of described calculating and described impact point, comprising:
Wherein, i is candidate point numbering, and n is candidate point number, d
ibe the distance of i-th candidate point and described impact point, Sd is the Distance geometry of each candidate point and described impact point.
5. method as claimed in claim 4, it is characterized in that, the distance weights of each candidate point of described calculating and described impact point, comprising:
w
i＝d
i/Sd，
Wherein, w
ibe the distance weights of i-th candidate point and described impact point.
6. method as claimed in claim 5, it is characterized in that, the height value of the described impact point of described calculating, comprising:
Wherein, h
ibe the height value of i-th candidate point, H is the height value of impact point.
7. method as claimed in claim 6, is characterized in that, after the height value calculating described impact point, described method also comprises:
Choose one other pixel point on described target figure as another impact point, and according to another impact point described, Euclidean Distance Transform is carried out to the pixel in Background, obtain another Euclidean Distance Transform data;
According to another Euclidean Distance Transform data described, searched for by contour, obtain the height value of another impact point described, another impact point described is mated;
The step of the one other pixel point on described target figure as another impact point is chosen, until complete the pixel coupling on all described target figure described in execution.
8., based on a cloud brim height matching system for Thermal Infrared Data, it is characterized in that, described system comprises:
Data capture unit, for obtaining visible ray cloud detection data and the thermal infrared cloud level data of same time zone;
Data extracting unit, for the visible ray cloud detection data that obtain according to described data capture unit and thermal infrared cloud level data, obtain target figure and Background, the pixel set that described target figure forms for visible ray cloud detection data edges pixel, the pixel set that described Background forms for thermal infrared cloud level data edges pixel;
Data matching unit, for the target figure that obtains according to data extracting unit and Background, chooses a pixel on described target figure as impact point, and obtains the height value of described impact point, mate described impact point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510885852.8A CN105528760B (en) | 2015-12-04 | 2015-12-04 | A kind of cloud brim height matching process and system based on Thermal Infrared Data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510885852.8A CN105528760B (en) | 2015-12-04 | 2015-12-04 | A kind of cloud brim height matching process and system based on Thermal Infrared Data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105528760A true CN105528760A (en) | 2016-04-27 |
| CN105528760B CN105528760B (en) | 2018-09-25 |
Family
ID=55770967
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510885852.8A Expired - Fee Related CN105528760B (en) | 2015-12-04 | 2015-12-04 | A kind of cloud brim height matching process and system based on Thermal Infrared Data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105528760B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109409186A (en) * | 2017-08-17 | 2019-03-01 | 哈曼国际工业有限公司 | Driver assistance system and method for object detection and notice |
| CN110634259A (en) * | 2019-09-25 | 2019-12-31 | 云南电网有限责任公司电力科学研究院 | Forest fire detection method |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101566692A (en) * | 2009-05-26 | 2009-10-28 | 吉林大学 | Method for detecting cloud height by utilizing cloud shadow information in satellite remote sensing data |
-
2015
- 2015-12-04 CN CN201510885852.8A patent/CN105528760B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101566692A (en) * | 2009-05-26 | 2009-10-28 | 吉林大学 | Method for detecting cloud height by utilizing cloud shadow information in satellite remote sensing data |
Non-Patent Citations (3)
| Title |
|---|
| CATHERINE NAUD ET AL: "Assessment of multispectral ATSR2 stereo cloud-top height retrievals", 《REMOTE SENSING OF ENVIRONMENT》 * |
| 章文星 等: "地基热红外亮温遥感云底高度可行性的模拟研究", 《地球物理学报》 * |
| 黄兴友 等: "云底高度的激光云高仪、红外测云仪以及云雷达观测比对分析", 《量子电子学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109409186A (en) * | 2017-08-17 | 2019-03-01 | 哈曼国际工业有限公司 | Driver assistance system and method for object detection and notice |
| CN109409186B (en) * | 2017-08-17 | 2023-09-26 | 哈曼国际工业有限公司 | Driver assistance system and method for object detection and notification |
| CN110634259A (en) * | 2019-09-25 | 2019-12-31 | 云南电网有限责任公司电力科学研究院 | Forest fire detection method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105528760B (en) | 2018-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20200226413A1 (en) | Fast and robust multimodal remote sensing images matching method and system | |
| CN104574347B (en) | Satellite in orbit image geometry positioning accuracy evaluation method based on multi- source Remote Sensing Data data | |
| CN103727930B (en) | A kind of laser range finder based on edge matching and camera relative pose scaling method | |
| Huang et al. | Cloud motion estimation for short term solar irradiation prediction | |
| US9396583B2 (en) | Method of modelling buildings on the basis of a georeferenced image | |
| US9582885B2 (en) | Zonal underground structure detection method based on sun shadow compensation | |
| CN110689622A (en) | Synchronous positioning and composition algorithm based on point cloud segmentation matching closed-loop correction | |
| CN112837315B (en) | Deep learning-based transmission line insulator defect detection method | |
| CN104751465A (en) | ORB (oriented brief) image feature registration method based on LK (Lucas-Kanade) optical flow constraint | |
| CN105046251A (en) | Automatic ortho-rectification method based on remote-sensing image of environmental No.1 satellite | |
| CN106096497B (en) | A kind of house vectorization method for polynary remotely-sensed data | |
| CN107909018B (en) | Stable multi-mode remote sensing image matching method and system | |
| Sheng et al. | Automated georeferencing and orthorectification of Amazon basin-wide SAR mosaics using SRTM DEM data | |
| CN103136525A (en) | High-precision positioning method for special-shaped extended target by utilizing generalized Hough transformation | |
| CN103198447A (en) | Wind arrow field real-time measuring method based on satellite cloud pictures | |
| CN103745453B (en) | Urban residential areas method based on Google Earth remote sensing image | |
| CN104766341A (en) | Detection method for construction change of remote-sensing image based on DSM and kernel density estimation | |
| CN109117851A (en) | A kind of video image matching process based on lattice statistical constraint | |
| CN107590782B (en) | High-resolution optical image thick cloud removing method based on full convolution network | |
| CN105528760A (en) | Cloud edge height matching method and system based on thermal infrared data | |
| CN114387446A (en) | Automatic water body extraction method for high-resolution remote sensing image | |
| Li et al. | An automated system for image-to-vector georeferencing | |
| CN103473782B (en) | Least square matching method based on object-space vertical double-surface element | |
| CN111862005A (en) | Method and system for accurately positioning tropical cyclone center by using synthetic radar image | |
| CN113592744B (en) | Geometric fine correction method suitable for high-resolution remote sensing image |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180925 Termination date: 20191204 |