CN109009226A - A kind of method of 3-D supersonic imaging - Google Patents
A kind of method of 3-D supersonic imaging Download PDFInfo
- Publication number
- CN109009226A CN109009226A CN201810827517.6A CN201810827517A CN109009226A CN 109009226 A CN109009226 A CN 109009226A CN 201810827517 A CN201810827517 A CN 201810827517A CN 109009226 A CN109009226 A CN 109009226A
- Authority
- CN
- China
- Prior art keywords
- dimensional data
- smoothing processing
- median sagittal
- plane
- voxel
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0866—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving foetal diagnosis; pre-natal or peri-natal diagnosis of the baby
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Gynecology & Obstetrics (AREA)
- Pregnancy & Childbirth (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The present invention provides a kind of methods of 3-D supersonic imaging, comprising: emits ultrasonic wave to fetus head；Ultrasonic echo is received, ultrasound echo signal is obtained；The three-dimensional data of fetus head is obtained according to the ultrasound echo signal；Multiple smoothing processing is carried out to three-dimensional data；According to the feature of fetus head median sagittal section, median sagittal section is detected from the three-dimensional data after smoothing processing；Show the median sagittal section.The present invention can carry out the three-dimensional data that ultrasonic scanning obtains fetus head to fetus, and the three-dimensional data of acquisition is smoothed, the median sagittal section of Fetal Brain is detected automatically according to the three-dimensional data after smoothing processing and is shown, solves the problems, such as that doctor is difficult to that median sagittal section is accurately positioned manually.
Description
Technical field
The present invention relates to medical ultrasound imaging technical fields, and in particular to a kind of method of 3-D supersonic imaging.
Background technique
Ultrasonic instrument is generally used for the internal organizational structure that doctor observes human body, and ultrasonic probe is placed on human body by doctor
Corresponding skin surface, the ultrasound image at the available position.Ultrasound due to it conveniently, safely, it is lossless, cheap the features such as,
Have become one of the main supplementary means of diagnosis.
In fetal central nervous system inspection, corpus callosum and vermis of cerebellum are two critically important inspection items, wherein callosity
Body is maximum commissural fibre in cerebral hemisphere, be responsible for brain hemispheres between communication, missing or depauperation will lead to epilepsy,
A series of complication such as feeblemindedness, dyskinesia.Vermian missing or depauperation are Dandy-walker
The performance of syndrome, 50% Dandy-walker patient has psychomotor development sluggish and feeblemindedness, and is often accompanied by dyeing
Body is abnormal and other deformities, poor prognosis, the death rate are high.
In recent years, the extensive use with three-D ultrasonic clinically, part doctor by using Double Tops radial longitudinal section as rise
Beginning plane carries out three-dimensional scanning to fetus, and then passing through the geometric transformations such as can be manually rotated, translate, the third plane in 3D ultrasound
Median sagittal plane is recalled, corpus callosum and vermis of cerebellum are checked under the section.As can be seen that using from various open selected works
Although the picture quality for the sagittal plane that this method obtains is slightly poorer than two dimensional image, corpus callosum and vermian display rate are non-
Chang Gao can quickly and accurately judge whether corpus callosum and vermis of cerebellum are abnormal by this method.However, doctor needs to three
Dimension space has very deep understanding, can adjust out median sagittal plane by can be manually rotated, translating geometric operation at 3D.
Summary of the invention
In view of the above-mentioned problems, the present invention provides a kind of method of 3-D supersonic imaging.
The purpose of the present invention is realized using following technical scheme:
A kind of method of 3-D supersonic imaging is provided, this method comprises:
Emit ultrasonic wave to fetus head；
Ultrasonic echo is received, ultrasound echo signal is obtained；
The three-dimensional data of fetus head is obtained according to the ultrasound echo signal；
Multiple smoothing processing is carried out to three-dimensional data；
According to the feature of fetus head median sagittal section, median sagittal is detected from the three-dimensional data after smoothing processing
Section；
Show the median sagittal section.
Preferably, according to the gray feature of fetus head median sagittal section, from the three-dimensional data after smoothing processing
Detect the median sagittal section.
Preferably, the median sagittal section is detected from the three-dimensional data after smoothing processing, comprising:
It extracts to represent in three-dimensional data after smoothing processing and meets the ash that the gray value in plane is greater than plane two sides
The sagittal plane characteristic area of the plane of the condition of angle value；
At least three characteristic points are selected from the sagittal plane characteristic area；
The plane where the median sagittal section is determined according at least three feature points.
In another preferred scheme, the median sagittal section is detected from the three-dimensional data after smoothing processing,
Include:
At least two sections are extracted in the three-dimensional data；
It is greater than the feature of the gray value of brain middle line two sides according to the gray value of brain middle line, is mentioned at least two section
Brain middle line is taken, at least two brain middle lines are obtained；
The plane where the median sagittal plane is determined according to the plane that at least two brain middle lines limit.
The invention has the benefit that the three-dimensional data that ultrasonic scanning obtains fetus head can be carried out to fetus, and
The three-dimensional data of acquisition is smoothed, Fetal Brain is being detected just according to the three-dimensional data after smoothing processing automatically
Middle sagittal section is simultaneously shown solve the problems, such as that doctor is difficult to that median sagittal section is accurately positioned manually, so that doctor can
The case where easily to observe Fetal Brain median sagittal section, can provide a large amount of important key messages for doctor.
Detailed description of the invention
The present invention will be further described with reference to the accompanying drawings, but the embodiment in attached drawing is not constituted to any limit of the invention
System, for those of ordinary skill in the art, without creative efforts, can also obtain according to the following drawings
Other attached drawings.
Fig. 1 is the flow chart of the 3-D supersonic imaging method of an illustrative embodiment of the invention；
Fig. 2 is the method flow that multiple smoothing processing is carried out to three-dimensional data of an illustrative embodiment of the invention
Figure.
Specific embodiment
The invention will be further described with the following Examples.
Referring to Fig. 1, the embodiment of the invention provides a kind of methods of 3-D supersonic imaging, this method comprises:
S01 emits ultrasonic wave to fetus head.
S02 receives ultrasonic echo, obtains ultrasound echo signal.
S03 obtains the three-dimensional data of fetus head according to the ultrasound echo signal.
S04 carries out multiple smoothing processing to three-dimensional data.
S05 detects center according to the feature of fetus head median sagittal section from the three-dimensional data after smoothing processing
Sagittal section.In a kind of mode that can be realized, according to the gray feature of fetus head median sagittal section, from smoothing processing
The median sagittal section is detected in three-dimensional data afterwards.
In one embodiment, the median sagittal section is detected from the three-dimensional data after smoothing processing, comprising:
It extracts to represent in three-dimensional data after smoothing processing and meets the ash that the gray value in plane is greater than plane two sides
The sagittal plane characteristic area of the plane of the condition of angle value；
At least three characteristic points are selected from the sagittal plane characteristic area；
The plane where the median sagittal section is determined according at least three feature points.
In another embodiment, the median sagittal section is detected from the three-dimensional data after smoothing processing, comprising:
At least two sections are extracted in the three-dimensional data；
It is greater than the feature of the gray value of brain middle line two sides according to the gray value of brain middle line, is mentioned at least two section
Brain middle line is taken, at least two brain middle lines are obtained；
The plane where the median sagittal plane is determined according to the plane that at least two brain middle lines limit.
The method of the 3-D supersonic imaging of the above embodiment of the present invention can carry out ultrasonic scanning to fetus and obtain fetal head
The three-dimensional data in portion, and the three-dimensional data of acquisition is smoothed, certainly according to the three-dimensional data after smoothing processing
The median sagittal section of dynamic detection Fetal Brain is simultaneously shown solve doctor and be difficult to that median sagittal section is accurately positioned manually
The problem of, so that the case where doctor observes Fetal Brain median sagittal section in which can be convenient, can provide a large amount of weights for doctor
The key message wanted.
As shown in Fig. 2, in step S04, it is described that multiple smoothing processing is carried out to three-dimensional data, comprising:
S10 uses the first L0Gradient minimisation model carries out first smoothing processing to three-dimensional data；
S11 carries out secondary smoothing processing to the three-dimensional data after first smoothing processing using Gaussian filter；
S12 is to the blurring three-dimensional data generated after secondary smoothing processing, using the 2nd L0Gradient minimisation model into
Row smoothing processing.
Noise is inevitably introduced into three-dimensional data acquisition process, this noise like would generally pollute said three-dimensional body
Important feature in data, to influence the detection to subsequent median sagittal section.In the prior art, L0Gradient minimisation model
Essence be done under an Optimization Framework sparse gradient statistics, L0Gradient minimisation model passes through control non-zero image gradient
Number enhance image prominent edge part, successively reach the global optimization of image.If only passing through L0Gradient minimisation method
Noise is handled, although can globally retain important feature and the region of volume data with equivalent form, cannot efficiently be gone
Except borderline noise.And gaussian filtering is most widely used smoothing filter, Gaussian convolution for flat site have compared with
Good effect, but in boundary or the place of detail textures, it will generate and obscure or cross smooth effect.Lacking based on the prior art
Fall into, the present embodiment improves filter method in the prior art, innovatively propose it is a kind of to collected three-dimensional data into
The mechanism of the multiple smoothing processing of row uses the first L in the mechanism first0Gradient minimisation model carries out three-dimensional data first
Smoothing processing globally can retain important feature and region with equivalent form, then be carried out using Gaussian filter secondary flat
It is sliding, the 2nd L is used to the secondary blurring three-dimensional data smoothly generated0Gradient minimisation model is smoothed, can
Effectively enhance boundary information, to be beneficial to the subsequent feature according to fetus head median sagittal section, after smoothing processing
Three-dimensional data in detect median sagittal section, improve the efficiency of detection.
In one embodiment, three-dimensional of the three-dimensional data as Q, after first smoothing processing before setting first smoothing processing
Volume data is Q ', sets the first L0The objective function of gradient minimisation model are as follows:
Wherein
E (Q ')=# i | | B 'i,x|+|B′i,y|+|B′i,z|≠0}
In formula, λ is the weight coefficient for controlling smooth item E (Q '), and by artificially being set, # { } is indicated the occurrence of λ
One statistical operation meets for counting | B 'i,x|+| B 'i,y|+|B′i,z| the number of the voxel of result non-zero, wherein B 'i,xFor
In three-dimensional data after first smoothing processing i-th of voxel in the x-direction with the scalar value difference of nearest voxel, B 'i,yTo put down for the first time
In sliding treated three-dimensional data i-th of voxel in the y-direction with the scalar value difference of nearest voxel, B 'i,zFor first smoothing processing
In three-dimensional data afterwards i-th of voxel in the z-direction with the scalar value difference of nearest voxel.
The present embodiment is to L0The objective function of gradient minimisation model improves, so that L0Gradient minimisation method can
Suitable for handling three-dimensional data, wherein the specific formula for calculation for setting smooth item E (Q ') can by the calculation formula
Preferably, more easily estimate the gradient information in three-dimensional data, improve and collected three-dimensional data is carried out for the first time
The efficiency of smoothing processing.
In one embodiment, if the blurring three-dimensional data before smoothing processing is P, the blurring three after smoothing processing
Dimension volume data is P ', sets the 2nd L0The objective function of gradient minimisation model are as follows:
Wherein
In formula, λ is the weight coefficient for controlling smooth item E (Q '), and by artificially being set, # { } is indicated the occurrence of λ
One statistical operation meets for countingAs a result the number of the voxel of non-zero, wherein B 'j,xFor
In blurring three-dimensional data after smoothing processing j-th of voxel in the x-direction with the scalar value difference of nearest voxel, B 'j,yIt is smooth
Treated blurring three-dimensional data in j-th of voxel in the y-direction with the scalar value difference of nearest voxel, B 'j,zFor smoothing processing
In blurring three-dimensional data afterwards j-th of voxel in the z-direction with the scalar value difference of nearest voxel.
The present embodiment is to L0The objective function of gradient minimisation model improves, so that L0Gradient minimisation method can
It is blurred three-dimensional data suitable for processing, wherein setting the specific formula for calculation of smooth item E (Q '), passes through calculating public affairs
Formula, can preferably, the more easily gradient information in ambiguous estimation three-dimensional data, improve to blurring said three-dimensional body number
According to the efficiency being smoothed.
Finally it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than the present invention is protected
The limitation of range is protected, although explaining in detail referring to preferred embodiment to the present invention, those skilled in the art are answered
Work as understanding, it can be with modification or equivalent replacement of the technical solution of the present invention are made, without departing from the reality of technical solution of the present invention
Matter and range.
Claims (6)
1. a kind of method of 3-D supersonic imaging, characterized in that include:
Emit ultrasonic wave to fetus head；
Ultrasonic echo is received, ultrasound echo signal is obtained；
The three-dimensional data of fetus head is obtained according to the ultrasound echo signal；
Multiple smoothing processing is carried out to three-dimensional data；
According to the feature of fetus head median sagittal section, median sagittal is detected from the three-dimensional data after smoothing processing and is cut
Face；
Show the median sagittal section.
2. a kind of method of 3-D supersonic imaging according to claim 1, characterized in that according to fetus head median sagittal
The gray feature of section detects the median sagittal section from the three-dimensional data after smoothing processing.
3. a kind of method of 3-D supersonic imaging according to claim 2, characterized in that from the said three-dimensional body after smoothing processing
The median sagittal section is detected in data, comprising:
It extracts to represent in three-dimensional data after smoothing processing and meets the gray value that the gray value in plane is greater than plane two sides
Condition plane sagittal plane characteristic area；
At least three characteristic points are selected from the sagittal plane characteristic area；
The plane where the median sagittal section is determined according at least three feature points.
4. a kind of method of 3-D supersonic imaging according to claim 2, characterized in that from the said three-dimensional body after smoothing processing
The median sagittal section is detected in data, comprising:
At least two sections are extracted in the three-dimensional data；
It is greater than the feature of the gray value of brain middle line two sides according to the gray value of brain middle line, extracts brain at least two section
Middle line obtains at least two brain middle lines；
The plane where the median sagittal plane is determined according to the plane that at least two brain middle lines limit.
5. a kind of method of 3-D supersonic imaging according to claim 1, characterized in that described to be carried out to three-dimensional data
Multiple smoothing processing, comprising:
(1) the first L is used0The three-dimensional data of gradient minimisation model pair carries out first smoothing processing；
(2) secondary smoothing processing is carried out to the three-dimensional data after first smoothing processing using Gaussian filter；
(3) to the blurring three-dimensional data generated after secondary smoothing processing, using the 2nd L0Gradient minimisation model carries out smooth
Processing.
6. a kind of method of 3-D supersonic imaging according to claim 1, characterized in that before setting first smoothing processing
Three-dimensional data is Q, and the three-dimensional data after first smoothing processing is Q ', sets the first L0The target letter of gradient minimisation model
Number are as follows:
Wherein
E (Q ')=# i | | B 'I, x|+|B′I, y|+|B′I, z|≠0}
In formula, λ is the weight coefficient for controlling smooth item E (Q '), and for the occurrence of λ by artificially being set, # { } indicates one
Statistical operation meets for counting | B 'I, x|+|B′I, y|+|B′I, z| the number of the voxel of result non-zero, wherein B 'I, xIt is first
In three-dimensional data after smoothing processing i-th of voxel in the x-direction with the scalar value difference of nearest voxel, B 'I, yFor first smooth place
In three-dimensional data after reason i-th of voxel in the y-direction with the scalar value difference of nearest voxel, B 'I, zAfter first smoothing processing
In three-dimensional data i-th of voxel in the z-direction with the scalar value difference of nearest voxel.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810827517.6A CN109009226A (en) | 2018-07-25 | 2018-07-25 | A kind of method of 3-D supersonic imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810827517.6A CN109009226A (en) | 2018-07-25 | 2018-07-25 | A kind of method of 3-D supersonic imaging |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109009226A true CN109009226A (en) | 2018-12-18 |
Family
ID=64645091
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810827517.6A Pending CN109009226A (en) | 2018-07-25 | 2018-07-25 | A kind of method of 3-D supersonic imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109009226A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111368586A (en) * | 2018-12-25 | 2020-07-03 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method and system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070285421A1 (en) * | 2006-06-12 | 2007-12-13 | Kabushiki Kaisha Toshiba | Ultrasonic imaging apparatus, an image-processing apparatus, and an ultrasonic image-processing method |
| CN101292883A (en) * | 2007-04-23 | 2008-10-29 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic three-dimensional quick imaging method and apparatus |
| CN103338707A (en) * | 2011-01-26 | 2013-10-02 | 株式会社日立医疗器械 | Ultrasonic diagnostic device and image processing method |
| CN104414680A (en) * | 2013-08-21 | 2015-03-18 | 深圳迈瑞生物医疗电子股份有限公司 | Three-dimensional ultrasonic imaging method and system |
| CN105389811A (en) * | 2015-10-30 | 2016-03-09 | 吉林大学 | Multi-modality medical image processing method based on multilevel threshold segmentation |
-
2018
- 2018-07-25 CN CN201810827517.6A patent/CN109009226A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070285421A1 (en) * | 2006-06-12 | 2007-12-13 | Kabushiki Kaisha Toshiba | Ultrasonic imaging apparatus, an image-processing apparatus, and an ultrasonic image-processing method |
| CN101292883A (en) * | 2007-04-23 | 2008-10-29 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic three-dimensional quick imaging method and apparatus |
| CN103338707A (en) * | 2011-01-26 | 2013-10-02 | 株式会社日立医疗器械 | Ultrasonic diagnostic device and image processing method |
| CN104414680A (en) * | 2013-08-21 | 2015-03-18 | 深圳迈瑞生物医疗电子股份有限公司 | Three-dimensional ultrasonic imaging method and system |
| CN105389811A (en) * | 2015-10-30 | 2016-03-09 | 吉林大学 | Multi-modality medical image processing method based on multilevel threshold segmentation |
Non-Patent Citations (1)
| Title |
|---|
| 王启超: "面向体数据可视化的数据预处理方法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111368586A (en) * | 2018-12-25 | 2020-07-03 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10561403B2 (en) | Sensor coordinate calibration in an ultrasound system | |
| US11229419B2 (en) | Method for processing 3D image data and 3D ultrasonic imaging method and system | |
| CN108053417B (en) | lung segmentation device of 3D U-Net network based on mixed rough segmentation characteristics | |
| KR101121353B1 (en) | System and method for providing 2-dimensional ct image corresponding to 2-dimensional ultrasound image | |
| KR101121396B1 (en) | System and method for providing 2-dimensional ct image corresponding to 2-dimensional ultrasound image | |
| US20110125016A1 (en) | Fetal rendering in medical diagnostic ultrasound | |
| DE102012108121A1 (en) | Method and system for ultrasound-assisted automatic detection, quantification and tracking of pathologies | |
| CN105266849B (en) | Real-time ultrasound elastograph imaging method and system | |
| KR20090127100A (en) | System and method for extracting anatomical feature | |
| NL2025814B1 (en) | Precise positioning system for three-dimensional position of tumor by multi-image fusion | |
| CN107680110B (en) | Inner ear three-dimensional level set segmentation method based on statistical shape model | |
| WO2005011470A2 (en) | System and method for assessing fetal abnormality based on landmarks | |
| CN108670301B (en) | Transverse process positioning method for vertebral column based on ultrasonic image | |
| CN104318565A (en) | Interactive method for retinal vessel segmentation based on bidirectional region growing of constant-gradient distance | |
| CN109009226A (en) | A kind of method of 3-D supersonic imaging | |
| Malathi et al. | Statistical Measurement of Ultrasound Placenta Images Complicated by Gestational Diabetes Mellitus Using Segmentation Approach. | |
| CN116994113A (en) | Automatic segmentation of liver and tumor in CT image based on residual UNet and efficient multi-scale attention method | |
| Lawankar et al. | Segmentation of liver using marker watershed transform algorithm for CT scan images | |
| Lazrag et al. | Combination of the Level‐Set Methods with the Contourlet Transform for the Segmentation of the IVUS Images | |
| Al-jaburi et al. | Medical ultrasound image quality enhancement and regions segmentation | |
| CN111242853B (en) | Medical CT image denoising method based on optical flow processing | |
| KR101024857B1 (en) | Ultrasound system and method for performing color modeling processing on three-dimensional ultrasound image | |
| Sriraam et al. | Performance evaluation of linear and nonlinear filters for despeckling B mode foetal heart ultrasound images | |
| Maitra et al. | Artifact Suppression and Homogenous Orientation of Digital Mammogram using Seeded Region Growing Algorithm | |
| Ruyuan et al. | 3-D target recognition based on sliced images |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181218 |
|
| RJ01 | Rejection of invention patent application after publication |