CN107945228B - Method for extracting height of oil tank based on single satellite image - Google Patents

Method for extracting height of oil tank based on single satellite image Download PDF

Info

Publication number
CN107945228B
CN107945228B CN201711341911.0A CN201711341911A CN107945228B CN 107945228 B CN107945228 B CN 107945228B CN 201711341911 A CN201711341911 A CN 201711341911A CN 107945228 B CN107945228 B CN 107945228B
Authority
CN
China
Prior art keywords
oil tank
arc
same name
image
satellite 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.)
Expired - Fee Related
Application number
CN201711341911.0A
Other languages
Chinese (zh)
Other versions
CN107945228A (en
Inventor
龙恩
孟钢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Remote Sensing and Digital Earth of CAS
Original Assignee
Institute of Remote Sensing and Digital Earth of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Remote Sensing and Digital Earth of CAS filed Critical Institute of Remote Sensing and Digital Earth of CAS
Priority to CN201711341911.0A priority Critical patent/CN107945228B/en
Publication of CN107945228A publication Critical patent/CN107945228A/en
Application granted granted Critical
Publication of CN107945228B publication Critical patent/CN107945228B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • 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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10004—Still image; Photographic image

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

The invention relates to a method for extracting the height of an oil tank based on a single satellite image, which comprises the following steps: s0, acquiring imaging parameters when the oil tank satellite images are shot; step S1, acquiring an oil tank satellite image, and obtaining a geometric fine correction image through geometric fine correction; step S2, identifying and positioning homonymous arc edges P 'and P' of the oil tank according to shadow boundary characteristics of the oil tank on the geometric precise correction satellite image; step S3, determining corresponding homonymous arc points T ' and T ' according to the homonymous arc edges P ' and P ' positioned in the step S2, and calculating to obtain homonymous arc distance T '; and S4, substituting the imaging parameters of the oil tank satellite images acquired in the S0 and the arc distances with the same name acquired in the S3 into an oil tank height measurement model to acquire the height of the oil tank. The method can accurately extract the height of the oil tank based on a single satellite image, effectively avoids calculation errors caused by artificial determination of the arc points with the same name, and greatly improves the extraction speed of the height of the oil tank.

Description

Method for extracting height of oil tank based on single satellite image
Technical Field
The invention relates to the technical field of remote sensing image processing and information extraction, in particular to a method for extracting the height of an oil tank based on a single satellite image.
Background
The oil tank is a core facility of a petroleum reserve base, belongs to the strategic target of national key construction, and the height information of the oil tank is important content and marks of oil depot identification, reservoir capacity estimation and quantitative analysis, and has important significance in the aspects of national economic construction, national defense safety capability and the like. However, in the conventional oil tank height extraction method based on satellite images, firstly, the laser ranging method and the three-dimensional image have the defects of high requirements on equipment and instruments and data sources, weak timeliness and complex calculation under the emergency task condition; secondly, the shadow length method is based on the arc edge characteristics of the shadow boundary of the oil tank, so that the shadow length is difficult to accurately measure, and the height of the oil tank cannot be accurately calculated.
Disclosure of Invention
In view of the above analysis, the present invention aims to provide a method for extracting the height of a tank based on a single satellite image, so as to solve some defects in the conventional technology for extracting the height of the tank based on the satellite image.
To achieve the above object, in one embodiment of the method according to the present invention, the method comprises the following steps:
s0, acquiring imaging parameters when the oil tank satellite images are shot;
step S1, acquiring an oil tank satellite image, and obtaining a geometric fine correction image through geometric fine correction;
step S2, identifying and positioning homonymous arc edges P ' and P ' of the oil tank according to the shadow boundary characteristics of the oil tank on the geometric precise correction satellite image, wherein the homonymous arc edges refer to an image P ' of the top edge of the oil tank on the remote sensing image and an image P of the shadow of the edge;
step S3, determining corresponding homonymous arc points T ' and T ' according to the homonymous arc edges P ' and P ' positioned in the step S2, and calculating to obtain homonymous arc distance T '; wherein the intersection point of the straight line passing through the centers of the two homonymous arc edges and having the slope in the solar direction and the two homonymous arc edges is a homonymous arc point, and the distance between the two homonymous arc points is a homonymous arc distance;
and S4, substituting the imaging parameters of the oil tank satellite images acquired in the S0 and the arc distances with the same name acquired in the S3 into an oil tank height measurement model to acquire the height of the oil tank.
in another embodiment of the method of the present invention, the imaging parameters in step S0 include a solar altitude α, a satellite altitude β, a solar azimuth γ, and a satellite azimuth θ.
In another embodiment of the method of the present invention, the step of geometrically fine correcting the satellite image in step S1 includes:
1) collecting ground control points: collecting more than six control points which have clear and obvious positioning marks and are distributed as uniformly as possible on the acquired oil tank satellite image;
2) selecting a geometric correction model: after the ground control points are collected, selecting a polynomial geometric correction model for establishing the relation between the image coordinates and the reference coordinates;
3) image resampling: and (4) carrying out re-interpolation calculation on the oil tank satellite image by adopting a resampling method to generate a geometric fine correction satellite image.
In another embodiment of the method of the present invention, the identification and positioning of the arc edges with the same name in step S2 adopts a human-computer interaction mode or an expert interpretation mode;
a man-machine interaction mode is adopted, specifically, shadow boundaries are automatically detected based on a computer, and then optimal positioning of homonymous arc edges is carried out according to the shadow boundaries of the oil tank in combination with manual work;
and adopting an expert mode, and directly carrying out interpretation positioning on the same-name arc edges based on expert interpretation knowledge.
In another embodiment of the method of the present invention, step S3 specifically includes the steps of:
s31, selecting key elements on the homonymous arc edges positioned in the S2;
step S32, determining corresponding homonymous arc points T 'and T' according to the key elements selected in the step S31;
and step S33, calculating to obtain the arc distance with the same name.
In another embodiment of the method of the present invention, step S31 specifically includes:
and (3) selecting three points on the homonymous arc edges P 'and P' of the positioned oil tank respectively based on the geometric fine correction image as input key elements, wherein:
coordinates of any three points on the same-name arc edge P': p'1(x’1，y’1)、P’2(x’2，y’2)、P’3(x’3，y’3)；
Coordinates of any three points on the same-name arc edge P': p'1(x"1，y"1)、P”2(x"2，y"2)、P”3(x"3，y"3)。
In another embodiment of the method of the present invention, step S32 specifically includes:
the intersection point coordinate model of the arc and the straight line which passes through the circle center of the arc and has the slope of the sun direction is calculated as follows:
Figure BDA0001508548810000031
wherein: the equation of a circle with the center coordinates M (A, B) and the radius R is (x-A)2+(y-B)2＝R2(ii) a A linear equation which passes through the circle center M and has the slope K is y which is Kx + B-KA;
Figure BDA0001508548810000032
Figure BDA0001508548810000033
Figure BDA0001508548810000034
K＝ctgγ
the coordinates of the arc points T '(x', y ') and T "(x", y ") of the same name can be obtained by substituting the coordinates of the center of the circle, the radius, and the slope ctg γ of the straight line of P' and P", respectively, into the intersection coordinate model.
In another embodiment of the method of the present invention, the distance between the homonymous arcs in step S33 is the distance between the homonymous arc points T ' (x ', y ') and T "(x", y "), which is expressed as:
Figure BDA0001508548810000041
in another embodiment of the method of the present invention, step S4 specifically includes:
step S41, constructing an oil tank height measurement model;
and S42, obtaining the height H of the oil tank by using the oil tank height measurement model constructed in the step S41 according to the oil tank satellite image imaging parameters acquired in the step S0 and the arc distance with the same name extracted in the step S33.
In another embodiment of the method of the present invention, the tank height measurement model of step S41 is
Figure BDA0001508548810000042
The invention has the characteristics and beneficial effects that:
the method can accurately extract the height of the oil tank based on a single satellite image, does not need stereopair or other special equipment, relaxes the requirements on a data source and an instrument, effectively avoids the calculation error caused by manually determining the homonymous arc point, greatly improves the extraction speed of the height of the oil tank, has stronger timeliness, and can be widely applied to the fields of oil depot identification, oil storage amount evaluation, emergency response and the like.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The drawings are only for purposes of illustrating particular embodiments and are not to be construed as limiting the invention, wherein like reference numerals are used to designate like parts throughout.
FIG. 1 is a flow chart of tank height extraction;
FIG. 2 is a schematic diagram of homonymous arc edges, homonymous arc points, and homonymous arc distances on a satellite image;
FIG. 3 is a diagram of the same arc pitch
FIG. 4 is a schematic diagram of the principle of construction of a model relating to the lateral height of an oil tank
Detailed Description
The preferred embodiments of the present invention will now be described in detail with reference to the accompanying drawings, which form a part hereof, and which together with the embodiments of the invention serve to explain the principles of the invention.
According to an embodiment of the invention, a method for extracting the height of an oil tank based on a single satellite image is disclosed, which comprises the following steps:
s0, acquiring imaging parameters when the oil tank satellite images are shot;
the optical satellite images are generally shot under relatively good illumination conditions during imaging, and each image corresponds to sun and satellite position information at the imaging moment, namely instantaneous imaging parameters. However, in the traditional model for extracting the ground feature height based on the shadow length method, only a single imaging parameter of the sun altitude angle is utilized, and imaging parameters such as the sun azimuth, the satellite height, the satellite azimuth and the like are ignored, so that the rigidness of the height measurement model is lacked, the error is large, and particularly with the common use of sub-meter high-resolution satellite images, the influence of the imaging parameters on the precision is also very obvious, and the imaging parameters need to be considered comprehensively.
The satellite image imaging parameters are generally recorded in a header file or metadata information or auxiliary data corresponding to the image product, and can be directly acquired. In one embodiment of the invention, the metadata is obtained by directly reading from the metadata file.
reading imaging parameters during satellite image shooting from metadata or header files or auxiliary data corresponding to the satellite images, wherein the imaging parameters comprise a solar altitude α, a satellite altitude angle β, a solar azimuth angle gamma and a satellite azimuth angle theta;
step S1, acquiring an oil tank satellite image, and obtaining a geometric fine correction image through geometric fine correction;
in the imaging process of the satellite image, due to images of various factors such as the sensor, the earth curvature, the topographic relief, the earth rotation, the orbit attitude, the projection mode and the like, the satellite image generally has various geometric distortions, so that the geometric shapes and the positions of the ground objects on the image are distorted, and the geometric shapes and the positions of the ground objects on the image are represented by displacement, rotation, scaling, affine, bending, higher-order distortion and the like. To eliminate these distortions and improve accuracy, the geometric distortion of the satellite image needs to be corrected to generate a new image with a certain map projection.
The step of geometrically fine correcting the satellite image comprises the following steps:
1) collecting ground control points: on the acquired oil tank satellite images, a plurality of control points which have clear and obvious positioning marks and are distributed as uniformly as possible are acquired, and the number of the control points is preferably 25.
2) Selecting a geometric correction model: and after the ground control points are acquired, selecting a polynomial geometric correction model for establishing the relation between the image coordinates and the reference coordinates.
3) Image resampling: and (3) carrying out re-interpolation calculation on the oil tank satellite image by adopting a resampling method such as bilinear interpolation, and generating a geometric precise correction satellite image which is high in positioning precision, has geographic information and can be calculated quantitatively.
Step S2, identifying and positioning homonymous arc edges P 'and P' of the oil tank according to shadow boundary characteristics of the oil tank on the geometric precise correction satellite image;
the method mainly extracts the height of the oil tank by means of the idea of homonymous characteristics, wherein the homonymous characteristics refer to homonymous arc edges of the oil tank on a single satellite image, as shown in figures 2 and 3, an oil tank top edge arc P 'and an oil tank shadow outer edge arc P' are arranged on the image along the shadow direction, the P 'and the P' are homonymous arc edges, the physical meaning of the P 'and the P' is the image of the oil tank top edge, the identification and positioning of the homonymous arc edges of the oil tank are mainly carried out according to the shadow boundary characteristics of oil tank facilities on a geometric correction image, and the identification and positioning of the homonymous arc edges are carried out;
the identification and positioning of the arc edges with the same name can adopt a man-machine interaction mode or an expert interpretation mode. In one embodiment of the invention, a man-machine interaction mode is adopted, firstly, a shadow boundary is automatically detected based on a computer, and then optimal positioning of homonymous arc edges is carried out according to the shadow boundary of the oil tank in combination with manual work; in another embodiment of the invention, an expert mode is adopted, and interpretation positioning of the arc edges with the same name is directly carried out based on expert interpretation knowledge.
Step S3, determining corresponding homonymous arc points T ' and T ' according to the homonymous arc edges P ' and P ' positioned in the step S2, and calculating to obtain homonymous arc distance T ';
the homonymous arc distance is an important parameter for extracting the height of the oil tank, and the definition of the homonymous arc distance refers to the distance between a circle where two homonymous arc sides are located and an intersection point of a straight line which passes through the center of the circle and has the slope of the sun direction in the shadow direction. The intersection point of a straight line passing through the center of the circle and having the slope of the sun direction and two homonymous arc edges is a homonymous arc point, and the important link in extracting the homonymous arc distance is the selection of key elements and the determination of the homonymous arc point.
Step S3 specifically includes the steps of:
s31, selecting key elements on the homonymous arc edges positioned in the S2;
and (3) selecting three points on the homonymous arc edges P 'and P' of the positioned oil tank respectively based on the geometric fine correction image as input key elements, wherein:
coordinates of any three points on the same-name arc edge P': p'1(x’1，y’1)、P’2(x’2，y’2)、P’3(x’3，y’3)；
Coordinates of any three points on the same-name arc edge P': p'1(x"1，y"1)、P”2(x"2，y"2)、P”3(x"3，y"3)；
Step S32, determining corresponding homonymous arc points T 'and T' according to the key elements selected in the step S31;
known as P'1、P’2、P’3Three points on the same arc edge P', P "1、P”2、P”3Three points on the same-name arc edge P' are provided, the coordinate values can be directly read based on the geometric fine correction image, and then the following steps exist:
the equation of the circle on which the arc sides P 'and P' lie can be expressed as:
(x-A)2+(y-B)2＝R2the circle center coordinates M (A, B) and the radius are R;
wherein, the circle center coordinate M ' (A ', B ') of the circle where the arc edge P ' is located has a radius of R '; the arc edge P ' is located at the center coordinates M ' (A ', B ') of the circle, and the radius is R '.
The equation of a straight line passing through the circle center M and having the slope K is as follows:
y＝Kx+B-KA；
wherein the slopes of the straight lines passing through the circle center M 'and the circle center M' are both ctg gamma.
The coordinate expressions of the arc points of the same name T ' (x ', y ') and T "(x", y ") are (i.e. the intersection of the straight line and the circle in the shadow direction):
Figure BDA0001508548810000081
wherein:
Figure BDA0001508548810000082
Figure BDA0001508548810000083
Figure BDA0001508548810000084
k' ═ ctg γ (γ sun azimuth)
Figure BDA0001508548810000085
Figure BDA0001508548810000086
Figure BDA0001508548810000087
K ═ ctg γ (γ sun azimuth)
The coordinates of the arc points T '(x', y ') and T' (x ', y') with the same name can be obtained after the substitution respectively.
Step S33, calculating to obtain the arc distance of the same name;
the homonymous arc distance is the distance between homonymous arc points T '(x', y ') and T' (x ', y'), and the expression is as follows:
Figure BDA0001508548810000091
step S4, substituting the oil tank satellite image imaging parameters obtained in the step S0 and the arc distances with the same name obtained in the step S33 into an oil tank height measurement model to obtain the height of the oil tank;
the oil tank height measurement model is important content of the invention, and the oil tank height extraction can be carried out by substituting the oil tank homonymous arc distance and imaging parameters based on the height measurement model.
Step S4 specifically includes the steps of:
step S41, constructing an oil tank height measurement model;
as shown in fig. 4, when TT "is a sun ray, TT 'is a satellite shooting sight line, the physical meanings of arc points T' and T" of the same name are image points of the edge T point at the top of the oil tank on the satellite image and shadow image points thereof, the height H of the oil tank is TO. the sun altitude angle α, the satellite altitude angle β, and the included angle between the sun azimuth and the satellite azimuth space is | γ - θ |, then
As can be seen from the spatial geometry, the,
T"O＝H*ctgα
T’O＝H*ctgβ
∠T’OT"＝|γ-θ|
by trigonometric relationship T'2＝T’O2+T"O2-2 ∠ T 'O ∠ T "O COS ∠ T' OT", knowing:
Figure BDA0001508548810000092
namely the height measuring model of the cylindrical facility.
And S42, obtaining the height H of the oil tank by using the oil tank height measurement model constructed in the step S41 according to the oil tank satellite image imaging parameters acquired in the step S0 and the arc distance with the same name extracted in the step S33.
The invention has the characteristics and beneficial effects that:
the method can accurately extract the height of the oil tank based on a single satellite image, does not need stereopair or other special equipment, relaxes the requirements on a data source and an instrument, effectively avoids the calculation error caused by manually determining the homonymous arc point, greatly improves the extraction speed of the height of the oil tank, has stronger timeliness, and can be widely applied to the fields of oil depot identification, oil storage amount evaluation, emergency response and the like.
Although illustrative embodiments of the present invention have been described above to facilitate the understanding of the present invention by those skilled in the art, it should be understood that the present invention is not limited to the scope of the embodiments, and various changes may be made apparent to those skilled in the art as long as they are within the spirit and scope of the present invention as defined and defined by the appended claims, and all matters of the invention which utilize the inventive concepts are protected.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention.

Claims (8)

1.一种基于单幅卫星图像提取油罐高度的方法，其特征在于，包括1. a method for extracting oil tank height based on a single satellite image, is characterized in that, comprising 步骤S0、获取油罐卫星图像拍摄时的成像参数；Step S0, acquiring imaging parameters when the oil tank satellite image is photographed; 步骤S1、获取油罐卫星图像，并经几何精校正得到几何精校正卫星图像；Step S1, obtaining the satellite image of the oil tank, and obtaining the satellite image of the geometric fine correction through the geometric fine correction; 步骤S2、根据油罐在几何精校正卫星图像上的阴影边界特征，识别和定位油罐同名弧边P’和P"，其中同名弧边是指图像上油罐顶部边缘的像P’和该边缘阴影的像P"；Step S2, according to the shadow boundary feature of the oil tank on the geometric fine-correction satellite image, identify and locate the arc edges P' and P" of the same name of the oil tank, wherein the arc edges of the same name refer to the image P' and the same name of the top edge of the oil tank on the image. Edge shadow like P"; 步骤S3、根据步骤S2中定位的油罐同名弧边P’和P"，确定相对应的同名弧点T’和T"，计算得到同名弧距T’T"；其中分别过两个同名弧边的圆心且斜率为太阳方位的直线与两个同名弧边的交点为同名弧点，两个同名弧点之间的距离为同名弧距；Step S3, according to the same name arc edge P' and P" of the oil tank positioned in step S2, determine the corresponding same name arc point T' and T", and calculate the same name arc distance T'T"; The intersection of the line whose center is the center of the side and whose slope is the azimuth of the sun and the two arc sides of the same name is the arc point of the same name, and the distance between the two arc points of the same name is the arc distance of the same name; 步骤S4、将步骤S0中获取的油罐卫星图像成像参数及步骤S3得到的同名弧距代入油罐测高模型，得到油罐高度；Step S4, substituting the oil tank satellite image imaging parameters obtained in step S0 and the arc distance of the same name obtained in step S3 into the oil tank altimetry model to obtain the oil tank height; 所述步骤S3具体包括：The step S3 specifically includes: 步骤S31、在步骤S2中定位的同名弧边上选取关键元素；Step S31, select key elements on the arc edge of the same name positioned in step S2; 步骤S32、根据步骤S31中选取的关键元素确定相对应的同名弧点T’和T"；Step S32, determine corresponding arc points T' and T" of the same name according to the key elements selected in step S31; 步骤S33、计算得到同名弧距；Step S33, calculate the arc distance of the same name; 所述步骤S32具体包括：The step S32 specifically includes: 计算圆弧和过圆弧圆心且斜率为太阳方位的直线的交点坐标模型为：The coordinate model for calculating the intersection point of the arc and the straight line passing through the center of the arc and the slope is the azimuth of the sun is:
Figure FDA0002374009120000011
Figure FDA0002374009120000011
 其中：圆心坐标M(A，B)、半径为R的圆的方程为(x-A)2+(y-B)2＝R2；过圆心M且斜率为K的直线方程为y＝Kx+B-KA；Wherein: the circle center coordinates M(A, B) and the equation of the circle with radius R is (xA) 2 +(yB) 2 =R 2 ; the equation of the straight line passing through the circle center M and the slope is K is y=Kx+B-KA ; 分别将P’和P"的圆心坐标、半径以及直线斜率ctgγ代入上述交点坐标模型，得到同名弧点T’(x'，y')和T"(x"，y")的坐标。Substitute the center coordinates, radius and straight line slope ctgγ of P' and P" into the above-mentioned intersection coordinate model to obtain the coordinates of the arc points T'(x', y') and T"(x", y") with the same name. 2.如权利要求1所述的方法，其特征在于，所述步骤S0中成像参数包括：太阳高度角α、卫星高度角β、太阳方位角γ和卫星方位角θ。2 . The method according to claim 1 , wherein the imaging parameters in the step S0 include: solar altitude α, satellite altitude β, solar azimuth γ and satellite azimuth θ. 3 . 3.如权利要求2所述的方法，其特征在于，所述步骤S1中卫星图像几何精校正的步骤包括：3. The method according to claim 2, wherein the step of fine-correcting the satellite image geometry in the step S1 comprises: 1)地面控制点采集：在获取的油罐卫星图像上，采集六个以上、具有清晰明显定位标志并尽可能均匀分布的控制点；1) Ground control point collection: On the acquired satellite image of the oil tank, collect more than six control points with clear and obvious positioning marks and as evenly distributed as possible; 2)几何校正模型选择：采集完地面控制点之后，选择多项式几何校正模型，用于建立图像坐标与参考坐标的关系；2) Geometric correction model selection: After collecting ground control points, select a polynomial geometric correction model to establish the relationship between image coordinates and reference coordinates; 3)图像重采样：采用重采样方法，对油罐卫星图像进行重新插值计算，生成几何精校正卫星图像。3) Image resampling: The resampling method is used to re-interpolate the satellite image of the oil tank to generate a geometrically fine-corrected satellite image. 4.如权利要求3所述的方法，其特征在于，所述步骤S2中同名弧边的识别与定位采取人机交互模式或专家判读模式；4. method as claimed in claim 3, is characterized in that, in described step S2, the identification and positioning of the arc edge of the same name adopts human-computer interaction mode or expert interpretation mode; 采用人机交互模式，具体为首先基于计算机自动检测阴影边界，然后根据油罐阴影边界结合人工进行同名弧边的优选定位；The human-computer interaction mode is adopted, specifically, the shadow boundary is automatically detected based on the computer, and then the optimal positioning of the arc edge of the same name is carried out manually according to the shadow boundary of the oil tank; 采用专家模式，基于专家判读知识，直接进行同名弧边的判读定位。Using the expert mode, based on the expert interpretation knowledge, the interpretation and positioning of the arc edge with the same name is directly carried out. 5.如权利要求1所述的方法，其特征在于，所述步骤S31具体包括：5. The method of claim 1, wherein the step S31 specifically comprises: 基于几何精校正图像，分别在定位的所述油罐同名弧边P’和P"上任选三点，作为输入关键元素，其中：Based on the geometrically fine-corrected image, three points are selected respectively on the positioned arc edges P' and P" of the same name of the oil tank as input key elements, wherein: 同名弧边P’上的任三点坐标：P’1(x’1，y’1)、P’2(x’2，y’2、P’3(x’3，y’3)；Coordinates of any three points on the arc edge P' of the same name: P' 1 (x' 1 , y' 1 ), P' 2 (x' 2 , y' 2 , P' 3 (x' 3 , y' 3 ); 同名弧边P"上的任三点坐标：P”1(x"1，y"1)、P”2(x"2，y"2)、P”3(x"3，y"3)。Coordinates of any three points on the arc edge P" with the same name: P" 1 (x" 1 , y" 1 ), P" 2 (x" 2 , y" 2 ), P" 3 (x" 3 , y" 3 ) . 6.如权利要求5所述的方法，其特征在于，所述步骤S33中同名弧距为同名弧点T’(x'，y')和T”(x"，y")之间的距离，其表达式为：6. method as claimed in claim 5 is characterized in that, in described step S33, the arc distance of the same name is the distance between the arc points T'(x', y') of the same name and T" (x", y") , whose expression is:
Figure FDA0002374009120000031
Figure FDA0002374009120000031
 7.如权利要求6所述的方法，其特征在于，所述步骤S4具体包括：7. The method of claim 6, wherein the step S4 specifically comprises: 步骤S41、构建油罐测高模型；Step S41, constructing an altimetry model of an oil tank; 步骤S42、将步骤S0中获取的油罐卫星图像成像参数及步骤S33提取的同名弧距，利用步骤S41中构建的油罐测高模型，得到油罐高度H。In step S42, the oil tank height H is obtained by using the oil tank height measurement model constructed in step S41 using the imaging parameters of the oil tank satellite image obtained in step S0 and the arc distance of the same name extracted in step S33. 8.如权利要求7所述的方法，其特征在于，所述步骤S41的油罐测高模型为8. The method of claim 7, wherein the oil tank altimetry model of the step S41 is
Figure FDA0002374009120000032
Figure FDA0002374009120000032
CN201711341911.0A 2017-12-14 2017-12-14 Method for extracting height of oil tank based on single satellite image Expired - Fee Related CN107945228B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711341911.0A CN107945228B (en) 2017-12-14 2017-12-14 Method for extracting height of oil tank based on single satellite image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711341911.0A CN107945228B (en) 2017-12-14 2017-12-14 Method for extracting height of oil tank based on single satellite image

Publications (2)

Publication Number Publication Date
CN107945228A CN107945228A (en) 2018-04-20
CN107945228B true CN107945228B (en) 2020-06-02

Family

ID=61944322

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711341911.0A Expired - Fee Related CN107945228B (en) 2017-12-14 2017-12-14 Method for extracting height of oil tank based on single satellite image

Country Status (1)

Country Link
CN (1) CN107945228B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110276796B (en) * 2019-06-25 2020-01-03 北京四象爱数科技有限公司 Floating roof oil tank reserve estimation method based on high-resolution SAR image
CN114036971B (en) * 2021-11-18 2022-05-17 中科星睿科技(北京)有限公司 Oil tank information generation method, oil tank information generation device, electronic device, and computer-readable medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101566692A (en) * 2009-05-26 2009-10-28 吉林大学 Method for detecting cloud height by utilizing cloud shadow information in satellite remote sensing data
CN104820984A (en) * 2015-04-24 2015-08-05 国家测绘地理信息局卫星测绘应用中心 Satellite remote sensing stereo image processing system and method
CN105758806A (en) * 2016-02-04 2016-07-13 中国农业科学院农业资源与农业区划研究所 Spectral characteristic based remote sensing monitoring method of plastic film mulched farmland

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10055885B2 (en) * 2015-09-16 2018-08-21 Raytheon Company Systems and methods for digital elevation map filters for three dimensional point clouds

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101566692A (en) * 2009-05-26 2009-10-28 吉林大学 Method for detecting cloud height by utilizing cloud shadow information in satellite remote sensing data
CN104820984A (en) * 2015-04-24 2015-08-05 国家测绘地理信息局卫星测绘应用中心 Satellite remote sensing stereo image processing system and method
CN105758806A (en) * 2016-02-04 2016-07-13 中国农业科学院农业资源与农业区划研究所 Spectral characteristic based remote sensing monitoring method of plastic film mulched farmland

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"TIN图像配准中的同名特征点的可靠性提高";崔晨耕;《信息技术》;20160921;第45卷(第8期);期刊第2节 *
"基于单张遥感影像的城市建筑物高度提取研究";王京卫等;《测绘通报》;20120425(第4期);期刊第2节 *
王京卫等."基于单张遥感影像的城市建筑物高度提取研究".《测绘通报》.2012,(第4期), *

Also Published As

Publication number Publication date
CN107945228A (en) 2018-04-20

Similar Documents

Publication Publication Date Title
CN108053416B (en) Maximum oil storage capacity extraction system based on single satellite image
US5309522A (en) Stereoscopic determination of terrain elevation
CN109544628B (en) Accurate reading identification system and method of a pointer meter
CN111368906B (en) A reading recognition method of pointer oil level gauge based on deep learning
CN110245654B (en) Pointer reading identification method and device, computer device and readable storage medium
CN105046251B (en) A kind of automatic ortho-rectification method based on environment No.1 satellite remote-sensing image
CN104657711B (en) A kind of readings of pointer type meters automatic identifying method of robust
CN102521560A (en) Instrument pointer image identification method of high-robustness rod
CN102855628B (en) Automatic matching method for multisource multi-temporal high-resolution satellite remote sensing image
CN110503144A (en) A kind of pointer instrument recognition methods for crusing robot
CN103177444A (en) Automatic SAR (synthetic-aperture radar) image rectification method
CN110906830B (en) Intelligent distance measuring wheel system and distance measuring method thereof
CN111738229B (en) Automatic reading method for scale of pointer dial
CN110852954B (en) Image inclination correction method and system for pointer instrument
CN110223355A (en) A kind of feature mark poiX matching process based on dual epipolar-line constraint
CN118999479B (en) Hydraulic engineering based on Beidou and equipment inclination settlement monitoring device and method thereof
CN107945228B (en) Method for extracting height of oil tank based on single satellite image
CN114399677A (en) Pointer instrument identification method based on text region reading
CN109188483B (en) Time-sequential high-precision automatic calibration method for exterior orientation elements
CN112464986B (en) Reading method and device for pointer type disc instrument
CN105571598A (en) Satellite laser altimeter footprint camera pose measuring method
CN114255051B (en) Authenticity inspection method of orthometric product based on stereo mapping satellite
CN116989825A (en) A roadside lidar-camera-UTM coordinate system joint calibration method and system
CN105825517B (en) A kind of image correcting method and system of navigation height error
CN103557849B (en) Be applied to the single wood in mountain area absolute fix and the coordinates compensation method fast of fast bird image

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200602