CN106384337A - Laser speckle blood flow imaging enhancement method - Google Patents
Laser speckle blood flow imaging enhancement method Download PDFInfo
- Publication number
- CN106384337A CN106384337A CN201610805929.0A CN201610805929A CN106384337A CN 106384337 A CN106384337 A CN 106384337A CN 201610805929 A CN201610805929 A CN 201610805929A CN 106384337 A CN106384337 A CN 106384337A
- Authority
- CN
- China
- Prior art keywords
- laser speckle
- blood flow
- algorithm
- speckle blood
- enhancement 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000003384 imaging method Methods 0.000 title claims abstract description 36
- 230000017531 blood circulation Effects 0.000 title abstract description 21
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 44
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 17
- 239000008280 blood Substances 0.000 claims description 26
- 210000004369 blood Anatomy 0.000 claims description 26
- 238000010586 diagram Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000011946 reduction process Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 12
- 238000001727 in vivo Methods 0.000 abstract description 4
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 210000004088 microvessel Anatomy 0.000 abstract 1
- 230000008092 positive effect Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 9
- 238000002474 experimental method Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 101100204393 Arabidopsis thaliana SUMO2 gene Proteins 0.000 description 1
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 101150112492 SUM-1 gene Proteins 0.000 description 1
- 101150096255 SUMO1 gene Proteins 0.000 description 1
- 101100311460 Schizosaccharomyces pombe (strain 972 / ATCC 24843) sum2 gene Proteins 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 210000003710 cerebral cortex Anatomy 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 230000010247 heart contraction Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Public Health (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Hematology (AREA)
- Cardiology (AREA)
- Image Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The present invention provides a laser speckle blood flow imaging enhancement method. The laser speckle blood flow imaging enhancement method is configured to perform noise reduction processing and the repairing processing of the laser speckle blood flow imaging to enhance the sharpness of the laser speckle blood flow imaging. Because the laser speckle blood flow imaging enhancement method employs the BEEMD processing algorithm for the laser speckle blood flow imaging, the microvessel image is more clearly displayed and more details are reflected so as to reduce the influence of the vibration noise on the blood flow dispersed map; because the employed BEEMD algorithm can effectively avoid the information loss at the aspect of the image rank correlation of the original EEMD algorithm to retain the completer image information as much as possible prior to the noise reduction operation; and finally, because the employed BEEMD algorithm is far better than each two-dimensional empirical mode decomposition (BEMD) on the computation speed, so that the laser speckle blood flow imaging enhancement method has a positive effect on the real-time monitoring in vivo.
Description
Technical field
The present invention relates to a kind of enhancing algorithm of laser speckle blood current imaging, it is based on two-dimensional ensemble experience particularly to a kind of
Mode decomposition (Bi-Dimensional Ensemble Empirical Mode Decomposition, BEEMD) algorithm swash
Light speckle blood flow imaging Enhancement Method, belongs to processing of biomedical signals technical field.
Background technology
Laser speckle contrast imaging technique (Laser Speckle Contrast Imaging, LSCI) is to be transported by analysis
Dynamic granule to obtain the technology of particle velocity to the scattering propertiess of coherent laser, and the blood flow of two dimension also can be provided to divide simultaneously
Cloth image.By laser speckle relative analyses, the data of record is processed, obtain contrast images, this image can reflect blood
The velocity information of stream is it is also possible to reflect microvascular distribution situation by image.To examine from experimental viewpoint and hardware device angle
Consider, laser speckle contrast imaging technique is all highstrung for any shake in imaging process.For example：In simulation test
Bionical flexible pipe can produce vibrations noise, the breathing of white mouse and human body, heart beating etc. because of the creeping effect of micropump and all can produce
Different degrees of shake.What the above was shaken directly affects the resolution loss being exactly to cause speckle contrast images.
Gathering empirical mode decomposition (Ensemble Empirical Mode Decomposition, EEMD) is
A kind of noise assistance data analysis method that the researcheres such as N.E.Huang propose on the basis of EMD, compared to Fourier decomposition and
Wavelet-decomposing method, its clear superiority is that and is suitable for analyzing non-linear, non-stationary signal sequence and having very high signal to noise ratio,
Solve the problems, such as modal overlap present in EMD algorithm simultaneously.EEMD algorithm decomposes in one-dimensional data, such as ocean, air, sky
The aspect such as body observational data and earthquake record analysis, mechanical fault detection, clinical medicine monitoring is applied widely, but two
Also there are problems in the data processing of dimension and multidimensional signal and signal analysis.For two-dimensional space data or image, I
Often it can be regarded as one-dimensional data sequence combined sequence in an x or y direction, initial algorithm is by 2-D data (figure
Picture) every a line (or row) individually carry out EEMD process, then all of row (column) result is combined, though do so
The image of different levels so can be decomposited, but due to it by every one-dimensional vector of 2D signal as an independent process,
Have ignored the dependency of 2D signal, therefore process in this way the result drawing be difficult to satisfactory.Additionally, also there being research
Person proposes, using different two-dimensional empirical mode decomposition (BEMD) methods, to be such as based on Delaunay Triangulation and curved surface differential technique
BEMD algorithm, the maximum defect of this kind of algorithm is exactly that logic is complicated, and the processing procedure used time is oversize.
Content of the invention
The present invention be directed to the pixel loss problem two-dimensional ensemble producing because of vibrations in laser speckle blood current imaging
The noise reduction process that empirical mode decomposition BEEMD algorithm is carried out, enables more clearly to show micro- blood-stream image, meanwhile, should
Algorithm passes through whole to natural mode of vibration component layers (Intrinsic Mode Function, IMF) the procession superposition decompositing
Close, the row (column) data dependence disappearance problem that original EEMD algorithm is occurred in 2-D data processing procedure is repaiied
Just.
What the present invention provided employs following technical scheme：
A kind of Enhancement Method of laser speckle blood current imaging, for carrying out noise reduction process and repairing to laser speckle blood current imaging
Process again and obtain speckle comparison diagram to strengthen the definition of laser speckle blood current imaging it is characterised in that comprising the following steps：
Step 1, gathers the gray-scale maps of speckle using laser speckle blood current imaging device, and the gray scale collecting pattern is stored in calculating
To treat to process further in machine, step 2, with EEMD algorithm, process is carried out on y (arranging) direction to gray scale pattern and obtains multilamellar IMF
Tomographic image；Step 3, carries out resolution process with EEMD algorithm to every layer of IMF image on x (OK) direction；Step 4, at decomposition
Result after reason is overlapped according to ranks rule, obtains final natural mode of vibration component map；Step 5, to natural mode of vibration component
Figure carries out threshold filter process, removes noise and residual, information, and the component map superposition calculation of reservation is synthesized an opening and closing one-tenth
Image；Step 6, carries out speckle comparing calculation to the image of synthesis, obtains laser speckle comparison diagram.
The Enhancement Method of the laser speckle blood current imaging that the present invention provides, can also have the feature that：Wherein, step
2 include following sub-step：
Step 2.1, the gray-scale maps by collection in step 1 are the 2-D data sequence of N row M row as formula (1)
Take m column data, the span (1, M) of m
Step 2.2, is decomposed to m column data with EEMD algorithm
In formula (3), J represents the IMF number of plies of decomposition, and j represents jth layer data；
Step 2.3, after the whole decomposition of M column data finishes, we are rearranged into the matrix of such as formula (4)
The Enhancement Method of the laser speckle blood current imaging that the present invention provides, can also have the feature that：Wherein, step
3 include following sub-step：
Step 3.1, in gjLine n data, the span (1, N) of n is taken in (m, n)：
gj(～, n)=(C1, n, jC2, n, j… CM, n, j) (5)；
Step 3.2, is decomposed to line n data with EEMD algorithm：
In formula (6), K represents the IMF number of plies of decomposition, and k represents kth layer data；
Step 3.3, each layer data decomposed is rearranged into the matrix of a formula (7) in step 6：
Invention effect and effect
According to the present invention provide laser speckle blood current imaging Enhancement Method, due to using BEEMD for laser speckle
The Processing Algorithm of blood flow imaging, so that blood capillary image shows relatively sharp, can reflect more details, reduce vibrations and make an uproar
The impact to blood flow speckle pattern for the sound；BEEMD algorithm due to adopting can be effectively prevented from former EEMD algorithm in image ranks phase again
The loss of learning that closing property aspect causes is so that remained complete image information as much as possible before carrying out denoising operation；?
Afterwards due also to adopt BEEMD algorithm all kinds of two-dimensional empirical mode decompositions (BEMD) algorithm is much better than on calculating speed so that
There is positive role to real time monitor in vivo blood flow.
Brief description
The image capturing system laser speckle blood current imaging system principle diagram (1 that Fig. 1 is related to for embodiments of the invention：
Laser instrument 2：Plane mirror 3：Beam expanding lens 4：Object stage 5：CCD camera 6：Computer)；
The BEEMD algorithm general diagram that Fig. 2 is related to for embodiments of the invention；
The BEEMD algorithm all directions computing block diagram that Fig. 3 is related to for embodiments of the invention, 3 (a) decomposes for first direction, and 3
B () decomposes for second direction；
The Simulation of Blood contrast experiment figure that Fig. 4 is related to for embodiments of the invention is schemed in body cortex contrast experiment with rat
(lower row is former EEMD algorithm effect, and upper row is BEEMD algorithm effect used in the present invention).
Specific embodiment
In order that technological means, creation characteristic, reached purpose and effect that the present invention realizes are easy to understand, real below
Apply that example combines the flow process of Enhancement Method of the laser speckle blood current imaging to the present invention for the accompanying drawing and concrete operation method, principle are made to have
Body illustrates.
The image capturing system laser speckle blood current imaging system principle that Fig. 1 is related to for embodiments of the invention
Block diagram (1：Laser instrument 2：Plane mirror 3：Beam expanding lens 4：Object stage 5：CCD camera 6：Computer).
1. gather speckle pattern f (x, y) using laser speckle blood current imaging device (as shown in Figure 1), picture size is
512(M)×512(N).Bionical flexible pipe has been used as simulation experiment material, experiments in vivo then selects rat brain cortex in experiment
It is sampled.The flexible pipe collecting (cerebral cortex) speckle gray scale pattern is stored in computer to treat to process further.
The BEEMD algorithm general diagram that Fig. 2 is related to for embodiments of the invention.
2., in units of each column data, value from left to right, using set empirical modal for the blood flow speckle pattern collecting
Decompose all row of EEMD algorithm Ergodic Matrices (as accompanying drawing 2), each column data can be broken down into multiple IMF layer cj(xi, y), this
Wherein both comprise our blood flow informations interested, also contain vibrations noise, other noise and residual, information.
3. the natural mode of vibration component IMF layer after every string being decomposed combines according to original data order, so
It is the formation of one group of image carrying out column processing with EEMD algorithm：
gj(x, y)=cj(x1, y)+cj(x1, y)+...+cj(xM, y)
This group image is the natural mode of vibration component IMF figure of original image.
The BEEMD algorithm all directions computing block diagram that Fig. 3 is related to for embodiments of the invention, 3 (a) decomposes for first direction, and 3
B () decomposes for second direction.
Process 2,3 refers to accompanying drawing 3 (a) row exploded block diagram.
4. pair one group of natural mode of vibration component map g obtainingj(x, y) enters every trade and processes, with behavior unit repeat step 2, step
Rapid 3 operations, have carried out an EEMD row again and have processed, each IMF has resolved into a series of again to each natural mode of vibration IMF layer
New IMF layer DJ, k(x, y), combines exploded view of embarking on journey,
hJ, k(x, y)=DJ, 1(x, y)+DJ, 2(x, y)+...+DJ, K(x, y)
Process 4 refers to accompanying drawing 3 (b) row exploded block diagram.
5. carry out IMF layer ranks superposition calculation.Each layer IMF figure after column processing is processed through new EEMD row, each layer
Row IMF figure is got back the row IMF figure of same number, all row IMF figure superposition of first row IMF figure, then with each row IMF figure
First row IMF figure superposition obtain final ground floor IMF figure.
SUM1=h1,1(x, y)+h1,2(x, y)+...+h1, K(x, y)+h2,1(x, y)+...+hJ, 1(x, y)
In the same manner, all row IMF superposition of secondary series IMF figure, then be superimposed with the second row IMF figure of every string IMF figure and obtain
Second layer IMF schemes SUM2, and the rest may be inferred, calculates the IMF layer of all speckle patterns.Said process refers to accompanying drawing 2.
6. pair each tomographic image analysis, carries out threshold filter denoising, retains useful information, remove noise and residual risk.Will
Useful information layer is overlapped calculating, and synthesizes an image.
7. the image of pair synthesis carries out speckle comparing calculation, obtains final laser speckle blood flow comparison diagram,
The Simulation of Blood contrast experiment figure that Fig. 4 is related to for embodiments of the invention is schemed in body cortex contrast experiment with rat
(lower row is former EEMD algorithm effect, and upper row is BEEMD algorithm effect used in the present invention).
Three groups of experimental paradigm, one group (1,2), two groups (3,4), three groups (5,6) are illustrated, first two groups is blood vessel in accompanying drawing 4
Human simulation is tested, Bottomhole pressure Fluid simulation blood flow, the 3rd group of cortex rheography shooting in body for rat.Upper row three figure
Using the BEEMD algorithm in the present invention, lower row three figure employs original EEMD algorithm, and 1 figure more clearly visible shows compared with 2 figures
Liquid in pipe, 3 figures decrease liquid in pipe compared with 4 figures and damage because of the pixel that vibration facter or processing procedure associated deletion lead to
Consumption, 5 figures show more microvascular details compared with 6 figures, and the definition of blood vessel is effectively improved.
The beneficial effect of embodiment and effect
According to the present embodiment provide laser speckle blood current imaging Enhancement Method, due to due to using BEEMD for swash
The Processing Algorithm of light speckle blood flow imaging, so that blood capillary image shows relatively sharp, can reflect more details, reduce
The vibrations impact to blood flow speckle pattern for the noise；
BEEMD algorithm due to adopting can be effectively prevented from what former EEMD algorithm caused in terms of image ranks dependency again
Loss of learning is so that remained complete image information as much as possible before carrying out denoising operation；
Finally due also to the BEEMD algorithm adopting is much better than all kinds of two-dimensional empirical mode decompositions (BEMD) in calculating speed
Algorithm is so that have positive role to real time monitor in vivo blood flow.
Claims (3)
1. a kind of Enhancement Method of laser speckle blood current imaging, for carrying out noise reduction process and reparation to laser speckle blood current imaging
Process and obtain speckle comparison diagram to strengthen the definition of described laser speckle blood current imaging it is characterised in that including following walking
Suddenly：
Step 1, gathers the gray-scale maps of speckle using laser speckle blood current imaging device, and the described gray scale pattern collecting is deposited
Storage is to treat to process further；
Step 2, carries out process with EEMD algorithm to described gray scale pattern on y (arranging) direction and obtains multilamellar IMF tomographic image；
Step 3, carries out resolution process with EEMD algorithm to IMF image every layer described on x (OK) direction；
Step 4, the result after described resolution process is overlapped according to ranks rule, obtains final natural mode of vibration component
Figure；
Step 5, carries out threshold filter process to described natural mode of vibration component map, remove noise and residual, information, and will retain
Component map superposition calculation synthesizes a composograph；
Step 6, carries out speckle comparing calculation to the image of described synthesis, obtains laser speckle comparison diagram.
2. laser speckle blood current imaging according to claim 1 Enhancement Method it is characterised in that：
Wherein, described step 2 includes following sub-step：
Step 2.1, the described gray-scale maps by collection in step 1 are the 2-D data sequence of N row M row as formula (1)
Take m column data, the span (1, M) of m
Step 2.2, is decomposed to described m column data with EEMD algorithm
In formula (3), J represents the IMF number of plies of decomposition, and j represents jth layer data；
Step 2.3, after the whole decomposition of M column data finishes, is rearranged into the matrix of such as formula (4)
3. laser speckle blood current imaging according to claim 1 Enhancement Method it is characterised in that：
Wherein, step 3 includes following sub-step：
Step 3.1, in gjLine n data, the span (1, N) of n is taken in (m, n)：
gj(～, n)=(c1, n, jc2, n, j… cM, n, j) (5)；
Step 3.2, is decomposed to line n data with EEMD algorithm：
In formula (6), K represents the IMF number of plies of decomposition, and k represents kth layer data；
Step 3.3, each layer data decomposed is rearranged into the matrix of a formula (7) in step 6：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610805929.0A CN106384337A (en) | 2016-09-07 | 2016-09-07 | Laser speckle blood flow imaging enhancement method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610805929.0A CN106384337A (en) | 2016-09-07 | 2016-09-07 | Laser speckle blood flow imaging enhancement method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106384337A true CN106384337A (en) | 2017-02-08 |
Family
ID=57938953
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610805929.0A Pending CN106384337A (en) | 2016-09-07 | 2016-09-07 | Laser speckle blood flow imaging enhancement method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106384337A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107862724A (en) * | 2017-12-01 | 2018-03-30 | 中国医学科学院生物医学工程研究所 | A kind of improved microvascular blood flow imaging method |
| CN109978782A (en) * | 2019-03-18 | 2019-07-05 | 新疆大学 | The minimizing technology and device of speckle noise in a kind of Reconstructed Image of Digital Holography |
| CN113160080A (en) * | 2021-04-16 | 2021-07-23 | 桂林市啄木鸟医疗器械有限公司 | CR image noise reduction method, device, equipment and medium |
| CN117409153A (en) * | 2023-12-15 | 2024-01-16 | 深圳大学 | Three-dimensional target transmission imaging method in turbid medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102357033A (en) * | 2011-09-27 | 2012-02-22 | 华中科技大学 | Laser speckle blood stream imaging processing system and method |
| CN104392444A (en) * | 2014-11-18 | 2015-03-04 | 陕西师范大学 | Method of extracting characteristics of medical MR (magnetic resonance) images based on ensemble empirical mode decomposition |
-
2016
- 2016-09-07 CN CN201610805929.0A patent/CN106384337A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102357033A (en) * | 2011-09-27 | 2012-02-22 | 华中科技大学 | Laser speckle blood stream imaging processing system and method |
| CN104392444A (en) * | 2014-11-18 | 2015-03-04 | 陕西师范大学 | Method of extracting characteristics of medical MR (magnetic resonance) images based on ensemble empirical mode decomposition |
Non-Patent Citations (4)
| Title |
|---|
| ANNE HUMEAU-HEURTIER,ET,AL.: "Multi-Dimensional Complete Ensemble Empirical Mode Decomposition With Adaptive Noise Applied to Laser Speckle Contrast Images", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 * |
| SHANG-CHING LIN,ET,AL.: "Automatic Contrast Enhancement Using Ensemble Empirical Mode Decomposition", 《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND 2680 FREQUENCY CONTROL》 * |
| ZHAOHUA WU,ET,AL.: "THE MULTI-DIMENSIONAL ENSEMBLE EMPIRICAL MODE DECOMPOSITION METHOD", 《ADVANCES IN ADAPTIVE DATA ANALYSIS》 * |
| 孔平等: "激光散斑血流成像技术研究新进展", 《光学技术》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107862724A (en) * | 2017-12-01 | 2018-03-30 | 中国医学科学院生物医学工程研究所 | A kind of improved microvascular blood flow imaging method |
| CN107862724B (en) * | 2017-12-01 | 2021-10-01 | 中国医学科学院生物医学工程研究所 | Improved microvascular blood flow imaging method |
| CN109978782A (en) * | 2019-03-18 | 2019-07-05 | 新疆大学 | The minimizing technology and device of speckle noise in a kind of Reconstructed Image of Digital Holography |
| CN113160080A (en) * | 2021-04-16 | 2021-07-23 | 桂林市啄木鸟医疗器械有限公司 | CR image noise reduction method, device, equipment and medium |
| CN113160080B (en) * | 2021-04-16 | 2023-09-22 | 桂林市啄木鸟医疗器械有限公司 | CR image noise reduction method, device, equipment and medium |
| CN117409153A (en) * | 2023-12-15 | 2024-01-16 | 深圳大学 | Three-dimensional target transmission imaging method in turbid medium |
| CN117409153B (en) * | 2023-12-15 | 2024-05-07 | 深圳大学 | Three-dimensional target transmission imaging method in turbid medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11120582B2 (en) | Unified dual-domain network for medical image formation, recovery, and analysis | |
| CN101879073B (en) | Ultrasound diagnosis apparatus | |
| US9307958B2 (en) | Ultrasonic diagnostic apparatus and an ultrasonic image processing apparatus | |
| CN106384337A (en) | Laser speckle blood flow imaging enhancement method | |
| Mitsuhashi et al. | A forward-adjoint operator pair based on the elastic wave equation for use in transcranial photoacoustic computed tomography | |
| JP2007209756A (en) | Noise reduction method in tomographic image data set | |
| Hesse et al. | Nonlinear simultaneous reconstruction of inhomogeneous compressibility and mass density distributions in unidirectional pulse-echo ultrasound imaging | |
| CN109978888A (en) | A kind of image partition method, device and computer readable storage medium | |
| Zhao et al. | Ultrasound transmission tomography image reconstruction with a fully convolutional neural network | |
| KR20190131800A (en) | An image processing apparatus using neural network and a method performed by the image processing apparatus | |
| Akkasaligar et al. | Diagnosis of renal calculus disease in medical ultrasound images | |
| Cheng et al. | Dual-IRT-GAN: A defect-aware deep adversarial network to perform super-resolution tasks in infrared thermographic inspection | |
| JP5451414B2 (en) | Subject information processing apparatus and subject information processing method | |
| Zhang et al. | Deep image translation for enhancing simulated ultrasound images | |
| Jiang | MFCT-GAN: multi-information network to reconstruct CT volumes for security screening | |
| CN117115452B (en) | Controllable medical ultrasonic image denoising method, system and computer storage medium | |
| Gao et al. | Ultrasound image super-resolution reconstruction based on semi-supervised CycleGAN | |
| CN112347828A (en) | Self-adaptive micro-motion highlighting system | |
| Pasha et al. | A computational framework for edge-preserving regularization in dynamic inverse problems | |
| CN113272869A (en) | Three-dimensional shape reconstruction from topograms in medical imaging | |
| Molinier et al. | Ultrasonic imaging using conditional generative adversarial networks | |
| Pasha et al. | Efficient edge-preserving methods for dynamic inverse problems | |
| Joshi et al. | Multi-Modality Medical Image Fusion Using SWT & Speckle Noise Reduction with Bidirectional Exact Pattern Matching Algorithm | |
| Tran | 3D image analysis with variational methods and wavelets: applications to medical image processing | |
| Ekmekci et al. | Model-based Bayesian deep learning architecture for linear inverse problems in computational imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170208 |