CN109410223A - A kind of SAR image segmentation method based on watershed algorithm and dictionary learning - Google Patents
A kind of SAR image segmentation method based on watershed algorithm and dictionary learning Download PDFInfo
- Publication number
- CN109410223A CN109410223A CN201811322086.4A CN201811322086A CN109410223A CN 109410223 A CN109410223 A CN 109410223A CN 201811322086 A CN201811322086 A CN 201811322086A CN 109410223 A CN109410223 A CN 109410223A
- Authority
- CN
- China
- Prior art keywords
- data
- image
- histogram
- algorithm
- dictionary
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/155—Segmentation; Edge detection involving morphological operators
Abstract
The invention discloses the SAR image segmentation methods of a kind of watershed algorithm and dictionary learning, this method comprises: pre-processing using Denoising Algorithm to original image, so that image is become more smooth, while extracting three layers of wavelet character;Initial partitioning is carried out to original image using K-means algorithm;The morphometric characters of training image are found out, calculate floating point activities image, then be split with watershed algorithm;Data in each pocket regard the data acquisition system of a kind of data as, by all data in the set respectively with the atom computing arest neighbors in each category dictionary, we can obtain a statistic histogram to corresponding each category dictionary, calculate the error between such histogram that the histogram for corresponding to every one kind and original training data obtain;Judge that the error on which kind of is minimum, which kind of the data in set are just divided into, successively all pockets are divided, obtain final segmentation result.Compared with traditional SAR image segmentation method, validity and accuracy of the present invention in identification process are higher, and algorithm complexity is lower.
Description
Technical field
The invention belongs to synthetic aperture radar (SAR) image application fields, are related to a kind of SAR image segmentation method, especially
It is a kind of SAR image segmentation method based on watershed algorithm and dictionary learning.
Background technique
Since SAR imaging technique has the advantage round-the-clock, the factors such as climate do not influence so that it in national product and
Very big effect has been played in construction.And also just become the hot spot and emphasis of people's research to the processing of SAR image, wherein
SAR image segmentation is as to the fundamental operation during image procossing, the also concern by more and more people.Spectral clustering
Algorithm is common a kind of algorithm in processing image segmentation problem, but since the algorithm is difficult to handle the data of magnanimity, so its
Using being restricted.
Dividing method based on mathematical morphology is method important always in the field of image segmentation.Based on Mathematical Morphology
A kind of mathematical tool of the method based on image aspects, wherein more classical algorithm, such as watershed algorithm
(watershed algorithm).Watershed algorithm has many good qualities, his scale is intuitive, and the width of cut-off rule is single pixel
And closure is continuous, so being paid attention in recent years by many people, the method for dictionary learning is successfully applied to many fields, such as:
Classification, segmentation, identification, super-resolution, denoising etc..
It is that everybody is common using the pixel in image as the method that a sample is handled in image segmentation problem
Method, and traditional clustering method has spectral clustering etc., but the algorithm is difficult to that a large amount of data are effectively treated.To understand
The problem of certainly above-mentioned SAR image is divided, the invention proposes the SAR image segmentation sides based on watershed algorithm and dictionary learning
Method.Initial segmentation is carried out to image first, obtains the training sample of some tape labels, and the data of every a kind of tape label, benefit
It obtains a statistic histogram under such obtained dictionary with arest neighbors, counts in dictionary atomic distance training data most
Close number.Then in order to divide an image into some pockets, make over-segmentation of watershed, then by inside with picture
Vegetarian refreshments is that unit is handled, and by using the method for statistic histogram, is sorted out, obtains final segmentation result.
Summary of the invention
In view of above-mentioned deficiencies of the prior art, the technical problem to be solved by the present invention is to how be directed to high resolution SAR
Image carries out the segmentation identification of efficiently and accurately using the principle of watershed algorithm and dictionary learning.
To achieve the above object, the present invention provides the SAR image segmentation methods of a kind of watershed algorithm and dictionary learning.
Its feature includes:
(1) SAR image pre-processes: since there are a large amount of speckle noises for SAR image itself, therefore first with Denoising Algorithm pair
Original image is pre-processed, and so that image is become more smooth, while extracting three layers of wavelet character;
(2) initial partitioning is carried out to original image using K-means algorithm: selects some lean on from the segmentation result of every one kind
The data of such nearly cluster centre;Every one kind data obtain such dictionary, the instruction of every one kind using K-means algorithm
The arest neighbors for practicing data calculating and the category dictionary atom, obtains a statistic histogram, and every one kind data all obtain a histogram
Figure;
(3) initial partitioning is carried out using watershed algorithm: to treated in (2) training image, finding out its morphocline
Image calculates floating point activities image, then is split with watershed algorithm;
(4) region division: to the image after segmentation in (3), the data in each pocket regard a kind of data as
Data acquisition system corresponds to every a kind of word by all data in the set respectively with the atom computing arest neighbors in each category dictionary
We can obtain a statistic histogram to allusion quotation, calculate and be somebody's turn to do corresponding to what the histogram and original training data of every one kind obtained
Error between the histogram of class;Judge that the error on which kind of is minimum, which kind of the data in set are just divided into, according to
It is secondary to divide all pockets, obtain final segmentation result.
Further, in described (1) using open form state reconstruction filter remove in image some inessential details and
Small noise.When calculating wavelet character, three layers of Stationary Wavelet Transform are carried out to original image, obtain coefficient matrix coefm1(i1,
j1), m1=1,10, as m1=1, coefm1(i1,j1) represent low frequency coefficient;As m1 > 1, coefm1(i1,j1) represent height
Frequency coefficient.10 dimension sub-belt energy feature e (i, j)=[e are extracted to each pixel1(i,j),...,e10(i,j)]T, as this
The wavelet character of pixel:
Wherein, w × w is the size of sliding window, coefm1(i1,j1) it is i-th in stationary wavelet subband1Row jth1Column correspond to
Coefficient value.
The classification number for needing to cluster in (2) is K, is selected in the segmentation result of every one kind some in such cluster
The data of the heart;Every one kind data obtain such dictionary using K-means algorithm, the training data of every one kind calculate and
The arest neighbors of the category dictionary atom, obtains a statistic histogram, and every one kind data all obtain a histogram;
Morphometric characters subtract corrosion transformation by dilation transformation and obtain in (3):
Wherein, b is structural element, and Θ indicates erosion operation,Dilation operation is indicated, if having grey scale change in image very
Big pocket, then the value of structural element will suitably reduce, and otherwise may be filtered.Floating point activities image definition are as follows:
Fimg (f)=grad (f) * grad (f)/255.0
The pocket for generating watershed algorithm over-segmentation in (4) acquires statistics with the training sample of tape label
Histogram seeks residual error by the statistic histogram as dictionary, to determine which kind of the region belongs to.Scheme with traditional SAR
As dividing method is compared, validity and accuracy of the present invention in identification process are higher, and algorithm complexity is lower.
The technical effect of design of the invention, concrete scheme and generation will be described further below, with fully
Solve the purpose of the present invention, feature and effect.
Specific embodiment
A specific embodiment of the invention addressed below
(1) remove some inessential details and small noise in image using open form state reconstruction filter.It is counting
When calculating wavelet character, three layers of Stationary Wavelet Transform are carried out to original image, obtain coefficient matrix coefm1(i1,j1), m1=1,10, when
When m1=1, coefm1(i1,j1) represent low frequency coefficient;As m1 > 1, coefm1(i1,j1) represent high frequency coefficient.To each pixel
Point extracts 10 dimension sub-belt energy feature e (i, j)=[e1(i,j),...,e10(i,j)]T, wavelet character as the pixel:
Wherein, w × w is the size of sliding window, coefm1(i1,j1) it is i-th in stationary wavelet subband1Row jth1Column correspond to
Coefficient value.
(2) initial partitioning is carried out to original image using K-means algorithm: selects some lean on from the segmentation result of every one kind
The data of such nearly cluster centre;Every one kind data obtain such dictionary, the instruction of every one kind using K-means algorithm
The arest neighbors for practicing data calculating and the category dictionary atom, obtains a statistic histogram, and every one kind data all obtain a histogram
Figure;
(3) corrosion transformation is subtracted by dilation transformation and obtains morphometric characters:
Wherein, b is structural element, and Θ indicates erosion operation,Dilation operation is indicated, if having grey scale change in image very
Big pocket, then the value of structural element will suitably reduce, and otherwise may be filtered.Floating point activities image definition are as follows:
Fimg (f)=grad (f) * grad (f)/255.0
(4) data in each pocket regard the data acquisition system of a kind of data as, by all data in the set
Respectively with the atom computing arest neighbors in each category dictionary, corresponding to each category dictionary, we can obtain a statistics histogram
Figure calculates the error between such histogram that the histogram for corresponding to every one kind and original training data obtain;Judge
Error on which kind of is minimum, which kind of the data in set are just divided into, successively divides all pockets,
Obtain final segmentation result.
The preferred embodiment of the present invention has been described in detail above.It should be appreciated that the ordinary skill of this field is without wound
The property made labour, which according to the present invention can conceive, makes many modifications and variations.Therefore, all technician in the art
Pass through the available technology of logical analysis, reasoning, or a limited experiment on the basis of existing technology under this invention's idea
Scheme, all should be within the scope of protection determined by the claims.
Claims (4)
1. the SAR image segmentation method of a kind of watershed algorithm and dictionary learning characterized by comprising
Step (1) SAR image pretreatment: since there are a large amount of speckle noises for SAR image itself, therefore first with Denoising Algorithm pair
Original image is pre-processed, and so that image is become more smooth, while extracting three layers of wavelet character;
Step (2) carries out initial partitioning to original image using K-means algorithm: selecting some lean on from the segmentation result of every one kind
The data of such nearly cluster centre;Every one kind data obtain such dictionary, the instruction of every one kind using K-means algorithm
The arest neighbors for practicing data calculating and the category dictionary atom, obtains a statistic histogram, and every one kind data all obtain a histogram
Figure;
Step (3) carries out initial partitioning using watershed algorithm: to treated in (2) training image, finding out its morphocline
Image calculates floating point activities image, then is split with watershed algorithm;
Step (4) region division: to the image after segmentation in (3), the data in each pocket regard a kind of data as
Data acquisition system corresponds to every a kind of word by all data in the set respectively with the atom computing arest neighbors in each category dictionary
We can obtain a statistic histogram to allusion quotation, calculate and be somebody's turn to do corresponding to what the histogram and original training data of every one kind obtained
Error between the histogram of class;Judge that the error on which kind of is minimum, which kind of the data in set are just divided into, according to
It is secondary to divide all pockets, obtain final segmentation result.
2. as removed in image some inessential details using open form state reconstruction filter in claim 1 and small making an uproar
Sound carries out three layers of Stationary Wavelet Transform to original image, obtains coefficient matrix coef when calculating wavelet characterm1(i1,j1), m1=
1,10, as m1=1, coefm1(i1,j1) represent low frequency coefficient;As m1 > 1, coefm1(i1,j1) high frequency coefficient is represented, it is right
Each pixel extracts 10 dimension sub-belt energy feature e (i, j)=[e1(i,j),...,e10(i,j)]T, as the small of the pixel
Wave characteristic:
Wherein, w × w is the size of sliding window, coefm1(i1,j1) it is i-th in stationary wavelet subband1Row jth1Arrange corresponding system
Numerical value.
3. being selected in the segmentation result of every one kind some close to such cluster if the classification number for needing to cluster in claim 1 is K
The data at center;Every one kind data obtain such dictionary, the training data calculating of every one kind using K-means algorithm
With the arest neighbors of the category dictionary atom, a statistic histogram is obtained, every one kind data all obtain a histogram;
Morphometric characters subtract corrosion transformation by dilation transformation and obtain in its described step (3):
Wherein, b is structural element, and Θ indicates erosion operation,Dilation operation is indicated, if there is grey scale change very big in image
Pocket, then the value of structural element will suitably reduce, and otherwise may be filtered, floating point activities image definition are as follows:
Fimg (f)=grad (f) * grad (f)/255.0
4. acquiring system with the training sample of tape label such as the pocket for generating watershed algorithm over-segmentation in claim 1
Histogram is counted, residual error is sought by the statistic histogram as dictionary, to determine which kind of the region belongs to.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811322086.4A CN109410223A (en) | 2018-11-07 | 2018-11-07 | A kind of SAR image segmentation method based on watershed algorithm and dictionary learning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811322086.4A CN109410223A (en) | 2018-11-07 | 2018-11-07 | A kind of SAR image segmentation method based on watershed algorithm and dictionary learning |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109410223A true CN109410223A (en) | 2019-03-01 |
Family
ID=65472246
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811322086.4A Pending CN109410223A (en) | 2018-11-07 | 2018-11-07 | A kind of SAR image segmentation method based on watershed algorithm and dictionary learning |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109410223A (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101510309A (en) * | 2009-03-30 | 2009-08-19 | 西安电子科技大学 | Segmentation method for improving water parting SAR image based on compound wavelet veins region merge |
CN102129573A (en) * | 2011-03-10 | 2011-07-20 | 西安电子科技大学 | SAR (Synthetic Aperture Radar) image segmentation method based on dictionary learning and sparse representation |
CN104166856A (en) * | 2014-07-30 | 2014-11-26 | 西安电子科技大学 | Polarization SAR image classification method based on neighbor propagation clustering and region growing |
CN104915676A (en) * | 2015-05-19 | 2015-09-16 | 西安电子科技大学 | Deep-level feature learning and watershed-based synthetic aperture radar (SAR) image classification method |
KR101629163B1 (en) * | 2015-01-12 | 2016-06-13 | 한밭대학교 산학협력단 | Signal Filling Method Using Watershed Algorithm for MRC-based Image Compression |
CN108460400A (en) * | 2018-01-02 | 2018-08-28 | 南京师范大学 | A kind of hyperspectral image classification method of combination various features information |
-
2018
- 2018-11-07 CN CN201811322086.4A patent/CN109410223A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101510309A (en) * | 2009-03-30 | 2009-08-19 | 西安电子科技大学 | Segmentation method for improving water parting SAR image based on compound wavelet veins region merge |
CN102129573A (en) * | 2011-03-10 | 2011-07-20 | 西安电子科技大学 | SAR (Synthetic Aperture Radar) image segmentation method based on dictionary learning and sparse representation |
CN104166856A (en) * | 2014-07-30 | 2014-11-26 | 西安电子科技大学 | Polarization SAR image classification method based on neighbor propagation clustering and region growing |
KR101629163B1 (en) * | 2015-01-12 | 2016-06-13 | 한밭대학교 산학협력단 | Signal Filling Method Using Watershed Algorithm for MRC-based Image Compression |
CN104915676A (en) * | 2015-05-19 | 2015-09-16 | 西安电子科技大学 | Deep-level feature learning and watershed-based synthetic aperture radar (SAR) image classification method |
CN108460400A (en) * | 2018-01-02 | 2018-08-28 | 南京师范大学 | A kind of hyperspectral image classification method of combination various features information |
Non-Patent Citations (2)
Title |
---|
郝阳阳: "基于字典学习的SAR图像分割", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 * |
马秀丽 等: "基于分水岭-谱聚类的SAR图像分割", 《红外与毫米波学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101520894B (en) | Method for extracting significant object based on region significance | |
CN105005765B (en) | A kind of facial expression recognizing method based on Gabor wavelet and gray level co-occurrence matrixes | |
CN104834922B (en) | Gesture identification method based on hybrid neural networks | |
CN105825183B (en) | Facial expression recognizing method based on partial occlusion image | |
Sutthiwichaiporn et al. | Adaptive boosted spectral filtering for progressive fingerprint enhancement | |
CN108122008A (en) | SAR image recognition methods based on rarefaction representation and multiple features decision level fusion | |
CN106815819B (en) | More strategy grain worm visible detection methods | |
Guedri et al. | Indexing and images retrieval by content | |
CN103778435A (en) | Pedestrian fast detection method based on videos | |
CN103268485A (en) | Sparse-regularization-based face recognition method capable of realizing multiband face image information fusion | |
CN105718866A (en) | Visual target detection and identification method | |
CN103336942A (en) | Traditional Chinese painting identification method based on Radon BEMD (bidimensional empirical mode decomposition) transformation | |
CN103279744A (en) | Multi-scale tri-mode texture feature-based method and system for detecting counterfeit fingerprints | |
CN110533100A (en) | A method of CME detection and tracking is carried out based on machine learning | |
Kosarevych et al. | Image segmentation based on the evaluation of the tendency of image elements to form clusters with the help of point field characteristics | |
Chen et al. | Segmentation of overlapping cervical cells with mask region convolutional neural network | |
CN106844785A (en) | A kind of CBIR method based on conspicuousness segmentation | |
CN109523484B (en) | Fractal feature-based finger vein network repair method | |
CN105930788A (en) | Non-downsampling contour wave and PCA (principal component analysis) combining human face recognition method | |
CN105528791B (en) | A kind of quality evaluation device and its evaluation method towards touch screen hand-drawing image | |
CN109460768B (en) | Text detection and removal method for histopathology microscopic image | |
Kaur et al. | Illumination invariant face recognition | |
CN109410223A (en) | A kind of SAR image segmentation method based on watershed algorithm and dictionary learning | |
CN108898052A (en) | The detection method and equipment of man-made features in remote sensing images | |
Liu et al. | A novel image segmentation algorithm based on visual saliency detection and integrated feature extraction |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190301 |