[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN105574855A - Method for detecting infrared small targets under cloud background based on temperate filtering and false alarm rejection - Google Patents

Method for detecting infrared small targets under cloud background based on temperate filtering and false alarm rejection Download PDF

Info

Publication number
CN105574855A
CN105574855A CN201510917030.3A CN201510917030A CN105574855A CN 105574855 A CN105574855 A CN 105574855A CN 201510917030 A CN201510917030 A CN 201510917030A CN 105574855 A CN105574855 A CN 105574855A
Authority
CN
China
Prior art keywords
template
cloud
max
pixel
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.)
Granted
Application number
CN201510917030.3A
Other languages
Chinese (zh)
Other versions
CN105574855B (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201510917030.3A priority Critical patent/CN105574855B/en
Publication of CN105574855A publication Critical patent/CN105574855A/en
Application granted granted Critical
Publication of CN105574855B publication Critical patent/CN105574855B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/10Image enhancement or restoration using non-spatial domain filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10048Infrared image

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

The invention diWscloses a method for detecting infrared small targets under a cloud background based on temperate filtering and false alarm rejection. The method comprises the following steps: firstly, obvious noises are removed from an image through maximum median filtering to complete image preprocessing; secondly, Robinson template filtering is used to suppress the background and highlight targets; thirdly, cloud partitioning of the original image is carried out, binarization processing of a result obtained after the Robinson template filtering is carried out through a low threshold in a cloud partition, and binarization processing of the result obtained after the Robinson template filtering is carried out through a high threshold in a non-cloud partition; finally, a plurality of ''false target points'' generated by a same target are further removed from the binarization results to complete ''coarse detection'', and time domain processing is continued for adjacent frame images that have gone through space domain processing to complete ''fine detection'', so that the detection of the infrared small targets is achieved. The method provided by the invention has the advantage that the constant false alarm rejection algorithm is added in interframe track correlation, so that the detection false alarm rate is greatly reduced.

Description

Based on the infrared small target detection method of template convolution and false alarm rejection under cloud background
Technical field
The invention belongs to infrared image detection and processing technology field, particularly under a kind of cloud background based on the infrared small target detection method of template convolution and false alarm rejection.
Background technology
Infrared small target detection is a core technology in the systems such as Infra-Red Search & Track System, Large visual angle object detection system, satellite remote sensing, disaster alarm, fire control and disaster rescue.Be subject to the factor impacts such as air, radiation from sea surface, operating distance and noise of detector due to infrared sensor, make remote target size on infrared image less, even present point-like; In addition, the signal to noise ratio (S/N ratio) of image is lower, adds background more complicated under normal circumstances, target to be easy to by noise and background clutter flood, make the detection of infrared small target become more difficult.
The detection probability and the reduction false alarm rate that how to improve target are mainly devoted to the research of infrared small target detection method.Traditional detection method is divided into two large classes: Time Domain Processing and spatial processing.Time Domain Processing, based on the multiple image sequence do not gathered in the same time, fully takes into account the correlativity between picture frame, can suppress between consecutive frame to false-alarm, keeps higher accuracy of detection.Conventional time-domain processing method has entropy difference method (Wang Guangjun, field inscription on ancient bronze objects, Liu Jian. based on the infrared image small target deteection [J] of local entropy. infrared and laser engineering, 2000,04:26-29.), sequence image detection method (Sun Jigang. sequence image infrared small target detection and track algorithm are studied [D]. Postgraduate School, Chinese Academy of Sciences (Changchun optical precision optical machinery and physics Institute), 2014.) etc.But must carry out registration operation between the method consecutive frame, intractability is comparatively large, and algorithm is comparatively complicated, and is not suitable for the detection of single-frame images.Spatial processing then utilizes object pixel to there is significant difference with neighborhood territory pixel in gray scale, and there is not the feature of correlativity with background, directly in its spatial domain, carry out Pixel-level process to single-frame images, owing to usually utilizing template operation, therefore algorithm is easy to be transplanted in hardware.Conventional method has morphology filter method (to spend the profit autumn, Zhang Ying, Lin Xiaochun. based on the infrared small target detection method [J] of shape filtering. laser and infrared, 2005,06:451-453.), high-pass filtering method (Dong Hongyan, Li Jicheng, Shen Zhenkang etc. based on the small target deteection [J] of high-pass filtering and Order Filtering. systems engineering and electronic technology, 2004,26 (5) .) etc.But the false alarm rate of the method is higher, and robustness is poor, and cannot be applied in the process of video sequence.
Summary of the invention
The invention provides the infrared small target detection method based on template convolution and false alarm rejection under a kind of cloud background, the infrared small target under cloud background accurately can be detected.
The present invention is the technical scheme solving prior art problem: based on the infrared small target detection method of template convolution and false alarm rejection under a kind of cloud background, for single-frame images, adopt the spatial domain operation of template convolution, first carry out max-medium filter to image to remove significant noise and complete Image semantic classification, secondly with Robinson's template convolution Background suppression, outstanding target, then cloud sector division is carried out to former figure, Low threshold is adopted to carry out binary conversion treatment in cloud sector part to the filtered result of Robinson, but not cloud sector part then adopts high threshold process, finally the result after binaryzation is rejected further multiple " the pseudo-impact points " of the generation of same target, thus complete " rough detection ", the consecutive frame image carrying out spatial processing is continued to take time domain operation process, namely track association is carried out in interframe, and the difference in gamma characteristic and kinetic characteristic carries out CFAR suppression operation for real goal and false-alarm point, complete " essence detects ", thus realize the detection of infrared small target.
Compared with prior art, its remarkable advantage is in the present invention: spatial processing and Time Domain Processing combine by (1), have concurrently simultaneously spatial processing algorithm complex low, be easy to Hardware and realize and Time Domain Processing accuracy of detection is high, false alarm rate is low feature.(2) preprocessing part, the basis of medium filtering proposes a kind of max-medium filter method, thus considers the gamma characteristic in image all directions.(3) target rough detection part: introduce Robinson's filter method Background suppression, outstanding target; Divide cloud sector and non-cloud sector, employing dual threshold carries out binary conversion treatment to different piece, ensures that Small object in cloud sector is not by blindness filtering.(4) target essence detecting portion, utilizes the difference in gamma characteristic and kinetic characteristic between real goal and false-alarm point, adds CFAR Restrainable algorithms, greatly reduce detection false alarm rate in interframe track association.
Below in conjunction with accompanying drawing, the present invention is described in further detail.
Accompanying drawing explanation
Fig. 1 is the process flow diagram based on the infrared small target detection method of template convolution and false alarm rejection under cloud background of the present invention.
Fig. 2 is the infrared image comprising Small object under the complicated cloud background of single frames.
Fig. 3 is the design sketch after max-medium filter process.
Fig. 4 is the design sketch after Robinson's filtering process.
Fig. 5 is the design sketch after cloud sector differentiates.
Fig. 6 is the design sketch after binaryzation operation.
Fig. 7 is the design sketch after adjacent point target combination.
Fig. 8 is the design sketch after interframe track association.
Fig. 9 is the design sketch after false alarm rejection.
Figure 10 is false alarm rate statistical graph.
Figure 11 is detectivity statistical graph.
Embodiment
Composition graphs 1, based on the infrared small target detection method of template convolution and false alarm rejection under cloud background of the present invention, for single-frame images, adopt the spatial domain operation of template convolution, first carry out max-medium filter to image to remove significant noise and complete Image semantic classification, secondly with Robinson's template convolution Background suppression, outstanding target, then cloud sector division is carried out to former figure, Low threshold is adopted to carry out binary conversion treatment in cloud sector part to the filtered result of Robinson, but not cloud sector part then adopts high threshold process, " neighborhood non-maxima suppression " principle is finally utilized to reject multiple " the pseudo-impact points " of same target generation further to the result after binaryzation, thus complete " rough detection ", the consecutive frame image carrying out spatial processing is continued to take time domain operation process, namely track association is carried out in interframe, and the difference in gamma characteristic and kinetic characteristic carries out CFAR suppression operation for real goal and false-alarm point, complete " essence detects ", thus realize the detection of infrared small target.Its concrete implementation step is as follows:
1. Image semantic classification.Max-medium filter is carried out to input picture as shown in Figure 2, remove the remarkable noise in image, i.e. max-medium filter water intaking flat, vertical, left 45 degree, 45 degree of four filtering directions, the right side, the maximal value of getting pixel grey scale intermediate value in all directions gives center pixel.Max-medium filter has fully taken into account the pixel grey scale distribution in multiple directions, and retains the energy of target in image better, and the basis not destroying the original gray feature of target realizes efficient denoising, and result as shown in Figure 3.
For the max-medium filter template of (2N+1) * (2N+1), its computing formula is:
f max-med(i,j)=max(z 1,z 2,z 3,z 4)(1)
In formula (1),
z 1=med[f(i,j-N),...,f(i,j),...,f(i,j+N)]
z 2=med[f(i-N,j),...,f(i,j),...,f(i+N,j)]
(2)
z 3=med[f(i+N,j-N),...,f(i,j),...,f(i-N,j+N)]
z 4=med[f(i-N,j-N),...,f(i,j),...,f(i+N,j+N)]
Wherein, pixel coordinate centered by (i, j), med is for getting median operation, and max is for getting maxima operation, and N characterizes template size size, and in the inventive method, for N=2, namely template size is 5*5.
2. rough detection, its process is:
(1) Robinson's template convolution is carried out to the image after max-medium filter, suppress cloud background, and outstanding target.This template utilizes this characteristic comparatively large of the gray difference between real goal and cloud background, the maximum gradation value of central pixel point gray-scale value and surrounding pixel is compared and arithmetical operation.If the gray scale of center pixel is comparatively strong, then can be retained, otherwise will be suppressed.Meanwhile, template-setup has isolation strip, can ensure that the gamma characteristic of real goal is not destroyed.Its result as shown in Figure 4.
The method of Robinson's template convolution is: for Robinson's Filtering Template of (2N+1) * (2N+1), its computing formula is:
f R o b i s o n ( i , j ) = f ( i , j ) - m a x ( z 1 , z 2 , z 3 , z 4 ) , f ( i , j ) > m a x ( z 1 , z 2 , z 3 , z 4 ) f ( i , j ) - m i n ( z 1 , z 2 , z 3 , z 4 ) , f ( i , j ) < m a x ( z 1 , z 2 , z 3 , z 4 ) 0 , e l s e - - - ( 3 )
In formula (3),
z 1=max[f(i-N:i+N,j-N)]
z 2=max[f(i-N:i+N,j+N)]
(4)
z 3=max[f(i-N,i-N:j+N)]
z 4=max[f(i+N,i-N:j+N)]
Wherein, pixel coordinate centered by (i, j), max is for getting maxima operation, and N characterizes in template size size the inventive method, and for N=4, namely template size is 9*9.
(2) cloud sector differentiation is carried out to original image.Because cloud sector usually presents wire or stratiform in infrared image, utilize this feature, can take the method for template matches filtering in cloud, cloud exterior domain distinguished.This template is provided with blank isolation strip on the one hand and is not damaged to protect the gamma characteristic of real goal; on the other hand by the gray scale sum of all for central row pixels and the pixel grey scale sum of the top a line or bottom a line are carried out doing difference operation and normalizing to [0; 255] interval, and utilize binaryzation to operate cloud sector segmentation.Particularly, this operation is divided into following 3 steps:
I () carries out matched filtering operation to original image.For the matched filtering template of (2M+1) * (2N+1) size, its computing formula is:
f c ( i , j ) = z 2 - z 3 , z 1 > z 3 z 2 - z 1 , e l s e - - - ( 5 )
In formula (5),
z 1=sum(f(i-N,j-M:j+M))
z 2=sum(f(i,j-M:j+M))(6)
z 3=sum(f(i+N,j-M:j+M))
Wherein, pixel coordinate centered by (i, j), sum is sum operation, and M, N characterize template size size, and in the present invention, M gets 5, N and gets 1, and namely template size is 11*3.
(ii) gray scale normalization process is carried out to the result after matched filtering.Normalization formula is:
f c &prime; ( i , j ) = 255 , f c ( i , j ) > 255 0 , f c ( i , j ) < 0 f c ( i , j ) , e l s e - - - ( 7 )
(iii) binarization segmentation is carried out to the result after normalization.The formula of binaryzation operation is:
f c l o u d ( i , j ) = 255 , f c &prime; ( i , j ) > f t h r e s h o l d 0 , e l s e - - - ( 8 )
F in formula (8) thresholdfor binary-state threshold.Herein, gray-scale value be 255 pixel represent cloud sector, and the pixel that gray-scale value is 0 represents non-cloud sector.
Through aforesaid operations, the cloud sector that can obtain as shown in Figure 5 differentiates result.
(3) in the filtered image of Robinson, different threshold value is adopted to carry out binarization segmentation to cloud sector and non-cloud sector part.Because the contrast difference between cloud sector target and cloud background is less, but not cloud sector part is then comparatively large, therefore takes Low threshold binarization segmentation in cloud sector, takes high threshold binarization segmentation in non-cloud sector.Concrete formula is:
Cloud sector part: f s e g ( i , j ) = 255 , f R o b i n s o n ( i , j ) > f L 0 , e l s e - - - ( 9 )
Non-cloud sector part: f s e g ( i , j ) = 255 , f R o b i n s o n ( i , j ) > f H 0 , e l s e - - - ( 10 )
F in formula (9) (10) land f hrepresent the binary-state threshold of two different sizes respectively.Wherein the former is Low threshold, and the latter is high threshold.As shown in Figure 6, comparatively complete extraction effect can be obtained after Double Thresholding Segmentation.
(4) continuous point target combination is carried out to the result after binarization segmentation, guarantees that target only stays next pixel size " isolated point ".Concrete operations are, are the pixel of 255 to gray-scale value after binaryzation, record its coordinate position, and be mapped in the filtered image of Robinson, pixel centered by it, offers the template (L refers to template length) of a L*L, as offered the template of a 5*5.If the gray-scale value of center pixel is the maximal value of all pixels in template, then retained in binary image, otherwise just rejected (namely utilizing " neighborhood non-maxima suppression " principle to reject multiple " the pseudo-impact points " of same target generation further to the result after binaryzation).As shown in Figure 7, each continuity point target is successfully merged into an isolated point to its result.
3. the step that essence detects is:
(1) to the image that continuous multiple frames operates through above-mentioned spatial domain, interframe track association is carried out.First set up flight path manager with list structure, from the second two field picture, successively the point coordinate in present frame and the point coordinate in previous frame are compared, if drop on certain some institute opened round Bo Mennei in previous frame, be then considered as being successfully associated.The ripple door of tradition track association is chosen and is usually needed the method utilizing prediction of speed dynamically to generate, but the number of pixels that infrared small target occupies in the picture is little, the amount of movement of interframe is also very little, and therefore ripple door radius can be taken as a definite value (a general desirable 1-2 pixel).As shown in Figure 8, the suspected target obtained after each rough detection all defines stable flight path (flight path is red annotate portions) in interframe.
(2) CFAR is carried out to the result after track association and suppress operation.The inventive method utilizes the difference that there is kinetic characteristic and gamma characteristic two aspect between real goal and false-alarm point, distinguishes both.Suppose that current time is k, object is v in the speed in k moment k, in the speed sum in front k-1 moment be order (wherein, alpha+beta=1).Because noise motion meets standardized normal distribution, so its Δ v levels off to 0 within a period of time, and the Δ v of real goal will level off to a positive number.Setting speed threshold value V t, as Δ v > V ttime, be real goal by this target-recognition, otherwise, be then false-alarm noise (kinetic characteristic);
Gamma characteristic aspect, because real goal is approach ours gradually in visual field, in the process that this draws near, its gray scale is inevitable is also ascending conversion, and the gray scale of noise remains unchanged substantially.Investigate gamma characteristic for " shade of gray " this physical quantity: suppose that current time is k, object is I in the gray scale in k moment k, be I in the gray scale in k-1 moment k-1, then the shade of gray J of current time object k=I k-I k-1, in the gradient sum in k-1 moment be before this order (wherein, alpha+beta=1).For noise, within a period of time, its Δ J levels off to 0, and the Δ J convergence of real goal is a positive number.Setting Grads threshold J t, as Δ J > J ttime, be real goal by this target-recognition, otherwise, be then considered as false-alarm noise.Based on above two criterions, the CFAR just completed in track association process suppresses operation.As shown in Figure 9, the flight path produced by false-alarm is not labeled as real goal by yellow box to concrete result.
For the advantage of the present invention in infrared small target detection reliability is described, use the inventive method to carry out simulation process to the mobile Small object video under the cloud background of 90 ° of Large visual angle infrared thermovision systems collection, and carry out statistical computation for false alarm rate and detectivity two indexs.Wherein, false alarm rate is defined as the detection number of times of interior false target per hour; Verification and measurement ratio is defined as the number percent of correct detection destination number relative to realistic objective quantity.
As shown in Figure 10, horizontal ordinate is hourage, and ordinate is wrong report number of times, sets up false alarm rate statistical graph.As calculated, under this scene, false alarm rate can control at 1.7 times/hour by the inventive method.
As shown in figure 11, horizontal ordinate is hourage, and ordinate is detectivity, sets up detectivity statistical graph.As calculated, under this scene, detectivity can control in 96% by the inventive method.

Claims (7)

1. under a cloud background based on the infrared small target detection method of template convolution and false alarm rejection, it is characterized in that for single-frame images, adopt the spatial domain operation of template convolution, first carry out max-medium filter to image to remove significant noise and complete Image semantic classification, secondly with Robinson's template convolution Background suppression, outstanding target, then cloud sector division is carried out to former figure, Low threshold is adopted to carry out binary conversion treatment in cloud sector part to the filtered result of Robinson, but not cloud sector part then adopts high threshold process, finally the result after binaryzation is rejected further multiple " the pseudo-impact points " of the generation of same target, thus complete " rough detection ", the consecutive frame image carrying out spatial processing is continued to take time domain operation process, namely track association is carried out in interframe, and the difference in gamma characteristic and kinetic characteristic carries out CFAR suppression operation for real goal and false-alarm point, complete " essence detects ", thus realize the detection of infrared small target.
2. according to the infrared small target detection method based on template convolution and false alarm rejection under the cloud background described in claim 1, it is characterized in that Image semantic classification process is: max-medium filter is carried out to input picture, remove the remarkable noise in image, i.e. max-medium filter water intaking flat, vertical, left 45 degree, 45 degree of four filtering directions, the right side, the maximal value of getting pixel grey scale intermediate value in all directions gives center pixel.
3., according to the infrared small target detection method based on template convolution and false alarm rejection under the cloud background described in claim 2, it is characterized in that the max-medium filter template for (2N+1) * (2N+1), its computing formula is:
f max-med(i,j)=max(z 1,z 2,z 3,z 4)(1)
In formula (1),
z 1=med[f(i,j-N),...,f(i,j),...,f(i,j+N)]
z 2=med[f(i-N,j),...,f(i,j),...,f(i+N,j)](2)
z 3=med[f(i+N,j-N),...,f(i,j),...,f(i-N,j+N)]
z 4=med[f(i-N,j-N),...,f(i,j),...,f(i+N,j+N)]
Wherein, pixel coordinate centered by (i, j), med is for getting median operation, and max is for getting maxima operation, and N characterizes template size size.
4., according to the infrared small target detection method based on template convolution and false alarm rejection under the cloud background described in claim 1, it is characterized in that the process of rough detection is:
(1) Robinson's template convolution is carried out to the image after max-medium filter, suppress cloud background, and outstanding target, this template utilizes this characteristic comparatively large of the gray difference between real goal and cloud background, the maximum gradation value of central pixel point gray-scale value and surrounding pixel is compared and arithmetical operation, if the gray scale of center pixel is comparatively strong, then can be retained, otherwise will be suppressed; Meanwhile, template-setup has isolation strip, ensures that the gamma characteristic of real goal is not destroyed;
(2) cloud sector differentiation is carried out to original image, take the method for template matches filtering in cloud, cloud exterior domain distinguished, this template is provided with blank isolation strip on the one hand and is not damaged to protect the gamma characteristic of real goal, on the other hand by the gray scale sum of all for central row pixels and the pixel grey scale sum of the top a line or bottom a line are carried out doing difference operation and normalizing to [0,255] interval, and utilize binaryzation to operate cloud sector segmentation;
(3) in the filtered image of Robinson, adopt different threshold value to carry out binarization segmentation to cloud sector and non-cloud sector part, take Low threshold binarization segmentation in cloud sector, take high threshold binarization segmentation in non-cloud sector, be specially:
Cloud sector part: f s e g ( i , j ) = 255 , f R o b i n s o n ( i , j ) > f L 0, e l s e - - - ( 9 )
Non-cloud sector part: f s e g ( i , j ) = 255 , f R o b i n s o n ( i , j ) > f H 0 , e l s e - - - ( 10 )
F in formula (9) (10) land f hrepresent the binary-state threshold of two different sizes respectively, wherein the former is Low threshold, and the latter is high threshold;
(4) continuous point target combination is carried out to the result after binarization segmentation, guarantee that target only stays next pixel size " isolated point ", namely be the pixel of 255 to gray-scale value after binaryzation, record its coordinate position, and be mapped in the filtered image of Robinson, pixel centered by it, offer the template of a L*L, if the gray-scale value of center pixel is the maximal value of all pixels in template, then retained in binary image, otherwise just rejected.
5. according to the infrared small target detection method based on template convolution and false alarm rejection under the cloud background described in claim 4, it is characterized in that the method for Robinson's template convolution is: for Robinson's Filtering Template of (2N+1) * (2N+1), its computing formula is:
f R o b i s o n ( i , j ) = f ( i , j ) - m a x ( z 1 , z 2 , z 3 , z 4 ) , f ( i , j ) > m a x ( z 1 , z 2 , z 3 , z 4 ) f ( i , j ) - m i n ( z 1 , z 2 , z 3 , z 4 ) , f ( i , j ) < m a x ( z 1 , z 2 , z 3 , z 4 ) 0 , e l s e - - - ( 3 )
In formula (3),
z 1=max[f(i-N:i+N,j-N)]
z 2=max[f(i-N:i+N,j+N)](4)
z 3=max[f(i-N,i-N:j+N)]
z 4=max[f(i+N,i-N:j+N)]
Wherein, pixel coordinate centered by (i, j), max is for getting maxima operation, and N characterizes template size size.
6., according to the infrared small target detection method based on template convolution and false alarm rejection under the cloud background described in claim 4, it is characterized in that the step that cloud sector differentiates is:
I () carries out matched filtering operation to original image, for the matched filtering template of (2M+1) * (2N+1) size, its computing formula is:
f c ( i , j ) = z 2 - z 3 , z 1 > z 3 z 2 - z 1 , e l s e - - - ( 5 )
In formula (5),
z 1=sum(f(i-N,j-M:j+M))
z 2=sum(f(i,j-M:j+M))(6)
z 3=sum(f(i+N,j-M:j+M))
Wherein, pixel coordinate centered by (i, j), sum is sum operation, and M, N characterize template size size;
(ii) carry out gray scale normalization process to the result after matched filtering, normalization formula is:
f c &prime; ( i , j ) = 255 , f c ( i , j ) > 255 0 , f c ( i , j ) < 0 f c ( i , j ) , e l s e - - - ( 7 )
(iii) carry out binarization segmentation to the result after normalization, the formula of binaryzation operation is:
f c l o u d ( i , j ) = 255 , f c &prime; ( i , j ) > f t h r e s h o l d 0 , e l s e - - - ( 8 )
F in formula (8) thresholdfor binary-state threshold, herein, gray-scale value be 255 pixel represent cloud sector, and the pixel that gray-scale value is 0 represents non-cloud sector.
7., according to the infrared small target detection method based on template convolution and false alarm rejection under the cloud background described in claim 4, it is characterized in that the step that essence detects is:
(1) to the image that continuous multiple frames operates through spatial domain, carry out interframe track association, first flight path manager is set up with list structure, from the second two field picture, successively the point coordinate in present frame and the point coordinate in previous frame are compared, if drop on certain some institute opened round Bo Mennei in previous frame, be then considered as being successfully associated; Ripple door radius is taken as a definite value;
(2) carry out CFAR to the result after track association and suppress operation, utilize the difference that there is kinetic characteristic and gamma characteristic two aspect between real goal and false-alarm point, distinguish both, namely suppose that current time is k, object is v in the speed in k moment k, in the speed sum in front k-1 moment be order (wherein, alpha+beta=1); Because noise motion meets standardized normal distribution, so its Δ v levels off to 0 within a period of time, and the Δ v of real goal will level off to a positive number; Setting speed threshold value V t, as Δ v > V ttime, be real goal by this target-recognition, otherwise, be then false-alarm noise;
For gamma characteristic aspect, investigate gamma characteristic for " shade of gray " this physical quantity: suppose that current time is k, object is I in the gray scale in k moment k, be I in the gray scale in k-1 moment k-1, then the shade of gray J of current time object k=I k-I k-1, in the gradient sum in k-1 moment be before this order (wherein, alpha+beta=1), for noise, within a period of time, its Δ J levels off to 0, and the Δ J convergence of real goal is a positive number; Setting Grads threshold J t, as Δ J > J ttime, be real goal by this target-recognition, otherwise, be then considered as false-alarm noise; Based on above two characteristics, just complete the false alarm rejection operation in track association process.
CN201510917030.3A 2015-12-10 2015-12-10 Infrared small target detection method under cloud background based on template convolution and false alarm rejection Expired - Fee Related CN105574855B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510917030.3A CN105574855B (en) 2015-12-10 2015-12-10 Infrared small target detection method under cloud background based on template convolution and false alarm rejection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510917030.3A CN105574855B (en) 2015-12-10 2015-12-10 Infrared small target detection method under cloud background based on template convolution and false alarm rejection

Publications (2)

Publication Number Publication Date
CN105574855A true CN105574855A (en) 2016-05-11
CN105574855B CN105574855B (en) 2019-01-22

Family

ID=55884946

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510917030.3A Expired - Fee Related CN105574855B (en) 2015-12-10 2015-12-10 Infrared small target detection method under cloud background based on template convolution and false alarm rejection

Country Status (1)

Country Link
CN (1) CN105574855B (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106056115A (en) * 2016-05-25 2016-10-26 西安科技大学 Infrared small target detection method under non-uniform background
CN106803260A (en) * 2016-12-28 2017-06-06 辽宁师范大学 Infrared ship activity of imagination contours segmentation method based on the convex optimization of local entropy
CN107633505A (en) * 2017-08-24 2018-01-26 南京理工大学 A kind of undercarriage detection method based on target gray distribution character
CN107644414A (en) * 2017-08-27 2018-01-30 南京理工大学 A kind of sea wake detection method based on constant statistics and Radon transform
CN109410137A (en) * 2018-10-11 2019-03-01 中国科学院上海技术物理研究所 A kind of detection method of dark weak signal target
CN109544535A (en) * 2018-11-26 2019-03-29 马杰 It is a kind of that camera detection method and system are pried through based on infrared cutoff filter optical filtration characteristic
CN110084231A (en) * 2019-04-29 2019-08-02 北京富吉瑞光电科技有限公司 A kind of method and apparatus of infrared panoramic system multi-target detection
CN111493853A (en) * 2020-04-24 2020-08-07 天津恒宇医疗科技有限公司 Blood vessel parameter evaluation method and system for angiodermic diseases
CN112967305A (en) * 2021-03-24 2021-06-15 南京莱斯电子设备有限公司 Image cloud background detection method under complex sky scene
CN113111883A (en) * 2021-03-23 2021-07-13 浙江大华技术股份有限公司 License plate detection method, electronic equipment and storage medium
CN113240664A (en) * 2021-06-03 2021-08-10 郑州航空工业管理学院 Infrared detection false alarm detection method based on scene space-time significance and application thereof
CN113776676A (en) * 2021-09-23 2021-12-10 中国航天科工集团八五一一研究所 Infrared small target detection method based on image curvature and gradient
CN115311460A (en) * 2022-08-16 2022-11-08 哈尔滨工业大学 Infrared small target detection method fusing time-space domain information under slow motion background
CN115631114A (en) * 2022-12-06 2023-01-20 北京九章星图科技有限公司 Dark and weak moving target indication enhancement on-track processing method based on time domain profile analysis
CN117315498A (en) * 2023-10-10 2023-12-29 中国人民解放军战略支援部队航天工程大学 False alarm discrimination method based on space target detection result

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070076917A1 (en) * 2003-03-21 2007-04-05 Lockheed Martin Corporation Target detection improvements using temporal integrations and spatial fusion
CN104125372A (en) * 2014-07-29 2014-10-29 北京机械设备研究所 Target photoelectric search and detection method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070076917A1 (en) * 2003-03-21 2007-04-05 Lockheed Martin Corporation Target detection improvements using temporal integrations and spatial fusion
CN104125372A (en) * 2014-07-29 2014-10-29 北京机械设备研究所 Target photoelectric search and detection method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
秦剑等: "基于光流估计和自适应背景抑制的弱小目标检测", 《光子学报》 *
陆恺立: "基于面阵探测器的红外搜索预警系统关键技术研究", 《万方学位论文数据库》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106056115B (en) * 2016-05-25 2019-01-22 西安科技大学 A kind of infrared small target detection method under non-homogeneous background
CN106056115A (en) * 2016-05-25 2016-10-26 西安科技大学 Infrared small target detection method under non-uniform background
CN106803260B (en) * 2016-12-28 2019-08-09 辽宁师范大学 Infrared ship activity of imagination contours segmentation method based on the convex optimization of local entropy
CN106803260A (en) * 2016-12-28 2017-06-06 辽宁师范大学 Infrared ship activity of imagination contours segmentation method based on the convex optimization of local entropy
CN107633505A (en) * 2017-08-24 2018-01-26 南京理工大学 A kind of undercarriage detection method based on target gray distribution character
CN107644414A (en) * 2017-08-27 2018-01-30 南京理工大学 A kind of sea wake detection method based on constant statistics and Radon transform
CN109410137A (en) * 2018-10-11 2019-03-01 中国科学院上海技术物理研究所 A kind of detection method of dark weak signal target
CN109410137B (en) * 2018-10-11 2021-10-01 中国科学院上海技术物理研究所 Method for detecting dim and weak target
CN109544535A (en) * 2018-11-26 2019-03-29 马杰 It is a kind of that camera detection method and system are pried through based on infrared cutoff filter optical filtration characteristic
CN110084231A (en) * 2019-04-29 2019-08-02 北京富吉瑞光电科技有限公司 A kind of method and apparatus of infrared panoramic system multi-target detection
CN111493853A (en) * 2020-04-24 2020-08-07 天津恒宇医疗科技有限公司 Blood vessel parameter evaluation method and system for angiodermic diseases
CN113111883A (en) * 2021-03-23 2021-07-13 浙江大华技术股份有限公司 License plate detection method, electronic equipment and storage medium
CN113111883B (en) * 2021-03-23 2023-06-06 浙江大华技术股份有限公司 License plate detection method, electronic device and storage medium
CN112967305A (en) * 2021-03-24 2021-06-15 南京莱斯电子设备有限公司 Image cloud background detection method under complex sky scene
CN112967305B (en) * 2021-03-24 2023-10-13 南京莱斯电子设备有限公司 Image cloud background detection method under complex sky scene
CN113240664A (en) * 2021-06-03 2021-08-10 郑州航空工业管理学院 Infrared detection false alarm detection method based on scene space-time significance and application thereof
CN113240664B (en) * 2021-06-03 2023-06-09 郑州航空工业管理学院 Infrared detection false alarm detection method based on scene space-time significance and application thereof
CN113776676A (en) * 2021-09-23 2021-12-10 中国航天科工集团八五一一研究所 Infrared small target detection method based on image curvature and gradient
CN115311460A (en) * 2022-08-16 2022-11-08 哈尔滨工业大学 Infrared small target detection method fusing time-space domain information under slow motion background
CN115631114A (en) * 2022-12-06 2023-01-20 北京九章星图科技有限公司 Dark and weak moving target indication enhancement on-track processing method based on time domain profile analysis
CN117315498A (en) * 2023-10-10 2023-12-29 中国人民解放军战略支援部队航天工程大学 False alarm discrimination method based on space target detection result
CN117315498B (en) * 2023-10-10 2024-05-24 中国人民解放军战略支援部队航天工程大学 False alarm discrimination method based on space target detection result

Also Published As

Publication number Publication date
CN105574855B (en) 2019-01-22

Similar Documents

Publication Publication Date Title
CN105574855A (en) Method for detecting infrared small targets under cloud background based on temperate filtering and false alarm rejection
CN103824070B (en) A kind of rapid pedestrian detection method based on computer vision
CN104899866B (en) A kind of intelligentized infrared small target detection method
CN108182690B (en) A kind of infrared Weak target detecting method based on prospect weighting local contrast
Fang et al. Infrared small UAV target detection based on depthwise separable residual dense network and multiscale feature fusion
CN103020606B (en) Pedestrian detection method based on spatio-temporal context information
CN103761731A (en) Small infrared aerial target detection method based on non-downsampling contourlet transformation
CN105869146B (en) SAR image change detection based on conspicuousness fusion
CN103985127B (en) The detection method of small target of a kind of intensive star background and device
CN105608691A (en) High-resolution SAR image individual building extraction method
Xia et al. A novel sea-land segmentation algorithm based on local binary patterns for ship detection
CN107507209B (en) Printogram extraction method of polarized SAR image
CN113743260B (en) Pedestrian tracking method under condition of dense pedestrian flow of subway platform
CN109711256B (en) Low-altitude complex background unmanned aerial vehicle target detection method
CN103593679A (en) Visual human-hand tracking method based on online machine learning
Zhang et al. A warning framework for avoiding vessel‐bridge and vessel‐vessel collisions based on generative adversarial and dual‐task networks
CN104268877A (en) Infrared image sea-sky-line self adaption detection method
CN103745216A (en) Radar image clutter suppression method based on airspace feature
CN102129559B (en) SAR (Synthetic Aperture Radar) image object detection method based on Primal Sketch algorithm
CN108010065A (en) Low target quick determination method and device, storage medium and electric terminal
Zhao et al. An adaptation of CNN for small target detection in the infrared
CN102663385A (en) Detection method for spot target on satellite
Sahba et al. 3D Object Detection Based on LiDAR Data
CN105022101B (en) Severe Convective Cloud Cluster method for tracing
Lian et al. A novel method on moving-objects detection based on background subtraction and three frames differencing

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190122

Termination date: 20201210

CF01 Termination of patent right due to non-payment of annual fee