CN104240190B - The processing method of headgear is automatically removed based on incidence CT images - Google Patents
The processing method of headgear is automatically removed based on incidence CT images Download PDFInfo
- Publication number
- CN104240190B CN104240190B CN201310242840.4A CN201310242840A CN104240190B CN 104240190 B CN104240190 B CN 104240190B CN 201310242840 A CN201310242840 A CN 201310242840A CN 104240190 B CN104240190 B CN 104240190B
- Authority
- CN
- China
- Prior art keywords
- headgear
- pixel
- image
- incidence
- processing method
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 14
- 238000006243 chemical reaction Methods 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 5
- 235000013399 edible fruits Nutrition 0.000 claims 2
- 238000003745 diagnosis Methods 0.000 abstract description 8
- 238000012800 visualization Methods 0.000 abstract description 8
- 238000011161 development Methods 0.000 abstract description 2
- 210000003128 head Anatomy 0.000 description 26
- 238000010586 diagram Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 241000272525 Anas platyrhynchos Species 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 208000028979 Skull fracture Diseases 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 206010020718 hyperplasia Diseases 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 210000001944 turbinate Anatomy 0.000 description 1
- 238000007794 visualization technique Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The present invention disclose it is a kind of the processing method of headgear is automatically removed based on incidence CT images, including carry out binary conversion treatment, obtain binary image sequence；Traversal binary image sequence, finds out largest connected domain, finds out layer where shoulder；Calculate maximal projection stop-layer；CT image sequences are carried out into binaryzation to maximal projection stop-layer MS from the crown according to default high threshold, projection superposition, generation projection masterplate is then carried out；Connected domain judgement is carried out to projection masterplate, the maximum connected domain of area is found out, as an initial head template image, initial head template image is extended up with maximum head circumference, obtain final head template image；Layer where the overhead of pending CT image sequences to shoulder is mapped with head template image, the CT images after removal headgear are obtained.The present invention solves the problems, such as headgear interference in three-dimensional visualization, facilitates diagnosis, has important meaning to the development of the three-dimensional visualization of medical science.
Description
Technical field
Head is automatically removed based on incidence CT images the invention belongs to technical field of medical image processing, more particularly to one kind
The processing method of set.
Background technology
In medical diagnosis, medical worker carries out three-dimensional in the brain by observing multigroup two dimension (2D) faultage image
The reconstruction of (3D) data, the space structure of lesion body is determined with this.This is just difficult to accurately to determine the locus of lesion body, big
Small, geometry and the relation between surrounding biological tissue.So three-dimensional visualization technique is for clinical medical essence
Make a definite diagnosis disconnected and surgical planning and play more and more important effect.
But simple three-dimensional visible technology can not meet the demand of medical personnel.The CT of incidence is such as in patient
When, in order to prevent head from moving, all can be to patient with individual headgear.So the image of headgear can be all carried in the CT slice, thin pieces taken,
If the image in two-dimentional tomography, these images with headgear do not interfere with diagnosis.Can be in three-dimensional visualization, headgear will
The rear of incidence can be intuitively shown, if patient has the disease that cervical vertebra hyperplasia or skull fracture etc. are blocked or disturbed by headgear
Become, doctor will be difficult to position area-of-interest.
Usual diagnostician solves the mode disturbed by headgear, is to allow medical personnel by analog scalpel come slowly labour contractor
Set hand cut removal, can so waste great effort.
The content of the invention
The technical problems to be solved by the invention are to propose a kind for the treatment of that headgear is automatically removed based on incidence CT images
Method, solves the interference that headgear in three-dimensional visualization brings, and facilitates diagnosis.
The present invention adopts the following technical scheme that realization：A kind for the treatment of side that headgear is automatically removed based on incidence CT images
Method, it includes step：
A, pending CT image sequences are carried out into binary conversion treatment by default Low threshold, obtain binary image sequence
Row；
B, traversal binary image sequence, find out largest connected domain, calculate largest connected domain two ends from image boundary away from
From, it is considered shoulder if distance is less than predeterminable range, obtain layer where shoulder；
C, maximal projection stop-layer MS is calculated using formula MS=(Size+2*S)/3, wherein Size is CT image sequences
The number of plies, S be shoulder where layer；
D, CT image sequences are carried out into binaryzation to maximal projection stop-layer MS from the crown according to default high threshold, then
Carry out projection superposition, generation projection masterplate；
E, connected domain judgement is carried out to projection masterplate, the maximum connected domain of area is found out, as an initial head masterplate
Image, is extended up initial head template image with maximum head circumference, obtains final head template image；
F, the overhead of pending CT image sequences is reflected to layer where shoulder with final head template image
Penetrate, obtain the CT images after removal headgear.
Wherein, the step of step A also includes binaryzation holes filling.
Wherein, the step of binaryzation holes filling is specifically included：
After by the treatment of binary image inverse, traversal is proceeded by from a certain pixel, when running into first picture of non-zero value
Vegetarian refreshments is off traversal；
Seed algorithm is proceeded by from the pixel of this non-zero value, all non-zero points being attached thereto are arranged to one
Block, forms a connected region；
The value of all pixels point in this connected region is set as 0 and inverse treatment.
Wherein, step A also disconnects the connection of headgear belt and head in CT images using opening operation.
Wherein, default Low threshold is 700HU.
Wherein, predeterminable range is set to 2 distances of pixel.
Wherein, default high threshold is 1224HU.
Wherein, the step of finding out largest connected domain in step B includes：
Four field flags are carried out to binary image；
The area of each mark is calculated, it is exactly largest connected domain to be eventually found the maximum domain of area.
Wherein, the step of carrying out four field flags to binary image includes：
From left to right, binary image is scanned from top to bottom, judge in each field of pixel four it is left, on either with or without
Point, is considered the beginning of a new region if not.；
If most left in this field of pixel four has a pixel and upper without pixel, this pixel is marked for most left
Value, if most left no pixel in this field of pixel four and on have pixel, mark this pixel be top point value；
If most left in four fields of this pixel and most go up pixel, it is this 2 pixels to mark this pixel
Minimum mark point in point, and change big mark and be.
Compared with prior art, the present invention has the advantages that：
The present invention proposes to automatically remove the implementation of headgear based on incidence CT images, combines the powerful meter of computer
Calculation ability, solves the problems, such as headgear interference in three-dimensional visualization, diagnosis is facilitated, to the hair of the three-dimensional visualization of medical science
Exhibition has important meaning.
Brief description of the drawings
Fig. 1 is that a preferred embodiment of the present invention realizes schematic flow sheet.
Fig. 2A and Fig. 2 B are respectively to carry out the headgear portion image schematic diagram before and after opening operation.
Fig. 3 is the schematic diagram of initial head template image.
Fig. 4 is the schematic diagram of the head template image of pre-production.
The schematic diagram of the CT images being respectively shown in Fig. 5 A and Fig. 5 B after pending CT images and removal headgear.
Specific embodiment
The present invention proposes a kind of processing method that headgear is automatically removed based on incidence CT images, is by pre-production head
Portion's template image, using pre-production head template image to incidence CT images remove headgear, so as to solve three-dimensional visible
The interference that headgear brings in change, facilitates diagnosis.
As shown in figure 1, in a preferred embodiment of the present invention, step S1~S6 is pre-production head template image
Flow is realized, step S7 is the handling process that headgear is automatically removed to CT images.Specifically, the present invention includes step is implemented as described below
Suddenly：
Step S1, pending CT images are pre-processed.
First, pending CT image sequences are carried out into binary conversion treatment by default Low threshold.Wherein, default low threshold
Value is taken as 700HU values, and (HU is CT units, can describe the brightness of image, and 700 are taken herein is counted by mass data
Arrive, this threshold value can be maximum comprising all of information of head skeleton).Pixel in CT image sequences less than 700H is considered
Background, all sets to 0；Pixel higher than 700H is considered, comprising all of information of head skeleton, to be all set to 1, that is, complete two-value
Change is processed, and obtains binary image.
Secondly, some low-density tissues can form cavity in binary image, it is necessary to carry out binaryzation holes filling.Cavity
It is exactly the pixel but with overall background UNICOM without overall background color.Here be that se ed filling algorithm is filled, it is false first
Certain point determined in occluding contour is known, then starts search point adjacent with seed point and in contour line, if phase
Adjoint point is not in contour line, then just reach the border of contour line, and this turns into new with you if consecutive points are located within contour line
Seed point, then proceed to search and go down, finally reach requirement.
Specifically, binaryzation holes filling is to be achieved by the steps of：After by the treatment of binary image inverse, from certain
One pixel proceeds by traversal, when the pixel for running into first non-zero value is off traversal；From the pixel of this non-zero value
Start, carry out seed algorithm, all non-zero points being attached thereto are arranged to a block, form a connected region；This is connected
The value of all pixels point deletes (being set as 0) and inverse is processed in logical region.
Finally, if the belt density ratio of headgear is larger, then headgear and head just have the situation of adhesion.In order to solve
This problem has used the masterplate that radius is 2 and use opening operation here, and such as Fig. 2A and Fig. 2 B are respectively before and after carrying out opening operation
Headgear portion image schematic diagram.
Opening operation belongs to morphological images treatment, is first to corrode to expand afterwards.The effect of opening operation is to disconnect narrow connection,
And keep size constant.The formula of opening operation is：
In above-mentioned opening operation formula, X is image, and B is masterplate.OPEN (X, B) and XBIt is the different representations of opening operation,
O in X o B is the symbol of opening operation.It is the symbol of Image erosion.It is the symbol of image expansion.
Wherein, masterplate is the circular masterplate (can regard an image for binaryzation as) that radius is 2, and it is expressed as：
Wherein B represents the circular masterplate (can regard the figure of binaryzation as) that radius is as 2, and 0 is to represent the value that pixel is as 0.1
It is to represent the value that pixel is as 1.
Step S2, each pretreated binary image sequence is traveled through, find out largest connected domain, calculate largest connected domain
Two ends are considered shoulder with a distance from image boundary if distance is less than predeterminable range, obtain layer S where shoulder.Here
Predeterminable range is set to 2 pixels, and the position of shoulder is then considered if less than 2 pixels.
Wherein, the method for finding out largest connected domain is as follows：First have to carry out four field flags to binary image, mark
Process includes (1) (2) and (3).(1) from left to right, scan image from top to bottom, judge the left side in each four field of point, on have
Do not have a little, the beginning of a new region is considered if not.(2) if in this four field of point it is most left a little, it is upper without point,
This point is then marked for most lvalue, if this four field of point is most left without point, on a little, then it is the value of top point to mark this point.
(3) if most left in four fields of this point, most go up all a little, it is the minimum mark point in this 2 to mark this point, and is changed
Mark is greatly.Made marks to binary image by (1) (2) and (3) these three processes, then calculate each mark
Area, it is exactly largest connected domain to be eventually found the maximum domain of area.
Step S3, maximal projection stop-layer MS is obtained by function.Empirical function is：MS=(Size+2*S)/3, wherein,
MS is maximal projection stop-layer, and Size is the number of plies of CT image sequences, and S is layer where shoulder.This function is by a large amount of
The rule that CT view data finds after being counted.
Step S4, CT image sequences are carried out high threshold binaryzation from the crown to maximal projection stop-layer MS, then carried out
Projection superposition, generation projection masterplate.That high threshold is selected here is 1224HU, the characteristics of this threshold value is headgear and head is certainly
It is kept completely separate, as multilevel projection is superimposed, head zone is increasingly more complete.The value of pixel in CT image sequences is more than
1224HU's puts 1, otherwise sets to 0.
Step S5, connected domain judgement is carried out to projection masterplate, find out the maximum connected domain of area, duck eye filling, as
Individual initial head template image, as shown in Figure 3.
Wherein, the implementation method for finding out the maximum connected domain of area is identical with the method in largest connected domain is looked in step S2,
No longer chat again.The implementation method of duck eye filling is identical the step of binaryzation holes filling with step S1, no longer chats again.
Step S6, in order to prevent head template image middle nasal concha lack, initial head template image with maximum head circumference to
Upper extension, obtains final head template image, as shown in Figure 4.
Step S7, the overhead the 1st layer of CT image sequences (overhead for) of pending CT image sequences are to shoulder institute
Mapped that (mapping is meant that with final head template image in layer S：Each layer of CT image sequence all with head masterplate figure
As carrying out and computing, the white portion in head template image retains, and the black region removal in head template image) and shoulder
Layer S where wing to sequence bottom layer then without modification, you can obtain the three-dimensional visualization effect of incidence CT images removal headgear.
By above step, the purpose of removal headgear is finally reached.Pending CT is respectively as fig. 5 a and fig. 5b
The schematic diagram of the CT images after image and removal headgear.
To sum up, the present invention proposes to automatically remove the implementation of headgear based on incidence CT images, combines computer strong
Big computing capability, solves the problems, such as headgear interference in three-dimensional visualization, diagnosis is facilitated, to the three-dimensional visible of medical science
The development of change has important meaning.
Presently preferred embodiments of the present invention is the foregoing is only, is not intended to limit the invention, it is all in essence of the invention
Any modification, equivalent and improvement made within god and principle etc., should be included within the scope of the present invention.
Claims (9)
1. a kind of processing method that headgear is automatically removed based on incidence CT images, it is characterised in that methods described includes step：
A, pending CT image sequences are carried out into binary conversion treatment by default Low threshold, obtain binary image sequence；
B, traversal binary image sequence, find out largest connected domain, calculate largest connected domain two ends with a distance from image boundary, such as
Fruit distance is then considered shoulder less than predeterminable range, obtains layer where shoulder；
C, maximal projection stop-layer MS is calculated using formula MS=(Size+2*S)/3, wherein Size is the layer of CT image sequences
Number, S is layer where shoulder；
D, CT image sequences are carried out into binaryzation to maximal projection stop-layer MS from the crown according to default high threshold, then carried out
Projection superposition, generation projection masterplate；
E, connected domain judgement is carried out to projection masterplate, find out the maximum connected domain of area, as an initial head template image,
Initial head template image is extended with maximum head circumference to orientation where nose, final head template image is obtained；
F, the overhead of pending CT image sequences is mapped to layer where shoulder with final head template image, obtained
The CT images after headgear must be removed.
2. the processing method of headgear is automatically removed based on incidence CT images according to claim 1, it is characterised in that step
The step of A also includes binaryzation holes filling.
3. the processing method of headgear is automatically removed based on incidence CT images according to claim 2, it is characterised in that described
The step of binaryzation holes filling, specifically includes：
After by the treatment of binary image inverse, traversal is proceeded by from a certain pixel, when running into first pixel of non-zero value
When stop traversal；
Seed algorithm is proceeded by from the pixel of this non-zero value, all non-zero points being attached thereto are arranged to a block, shape
Into a connected region；
The value of all pixels point in this connected region is set as 0 and inverse treatment.
4. the processing method of headgear is automatically removed based on incidence CT images according to claim 1, it is characterised in that step
A also disconnects the connection of headgear belt and head in CT images using opening operation.
5. the processing method of headgear is automatically removed based on incidence CT images according to claim 1, it is characterised in that default
Low threshold be 700HU.
6. the processing method of headgear is automatically removed based on incidence CT images according to claim 1, it is characterised in that default
Distance is set to 2 distances of pixel.
7. the processing method of headgear is automatically removed based on incidence CT images according to claim 1, it is characterised in that default
High threshold be 1224HU.
8. the processing method of headgear is automatically removed based on incidence CT images according to claim 1, it is characterised in that step
The step of largest connected domain is found out in B includes：
Four field flags are carried out to binary image；
The area of each mark is calculated, it is exactly largest connected domain to be eventually found the maximum domain of area.
9. the processing method of headgear is automatically removed based on incidence CT images according to claim 8, it is characterised in that to two
The step of value image carries out four field flags includes：
From left to right, binary image is scanned from top to bottom, judge in each field of pixel four it is left, on either with or without point, such as
Fruit is not considered the beginning of a new region then；
If most left in this field of pixel four has a pixel and upper without pixel, this pixel is marked for most lvalue,
If most left no pixel in this field of pixel four and on have pixel, mark this pixel be top point value；
If most left in four fields of this pixel and most go up pixel, this pixel is marked in this 2 pixels
Minimum mark point, and change big mark and be.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310242840.4A CN104240190B (en) | 2013-06-18 | 2013-06-18 | The processing method of headgear is automatically removed based on incidence CT images |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310242840.4A CN104240190B (en) | 2013-06-18 | 2013-06-18 | The processing method of headgear is automatically removed based on incidence CT images |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104240190A CN104240190A (en) | 2014-12-24 |
| CN104240190B true CN104240190B (en) | 2017-06-16 |
Family
ID=52228199
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310242840.4A Active CN104240190B (en) | 2013-06-18 | 2013-06-18 | The processing method of headgear is automatically removed based on incidence CT images |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104240190B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114972566B (en) * | 2022-05-27 | 2025-05-09 | 西安盈谷网络科技有限公司 | A method, device and storage medium for removing bed board in CT sequence images |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202446104U (en) * | 2012-02-23 | 2012-09-26 | 徐春林 | CT (Computed Tomography) examination fixing head cover |
| CN103020968A (en) * | 2012-12-21 | 2013-04-03 | 东软集团股份有限公司 | Head and neck CTA (computed tomography angiography) image layering method and device |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7929660B2 (en) * | 2008-07-18 | 2011-04-19 | Xcounter Ab | Support structure for mammography |
| US8861891B2 (en) * | 2010-03-05 | 2014-10-14 | Siemens Aktiengesellschaft | Hierarchical atlas-based segmentation |
-
2013
- 2013-06-18 CN CN201310242840.4A patent/CN104240190B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202446104U (en) * | 2012-02-23 | 2012-09-26 | 徐春林 | CT (Computed Tomography) examination fixing head cover |
| CN103020968A (en) * | 2012-12-21 | 2013-04-03 | 东软集团股份有限公司 | Head and neck CTA (computed tomography angiography) image layering method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104240190A (en) | 2014-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Saha et al. | Digital topology and geometry in medical imaging: a survey | |
| CN104462650B (en) | A kind of hypostazation heart 3D model production methods of achievable external and internal compositionses | |
| CN103440665B (en) | Automatic segmentation method of knee joint cartilage image | |
| KR101120250B1 (en) | Vascular image extraction and labeling system and method | |
| CN101452577B (en) | Rib auto-demarcating method and device | |
| CN100573581C (en) | Semi-automatic segmentation method for lung CT image lesions | |
| CN106355582B (en) | A method for liver medical image segmentation based on shape prior | |
| CN107045721B (en) | Method and device for extracting pulmonary blood vessels from chest CT (computed tomography) image | |
| CN107808156A (en) | ROI Extraction Method | |
| WO2015161728A1 (en) | Three-dimensional model construction method and device, and image monitoring method and device | |
| CN103489198A (en) | Method for partitioning brainstem areas automatically from MR (magnetic resonance) sequence images | |
| CN103985122B (en) | The dirtiest extracting method of based on heart CT image | |
| CN107895364B (en) | A 3D reconstruction system for preoperative planning of virtual surgery | |
| JP2006198411A (en) | How to visualize myocardial damage | |
| CN102622750A (en) | Stomach computed tomography (CT) sequence image segmentation method based on interactive region growth | |
| CN107595390B (en) | A Real-time Matching Fusion Method of Ultrasound Image and CT Image | |
| CN107067398A (en) | Complementing method and device for lacking blood vessel in 3 D medical model | |
| CN107945169B (en) | A Coronary Artery Imaging Analysis Method | |
| CN103793905B (en) | Determine the method and device in brain image Coronal direction | |
| CN102982531A (en) | Bronchial partition method and system thereof | |
| CN103077557A (en) | Adaptive hierarchical chest large data display implementation method | |
| CN110310726B (en) | Liver segment resection auxiliary system based on medical image three-dimensional reconstruction and holographic display | |
| Gunz et al. | Skull reconstruction of the late Miocene ape Rudapithecus hungaricus from Rudabánya, Hungary | |
| CN108961278B (en) | Method and system for abdominal wall muscle segmentation based on image data | |
| Zhao et al. | Research on 3D reconstruction algorithm of medical CT image based on parallel contour |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: Nanshan District Xueyuan Road in Shenzhen city of Guangdong province 518000 No. 1001 Nanshan Chi Park A3 building 6 floor Applicant after: LANWON TECHNOLOGY CO., LTD. Address before: Nanshan District Xueyuan Road in Shenzhen city of Guangdong province 518000 No. 1001 Nanshan Chi Park A3 building 6 floor Applicant before: LANWON SCIENCE & TECHNOLOGY CO., LTD. |
|
| COR | Change of bibliographic data | ||
| GR01 | Patent grant |