CN105654111B - Pathological image swept-focus classification method based on k-means cluster - Google Patents
Pathological image swept-focus classification method based on k-means cluster Download PDFInfo
- Publication number
- CN105654111B CN105654111B CN201510970704.6A CN201510970704A CN105654111B CN 105654111 B CN105654111 B CN 105654111B CN 201510970704 A CN201510970704 A CN 201510970704A CN 105654111 B CN105654111 B CN 105654111B
- Authority
- CN
- China
- Prior art keywords
- center
- focus
- mass point
- distance
- class
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000001575 pathological effect Effects 0.000 title claims abstract description 22
- 238000007405 data analysis Methods 0.000 claims abstract description 8
- 239000011521 glass Substances 0.000 description 2
- 238000000386 microscopy Methods 0.000 description 1
- 238000010827 pathological analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2413—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on distances to training or reference patterns
- G06F18/24133—Distances to prototypes
- G06F18/24137—Distances to cluster centroïds
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Image Analysis (AREA)
Abstract
The invention discloses the pathological image swept-focus classification method clustered based on k-means, this method provides a pathological section, includes the following steps: that S1. obtains K focus from the pathological section；S2. the K focuses are obtained into N number of center of mass point using iteration self-organizing data analysis method one by one；S3. classified according to N number of center of mass point of acquisition to the K focuses, to divide class belonging to each focus, each center of mass point is corresponding a kind of；S4. the position of the center of mass point of every class is adjusted according to the corresponding focus of every class after division, is terminated.The present invention sorts to classification of focal spot to be calculated using iteration self-organizing data analysis method, effectively accelerates focusing speed, improves focusing precision, improves the applicability of focal plane.
Description
Technical field
The present invention relates to technical field of microscopy more particularly to a kind of pathological image swept-focus based on k-means cluster
Classification method.
Background technique
Digital slices scanner is the Medical Devices quickly grown in recent years, and digital slices scanning system can be complete by glass slide
Information, comprehensive quick scanning make the glass slide of conventional matter become new-generation digital pathological section, digital slices scanning
Instrument is to realize epoch-making change to pathological diagnosis technology.A digital slices scanner essential ring during the scanning process
Section is exactly to focus, i.e., is cooperated by moving up and down for Z axis with camera grabgraf, finds the slice a certain number of focuses of different location
Afterwards, the focal plane for calculating slice, for preparing for slice scanning.There are mainly two types of existing focusing methods, and one is one
Secondary property is to complete all focuses, but the error focused in this way is big, and scan image is fuzzy；Another method is to complete one
Part focus scanning a part, i.e. dynamic focusing, this method there occur a problem that when being sliced larger or tissue be island
When, the focal plane calculated for the first time may be not suitable for farther away focus, cause farther away focus focusing difficult.
Summary of the invention
The above problem existing for focusing method for existing digital slices scanner, now providing that one kind aims at can
It calculates the first time focal plane of dynamic focusing and is suitable for all focuses to be calculated, accelerate focusing speed, improve focusing precision
Pathological image swept-focus classification method based on k-means cluster.
Specific technical solution is as follows:
Based on the pathological image swept-focus classification method of k-means cluster, a pathological section, including following steps are provided
It is rapid:
S1. K focus is obtained from the pathological section；
S2. the K focuses are obtained into N number of center of mass point using iteration self-organizing data analysis method one by one；
S3. classified according to N number of center of mass point of acquisition to the K focuses, to divide each focus institute
The class of category, each center of mass point are corresponding a kind of；
S4. the position of the center of mass point of every class is adjusted according to the corresponding focus of every class after division, is terminated.
Preferably, the detailed process of the step S2 are as follows:
S21. choose a focus as the i-th=1 center of mass point, obtain respectively the K-1 focuses with it is described
The distance D of i-th of center of mass pointi(j), and by K-1 distance Di(j) it adds up to obtain total distance sum (Di),
Wherein, i, j are positive integer, Di(j) the distance between j-th of focus and i-th of center of mass point are indicated,
1≤i≤K,1≤j≤K,i≠j；
S22. i+1 center of mass point is obtained according to formula (1):
Random-Di(j) < 0 (1),
Wherein, Random is less than or equal to sum (Di) integer, when Random meets formula (1), the distance Di
(j) corresponding j-th focus is i-th=i+1 center of mass point；
S23. judge whether i is equal to N, if so, executing step S24；If it is not, executing step S3；
S24. the distance D of K-1 the focuses and described i-th center of mass point is obtained respectivelyi(j), and by K-1
A distance Di(j) it adds up to obtain the total distance sum (Di), it returns, executes step S22.
Preferably, the distance that each focus arrives each center of mass point respectively is obtained one by one in the step S3,
The focus is classified as shortest apart from the corresponding class of the corresponding center of mass point in corresponding distance.
Preferably, each focus is obtained according to formula (2) in the step S3 arrive each center of mass point respectively
Distance:
Cj=argmin ‖ Xj-Yi‖2(2),
Wherein, XjIndicate the position coordinates of j-th of focus, YiIndicate the position coordinates of i-th of center of mass point, Cj
Indicate described j-th focus to i-th of center of mass point distance.
Preferably, being averaged for the position coordinates for all focuses for belonging to every class is obtained one by one in the step S4
Coordinate, using the average coordinates as the position coordinates of the center of mass point corresponding with the class.
Preferably, it is sat in the step S4 according to the position that formula (3) obtain all focuses for belonging to every class
Target average coordinates:
Wherein, M indicates all focuses belonged in the corresponding class of i-th of center of mass point, QiIndicate i-th of institute
State the average coordinates of the center of mass point of the corresponding class of center of mass point.
Above-mentioned technical proposal the utility model has the advantages that
In the technical scheme, it is sorted using iteration self-organizing data analysis method to classification of focal spot to be calculated, it is effective to add
Fast focusing speed, improves focusing precision, improves the applicability of focal plane.
Detailed description of the invention
Fig. 1 is a kind of implementation of the pathological image swept-focus classification method of the present invention based on k-means cluster
The flow chart of example.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art without creative labor it is obtained it is all its
His embodiment, shall fall within the protection scope of the present invention.
It should be noted that in the absence of conflict, the feature in embodiment and embodiment in the present invention can phase
Mutually combination.
The present invention will be further explained below with reference to the attached drawings and specific examples, but not as the limitation of the invention.
As shown in Figure 1, the pathological image swept-focus classification method based on k-means cluster, provides a pathological section, wrap
Include following step:
S1. K focus is obtained from pathological section；
S2. K focus is obtained into N number of center of mass point using iteration self-organizing data analysis method one by one；
S3. classified according to N number of center of mass point of acquisition to K focus, to divide class belonging to each focus, Mei Gezhi
Heart point is corresponding a kind of；
S4. the position of the center of mass point of every class is adjusted according to the corresponding focus of every class after division, is terminated.
In the present embodiment, it is sorted using iteration self-organizing data analysis method to classification of focal spot to be calculated, it is effective to accelerate
Focusing speed improves focusing precision, improves the applicability of focal plane.
In a preferred embodiment, the detailed process of step S2 are as follows:
S21. a focus is chosen as the i-th=1 center of mass point, obtains K-1 focus respectively at a distance from i-th of center of mass point
Di(j), and by K-1 distance Di(j) it adds up to obtain total distance sum (Di),
Wherein, i, j are positive integer, Di(j) expression the distance between j-th of focus and i-th of center of mass point, 1≤i≤K,
1≤j≤K,i≠j；
S22. i+1 center of mass point is obtained according to formula (1):
Random-Di(j) < 0 (1),
Wherein, Random is less than or equal to sum (Di) integer, when Random meets formula (1), distance Di(j) right
J-th of the focus answered is i-th=i+1 center of mass point；
S23. judge whether i is equal to N, if so, executing step S24；If it is not, executing step S3；
S24. K-1 focus and i-th of center of mass point distance D are obtained respectivelyi(j), and by K-1 distance Di(j) add up with
Obtain total distance sum (Di), it returns, executes step S22.
In the present embodiment, the selection number of center of mass point and position will have a direct impact on the quality of cluster result, therefore use
Iteration self-organizing data analysis method algorithm calculates center of mass point number, and focusing precision can be improved.
In a preferred embodiment, the distance that each focus arrives each center of mass point respectively is obtained one by one in step s3, it will
Focus is classified as shortest apart from the corresponding class of corresponding center of mass point in corresponding distance.
Further, the distance that each focus arrives each center of mass point respectively is obtained according to formula (2) in step s3:
Cj=argmin ‖ Xj-Yi‖2(2),
Wherein, XjIndicate the position coordinates of j-th of focus, YiIndicate the position coordinates of i-th of center of mass point, CjIt indicates j-th
Distance of the focus to i-th of center of mass point.
In a preferred embodiment, in step s 4 one by one obtain belong to every class institute focal position coordinates put down
Equal coordinate, using average coordinates as the position coordinates of center of mass point corresponding with class.
Further, in step s 4 according to formula (3) obtain belong to every class institute focal position coordinates put down
Equal coordinate:
Wherein, M indicates all focuses belonged in the corresponding class of i-th of center of mass point, QiIndicate that i-th of center of mass point is corresponding
Class center of mass point average coordinates.
In the present embodiment, using average coordinates as the position coordinates of center of mass point corresponding with class, after completing focus cluster,
By the sequence to inhomogeneity focus, efficient focus identification can be realized.
The foregoing is merely preferred embodiments of the present invention, are not intended to limit embodiments of the present invention and protection model
It encloses, to those skilled in the art, should can appreciate that all with made by description of the invention and diagramatic content
Equivalent replacement and obviously change obtained scheme, should all be included within the scope of the present invention.
Claims (5)
1. the pathological image swept-focus classification method based on k-means cluster, which is characterized in that provide a pathological section, packet
Include following step:
S1. K focus is obtained from the pathological section；
S2. the K focuses are obtained into N number of center of mass point using iteration self-organizing data analysis method one by one；
S3. classified according to N number of center of mass point of acquisition to the K focuses, to divide belonging to each focus
Class, each center of mass point are corresponding a kind of；
S4. the position of the center of mass point of every class is adjusted according to the corresponding focus of every class after division, is terminated；
The detailed process of the step S2 are as follows:
S21. a focus is chosen as the i-th=1 center of mass point, obtains the K-1 focuses and described i-th respectively
The distance D of the center of mass pointi(j), and by K-1 distance Di(j) it adds up to obtain total distance sum (Di),
Wherein, i, j are positive integer, Di(j) the distance between j-th of focus and i-th of center of mass point, 1≤i are indicated
≤K,1≤j≤K,i≠j；
S22. i+1 center of mass point is obtained according to formula (1):
Random-Di(j) < 0 (1),
Wherein, Random is less than or equal to sum (Di) integer, when Random meets formula (1), the distance Di(j) right
J-th answered the focus is i-th=i+1 center of mass point；
S23. judge whether i is equal to N, if so, executing step S24；If it is not, executing step S3；
S24. the distance D of K-1 the focuses and described i-th center of mass point is obtained respectivelyi(j), and by K-1 institute
State distance Di(j) it adds up to obtain the total distance sum (Di), it returns, executes step S22.
2. the pathological image swept-focus classification method as described in claim 1 based on k-means cluster, which is characterized in that
It obtains the distance that each focus arrives each center of mass point respectively one by one in the step S3, the focus is classified as phase
It is shortest apart from the corresponding class of the corresponding center of mass point in the distance answered.
3. the pathological image swept-focus classification method as claimed in claim 2 based on k-means cluster, which is characterized in that
The distance that each focus arrives each center of mass point respectively is obtained according to formula (2) in the step S3:
Cj=argmin | | Xj-Yi||2(2),
Wherein, XjIndicate the position coordinates of j-th of focus, YiIndicate the position coordinates of i-th of center of mass point, CjIt indicates
Distance of j-th of the focus to i-th of center of mass point.
4. the pathological image swept-focus classification method as claimed in claim 3 based on k-means cluster, which is characterized in that
Obtain the average coordinates for belonging to the position coordinates of all focuses of every class one by one in the step S4, it will be described average
Position coordinates of the coordinate as the center of mass point corresponding with the class.
5. the pathological image swept-focus classification method as claimed in claim 4 based on k-means cluster, which is characterized in that
The average coordinates for belonging to the position coordinates of all focuses of every class are obtained according to formula (3) in the step S4:
Wherein, M indicates all focuses belonged in the corresponding class of i-th of center of mass point, QiIndicate i-th of the matter
The average coordinates of the center of mass point of the corresponding class of heart point, CjIndicate described j-th focus to i-th center of mass point
Distance.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510970704.6A CN105654111B (en) | 2015-12-21 | 2015-12-21 | Pathological image swept-focus classification method based on k-means cluster |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510970704.6A CN105654111B (en) | 2015-12-21 | 2015-12-21 | Pathological image swept-focus classification method based on k-means cluster |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105654111A CN105654111A (en) | 2016-06-08 |
| CN105654111B true CN105654111B (en) | 2019-03-29 |
Family
ID=56477539
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510970704.6A Active CN105654111B (en) | 2015-12-21 | 2015-12-21 | Pathological image swept-focus classification method based on k-means cluster |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105654111B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109374621A (en) * | 2018-11-07 | 2019-02-22 | 杭州迪英加科技有限公司 | Focusing method, system and the device of slice scanner |
| CN113866974A (en) * | 2020-06-29 | 2021-12-31 | 深圳辉煌耀强科技有限公司 | Focusing method, device and equipment applied to cervical cell image and storage medium |
| CN115511831B (en) * | 2022-09-27 | 2023-04-04 | 佳木斯大学 | Data analysis processing system and method for tissue embryo pathological section |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104390981A (en) * | 2014-12-12 | 2015-03-04 | 宁波江丰生物信息技术有限公司 | Method for rapidly and accurately focusing and scanning pathological section tissue on basis of image collection device |
| CN105118088A (en) * | 2015-08-06 | 2015-12-02 | 曲阜裕隆生物科技有限公司 | 3D imaging and fusion method based on pathological slice scanning device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070031056A1 (en) * | 2005-08-02 | 2007-02-08 | Perz Cynthia B | System for and method of focusing in automated microscope systems |
-
2015
- 2015-12-21 CN CN201510970704.6A patent/CN105654111B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104390981A (en) * | 2014-12-12 | 2015-03-04 | 宁波江丰生物信息技术有限公司 | Method for rapidly and accurately focusing and scanning pathological section tissue on basis of image collection device |
| CN105118088A (en) * | 2015-08-06 | 2015-12-02 | 曲阜裕隆生物科技有限公司 | 3D imaging and fusion method based on pathological slice scanning device |
Non-Patent Citations (2)
| Title |
|---|
| 基于分层聚类的k-means算法;金微等;《河海大学常州分校学报》;20070331;第21卷(第1期);第8页 |
| 杨坤等.基于视场边界最大最小聚类分析的快速扫描方法.《光 电 子 •激 光》.2014,第25卷(第4期), |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105654111A (en) | 2016-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103631928B (en) | LSH (Locality Sensitive Hashing)-based clustering and indexing method and LSH-based clustering and indexing system | |
| CN105654111B (en) | Pathological image swept-focus classification method based on k-means cluster | |
| Ferguson | Population studies in groups and clusters of galaxies. II-A catalog of galaxies in the central 3.5 deg of the Fornax Cluster | |
| CN103874200B (en) | A kind of floor recognition methods and system | |
| JP6461128B2 (en) | Cell evaluation apparatus, method and program | |
| CN104390981B (en) | The method of quick and precisely focusing scanning pathological section tissue based on image collecting device | |
| CN106535134B (en) | Multi-room positioning method based on wifi and server | |
| CN101923712A (en) | Particle swarm optimization-based gene chip image segmenting method of K-means clustering algorithm | |
| CN106851571A (en) | WiFi localization methods in a kind of quick KNN rooms based on decision tree | |
| RU2018133474A (en) | SYSTEM FOR FORMING A SYNTHESIZED TWO-DIMENSIONAL IMAGE OF A BIOLOGICAL SAMPLE WITH AN INCREASED DEPTH OF SHARPNESS | |
| JP2014533852A (en) | Method and server for searching for nearby users in a social network service | |
| US20210286972A1 (en) | Determination method, elimination method and apparatus for electron microscope aberration | |
| CN109088715A (en) | Pilot distribution method and device | |
| WO2013142176A4 (en) | Scoring and adjusting pixels based on neighborhood relationships for revealing data in images | |
| EP2765557A2 (en) | Coregistering images of needle biopsies using multiple weighted landmarks | |
| US9612313B2 (en) | Supporting coverage area modeling | |
| CN105044699B (en) | Radar plot centroid method based on Radon-Fourier transformation | |
| WO2019225505A1 (en) | Biological tissue image processing system, and machine learning method | |
| EP3092614A1 (en) | Method and system for image correction using a quasiperiodic grid | |
| JP2014527640A (en) | Method for automatically adjusting the focal plane of a digital pathological image | |
| CN108459406A (en) | Microscope auto-focusing window selection method based on artificial fish-swarm algorithm | |
| JP2016087433A (en) | Medical image processing apparatus and medical image processing method | |
| CN101980251A (en) | Remote sensing classification method for binary tree multi-category support vector machines | |
| CN109567839A (en) | Hip joint x-ray image automatic analysis method | |
| CN107316306B (en) | A kind of diameter radar image fast partition method based on Markov model |
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 | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Scanning focus classification of pathological images based on K-means clustering Effective date of registration: 20210302 Granted publication date: 20190329 Pledgee: China Co. truction Bank Corp Yuyao branch Pledgor: KONFOONG BIOTECH INTERNATIONAL Co.,Ltd. Registration number: Y2021330000164 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Granted publication date: 20190329 Pledgee: China Co. truction Bank Corp Yuyao branch Pledgor: KONFOONG BIOTECH INTERNATIONAL Co.,Ltd. Registration number: Y2021330000164 |