Method for detecting zonal underground target in mountain land based on ridge heat radiation correction

Patent metadata

Field	Value
Publication number	US-9625611-B2
Application number	US-201515114651-A
Country	US
Kind code	B2
Filing date	Feb 10, 2015
Priority date	Dec 30, 2014
Publication date	Apr 18, 2017
Grant date	Apr 18, 2017

Abstract

The present invention discloses a method for detecting, recognizing, and positioning a zonal underground target in a mountain land environment by detecting a ridge position in the mountain land environment and carrying out energy correction. The method belongs to the interdisciplinary field of pattern recognition, remote sensing technology and terrain analysis. The zonal underground target can cause energy abnormity when the heat field thereof is different from that of a mountain mass, and the heat island effect of the ridge can also cause the energy of the mountain mass to be abnormal. However, the energy abnormity caused by the heat island effect is essentially different from the energy abnormity caused by the zonal underground target in the aspect of mode. Therefore, the present invention aims to achieve an effect of reducing a false alarm rate of detecting and recognizing a zonal underground target in the mountain land environment by eliminating the influence of the heat body effect generated by the ridge in the terrain on the weak energy abnormity mode presented by the zonal underground target. The present invention comprises steps of acquiring digital elevation information of terrain, performing de-noising pretreatment on the digital elevation information, detecting a ridge line, correcting energy at the ridge position, and detecting the zonal underground target.

First claim

The invention claimed is:

1. A method for detecting a zonal underground target in mountain land based on ridge heat radiation correction, the method comprising:
receiving an infrared energy image of terrain, the infrared energy image including the zonal underground target in mountain land;
acquiring digital elevation information of terrain, comprising:
determining longitudes and latitudes of coverage of the digital elevation information;
calculating a longitude-latitude array in the coverage of the digital elevation information according to the longitudes and latitudes; and
calculating the elevation of each point in the longitude-latitude array according to the longitude-latitude array, to obtain a digital elevation information array;

performing de-noising pretreatment on the digital elevation information array as obtained in the acquiring digital elevation information of terrain;
detecting a ridge line according to the digital elevation information array on which the de-noising treatment has been performed, comprising:
detecting a ridge line along a horizontal direction;
detecting a ridge line along a vertical direction; and
extracting a continuous ridge line according to the ridge line detected along the horizontal direction and the ridge line detected along the vertical direction;

determining an energy correction at a ridge position according to an energy distribution feature at the ridge position, and correcting the infrared energy image based on the energy correction as determined, the ridge position being determined based on the continuous ridge line as extracted; and
detecting the zonal underground target in the digital elevation information array, comprising:
setting parameters used for detecting, by means of traversing, the zonal underground target; and
traversing the digital elevation information array according to the parameters to detect the zonal underground target, and outputting the position of the detected zonal underground target.

2. The method of claim 1 , wherein the calculating the longitude-latitude array in the coverage of the digital elevation information according to the longitudes and latitudes includes:
calculating the width and height of a terrain range:
calculating, according to the longitudes and latitudes as determined of four vertices of a rectangle, the width and height of the rectangular coverage by using a distance measuring tool, the calculated width and height being marked as width and height in meters, and calculating the range of the longitudes and latitudes;

calculating the longitude-latitude array:
sampling at intervals of step meters according to a sampling interval of step meters, to respectively calculate the number of sampling points in the vertical direction and the number of sampling points in the horizontal direction, that is: height/step and width/step, and calculating a longitude step long_step and a latitude step lati_step between every two neighboring sampling points in the longitude-latitude array, thereby calculating the longitude and latitude of each sampling point in the longitude-latitude array.

3. The method of claim 1 , wherein the calculating the elevation of each point in the longitude-latitude array according to the longitude-latitude array, to obtain the digital elevation information array includes:
inputting the longitude-latitude array as calculated to obtain elevation data of each sampling point, and generating the digital elevation information array of the terrain according to the elevation data and outputting the digital elevation information array of the terrain.

4. The method of claim 1 , wherein the performing de-noising pretreatment on the digital elevation information array obtained in the acquiring digital elevation information of terrain includes:
performing de-noising pretreatment on the digital elevation information by using a mean filtering method.

5. The method of claim 1 , wherein the detecting the ridge line along the horizontal direction includes:
traversing along an x direction, to compare an elevation value of each sampling point with elevation values of sampling points in a preset range in a y direction, wherein if the sampling point has a maximum value in the y direction, the point is a candidate point of the ridge line along the horizontal direction.

6. The method of claim 1 , wherein the detecting the ridge line along the vertical direction includes:
traversing along a y direction, to compare an elevation value of each sampling point with elevation values of sampling points in a preset range in an x direction, wherein if the sampling point has a maximum value in the x direction, the point is a candidate point of the ridge line along the vertical direction.

7. The method of claim 1 , wherein the extracting the continuous ridge line according to the ridge line detected along the horizontal direction and the ridge line detected along the vertical direction includes:
determining the continuity of each ridge line candidate point as extracted, wherein a criterion for the determining is as follows: if there are different ridge line candidate points in a region defined by t in an x-direction and by t in a y direction from the ridge line candidate point, the ridge line candidate point is a final ridge line point; otherwise, it is judged that the isolate ridge line candidate point is a non-ridge line point, to finally obtain a continuous ridge line mark map, wherein t is a preset value.

8. The method of claim 1 , wherein the correcting energy at the ridge position: correcting energy at the ridge position according to an energy distribution feature at the ridge position includes:
correcting energy at the ridge position, and a specific correction method is: replacing energy of a sampling point on the ridge line with a mean value of energy of sampling points on two sides of the ridge line.

9. The method of claim 1 , wherein the parameters set in the setting parameters used for detecting, by means of traversing, the zonal underground target include:
size of a sampling section and spacing from a comparison section to a central section: suppose that a sampling section of a mountain land surface under which the zonal underground target is assumed to exist has a length of z pixels and a width of z pixels; sampling sections of mountain land surfaces on two sides of the zonal underground target also have a length of z pixels and a width of z pixels; and a distance from the center of either of the sampling sections of the mountain land surfaces on two sides to the center of the sampling section above the zonal underground target is s pixels, wherein z and s are preset values;
search direction: estimating the direction of the zonal underground target in this region, and setting, according to degree of the estimated direction, central head and tail coordinates P0 (x0, y0) and P1 (x1, y1) for traversing and searching, wherein a line connecting the two points passes through a position near the midpoint of this region;
distance for extension each time: each time after searching one group of images of the mountain land surface under which the zonal underground target is assumed to exist and mountain land surfaces on two sides of the zonal underground target, shifting a distance of d towards two sides to continue to search a new group, wherein searching is stopped automatically when image boundaries are reached; and
pulse threshold: only when absolute values of differences between an average gray value of each image sampling section, in the middle, of the mountain land surface under which the zonal underground target is assumed to exist and average gray values of image sampling sections of the mountain land surfaces on two sides of the zonal und

Assignees

• Univ Huazhong Science Tech

Inventors

• Zhang Tianxu
• Lu Cen
• Wang Yuehuan
• Yang Weidong
• Sang Nong
• Ma Wenxuan
• Hao Longwei

Classifications

• G06V20/194

  using hyperspectral data, i.e. more or other wavelengths than RGB · CPC title

• G06V20/13

  Satellite images · CPC title

• G01C21/005

  with correlation of navigation data from several sources, e.g. map or contour matching (G01C21/30 takes precedence) · CPC title

• G01V9/005 Primary

  by thermal methods, e.g. after generation of heat by chemical reactions · CPC title

• G06K9/00

  Physics · mapped topic