CN110210453A - A kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images determines method and system - Google Patents
A kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images determines method and system Download PDFInfo
- Publication number
- CN110210453A CN110210453A CN201910516536.1A CN201910516536A CN110210453A CN 110210453 A CN110210453 A CN 110210453A CN 201910516536 A CN201910516536 A CN 201910516536A CN 110210453 A CN110210453 A CN 110210453A
- Authority
- CN
- China
- Prior art keywords
- oil tank
- elevation
- solar
- remote sensing
- sensing images
- 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 32
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 claims abstract description 57
- 238000003384 imaging method Methods 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims description 26
- 101150095230 SLC7A8 gene Proteins 0.000 claims description 6
- 101150085091 lat-2 gene Proteins 0.000 claims description 6
- 101100511466 Caenorhabditis elegans lon-1 gene Proteins 0.000 claims description 3
- 101100021570 Caenorhabditis elegans lon-3 gene Proteins 0.000 claims description 3
- 101150044140 Slc7a5 gene Proteins 0.000 claims description 3
- 239000003921 oil Substances 0.000 description 112
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 241000127225 Enceliopsis nudicaulis Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 239000003502 gasoline Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000003703 image analysis method Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Life Sciences & Earth Sciences (AREA)
- Astronomy & Astrophysics (AREA)
- Animal Husbandry (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Agronomy & Crop Science (AREA)
- Primary Health Care (AREA)
- Remote Sensing (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Marine Sciences & Fisheries (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
A kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images determines method and system.Include: (1) in remote sensing images, is reference with the sharp atural object near oil tank, calculates solar azimuth；(2) the remote sensing images imaging date and solar elevation are determined according to solar azimuth；(3) in remote sensing images, it is reference with the sharp atural object near oil tank, calculates satellite aximuth；(4) according to oil tank, nearby sharply the shadow length of atural object, projected length and solar elevation calculate elevation of satellite；(5) selection needs to calculate the oil tank of oil reserve, calculates tank elevation using oil tank top shadow length, projected length, solar elevation, solar azimuth, elevation of satellite and satellite aximuth；(6) tank(storage)capacity is calculated according to tank elevation and oil tank radius；(7) oil tank filling rate is calculated using the shade of inner wall of oil tank and oil tank outer wall；(8) oil tank amount of storage is calculated using tank(storage)capacity and oil tank filling rate.
Description
Technical field
The invention belongs to field of remote sensing image processing, and in particular to optical strip when being imaged according to domestic Optical remote satellite
Part and the attitude of satellite, in conjunction with the method for geometrical principle estimation oil tank amount of storage.
Background technique
Remote sensing technology one of Spatial Information Technology with the fastest developing speed since being 21st century, it is a wide range of, quick
The ability of effective acquisition earth surface information is widely applied in different industries.Remote sensing images obtain earth's surface information compared to
The means of ground survey are at low cost, obtain information quantification, accuracy height, therefore are a kind of high-tech hands for substituting conventional method
Section.Petroleum is tellurian non-renewable resources, is mainly used as fuel oil and gasoline, referred to as " industrial blood ".One state
The oil reserve in family and area has very important influence, the method that tradition obtains petroleum storage amount to politics, economic, military affairs
For manual research, such method labor intensive, and time-effectiveness and regional limitation obtain the petroleum storage in the whole world in the short time
It is standby almost impossible, therefore remote sensing technology is a kind of economic, effective, quick and accurate means for obtaining oil reserve.
The method that remote sensing images are applied to oil tank detection has very much, such as Canny detection, Hough transform, symmetrical feature etc.
Algorithm has in document [1] general introduction to introduce.And the method and few calculated about oil tank amount of storage, document [1] also proposed oneself
Innovative oil tank detection method and reserve analysis method, document [2] propose oil tank volume acquiring method.Oil tank volume is sought
The basic principle of method is the diameter and height for seeking out oil tank, can be calculated using quadrature formula.For existing remote sensing images
The method for detecting oil tank, the prior art have the following problems: (1) most methods are for the purpose of detecting oil tank, and disregard
Calculate oil tank volume and oil tank amount of storage；(2) method for seeking volume and oil tank amount of storage assumes solar azimuth and satellite position is
Given value.
The following are bibliography:
[1] Zhang Jinfang, Wang Qingquan, Wang Siyu etc., oil tank detection and reserves point based on high-resolution optical remote sensing image
Analysis method [P].The patent No.: 201710727161
[2] Yu Shengtao, the oil tank volume based on high-resolution remote sensing image seek research [D].Dalian Maritime University,
2018。
Summary of the invention
Technology of the invention solves the problems, such as: in view of the deficienciess of the prior art, proposing a kind of special based on remote sensing images
The oil tank of sign stores analysis method, can not obtain solar elevation, solar azimuth, satellite altitude at imaging moment target
Solar elevation, solar azimuth, elevation of satellite and satellite aximuth at angle and satellite aximuth or imaging moment target
When inaccurate, by the shade and projection properties in image taking date and image, oil tank is calculated using geometrical relationship and is stored
Amount
The technical solution of the invention is as follows:
A kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images determines method, and steps are as follows:
(1) in remote sensing images, it is reference with the sharp atural object near oil tank, calculates solar azimuth；
(2) the remote sensing images imaging date and solar elevation are determined according to solar azimuth；
(3) in remote sensing images, it is reference with the sharp atural object near oil tank, calculates satellite aximuth；
(4) according to oil tank, nearby sharply the shadow length of atural object, projected length and solar elevation calculate elevation of satellite；
(5) selection needs to calculate the oil tank of oil reserve, utilizes oil tank top shadow length, projected length, solar elevation, too
Positive azimuth, elevation of satellite and satellite aximuth calculate tank elevation；
(6) tank(storage)capacity is calculated according to tank elevation and oil tank radius；
(7) oil tank filling rate is calculated using the shade of inner wall of oil tank and oil tank outer wall；
(8) oil tank amount of storage is calculated using tank(storage)capacity and oil tank filling rate.
The oil tank is floating roof formula oil tank.
The calculating of solar azimuth SoA is as follows:
Wherein, shade longitude and latitude of the sharp atural object vertex on ground is (Lon1, Lat1), contact of the sharp atural object on ground
Point longitude and latitude is (Lon2, Lat2).
The remote sensing images imaging date and solar elevation are determined according to solar azimuth, specifically:
(2.1) the solar hour angle T of imaging date atural object imaging is calculatedangle:
Tangle=(12-Tlocal) × 15,
Wherein, TlocalIt is satellite at the place that fixed latitude is passed by；
(2.2) solar declination δ is calculated:
δ=arcsin (0.39795cos (0.98563 (N-173)))；
Wherein, N is the day product that the date is imaged in image；
(2.3) solar elevation that it's the atural object moment pasts satellite is calculated using solar hour angle, solar declination and atural object latitude
SoH, calculation method use following formula:
SoH=arcsin (sinLatsin δ+cosLatcos δ cosTangle)；
Wherein, sharp atural object geographic latitude and longitude is (Lon, Lat) in remote sensing images；
(2.4) formula for calculating solar azimuth using solar elevation, atural object latitude and solar declination is as follows:
(2.5) the solar azimuth angle value being calculated in step (1) is substituted into step (2.4), it is counter to inquire into obtain imaging day
Phase and solar elevation SoH.
Satellite aximuth SaA calculation formula is as follows:
Wherein, contact point longitude and latitude of the sharp atural object on ground is (Lon2, Lat2), and sharp atural object vertex is in remote sensing images
In projection longitude and latitude be (Lon3, Lat3).
According to oil tank, nearby sharply the shadow length of atural object, projected length and solar elevation calculate elevation of satellite, tool
Body are as follows:
(4.1) shadow length calculation formula is as follows:
(4.2) projected length calculation formula is as follows:
(4.3) elevation of satellite calculation formula is as follows:
Selection needs to calculate the oil tank of oil reserve, utilizes oil tank top shadow length, projected length, solar elevation, the sun
Azimuth, elevation of satellite and satellite aximuth calculate tank elevation, specifically:
Tank elevation H calculation formula is as follows:
Wherein, θ=SaA-SoA, BC are shadow spots and sharp atural object top of the sharp atural object vertex on ground in remote sensing images
The distance between the subpoint of point.
Tank(storage)capacity calculates as follows:
Wherein, D is the oil tank abhiseca diameter measured on remote sensing images.
Oil tank filling rate R calculates as follows:Wherein, SoutFor oil tank outer wall shade, SinFor inner wall of oil tank shade；
Oil tank amount of storage VrCalculation method is as follows: Vr=V × R.
A kind of oil tank amount of storage that the oil tank amount of storage based on Characteristics of The Remote Sensing Images according to determines that method is realized is true
Determine system,
Solar azimuth computing module: for being reference with the sharp atural object near oil tank, calculating too in remote sensing images
Positive azimuth；
Solar elevation computing module: the remote sensing images imaging date and altitude of the sun are determined according to solar azimuth
Angle；
Satellite aximuth computing module: in remote sensing images, it is reference with the sharp atural object near oil tank, calculates satellite side
Parallactic angle；
Elevation of satellite computing module: according to shadow length, projected length and the altitude of the sun of atural object sharp near oil tank
Angle calculates elevation of satellite；
Tank elevation computing module: selection needs to calculate the oil tank of oil reserve, is grown using oil tank top shadow length, projection
Degree, solar elevation, solar azimuth, elevation of satellite and satellite aximuth calculate tank elevation；
Oil tank amount of storage computing module: tank(storage)capacity is calculated according to tank elevation and oil tank radius；Using inner wall of oil tank and
The shade of oil tank outer wall calculates oil tank filling rate；Oil tank amount of storage is calculated using tank(storage)capacity and oil tank filling rate.
Compared with the prior art, the invention has the advantages that:
(1) transit time is calculated using the longitude and latitude and Sun synchronization repeating orbit characteristic of the position of target, relative to it
He directly utilizes the time of image source file more acurrate at method, because image source file provides the north that sub-satellite point passes by
Capital time or UTC time, rather than when the place of imageable target.
(2) solar elevation and solar azimuth are calculated when the place on the date and target passed by using satellite, other oil
Tankage size calculation method assume the two parameters be given value, actually may not all satellite transit commercial cities the two ginsengs are provided
Number, and provide parameter be substar solar elevation and solar azimuth, rather than the solar elevation of imageable target and
Solar azimuth.
(3) elevation of satellite and satellite aximuth when being shot using box counting algorithm, other same tankometers
Calculation method assumes that the two parameters are given value, actually may not all satellite transit commercial cities the two parameters are provided, and scheme
As the elevation of satellite and satellite aximuth of upper different location cannot be all with the same values.
(4) capacity of oil tank has not only been calculated, but also has proposed the acquisition methods of the filling rate of oil tank, to obtain oil
The amount of storage of tank.Ningbo harbour oil tank reserve analysis is carried out using this method, with the true difference range of practical oil reserve 5%
Within.
Detailed description of the invention
Fig. 1 is calculation flow chart of the invention；
Fig. 2 is that the sun of the present invention and the attitude of satellite calculate schematic diagram.
Specific embodiment
As shown in Figure 1, a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images proposed by the present invention determines method, step is such as
Under:
(1) in remote sensing images, it is reference with the sharp atural object near oil tank, calculates solar azimuth；It is oily in the present invention
Tank is floating roof formula oil tank.
(2) the remote sensing images imaging date and solar elevation are determined according to solar azimuth；
(3) in remote sensing images, it is reference with the sharp atural object near oil tank, calculates satellite aximuth；
(4) according to oil tank, nearby sharply the shadow length of atural object, projected length and solar elevation calculate elevation of satellite；
(5) selection needs to calculate the oil tank of oil reserve, utilizes oil tank top shadow length, projected length, solar elevation, too
Positive azimuth, elevation of satellite and satellite aximuth calculate tank elevation；
(6) tank(storage)capacity is calculated according to tank elevation and oil tank radius；
(7) oil tank filling rate is calculated using the shade of inner wall of oil tank and oil tank outer wall；
(8) oil tank amount of storage is calculated using tank(storage)capacity and oil tank filling rate.
Specifically, the embodiment of each step is as follows:
Step 1: nearby there is the sharp atural object of certain altitude for reference with oil tank, calculate solar azimuth.
Solar azimuth is to be rotated clockwise to the angle of sunray using direct north as starting point.In remote sensing images
Can be with the discriminating direction solar azimuth of shade, but oil tank should not be elected to be the object of reference of solar azimuth calculating, because of satellite
The shade of oil tank may be blocked by oil tank itself when shooting, cause solar azimuth to calculate mistake, and oil tank is circle, is not present
Obvious characteristic point.Therefore, the sharp atural object of certain altitude nearby is advisable as solar azimuth calculating using oil tank, common atural object
Including electric pole, flagpole, lamppost etc..
By taking figure in attached drawing 2 as an example, it is assumed that the sharp atural object for having certain altitude in image is AO, and the height of AO is H, image
Middle atural object geographic latitude and longitude is (Lon, Lat), and shade BO, then using B as starting point, O is that the direction of terminal is solar azimuth.
The longitude and latitude of B point is (Lon1, Lat1) in image, and the longitude and latitude of O point is (Lon2, Lat2), then solar azimuth SoA calculates public
Formula is as follows:
Step 2: image imaging date and solar elevation are determined from solar azimuth.
Solar azimuth is one-to-one with the imaging date, it is known that solar azimuth can then calculate the imaging date.
In calculating with day product indicate imaging the date, such as January 1 day product be 1 day 1,2 months day product be 31.Leap year and non-
Leap year differs 1 day, little to the calculating image of solar azimuth, so with 1 year 365 for day product.Assuming that image is imaged
The day product on date is N, and satellite is T at the place that fixed latitude is passed bylocal, then imaging date atural object can be calculated
The solar hour angle T of imagingangle:
Tangle=(12-Tlocal)×15 (2)
Then solar declination is calculated, solar declination is the line of the sun and the earth's core and the angle of equatorial plane.Solar declination
Using year as mechanical periodicity, according to documents and materials, the variation of solar declination can be used following skeleton symbol and calculate δ:
δ=arcsin (0.39795cos (0.98563 (N-173))) (3)
Then the solar elevation that it's the atural object moment pasts satellite can be calculated using solar hour angle, solar declination and atural object latitude
SoH, calculation method use following formula:
SoH=arcsin (sinLatsin δ+cosLatcos δ cosTangle) (4)
The formula for calculating solar azimuth using solar elevation, atural object latitude and solar declination is as follows:
The unknown number that solar azimuth is calculated in this step is N, TlocalIt is the fixed value of satellite, therefore using in step 1
Calculated solar azimuth, which substitutes into formula, can then find out the imaging date, so as to find out solar elevation.
Step 3: nearby there is the sharp atural object of certain altitude for reference with oil tank, calculate satellite aximuth.
Satellite aximuth is to be rotated clockwise to the angle of satellite shooting using direct north as starting point.In remote sensing images
Can be with the discriminating direction satellite position of height displacement, but oil tank should not be elected to be the object of reference of satellite aximuth calculating, because defending
The projected base of oil tank may be blocked by oil tank itself when star is shot, and cause satellite position to calculate mistake, and oil tank is circle,
There is no obvious characteristic points.Therefore, the sharp atural object of certain altitude nearby is advisable as satellite aximuth calculating using oil tank, often
Insight object includes electric pole, flagpole, lamppost etc..
By taking figure in attached drawing 2 as an example, it is assumed that the sharp atural object for having certain altitude in image is AO, is projected as CO, is then with C
Starting point, O are that the direction of terminal is satellite aximuth.The longitude and latitude of C point is (Lon3, Lat3) in image, and the longitude and latitude of O point is
(Lon2, Lat2), then solar azimuth SaA calculation formula is as follows:
Step 4: shade, height displacement and the solar elevation of the sharp atural object of certain altitude nearby calculate using oil tank
Elevation of satellite.
The height of atural object can be calculated, according to height displacement and elevation of satellite according to the length of shade and solar elevation
The height of atural object can be calculated, therefore, elevation of satellite can be calculated using the length of shade, height displacement and solar elevation.
Shadow length calculation formula is as follows:
Projected length calculation formula is as follows:
Elevation of satellite calculation formula is as follows:
Step 5: selection needs to calculate the oil tank of oil reserve, utilizes the distance of oil tank top shade and height displacement, altitude of the sun
Angle, solar azimuth, elevation of satellite and satellite aximuth calculate tank elevation.
With reference to the accompanying drawings 2, in triangle OBC, there are following relationships:
Wherein, θ=SaA-SoA
Projection and shade bottom due to oil tank can be blocked by oil tank itself, so LshadowAnd LprojectOn the image without
Method is accurately measured, it is assumed that tank elevation H, it is known that solar elevation and elevation of satellite, then tank elevation H and Lshadow、
LprojectRelationship there are following relationships:
Lshadow=HtanSoH (11)
Lproject=HtanSaH (12)
Formula (11) and formula (12) are substituted into formula (10), then it can be high by solar elevation, solar azimuth, satellite
Degree angle and satellite aximuth calculate tank elevation H, and calculation formula is as follows:
Step 6: tank(storage)capacity is calculated by volume equation using tank elevation and oil tank radius.
It has found out tank elevation using satellite image, then has measured the diameter D of oil tank abhiseca on the image, then it can be by quadrature
Formula calculates tank(storage)capacity V, and calculation formula is as follows:
Step 7: oil tank filling rate is calculated using the shade of inner wall of oil tank and oil tank outer wall.
The method that oil tank filling rate can refer to bibliography [1] and [2] utilizes oil tank outer wall shade SoutAnd inner wall of oil tank
Shade SinRatio determine, SoutAnd SinIt can be measured on the image, oil tank filling rate R calculation formula is as follows:
Step 8: oil tank amount of storage is calculated using tank(storage)capacity and oil tank filling rate.
Oil tank amount of storage VrCalculation method is as follows:
Vr=V × R (16)
Embodiment:
This patent carries out 3 phases prison by a definite date to Ningbo Daxie harbour respectively with No. two satellites of high score and high scape No.1 satellite
It surveys, obtains the oil storage rate of 12 floating top oil tanks in harbour, as shown in Table 1, and compared with ground investigation, it can be deduced that
Conclusion, error is within 5%.
The oil storage rate data comparison table of one floating top oil tank of table
Claims (10)
1. a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images determines method, it is characterised in that steps are as follows:
(1) in remote sensing images, it is reference with the sharp atural object near oil tank, calculates solar azimuth；
(2) the remote sensing images imaging date and solar elevation are determined according to solar azimuth；
(3) in remote sensing images, it is reference with the sharp atural object near oil tank, calculates satellite aximuth；
(4) according to oil tank, nearby sharply the shadow length of atural object, projected length and solar elevation calculate elevation of satellite；
(5) selection needs to calculate the oil tank of oil reserve, utilizes oil tank top shadow length, projected length, solar elevation, sun side
Parallactic angle, elevation of satellite and satellite aximuth calculate tank elevation；
(6) tank(storage)capacity is calculated according to tank elevation and oil tank radius；
(7) oil tank filling rate is calculated using the shade of inner wall of oil tank and oil tank outer wall；
(8) oil tank amount of storage is calculated using tank(storage)capacity and oil tank filling rate.
2. a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images according to claim 1 determines method, it is characterised in that:
The oil tank is floating roof formula oil tank.
3. a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images according to claim 1 determines method, it is characterised in that:
The calculating of solar azimuth SoA is as follows:
Wherein, shade longitude and latitude of the sharp atural object vertex on ground is (Lon1, Lat1), and sharp atural object is passed through in the contact point on ground
Latitude is (Lon2, Lat2).
4. a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images according to claim 3 determines method, it is characterised in that:
The remote sensing images imaging date and solar elevation are determined according to solar azimuth, specifically:
(2.1) the solar hour angle T of imaging date atural object imaging is calculatedangle:
Tangle=(12-Tlocal) × 15,
Wherein, TlocalIt is satellite at the place that fixed latitude is passed by；
(2.2) solar declination δ is calculated:
δ=arcsin (0.39795cos (0.98563 (N-173)))；
Wherein, N is the day product that the date is imaged in image；
(2.3) the solar elevation SoH that it's the atural object moment pasts satellite, meter are calculated using solar hour angle, solar declination and atural object latitude
Calculation method uses following formula:
SoH=arcsin (sinLatsin δ+cosLatcos δ cosTangle)；
Wherein, sharp atural object geographic latitude and longitude is (Lon, Lat) in remote sensing images；
(2.4) formula for calculating solar azimuth using solar elevation, atural object latitude and solar declination is as follows:
(2.5) by the solar azimuth angle value being calculated in step (1) substitute into step (2.4) in, it is counter inquire into imaging the date and
Solar elevation SoH.
5. a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images according to claim 1 determines method, it is characterised in that:
Satellite aximuth SaA calculation formula is as follows:
Wherein, contact point longitude and latitude of the sharp atural object on ground is (Lon2, Lat2), and sharp atural object vertex is in remote sensing images
Projecting longitude and latitude is (Lon3, Lat3).
6. a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images according to claim 5 determines method, it is characterised in that:
According to oil tank, nearby sharply the shadow length of atural object, projected length and solar elevation calculate elevation of satellite, specifically:
(4.1) shadow length calculation formula is as follows:
(4.2) projected length calculation formula is as follows:
(4.3) elevation of satellite calculation formula is as follows:
7. a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images according to claim 6 determines method, it is characterised in that:
Selection needs to calculate the oil tank of oil reserve, using oil tank top shadow length, projected length, solar elevation, solar azimuth, defends
Elevation angle and satellite aximuth calculate tank elevation, specifically:
Tank elevation H calculation formula is as follows:
Wherein, θ=SaA-SoA, BC are shadow spots and sharp atural object vertex of the sharp atural object vertex on ground in remote sensing images
The distance between subpoint.
8. a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images according to claim 7 determines method, it is characterised in that:
Tank(storage)capacity calculates as follows:
Wherein, D is the oil tank abhiseca diameter measured on remote sensing images.
9. a kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images according to claim 8 determines method, it is characterised in that:
Oil tank filling rate R calculates as follows:Wherein, SoutFor oil tank outer wall shade, SinFor inner wall of oil tank shade；
Oil tank amount of storage VrCalculation method is as follows: Vr=V × R.
10. it is a kind of according to claim 1~any one of 9 described in method determined based on the oil tank amount of storage of Characteristics of The Remote Sensing Images
The oil tank amount of storage of realization determines system, it is characterised in that:
Solar azimuth computing module: for being reference with the sharp atural object near oil tank, calculating sun side in remote sensing images
Parallactic angle；
Solar elevation computing module: the remote sensing images imaging date and solar elevation are determined according to solar azimuth；
Satellite aximuth computing module: in remote sensing images, it is reference with the sharp atural object near oil tank, calculates satellite position
Angle；
Elevation of satellite computing module: according to shadow length, projected length and the solar elevation meter of atural object sharp near oil tank
Calculate elevation of satellite；
Tank elevation computing module: selection needs to calculate the oil tank of oil reserve, utilizes oil tank top shadow length, projected length, too
Positive elevation angle, solar azimuth, elevation of satellite and satellite aximuth calculate tank elevation；
Oil tank amount of storage computing module: tank(storage)capacity is calculated according to tank elevation and oil tank radius；Utilize inner wall of oil tank and oil tank
The shade of outer wall calculates oil tank filling rate；Oil tank amount of storage is calculated using tank(storage)capacity and oil tank filling rate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910516536.1A CN110210453B (en) | 2019-06-14 | 2019-06-14 | Remote sensing image feature-based oil tank storage capacity determination method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910516536.1A CN110210453B (en) | 2019-06-14 | 2019-06-14 | Remote sensing image feature-based oil tank storage capacity determination method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110210453A true CN110210453A (en) | 2019-09-06 |
| CN110210453B CN110210453B (en) | 2021-06-29 |
Family
ID=67792804
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910516536.1A Active CN110210453B (en) | 2019-06-14 | 2019-06-14 | Remote sensing image feature-based oil tank storage capacity determination method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110210453B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111144742A (en) * | 2019-12-24 | 2020-05-12 | 成都星时代宇航科技有限公司 | Satellite control method and device |
| CN111462222A (en) * | 2020-04-03 | 2020-07-28 | 深圳前海微众银行股份有限公司 | Method, device, equipment and medium for determining reserve of object to be detected |
| CN111540007A (en) * | 2020-04-10 | 2020-08-14 | 中国资源卫星应用中心 | Method for estimating storage capacity of oil tank by utilizing SAR (synthetic aperture radar) image |
| CN112200858A (en) * | 2020-10-10 | 2021-01-08 | 长光卫星技术有限公司 | External floating roof oil tank reserve analysis method based on high-resolution optical remote sensing image |
| CN112581595A (en) * | 2020-12-02 | 2021-03-30 | 中国人民解放军战略支援部队航天工程大学 | Multi-view satellite image consistency analysis method |
| CN113781478A (en) * | 2021-11-09 | 2021-12-10 | 中科星睿科技(北京)有限公司 | Oil tank image detection method, oil tank image detection device, electronic equipment and computer readable medium |
| CN114708512A (en) * | 2022-02-28 | 2022-07-05 | 尔特数据科技(江苏)有限公司 | Method and system for estimating crude oil reserve of floating-roof oil tank through satellite image |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103791885A (en) * | 2014-02-18 | 2014-05-14 | 西北工业大学 | Single satellite image-oriented self-acquired triangular element height calculation method |
| US20140219514A1 (en) * | 2013-02-07 | 2014-08-07 | Digitalglobe, Inc. | Automated metric information network |
| US20140363084A1 (en) * | 2013-06-05 | 2014-12-11 | Digitalglobe, Inc. | Oil tank farm storage monitoring |
| CN104637073A (en) * | 2014-12-30 | 2015-05-20 | 华中科技大学 | Zonal underground structure detection method based on sun shade compensation |
| CN104833335A (en) * | 2015-04-27 | 2015-08-12 | 中国资源卫星应用中心 | Method for obtaining satellite solar angle and time on basis of satellite orbit characteristics |
| CN104833336A (en) * | 2015-04-27 | 2015-08-12 | 中国资源卫星应用中心 | Satellite side-sway angle obtaining method based on image characteristics |
| CN107688782A (en) * | 2017-08-23 | 2018-02-13 | 中国科学院软件研究所 | Oil tank detection and reserve analysis method based on high-resolution optical remote sensing image |
| CN108053416A (en) * | 2017-12-14 | 2018-05-18 | 北京市遥感信息研究所 | A kind of maximum oil reserve extraction system based on single width satellite image |
| CN108363958A (en) * | 2018-01-19 | 2018-08-03 | 中国公路工程咨询集团有限公司 | A kind of oil tank detection method based on high-resolution optical remote sensing image |
| CN108427919A (en) * | 2018-02-22 | 2018-08-21 | 北京航空航天大学 | A kind of unsupervised oil tank object detection method guiding conspicuousness model based on shape |
-
2019
- 2019-06-14 CN CN201910516536.1A patent/CN110210453B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140219514A1 (en) * | 2013-02-07 | 2014-08-07 | Digitalglobe, Inc. | Automated metric information network |
| US20140363084A1 (en) * | 2013-06-05 | 2014-12-11 | Digitalglobe, Inc. | Oil tank farm storage monitoring |
| CN103791885A (en) * | 2014-02-18 | 2014-05-14 | 西北工业大学 | Single satellite image-oriented self-acquired triangular element height calculation method |
| CN104637073A (en) * | 2014-12-30 | 2015-05-20 | 华中科技大学 | Zonal underground structure detection method based on sun shade compensation |
| CN104833335A (en) * | 2015-04-27 | 2015-08-12 | 中国资源卫星应用中心 | Method for obtaining satellite solar angle and time on basis of satellite orbit characteristics |
| CN104833336A (en) * | 2015-04-27 | 2015-08-12 | 中国资源卫星应用中心 | Satellite side-sway angle obtaining method based on image characteristics |
| CN107688782A (en) * | 2017-08-23 | 2018-02-13 | 中国科学院软件研究所 | Oil tank detection and reserve analysis method based on high-resolution optical remote sensing image |
| CN108053416A (en) * | 2017-12-14 | 2018-05-18 | 北京市遥感信息研究所 | A kind of maximum oil reserve extraction system based on single width satellite image |
| CN108363958A (en) * | 2018-01-19 | 2018-08-03 | 中国公路工程咨询集团有限公司 | A kind of oil tank detection method based on high-resolution optical remote sensing image |
| CN108427919A (en) * | 2018-02-22 | 2018-08-21 | 北京航空航天大学 | A kind of unsupervised oil tank object detection method guiding conspicuousness model based on shape |
Non-Patent Citations (3)
| Title |
|---|
| TONG WANG 等: "Estimating the Volume of Oil Tanks Based on High-Resolution Remote Sensing Images", 《REMOTE SENSING》 * |
| 于圣涛: "基于高分辨率遥感影像的油罐体积求取研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 * |
| 李轩等: "基于似圆阴影的光学遥感图像油罐检测", 《电子与信息学报》 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111144742A (en) * | 2019-12-24 | 2020-05-12 | 成都星时代宇航科技有限公司 | Satellite control method and device |
| US11316584B2 (en) | 2019-12-24 | 2022-04-26 | Chengdu star era Aerospace Technology Co., Ltd. | Satellite control method and apparatus |
| CN111462222A (en) * | 2020-04-03 | 2020-07-28 | 深圳前海微众银行股份有限公司 | Method, device, equipment and medium for determining reserve of object to be detected |
| WO2021196698A1 (en) * | 2020-04-03 | 2021-10-07 | 深圳前海微众银行股份有限公司 | Method, apparatus and device for determining reserve of object to be detected, and medium |
| CN111462222B (en) * | 2020-04-03 | 2024-05-14 | 深圳前海微众银行股份有限公司 | Method, device, equipment and medium for determining reserves of objects to be detected |
| CN111540007B (en) * | 2020-04-10 | 2023-04-28 | 中国资源卫星应用中心 | Method for estimating oil tank storage capacity by SAR image |
| CN111540007A (en) * | 2020-04-10 | 2020-08-14 | 中国资源卫星应用中心 | Method for estimating storage capacity of oil tank by utilizing SAR (synthetic aperture radar) image |
| CN112200858A (en) * | 2020-10-10 | 2021-01-08 | 长光卫星技术有限公司 | External floating roof oil tank reserve analysis method based on high-resolution optical remote sensing image |
| CN112581595A (en) * | 2020-12-02 | 2021-03-30 | 中国人民解放军战略支援部队航天工程大学 | Multi-view satellite image consistency analysis method |
| CN112581595B (en) * | 2020-12-02 | 2023-12-19 | 中国人民解放军战略支援部队航天工程大学 | Multi-view satellite image consistency analysis method |
| CN113781478A (en) * | 2021-11-09 | 2021-12-10 | 中科星睿科技(北京)有限公司 | Oil tank image detection method, oil tank image detection device, electronic equipment and computer readable medium |
| CN113781478B (en) * | 2021-11-09 | 2022-05-24 | 中科星睿科技(北京)有限公司 | Oil tank image detection method, oil tank image detection device, electronic equipment and computer readable medium |
| CN114708512A (en) * | 2022-02-28 | 2022-07-05 | 尔特数据科技(江苏)有限公司 | Method and system for estimating crude oil reserve of floating-roof oil tank through satellite image |
| CN114708512B (en) * | 2022-02-28 | 2024-04-26 | 尔特数据科技(江苏)有限公司 | Method and system for estimating crude oil reserves of floating roof oil tank through satellite images |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110210453B (en) | 2021-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110210453A (en) | A kind of oil tank amount of storage based on Characteristics of The Remote Sensing Images determines method and system | |
| CN107367716B (en) | A kind of high-precision satellite-borne SAR geometric calibration method | |
| Chulliat et al. | The US/UK world magnetic model for 2015-2020 | |
| McLean et al. | The US/UK world magnetic model for 2005-2010 | |
| Maus et al. | The US/UK world magnetic model for 2010-2015 | |
| CN101354253B (en) | Geomagnetic auxiliary navigation algorithm based on matching degree | |
| CN109459027B (en) | Navigation method based on polarization-geomagnetic vector tight combination | |
| CN110057354B (en) | Geomagnetic matching navigation method based on declination correction | |
| US20120290199A1 (en) | Apparatus, system and method for self orientation | |
| CN104833336B (en) | A kind of satellite side-sway angle acquisition methods based on characteristics of image | |
| CN110187375A (en) | A kind of method and device improving positioning accuracy based on SLAM positioning result | |
| CN110017812B (en) | Measuring device and method for vertical line deviation, gravitational acceleration and gravity gradient | |
| CN103017762A (en) | Fast acquisition positioning method for space target of ground-based photoelectric telescope | |
| Stanchev et al. | Determination of the Black Sea area and coastline length using GIS methods and Landsat 7 satellite images | |
| CN109118520A (en) | A kind of minery power grid shaft tower displacement monitoring method and system | |
| CN101324435A (en) | Method and apparatus for obtaining direction information | |
| CN115291215B (en) | Long-time-sequence two-dimensional deformation rapid resolving method based on lifting orbit SAR satellite | |
| CN103630109A (en) | Method for determining geocentric vector based on stellar refraction | |
| CN104851130B (en) | A kind of three-dimensional generation method of satellite remote-sensing image | |
| CN115164871A (en) | Two-step autonomous positioning method based on polarized light field time difference | |
| Šugar et al. | Comparison of the reference mark azimuth determination methods | |
| CN115540856B (en) | Star sensor autonomous positioning method and system based on starlight three-dimensional refraction reconstruction | |
| Patat | Horizon synthesis for archaeo‐astronomical purposes | |
| Mohammad et al. | Evaluating the Relationship Between GPS Output and Planimetry System | |
| Park et al. | The Data Alignment Method Between GPS and IMU Based on ICP for Indoor Positioning |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |