CN101295401A - Infrared point target detecting method based on linear PCA - Google Patents
Infrared point target detecting method based on linear PCA Download PDFInfo
- Publication number
- CN101295401A CN101295401A CNA200810038548XA CN200810038548A CN101295401A CN 101295401 A CN101295401 A CN 101295401A CN A200810038548X A CNA200810038548X A CN A200810038548XA CN 200810038548 A CN200810038548 A CN 200810038548A CN 101295401 A CN101295401 A CN 101295401A
- Authority
- CN
- China
- Prior art keywords
- image
- target
- infrared
- subimage
- reconstructed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
The invention relates to a detection method based on a linear PCA for infrared point targets, comprising the following steps of: first, creating a training sample of the infrared point target through an improved gaussian gray model, and carrying out trainings for the PCA through the samples, to produce a group of principal components; then intercepting and acquiring sub-images from an infrared image to be detected, and respectively projecting the sub-images to each principal component, to reconstruct the sub-images through obtained projection coefficients and the principal components, and calculate the reconstruction error; finally, substituting a measure function with the reconstruction error, to produce a detection image and improve the testability of the target. The detection method based on the linear PCA for infrared point targets can greatly improve the testability of the infrared point targets. Different from the filtration-based detection method, the detection method based on the linear PCA for infrared point targets does not require a preprocessing for the infrared images to be detected.
Description
Technical field
The present invention relates to a kind of infrared target detection method of technical field of image processing, specifically is a kind of infrared point target detecting method based on linear PCA.
Background technology
In video monitoring and search system, the infrared photography The Application of Technology has very big singularity, can all weather operations, and be a kind of passive passive detection technology, can be hidden well when finding target oneself.Wherein, how detecting target under remote condition, is one of gordian technique of infrared monitoring and search system.If can detect target on farther distance, system just has the time of more target being reacted, for the system status that has the initiative in antagonism creates conditions.In addition, target detection is the first step that system reacts to target, and the performance of target detection directly has influence on the effect of subsequent treatment (such as target following and Target Recognition).Therefore, the infrared target detection technique is for the infrared video system, and is significant.
When the target range ir imaging system was far away, the area of target on imager was very little, is the little target of point-like (usually all less than 100 pixels).Because the influence of factors such as the contrast of infrared point target is low, area is little, the interference of texture and edge fog, noise and clutter is stronger makes the detection of infrared point target very difficult.Traditional infrared point target detecting method great majority are based on wave filter, have based on the method for airspace filter device: the filtering of high pass template, medium filtering, mathematical morphology filter, the filtering of local standard difference etc.; Method based on frequency domain filter has: desirable high-pass filtering, Butterworth high-pass filtering, Gauss's high-pass filtering method etc.If view data (in the image be that the gray-scale value of all pixels in the subimage at center arrange the vector that form with each pixel) is considered as the data set is made up of target class data and background classes data, the problem of detection target just converts to target data is discerned the problem that (classification) comes out from view data from image.Therefore, some mode identification methods can be used for realizing target detection.PCA (principal component analysis (PCA)) is an important method in the pattern recognition theory, often is used to realize data compression.Obtain describing the proper vector (principal component, major component) of training data by training, promptly Zui Da several characteristic is worth pairing proper vector.PCA has also extracted the feature of training data when realizing data compression.
Find by prior art documents, China's application number: 200410068024.7, patent name is " morphologic filter automatic destination detecting method ", this patent detects target based on the morphologic filtering method of Top-hat operator, concrete principle is: little target is in high band, and background is in low-frequency range, like this by High frequency filter or deduct the image of low-pass filtering gained with original image, the purpose that can realize giving prominence to target and suppress background, thus target detection finished.Obviously, this method is when detecting target, and some noise spots (also being in high band) also have bigger output.Like this, noise spot can be a target by flase drop also, and has improved false alarm rate, and this is based on the inherent shortcoming of the infrared point target detecting method of wave filter.
Summary of the invention
Of the present inventionly propose a kind of infrared point target detecting method, can realize detection the infrared point target of various wave bands based on linear PCA in order to solve the above-mentioned problems in the prior art.
The present invention is achieved by the following technical solutions, the present invention includes following steps:
At first,, PCA is trained, produce one group of principal component by these samples with the training sample of improved Gauss's gray level model generation infrared point target;
Then, from infrared image to be detected, intercept subimage, and it is projected to respectively on each principal component, this subimage is reconstructed, and calculate reconstructed error with the projection coefficient and the principal component that obtain;
At last,, produce detected image, improved the detectability of target reconstructed error substitution detection function.
Described improvement Gauss gray level model produces the infrared point target training sample, is specially: at traditional Gauss's gray level model
The basis on, add two constraint conditions
With
Make it to become improved Gauss's gray level model:
So just finished the improvement of Gauss's gray level model.Wherein, I
MaxBe target's center's pixel value (gray scale peak value); σ
xBe horizontal dispersion parameter, σ
yBe the vertical dispersion parameter, controlling the stroll characteristic of object pixel; (x
0, y
0) be the centre coordinate of target image; (i j) is other pixel coordinate of target image.I
Max, σ
xAnd σ
yThree parameters have been determined the characteristic of Gauss's gray level model.The infrared point target that traditional Gauss's gray level model produces can produce some and " point target " can not occur in actual infrared image, (works as I such as producing some point targets that are similar to noise spot
MaxBigger, and σ
xAnd σ
yHour), if train PCA, can reduce its detection performance greatly with such training sample (sample exterior point Outlier).
Describedly with sample PCA is trained, produce one group of principal component, be specially: produce some infrared point target images with above-mentioned improved Gauss's gray level model, these images are as training sample.With the grey scale pixel value of the image pattern of each infrared point target of producing, joining by the row first place is arranged in vector, with these vectors PCA is trained, and the generation principal component is got n proper vector of character pair value maximum, finishes target detection.The n value detects effect by compromise consideration and computation complexity (amplitude of dimensionality reduction) is decided, and is preferably n=6.
Describedly intercept subimage from infrared image to be detected, and it is projected to respectively on each principal component, with the projection coefficient and the principal component that obtain this subimage is reconstructed, and calculate reconstructed error, concrete steps are:
Step 1, (x, y) position intercepting subimage P in detected image
s, and generate vectorial P
Sc, with P
ScProject to respectively on n the principal component, obtain n projection coordinate, formula is
Wherein, ω
kBe P
ScK proper vector
On projection coordinate, Γ is an average criterion;
Step 2 is reconstructed subimage with n principal component and n projection coordinate, obtains reconstructed image, and formula is
Wherein,
Be reconstructed image;
Step 3, the reconstructed error of calculating subimage, formula is
Wherein, ε
RCBe reconstructed error.If subimage P
sBe target image, so reconstructed image
With P
sSimilarity just very big, otherwise then similarity is very little, so reconstructed error be exactly this subimage whether be that a kind of of target image estimates, the smaller value correspondence of reconstructed error the existence of target.
Described with reconstructed error substitution detection function, produce detected image, be specially: with reconstructed error ε
RCThe substitution detection function
Different with reconstructed error is ζ
D(x, higher value correspondence y) the existence of target, and this is consistent with human intuition; And ζ
D(x, value y) has saved the normalization link between 0~1, be more convenient for calculating detection indexs such as signal to noise ratio (snr).The value of parameter δ wherein is chosen for δ=0.25 by experience.With detection function value ζ corresponding on each location of pixels
D(x y), as the gray-scale value of this pixel, has just generated detected image.
Traditional infrared point target detection great majority are based on filtering method, when detecting target, also the noise spot erroneous judgement are impact point, have improved false alarm rate, have increased the load of system.The present invention is that the method with pattern-recognition detects infrared point target, PCA is an important branch of pattern-recognition, often be used to realize data compression, but in packed data, PCA has also extracted clarification of objective, utilize this characteristic, PCA is used to realize the detection of infrared point target, its essence is with PCA in detected image point target " identification " to be come out.For recognition methods, training sample is very big to the influence of its discrimination, and the present invention at first improves Gauss's gray level model, produces the training sample of infrared point target then with improved Gauss's gray level model.Through PCA is trained, produce principal component, detected subimage (intercepting from detected image) is projected to respectively on each principal component, obtain the respective projection coordinate, by the reconstruct of projection coordinate and principal component to subimage, calculate reconstructed error, less reconstructed error correspondence the existence of target.A detection function is devised, and is used for producing detected image, and in detected image, the detectability of infrared point target is improved greatly.
The present invention can improve the detectability of infrared point target largely, and different with method based on filtering, this method does not need infrared image to be detected is carried out pre-service.
Description of drawings
Fig. 1 is the training sample by the infrared point target of improved Gauss's gray level model generation
Fig. 2 is the synoptic diagram of the embodiment of the invention
Wherein: (a1)、(b
1) be respectively and comprise two original infrared images of 10 point targets, figure (b1) in target 6 is a real Ship Target, and other target in two width of cloth images all is to be produced by improved Gauss's gray level model Simulation objectives embed and to obtain in the infrared image, all targets among the figure all are Weak targets, detection difficulty is bigger; (a2)、(b
2) be respectively PCA to (a1)、(b
1) result that detects; (a3)、(b
3) be respectively (a2)、(b
2) 3D Structure chart.
Embodiment
Below in conjunction with specific embodiment technical scheme of the present invention is elaborated: present embodiment is being to implement under the prerequisite with the technical solution of the present invention; provided detailed embodiment and concrete operating process, but protection scope of the present invention is not limited to following embodiment.
Present embodiment has mainly comprised the content of two aspects, i.e. the production method of infrared point target training sample and based on the infrared point target detecting method of PCA.
As shown in Figure 1, be the training sample of the infrared point target that produced by improved Gauss's gray level model, wherein: the size of training sample is 11 * 11, and the image that is produced by improved Gauss's gray level model rotates at any angle and obtains.
(1) training
1) grey scale pixel value of 121 target training sample subimages (11 * 11) is arranged by end to end order, generated column vector (121 dimension);
2) column vector is arranged generator matrix, and matrix is carried out characteristic value decomposition, obtain optimum description clarification of objective vector, and proper vector is arranged from big to small by eigenwert; Present embodiment is got preceding 6 proper vectors
These 6 proper vectors are to describe an eigenvectors of target.
(2) detect
3) each location of pixels in detected image is used the sliding window with target sample picture size size identical (11 * 11), intercepts subimage to be detected, and converts thereof into column vector (121 dimension);
4) column vector that subimage to be detected generated is projected on 6 proper vectors, correspondingly obtain 6 projection coordinate { ω
1, ω
2... ω
n;
5) with these 6 proper vectors and projection coordinate subimage to be detected is reconstructed, and calculates reconstructed error ε
RC
6) with ε
RCSubstitution detection function ζ
D(x y), calculates the detection function value on each location of pixels, with the gray-scale value of these functional values as respective pixel, generates detected image.
Fig. 2 is the testing result of embodiment, (a among the figure
2), (b
2) be respectively PCA to (a
1), (b
1) result that detects; (a
3), (b
3) be respectively (a
2), (b
2) the 3D structural drawing.Table 1 has provided the detection index of Fig. 2, is respectively snr gain (SNRG) and background inhibiting factor (BSF).Snr gain is the tolerance to the signal to noise ratio (S/N ratio) load-carrying capacity, and computing formula is SNRG=SNR
Out/ SNR
In, SNR wherein
OutBe the signal to noise ratio (S/N ratio) of detected image ((a2) among Fig. 2 and (b2)), SNR
InSignal to noise ratio (S/N ratio) for original image ((a1) among Fig. 2 and (b1)); The background inhibiting factor is background to be suppressed the tolerance of ability, and computing formula is BSF=C
In/ C
Out, C wherein
InBe the standard deviation of original image background, C
OutStandard deviation for the detected image background.Two detection indexs are big more, and the detection performance of detection method is good more.Can find out that from table 1 a few target in Fig. 2 (b1), the snr gain of the inventive method gained illustrates that all greater than 1 this method can improve the detectability of infrared target.In addition, to two width of cloth images of embodiment, the background inhibiting factor illustrates that all much larger than 1 this method has very strong background to suppress ability.
Table 1 three seed space arithmetics compare the performance index of target detection in the accompanying drawing 2
Claims (5)
1. the infrared point target detecting method based on linear PCA is characterized in that, comprises the steps:
At first,, PCA is trained, produce one group of principal component by these samples with the training sample of improved Gauss's gray level model generation infrared point target;
Then, from infrared image to be detected, intercept subimage, and it is projected to respectively on each principal component, this subimage is reconstructed, and calculate reconstructed error with the projection coefficient and the principal component that obtain;
At last, with reconstructed error substitution detection function, produce detected image;
Described improved Gauss's gray level model is specially:
Wherein, I
MaxIt is target's center's pixel value; σ
xBe horizontal dispersion parameter, σ
yBe the vertical dispersion parameter, controlling the stroll characteristic of object pixel; (x
0, y
0) be the centre coordinate of target image; (i j) is other pixel coordinate of target image, I
Max, σ
xAnd σ
yThree parameters have been determined the characteristic of Gauss's gray level model.
2. the infrared point target detecting method based on linear PCA according to claim 1, it is characterized in that, describedly PCA is trained with sample, produce one group of principal component, be specially: produce some infrared point target images with improved Gauss's gray level model, these images are as training sample, grey scale pixel value with the image pattern of each infrared point target of producing, join and be arranged in vector by the row first place, with these vectors PCA is trained, produce principal component, get n proper vector of character pair value maximum, finish target detection.
3. the infrared point target detecting method based on linear PCA according to claim 2 is characterized in that, described n value is 6.
4. the infrared point target detecting method based on linear PCA according to claim 1, it is characterized in that, describedly from infrared image to be detected, intercept subimage, and it is projected to respectively on each principal component, with projection coefficient that obtains and principal component this subimage is reconstructed, and the calculating reconstructed error, concrete steps are:
Step 1, (x, y) position intercepting subimage P in detected image
s, and generate vectorial P
Sc, with P
ScProject to respectively on n the principal component, obtain n projection coordinate, formula is
Wherein, ω
kBe P
ScK proper vector
On projection coordinate, Γ is an average criterion;
Step 2 is reconstructed subimage with n principal component and n projection coordinate, obtains reconstructed image, and formula is
Wherein,
Be reconstructed image;
Step 3, the reconstructed error of calculating subimage, formula is
Wherein, ε
RCBe reconstructed error, if subimage P
sBe target image, so reconstructed image
With P
sSimilarity just very big, otherwise then similarity is very little, so reconstructed error be exactly this subimage whether be that a kind of of target image estimates.
5. the infrared point target detecting method based on linear PCA according to claim 1 is characterized in that, and is described with reconstructed error substitution detection function, produces detected image, is meant: with reconstructed error ε
RCThe substitution detection function
The value of parameter δ wherein is taken as δ=0.25, with detection function value ζ corresponding on each location of pixels
D(x y), as the gray-scale value of this pixel, has just generated detected image.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA200810038548XA CN101295401A (en) | 2008-06-05 | 2008-06-05 | Infrared point target detecting method based on linear PCA |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA200810038548XA CN101295401A (en) | 2008-06-05 | 2008-06-05 | Infrared point target detecting method based on linear PCA |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101295401A true CN101295401A (en) | 2008-10-29 |
Family
ID=40065663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA200810038548XA Pending CN101295401A (en) | 2008-06-05 | 2008-06-05 | Infrared point target detecting method based on linear PCA |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101295401A (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799875A (en) * | 2010-02-10 | 2010-08-11 | 华中科技大学 | Target detection method |
CN101882314A (en) * | 2010-07-20 | 2010-11-10 | 上海交通大学 | Infrared small target detection method based on overcomplete sparse representation |
CN101957993A (en) * | 2010-10-11 | 2011-01-26 | 上海交通大学 | Adaptive infrared small object detection method |
CN105023255A (en) * | 2015-08-05 | 2015-11-04 | 西安电子科技大学 | Infrared weak small target image sequence simulation method for jittering of Gaussian model area-array camera |
CN106643731A (en) * | 2016-12-29 | 2017-05-10 | 凌云光技术集团有限责任公司 | System and method for tracking and measuring point target |
CN108320024A (en) * | 2017-01-16 | 2018-07-24 | 佳能株式会社 | Dictionary creating apparatus and method, apparatus for evaluating and method and storage medium |
CN109511630A (en) * | 2017-09-20 | 2019-03-26 | 南京理工大学 | Intelligent sprayer with target following function |
CN109741396A (en) * | 2018-12-12 | 2019-05-10 | 天津津航技术物理研究所 | A kind of extremely small and weak infrared target detection method |
-
2008
- 2008-06-05 CN CNA200810038548XA patent/CN101295401A/en active Pending
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799875B (en) * | 2010-02-10 | 2011-11-30 | 华中科技大学 | Target detection method |
CN101799875A (en) * | 2010-02-10 | 2010-08-11 | 华中科技大学 | Target detection method |
CN101882314A (en) * | 2010-07-20 | 2010-11-10 | 上海交通大学 | Infrared small target detection method based on overcomplete sparse representation |
CN101882314B (en) * | 2010-07-20 | 2012-06-20 | 上海交通大学 | Infrared small target detection method based on overcomplete sparse representation |
CN101957993A (en) * | 2010-10-11 | 2011-01-26 | 上海交通大学 | Adaptive infrared small object detection method |
CN105023255B (en) * | 2015-08-05 | 2018-01-30 | 西安电子科技大学 | The Infrared DIM-small Target Image sequence emulation mode of Gauss model area array cameras shake |
CN105023255A (en) * | 2015-08-05 | 2015-11-04 | 西安电子科技大学 | Infrared weak small target image sequence simulation method for jittering of Gaussian model area-array camera |
CN106643731A (en) * | 2016-12-29 | 2017-05-10 | 凌云光技术集团有限责任公司 | System and method for tracking and measuring point target |
CN106643731B (en) * | 2016-12-29 | 2019-09-27 | 凌云光技术集团有限责任公司 | A kind of pair of point target carries out the system and method for tracking measurement |
CN108320024A (en) * | 2017-01-16 | 2018-07-24 | 佳能株式会社 | Dictionary creating apparatus and method, apparatus for evaluating and method and storage medium |
CN108320024B (en) * | 2017-01-16 | 2022-05-31 | 佳能株式会社 | Dictionary creation device and method, evaluation device and method, and storage medium |
US11521099B2 (en) | 2017-01-16 | 2022-12-06 | Canon Kabushiki Kaisha | Dictionary generation apparatus, evaluation apparatus, dictionary generation method, evaluation method, and storage medium for selecting data and generating a dictionary using the data |
CN109511630A (en) * | 2017-09-20 | 2019-03-26 | 南京理工大学 | Intelligent sprayer with target following function |
CN109511630B (en) * | 2017-09-20 | 2021-09-21 | 南京理工大学 | Intelligent sprayer with target tracking function |
CN109741396A (en) * | 2018-12-12 | 2019-05-10 | 天津津航技术物理研究所 | A kind of extremely small and weak infrared target detection method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101295401A (en) | Infrared point target detecting method based on linear PCA | |
US10852420B2 (en) | Object detection system, autonomous vehicle using the same, and object detection method thereof | |
CN109636742B (en) | Mode conversion method of SAR image and visible light image based on countermeasure generation network | |
CN103914847A (en) | SAR image registration method based on phase congruency and SIFT | |
CN101882314B (en) | Infrared small target detection method based on overcomplete sparse representation | |
Chen et al. | Feature-enhanced speckle reduction via low-rank and space-angle continuity for circular SAR target recognition | |
CN108090888A (en) | The infrared image of view-based access control model attention model and the fusion detection method of visible images | |
CN105809693A (en) | SAR image registration method based on deep neural networks | |
CN103761731A (en) | Small infrared aerial target detection method based on non-downsampling contourlet transformation | |
CN104834915B (en) | A kind of small infrared target detection method under complicated skies background | |
CN102494675A (en) | High-speed visual capturing method of moving target features | |
CN109740522B (en) | Personnel detection method, device, equipment and medium | |
CN103295221B (en) | The waterborne target method for testing motion of simulation compound eye visual mechanism and polarization imaging | |
CN102231844B (en) | Video image fusion performance evaluation method based on structure similarity and human vision | |
CN101957993A (en) | Adaptive infrared small object detection method | |
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 | |
CN103324921B (en) | A kind of mobile identification method based on interior finger band and mobile identification equipment thereof | |
Jing et al. | SAR image de-noising via grouping-based PCA and guided filter | |
CN107886526A (en) | Sequence image weak and small target detection method based on time domain filtering | |
CN103065320A (en) | Synthetic aperture radar (SAR) image change detection method based on constant false alarm threshold value | |
CN117665807A (en) | Face recognition method based on millimeter wave multi-person zero sample | |
CN107479052A (en) | Ground concealed target detection method based on Generalized Gaussian Distribution Model | |
CN111126508A (en) | Hopc-based improved heterogeneous image matching method | |
CN101847254A (en) | Pre-processing method of infrared small point target detection based on optimal design of matched filter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Open date: 20081029 |