CN108364264A - A kind of ocean temperature monitoring method and system based on unmanned plane infrared remote sensing technology - Google Patents
A kind of ocean temperature monitoring method and system based on unmanned plane infrared remote sensing technology Download PDFInfo
- Publication number
- CN108364264A CN108364264A CN201810122385.7A CN201810122385A CN108364264A CN 108364264 A CN108364264 A CN 108364264A CN 201810122385 A CN201810122385 A CN 201810122385A CN 108364264 A CN108364264 A CN 108364264A
- Authority
- CN
- China
- Prior art keywords
- unmanned plane
- temperature
- infrared
- remote sensing
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012544 monitoring process Methods 0.000 title claims abstract description 26
- 238000005516 engineering process Methods 0.000 title claims abstract description 17
- 238000012937 correction Methods 0.000 claims abstract description 32
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 16
- 238000012417 linear regression Methods 0.000 claims abstract description 5
- 230000005540 biological transmission Effects 0.000 claims abstract description 4
- 230000005855 radiation Effects 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 3
- 239000013535 sea water Substances 0.000 claims description 3
- 238000013210 evaluation model Methods 0.000 claims description 2
- 230000009467 reduction Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 12
- 230000001360 synchronised effect Effects 0.000 abstract description 5
- 238000009529 body temperature measurement Methods 0.000 abstract description 4
- 238000013461 design Methods 0.000 abstract description 4
- 230000007613 environmental effect Effects 0.000 abstract description 3
- 241001269238 Data Species 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000001931 thermography Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 239000000498 cooling water Substances 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 2
- 239000010865 sewage Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 230000005457 Black-body radiation Effects 0.000 description 1
- 208000028257 Joubert syndrome with oculorenal defect Diseases 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 240000005373 Panax quinquefolius Species 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 239000003637 basic solution Substances 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000009432 framing Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000008235 industrial water Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000004861 thermometry Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/0037—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the heat emitted by liquids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/48—Thermography; Techniques using wholly visual means
- G01J5/485—Temperature profile
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0077—Imaging
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radiation Pyrometers (AREA)
Abstract
The invention discloses a kind of ocean temperature monitoring methods and system based on unmanned plane infrared remote sensing technology, flying platform using unmanned plane as thermal infrared remote sensing, Differential positioning equipment after integrated PPK, design realizes the equidistant concentric circles fitting parameter correcting algorithm of wide-angle infrared lens barrel distortion, the improvement without photo control point geometric correction algorithm is carried out on the basis of traditional polynomial method, coordinate the high-precision location datas of PPK, realize the promotion of the image position precision of unmanned plane infrared remote sensing most critical, after finally carrying out temperature retrieval using Heat transmission equation, one layer of actual measurement aligning step is set again, utilize the sea measured data of acquisition synchronous with unmanned plane, equation of linear regression is established using ENVI analysis softwares, obtain calibration model, improve temperature measurement accuracy.Invention significantly improves the monitoring efficiency of thermal discharge waterfrom nuclear power plant, ensure that the data needed for environmental protection administration are normal, and greatly reduce the monitoring cost of warm water discharge.
Description
Technical field
The invention belongs to marine monitoring technology field, especially a kind of ocean temperature based on unmanned plane infrared remote sensing technology
Monitoring method and system.
Background technology
21st century is the century of ocean, and during developing and utilizing marine environment and resource, the activity of the mankind is not
It is evitable to bring pollution and destruction to marine environment.With the raising of environmental consciousness, it has been recognized that is " first dirty
Administered after dye " development model be not response environment problem correct thinking, it is believed that prevent environment by prior precautionary measures
The generation of problem is only basic solution.Pollution flux is being weighed in the index of marine pollution degree, temperature is then
One important indicator, its estuary can be promoted by coming from cities and towns sewage draining exit sanitary sewage, industrial water, power-plant cooling water
The temperature in marine site, the nuclear power station recirculated cooling water that especially borders on the sea discharge, since displacement is huge, can result in large stretch of marine site Wen Sheng,
Fish production and marine ecology are generated and greatly destroyed.
At present for it is all kinds of enter sea pollution source temperature monitoring method have using ship carry out the water surface acquisition monitoring, satellite it is distant
Sense monitoring and manned aircraft Aerial remote sensing monitoring.Water surface monitoring can provide zonal real-time, discrete, more accurate
Water surface temperature, the disadvantage is that can only reflect that dotted and linear sea and the temperature along depth of water direction index.Space satellite
Remote sensing thermometry be limited in that can only according to satellite transit time carry out image capturing, and Satellite Image Spatial Resolution and
Temperature measurement accuracy is low, is suitable for the monitoring of the not high big-and-middle regional extent of required precision.Air remote sensing monitoring system has speed fast
The characteristics of high resolution, is conducive to the source and the whereabouts that differentiate pollutant.But the shortcomings that manned aircraft air remote sensing is to landing
There is particular/special requirement in place, and time is long and implementation process is complicated, and flexibility is poor, and cost is high.
Invention content
In view of the deficiencies of the prior art, this application provides a kind of landings conveniently, implements flexible, temperature measurement accuracy and positioning accurate
Degree height and unmanned plane ocean temperature monitoring method and system at low cost, in addition while unmanned plane monitors, mistake during acquisition
The infrared remote sensing data of border satellite are used for the overall trend of assistant analysis distribution of water temperature.
To achieve the above object, the first technical solution of the application use is:One kind being based on unmanned plane infrared remote sensing skill
The ocean temperature monitoring method of art, is as follows:
Step 1:Unmanned Aerial Vehicle Data acquisition platform is built, using unmanned plane as platform, carries thermal infrared sensor equipment, union
At PPK reception antennas, to the precise positioning of unmanned plane;
Step 2:Nuclear power operating mode, meteorological data are collected, area course line is surveyed in planning, determines, fine situation smaller in wind-force
Lower progress field data acquisition, it is synchronous to carry out ship, the actual measurement of buoy temperature;
Step 3:Data processing is screened effective infrared image according to unmanned plane POS data in ArcGIS, and is designed
A kind of infrared image spatial information representation based on ArcEngine is realized, by establishing image vector file, is assigned each
Pixel position and temperature property.According to unmanned plane POS data, space ginseng of infrared image under the conditions of unmanned plane is just penetrated is determined
Examine information:Projected coordinate system, position, range ultimately produce under the conditions of unmanned plane just penetrates, the corresponding spatial data of infrared image;
Step 4:The camera lens of distortion correction, the infrared thermography that the present invention uses is wide-angle lens, and it is abnormal that there are wide-angles
Become, Distortion Law and parameter can determine, and have uniformity.For the technical parameter of this camera lens, design realizes one kind
Barrel distortion model of fit based on equidistant concentric circles, obtains the deformation parameter of the corresponding position of discrete radius, is corrected
Spatial data afterwards.In order to ensure the accuracy of temperature, the 10% of edge is dismissed so that other remaining to Duplication between image
10%;
Step 5:Geometric correction, the key that seawater thermometric is carried out using unmanned plane are to ensure the position essence of orthography
Degree, the present invention design realize it is a kind of on the basis of traditional polynomial method be suitable for geometric correction of the sea without photo control point situation,
The accurate angle information provided using elements of exterior orientation calculates deviation post;
Step 6:Temperature retrieval, it is contemplated that in shooting process under different height air, steam, incident angle difference, temperature
Atmospheric radiation correction, incident angle amendment are carried out when spending inverting, to improve temperature measurement accuracy, design is obtained using sea ship, buoy
Measured value is taken, realizes using temperature difference rule between ENVI analysis software different zones, school is established using linear regression fit
The method of positive model and precision evaluation model finally obtains ocean temperature by correction;
Step 7:The temperature value of the calculated infrared image of step 6 is assigned to by image reduction according to position correspondence
Then spatial data is converted into raster data in such a way that temperature value is pixel value, obtained by the spatial data of infrared image
Existing spatial referencing information has the image of temperature information again.
The application use second of technical solution be:A kind of ocean temperature monitoring based on unmanned plane infrared remote sensing technology
System, including unmanned plane and the infrared sensing equipment for measuring ocean temperature, the infrared sensor are built in unmanned plane
The camera lens of cabin bottom, infrared sensor leans out to come from unmanned plane cabin bottom.
Further, the unmanned plane is equipped with integrated PPK reception antennas.
Further, above system further includes:Buoy carries temperature sensor at the universal sensor interface of buoy,
Water temperature is monitored.
Further, above system further includes:Water quality instrument, for measuring seawater water temperature.
The present invention uses above technical scheme, can obtain following technique effect:Air remote sensing and sea actual measurement are mutually tied
It closes, is suitable for the temperature monitoring in large area marine site, has the features such as inexpensive, convenient, flexible, speed is fast.Choose UAV flight
Infrared sensor reduces the cost of similar temperature measuring application as sampling instrument, is improved using the geometric correction of no photo control point
The positioning accuracy of infrared image;Means are surveyed using synchronous sea to be corrected infrared data, can significantly improve temperature survey
Accuracy of measurement.This method causes temperature to rise environmental disruption and provides efficiently and effectively monitoring means for the pollution sources that border on the sea.
Description of the drawings
Fig. 1 is the ocean temperature monitoring method flow chart based on unmanned plane infrared remote sensing technology;
Fig. 2 is unmanned plane schematic diagram;
Fig. 3 is infrared sensing equipment schematic diagram;
Fig. 4 is buoy usage state diagram.
Specific implementation mode
To make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with the accompanying drawings with specific implementation case
The present invention will be described in detail.The application is using low profile thermal infrared thermoviewer as load, using fixed-wing unmanned plane as flying platform,
Synchronous progress sea ship, buoy actual measurement are corrected and are verified for infrared data inverting, i.e. unmanned plane+infrared sensor+sea
Survey alignment technique method.
The present embodiment provides a kind of ocean temperature monitoring methods based on unmanned plane infrared remote sensing technology, specifically:
1, in this example, the understanding of the factors such as environment, landform is carried out to measuring marine site first, selectes ward as nothing
Man-machine landing place of taking off, stands in local built-in on the spot face CORS, for Differential positioning after the PPK of unmanned plane.To seashore wind-force
Situation is understood, and determines the wind-force period less than normal.
2, the flight course planning of unmanned plane is carried out to measuring marine site, in this example, the range area to be measured is 20 × 6km long
120 square kilometres of wide marine sites are arranged Duplication to prevent leak-stopping bat situation to aerial photograph in unmanned plane planning, determine flight
Height is set as 600m, Duplication 20%.In the actual measurement of sea, discharge outlet and intake and several key areas is selected to lay floating
Mark, ship actual measurement divide 3 pieces of regions since region is larger, and every piece of region is measured using two ships.
3, the data processing based on ArcGIS, after obtaining infrared image, in conjunction with system for flight control computer POS data, camera
Real-time attitude angular data, infrared thermography photographed data carry out data prediction process it is as described below.
1) data screening
For unmanned plane in flight course, corner has larger shaking, this partial data is weeded out, according to overlapping
Rate requires to select available POS points.POS and infrared video time are all made of GPS time, are accurate to second, stringent synchronization, and reference has been selected
It carries out cutting frame in the time to infrared video file of POS points, obtains POS points and the infrared image with radiation value in the same time.
2) spatial resolution is determined
The resolution ratio of infrared image is 640*480, i.e., has 640*480 in the field range of infrared thermography camera lens
The radiation value of pixel.The actual size that pixel projects to sea is directly proportional to unmanned plane shooting height, 600m height unmanned planes
Visual field scope is 692.8m × 497.1m, size 1.1m × 1.0m of each pixel.
3) framing
According to condition is just penetrated, infrared image, which is navigated to the parameter needed on correct geographical location, to be had:Image center
Coordinate, shooting height, image actual size, course angle.
Picture centre point coordinates:In unmanned plane scanning process, infrared thermography is equipped with correction holder, ensures camera lens
Scanning sea is just penetrated, unmanned plane POS coordinate datas can be used as the center of the sea field range captured by infrared thermography
Point coordinates.
Shooting height:That is the difference of aircraft altitude height and the height above sea level on sea level, unmanned plane fly higher, and shooting regards
Field range is bigger.
Image actual size:According to field of view angle and unmanned plane shooting height, the field range of 600m height is as infrared
Sea actual range 692.8m × 497.1m of thermal imaging system shooting.
Course angle:With the clockwise angle of due north.
Vector file is established using ArcGIS Desktop, condition is just penetrated according to the length of side, angle drafting since central point
Under image capturing range.
4) pixel figure
Grid, lattice are generated according to picture position, Pixel size, pixel ranks number in the corresponding actual range of infrared image
Subnumber is 640*480, and each grid represents a pixel of infrared image, the pixel figure arranged in order.
4, lens distortion calibration, since damp state influences thermo parameters method situation, the crawl large area marine site in a damp state
Temperature field information certainly will use wide-angle lens, letter is carried out using wide-angle lens in this example in view of nobody current flying speed
Breath acquisition, and to be this example will solve the problems, such as barrel distortion situation existing for wide-angle lens.
According to this infrared camera wide-angle feature, the non standard template of an equidistant concentric circles is selected, one kind is used
Modification method is as follows:
It a., can be with if equidistantly the density of concentric circles is bigger close to the minimum concentric circles deformation at deformation center minimum
Think that minimum concentric circles does not have deformation.
B. after equidistant concentric circles is projected to by ideal camera lens on CCD or equidistant concentric circles.
C. barrel distortion is continuous, i.e. the deformation of the connected point of image is not mutated.According to this theoretical foundation, with making
Non-standard template alignment relative Calibration is carried out to wide-angle lens.
To the non-equidistant concentric circles on the fault image that actually obtains, equidistantly adjusted according to minimum concentric circles,
The deformation parameter of discrete radius corresponding position is obtained in this process.Assuming that lens centre and image plane center are straight in same
On line, then the different pixels on image are due to the radius r differences to picture centre, distortion situation is just different.
According to the deformation principle of continuity, use respectively linear regression and quartic polynomial respectively to the barrel distortion of wide-angle lens into
Row fitting, obtains the Deformation Law of wide-angle lens.
It is fitted barrel distortion rule with primary recurrence:
st(r)=ar+b
Wherein a, b are distortion factors.
It is fitted barrel distortion rule with quartic polynomial:
st(r)=ar4+br3+cr2+dr+e
Wherein a, b, c, d, e are distortion factors.
Scene carries out image distortion experiment by the sea for selection.Aircraft is acquired more using more rotary wind types under different height
Width image selects shooting height for 600 meters of a width thermal infrared imagery.Thermal infrared images imaging clearly, different atural objects (including sea
Face, harbour, shore protection, road, square) heat radiation imaging resolution is higher.And it is above-mentioned apparent in ground actual acquisition using GPS
Culture point and as more modified foundation,
With primary return barrel distortion rule is fitted with quartic polynomial:
st(r)=- 0.005190740r-0.9432991
It is fitted barrel distortion rule with quartic polynomial:
st(r)=0.017430r4-0.084409r3+0.085233r2-0.047001r+1.003114
5, geometric correction is done in geometric distortion correction, image geometry distortion caused by changing for sensor elements of exterior orientation,
Wherein elements of exterior orientation includes coordinate Xs of the photo centre point S in selected ground space coordinate systems、Ys、Zs, and photography moment picture
The posture of planeω, k, i.e. pitch angle, roll angle, course angle.
The collinearity equation for describing image object-image relation is as follows,
Wherein f is the focal length of photographic camera, and (x, y) and (X, Y, Z) is the image coordinate and ground coordinate of culture point, (Xs,
Ys, Zs) be video camera ground coordinate, aiFor elements of exterior orientation,Identified spin matrix element.
Using elements of exterior orientation as variable, differential process is carried out to above formula, while considering under approximately perpendicular photography conditions,ω, k are the angles of very little, can be approximately zero, and the first order public affairs that elements of exterior orientation variation influences picpointed coordinate can be obtained
Formula
DX as can be seen from the above equations、dYs、dZsInfluence with dk to whole picture image is to bring it about translation, rotation and scaling
Etc. linear changes, and only haveImage can be just set to generate nonlinear distortion with d ω.
When camera lens tilt, picture point offset and picture point to image center point apart from directly proportional, image center point it is inclined
Shifting amount is minimum, and edge offset is maximum.High-performance holder is carried during unmanned plane during flying, high-frequency adjusts the appearance of camera lens in real time
State ensures that camera lens is vertically shot so thatω (pitching, roll) angle reaches minimum, in order to ensure accuracy, infrared heat at
As fixing high-precision dual-axis inclinator on instrument, to obtain the accurate attitude angle parameter of camera lens, because infrared thermography is per second
32 frames are imaged, setting inclinator angle acquisition frequency is also 32Hz, the pitch angle at synchronous recording infrared thermography each moment
And roll angle.In 600m height, when the roll of camera lens, pitch angle are 0.1 °, caused image is 1m with ground offset, complete
The requirement of full up 5 meters of foot ground space resolution ratio, so, to improve data-handling efficiency as far as possible, not to pitching, roll angle
Image less than 0.1 ° does geometric correction, and the image primarily directed to pitching, roll angle more than 0.1 ° does geometric correction.
In geometric correction, multinomial correction method has very strong versatility, is suitable for the geometric correction of landform flat country.
And this example measures on large area sea, landform is flatter, so geometric distortion correction is carried out using multinomial correction method, and herein
On the basis of, realize the geometric correction process without photo control point.
Several steps crucial in the multinomial geometric correction method of no photo control point are described in detail below;
(1) image control point is chosen, on the image with spatial referencing information, according to the principle being evenly distributed, simultaneously
Since from image center to edge, distortion, so when selecting control point, leans on proximal edge more in rule is become larger.
(2) each control point calibration recoil mark is calculated, is calculated according to flying height, the angle of view, course and just penetrates condition
Under image position, several points are chosen, using unmanned plane flying height, the angle of view, attitude angle (pitching, roll) according to triangle letter
Number calculates the coordinate after each control point geometric distortion, is just under control a little pair.
(3) determine fine correction model, in order to enable positional accuracy higher after correction, using 15 pairs of control points, to image into
Row cubic polynomial corrects.
6, temperature retrieval and correction, temperature retrieval use radiation transfer equation method.Radiation transfer equation method is red by a heat
The radiation energy that outer tunnel (in atmospheric window) obtains calculates temperature.The energy that thermal infrared radiation sensor receives includes mainly three
Part:Earth's surface heat radiation after atmospheric extinction;The earth surface reflection of Downward atmospheric long-wave radiation;Air uplink radiation, such as following formula:
T in formulaSFor atural object temperature;Lai ↑For air uplink wave spectrum heat radiation;Lai ↓For downward atmospheric heat radiation;εiFor water
Table emissivity;τiFor atmospheric transmittance;Under indicate i parameter be corresponding wave spectrum channel i weighted mean.
Parameter in formula can pass through radiation transmission in addition to water meter emissivity (being determined according to the look-up table of foundation)
It simulates to seek.In actual application, surveyed with MODTRAN atmospheric radiation transmissions and sounding balloon and surface weather station
Data come assist calculate atmospheric parameter.
Under the flying height of hundreds of meters of unmanned plane, atmospheric effect with the variation of Vertical Temperature and moisture profile it is progressive gradually
Become.Above-mentioned radiation transfer equation can be expressed as following form:
B(T0)=τ [ε B (T1)+(1-ε)I↓]+I↑
In formula:T0For the brightness temperature of thermal infrared sensor detection;T1For temperature of surface water;τ is atmospheric transmissivity;ε is water meter
Radiance;B(T0) be thermal infrared sensor received by caloradiance;B(T1) be water meter in sensor wave band section
Blackbody radiation intensity;I↑And I↓Respectively up and down caloradiance of the air in detecting band section.
In order to reduce the Systematic Errors of infrared temperature inverting, in conjunction with boat measurement, buoy measurement data, with data acquisition time
For sequence, foundation takes photo by plane infrared temperature angle value and the one-to-one data of sea observed temperature value to set.It is established using ENVI softwares
In the equation of linear regression taken photo by plane between infrared temperature and sea observed temperature value by atmospheric correction, and then applied to infrared
The correction process of load data, inverting obtain true sea surface temperature value.
Present invention likewise provides the precision evaluation of the temperature value for acquisition and verification methods, and the time is used using non-modeling
Match point boat-carrying walks boat and fixed point buoy observation data make independent sample inspection to model algorithm result precision.Precision evaluation uses
Model inversion end value obtains precision evaluation and error statistics with the warm comparing mode in field measurement sea, such as formula:
In formula:Tci, TjiThe field measurement sea temperature value and model inversion of i-th of matched data pair respectively in verify data
End value;For the mean error of inversion algorithm.
Verify 30% of sample number not less than actual measurement sample number.
Boat is walked with time match point boat-carrying using non-modeling and fixed point buoy observation data make model algorithm result precision
Independent sample is examined, and total matched data collection is thus divided into two parts, most for correcting, least a portion of for examining
It tests, into the interference that sample may be not present, realizes that independent sample is examined.
7, temperature value after above-mentioned steps is assigned to the spatial data of infrared image according to position correspondence, then
Spatial data is converted into raster data in such a way that temperature value is pixel value, obtain existing spatial referencing information has temperature again
The image of information, finally complete complete thermo parameters method figure.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
Any one skilled in the art in the technical scope of present disclosure, according to the technique and scheme of the present invention and its
Inventive concept is subject to equivalent substitution or change, should be covered by the protection scope of the present invention.
Claims (8)
1. a kind of ocean temperature monitoring method based on unmanned plane infrared remote sensing technology, which is characterized in that be as follows:
Step 1:Build Unmanned Aerial Vehicle Data acquisition platform;
Step 2:Nuclear power operating mode, meteorological data are collected, area course line is surveyed in planning, then carries out field data acquisition;
Step 3:Data processing screens effective infrared image according to unmanned plane POS data, by establishing image vector file,
Assign each pixel position and temperature property;Determine infrared image unmanned plane spatial referencing information;
Step 4:Distortion correction obtains the correspondence of discrete radius using the barrel distortion model of fit based on equidistant concentric circles
The deformation parameter of position obtains revised spatial data;
Step 5:Geometric correction is realized the geometric correction of unmanned plane image using the multinomial algorithm of no photo control point, utilizes foreign side
The precision information that bit element provides calculates deviation post;
Step 6:Temperature retrieval and correction carry out Atmospheric radiation correction, incident angle amendment using Heat transmission equation, pass through line
Property regression fit establishes calibration model and precision evaluation model, the ocean temperature after being corrected;
Step 7:The temperature value of the calculated infrared image of step 6 is assigned to according to position correspondence infrared by image reduction
Then spatial data is converted into raster data in such a way that temperature value is pixel value, obtained existing by the spatial data of image
Spatial referencing information has the image of temperature information again.
2. a kind of ocean temperature monitoring method based on unmanned plane infrared remote sensing technology, feature exist according to claim 1
In using the barrel distortion model of fit based on equidistant concentric circles, specially:According to the deformation principle of continuity, use respectively linear
It returns and quartic polynomial is respectively fitted the barrel distortion of wide-angle lens, obtain the Deformation Law of wide-angle lens;
It is fitted barrel distortion size with primary recurrence:
st(r)=ar+b
Wherein, r is radius, and a, b are distortion factors;
It is fitted barrel distortion size with quartic polynomial:
st(r)=ar4+br3+cr2+dr+e
Wherein, r is radius, and a, b, c, d, e is distortion factor.
3. a kind of ocean temperature monitoring method based on unmanned plane infrared remote sensing technology, feature exist according to claim 1
In the multinomial algorithm of the no photo control point realizes the geometric correction of unmanned plane image, specially:
Using elements of exterior orientation as variable, differential process is carried out to the collinearity equation of geometric correction, obtains elements of exterior orientation variation pair
The first order formula that picpointed coordinate influences;When choosing control point in infrared image, marginal portion distortion is big, and choosing is by proximal edge
Point.
4. a kind of ocean temperature monitoring method based on unmanned plane infrared remote sensing technology, feature exist according to claim 1
In the temperature retrieval and correction are specially monitored while being surveyed using ship, buoy temperature, school in unmanned aerial vehicle remote sensing
Timing utilizes one-to-one measured value and infrared temperature angle value, establishes equation of linear regression, one-to-one temperature value time difference
It is required that control is within a certain period of time.
5. a kind of ocean temperature based on unmanned plane infrared remote sensing technology monitors system, which is characterized in that including unmanned plane and use
In the infrared sensing equipment for measuring ocean temperature, the infrared sensor is built in the cabin bottom of unmanned plane, infrared sensor
Camera lens from unmanned plane cabin bottom lean out come.
6. a kind of ocean temperature based on unmanned plane infrared remote sensing technology monitors system according to claim 5, feature exists
In the unmanned plane is equipped with integrated PPK reception antennas.
7. according to a kind of ocean temperature monitoring system based on unmanned plane infrared remote sensing technology of claim 5 or 6, feature
It is, further includes:Buoy carries temperature sensor at the universal sensor interface of buoy, is monitored to water temperature.
8. a kind of ocean temperature based on unmanned plane infrared remote sensing technology monitors system according to claim 7, it is characterized in that,
Further include:Water quality instrument, for measuring seawater water temperature.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810122385.7A CN108364264A (en) | 2018-02-07 | 2018-02-07 | A kind of ocean temperature monitoring method and system based on unmanned plane infrared remote sensing technology |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810122385.7A CN108364264A (en) | 2018-02-07 | 2018-02-07 | A kind of ocean temperature monitoring method and system based on unmanned plane infrared remote sensing technology |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108364264A true CN108364264A (en) | 2018-08-03 |
Family
ID=63005074
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810122385.7A Pending CN108364264A (en) | 2018-02-07 | 2018-02-07 | A kind of ocean temperature monitoring method and system based on unmanned plane infrared remote sensing technology |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108364264A (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110161533A (en) * | 2019-06-17 | 2019-08-23 | 国家卫星气象中心(国家空间天气监测预警中心) | Unattended satellite optical load intelligence spectral radiometric calibration system |
CN110375860A (en) * | 2019-07-09 | 2019-10-25 | 浙江大学 | A kind of inshore power station warm water discharge environmental detection set and method based on infrared image processing |
CN110599412A (en) * | 2019-08-15 | 2019-12-20 | 中国科学院遥感与数字地球研究所 | Remote sensing data processing method and system based on unmanned aerial vehicle |
CN111245945A (en) * | 2020-01-15 | 2020-06-05 | 北京工业大学 | Marine oil spilling supervisory systems based on buoy and unmanned aerial vehicle control |
CN111310309A (en) * | 2020-01-20 | 2020-06-19 | 中国矿业大学 | Unmanned aerial vehicle-based thermal infrared image temperature inversion correction method |
CN111784678A (en) * | 2020-07-03 | 2020-10-16 | 中煤航测遥感集团有限公司 | Temperature acquisition method, device and equipment based on thermal infrared image and storage medium |
CN112880842A (en) * | 2020-12-29 | 2021-06-01 | 中国海洋大学 | Water surface temperature infrared thermometer without sky temperature measurement |
CN113834572A (en) * | 2021-08-26 | 2021-12-24 | 电子科技大学 | Method for removing drift of temperature measurement result of non-refrigeration thermal imager of unmanned aerial vehicle |
CN114442667A (en) * | 2022-01-29 | 2022-05-06 | 国家海洋环境预报中心 | Sea skin temperature observation method and system based on composite wing unmanned aerial vehicle |
US11441951B2 (en) * | 2019-08-21 | 2022-09-13 | Fredrick S. Solheim | Apparatus and methods for remote measurement of sea surface temperature |
CN115144095A (en) * | 2022-09-06 | 2022-10-04 | 自然资源部第一海洋研究所 | Method and system for determining background water temperature of warm water drainage of operated nuclear power plant based on remote sensing |
CN111310309B (en) * | 2020-01-20 | 2024-06-28 | 中国矿业大学 | Temperature inversion correction method based on thermal infrared image of unmanned aerial vehicle |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551275A (en) * | 2009-04-30 | 2009-10-07 | 上海航遥信息技术有限公司 | Technical method of vehicular multispectral scanner for monitoring industrial warm discharge water |
CN105241429A (en) * | 2015-09-22 | 2016-01-13 | 中国科学院上海技术物理研究所 | Extraction method for offshore industrial warm discharge water based on aerial remote sensing |
CN105527657A (en) * | 2016-01-14 | 2016-04-27 | 河海大学 | Large-area farmland crop water status monitoring method and system based on unmanned aerial vehicle infrared thermal image acquisition |
CN107193025A (en) * | 2017-07-14 | 2017-09-22 | 陕西舜洋电子科技有限公司 | A kind of orientable unmanned environment measuring station |
-
2018
- 2018-02-07 CN CN201810122385.7A patent/CN108364264A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551275A (en) * | 2009-04-30 | 2009-10-07 | 上海航遥信息技术有限公司 | Technical method of vehicular multispectral scanner for monitoring industrial warm discharge water |
CN105241429A (en) * | 2015-09-22 | 2016-01-13 | 中国科学院上海技术物理研究所 | Extraction method for offshore industrial warm discharge water based on aerial remote sensing |
CN105527657A (en) * | 2016-01-14 | 2016-04-27 | 河海大学 | Large-area farmland crop water status monitoring method and system based on unmanned aerial vehicle infrared thermal image acquisition |
CN107193025A (en) * | 2017-07-14 | 2017-09-22 | 陕西舜洋电子科技有限公司 | A kind of orientable unmanned environment measuring station |
Non-Patent Citations (4)
Title |
---|
伊晓东 等: "小波神经网络在遥感测温数据拟合中的应用研究", 《测绘科学技术》 * |
武文波 主编: "《现代测绘科技丛书 遥感技术基础》", 30 May 2000 * |
谢颖 等: "广角镜头桶形畸变的多项式修正方法", 《合肥工业大学学报》 * |
贺佳惠 等: "核电站近岸温度场航空热红外遥感测量数据处理研究", 《国土资源遥感》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110161533A (en) * | 2019-06-17 | 2019-08-23 | 国家卫星气象中心(国家空间天气监测预警中心) | Unattended satellite optical load intelligence spectral radiometric calibration system |
CN110161533B (en) * | 2019-06-17 | 2024-05-03 | 国家卫星气象中心(国家空间天气监测预警中心) | Unattended satellite optical load intelligent radiation calibration system |
CN110375860A (en) * | 2019-07-09 | 2019-10-25 | 浙江大学 | A kind of inshore power station warm water discharge environmental detection set and method based on infrared image processing |
CN110599412A (en) * | 2019-08-15 | 2019-12-20 | 中国科学院遥感与数字地球研究所 | Remote sensing data processing method and system based on unmanned aerial vehicle |
US11441951B2 (en) * | 2019-08-21 | 2022-09-13 | Fredrick S. Solheim | Apparatus and methods for remote measurement of sea surface temperature |
CN111245945A (en) * | 2020-01-15 | 2020-06-05 | 北京工业大学 | Marine oil spilling supervisory systems based on buoy and unmanned aerial vehicle control |
CN111310309A (en) * | 2020-01-20 | 2020-06-19 | 中国矿业大学 | Unmanned aerial vehicle-based thermal infrared image temperature inversion correction method |
CN111310309B (en) * | 2020-01-20 | 2024-06-28 | 中国矿业大学 | Temperature inversion correction method based on thermal infrared image of unmanned aerial vehicle |
CN111784678A (en) * | 2020-07-03 | 2020-10-16 | 中煤航测遥感集团有限公司 | Temperature acquisition method, device and equipment based on thermal infrared image and storage medium |
CN112880842A (en) * | 2020-12-29 | 2021-06-01 | 中国海洋大学 | Water surface temperature infrared thermometer without sky temperature measurement |
CN113834572B (en) * | 2021-08-26 | 2023-05-12 | 电子科技大学 | Unmanned aerial vehicle non-refrigeration type thermal imager temperature measurement result drift removal method |
CN113834572A (en) * | 2021-08-26 | 2021-12-24 | 电子科技大学 | Method for removing drift of temperature measurement result of non-refrigeration thermal imager of unmanned aerial vehicle |
CN114442667A (en) * | 2022-01-29 | 2022-05-06 | 国家海洋环境预报中心 | Sea skin temperature observation method and system based on composite wing unmanned aerial vehicle |
CN115144095B (en) * | 2022-09-06 | 2022-11-11 | 自然资源部第一海洋研究所 | Method and system for determining background water temperature of warm water drainage of operated nuclear power plant based on remote sensing |
CN115144095A (en) * | 2022-09-06 | 2022-10-04 | 自然资源部第一海洋研究所 | Method and system for determining background water temperature of warm water drainage of operated nuclear power plant based on remote sensing |
US11830635B1 (en) | 2022-09-06 | 2023-11-28 | First Institute of Oceanography, Ministry of Natural Resources | Method and system for determining background water temperature of thermal discharge from operating nuclear power plants based on remote sensing |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108364264A (en) | A kind of ocean temperature monitoring method and system based on unmanned plane infrared remote sensing technology | |
CN103413272B (en) | Low spatial resolution multi-source Remote Sensing Images Space Consistency bearing calibration | |
CN106767714B (en) | Improve the equivalent mismatch model multistage Calibration Method of satellite image positioning accuracy | |
CN103345737B (en) | A kind of UAV high resolution image geometric correction method based on error compensation | |
US20220276332A1 (en) | Stellar atmospheric refraction measurement correction method based on collinearity of refraction surfaces | |
CN108603790A (en) | Thermal imaging system based on UAV system and method | |
CN104463470B (en) | Intertidal Habitats investigation method based on unmanned plane | |
CN106871787A (en) | Large space line scanning imagery method for three-dimensional measurement | |
CN106679676B (en) | A kind of monoscopic multifunctional optical sensor and implementation method | |
CN106873004B (en) | Based on the adaptive in-orbit geometry calibration method of high rail level array camera of solar elevation | |
CN112710311B (en) | Automatic planning method for three-dimensional live-action reconstruction aerial camera points of terrain adaptive unmanned aerial vehicle | |
CN105526906B (en) | Wide-angle dynamic high precision laser angular measurement method | |
CN106971408A (en) | A kind of camera marking method based on space-time conversion thought | |
CN105444778B (en) | A kind of star sensor based on imaging geometry inverting is in-orbit to determine appearance error acquisition methods | |
CN104361563B (en) | GPS-based (global positioning system based) geometric precision correction method of hyperspectral remote sensing images | |
CN110503687A (en) | A kind of aerial photoelectric measurement platform object localization method | |
CN108761271A (en) | A kind of power grid screen of trees detection method and system | |
CN112833861A (en) | Surveying and mapping method and surveying and mapping system based on oblique photography large-scale topographic map | |
CN115511956A (en) | Unmanned aerial vehicle imaging positioning method | |
CN108225563B (en) | Field environment sky polarization modeling alignment measuring device | |
Xiang et al. | Monitoring the thermal discharge of hongyanhe nuclear power plant with aerial remote sensing technology using a UAV platform | |
CN107784633A (en) | Suitable for the unmanned plane image calibrating method of plane survey | |
CN109883400A (en) | Fixed station Automatic Targets and space-location method based on YOLO-SITCOL | |
CN116824079A (en) | Three-dimensional entity model construction method and device based on full-information photogrammetry | |
CN107423858A (en) | A kind of urban planning method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180803 |