CN112509095A - Oct image dislocation correction method - Google Patents

Oct image dislocation correction method Download PDF

Info

Publication number
CN112509095A
CN112509095A CN202110165934.0A CN202110165934A CN112509095A CN 112509095 A CN112509095 A CN 112509095A CN 202110165934 A CN202110165934 A CN 202110165934A CN 112509095 A CN112509095 A CN 112509095A
Authority
CN
China
Prior art keywords
dislocation
oct image
point
angle
correction 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.)
Granted
Application number
CN202110165934.0A
Other languages
Chinese (zh)
Other versions
CN112509095B (en
Inventor
张真铨
滕忠照
沈金花
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Jingsan Medical Technology Co ltd
Original Assignee
Nanjing Jingsan Medical Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Jingsan Medical Technology Co ltd filed Critical Nanjing Jingsan Medical Technology Co ltd
Priority to CN202110165934.0A priority Critical patent/CN112509095B/en
Publication of CN112509095A publication Critical patent/CN112509095A/en
Application granted granted Critical
Publication of CN112509095B publication Critical patent/CN112509095B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • G06F18/232Non-hierarchical techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/44Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Multimedia (AREA)
  • Endoscopes (AREA)

Abstract

The invention discloses an oct image dislocation correction method. Reading an oct image in a blood vessel, and acquiring lumen contour information of the blood vessel from the oct image; the oct image and the lumen contour information are taken as the center (x) of the oct image0,y0) Performing polar coordinate conversion on the original point, and obtaining the distance from the lumen contour to the central point at each angle theta within the range of 0-360 degrees according to the lumen contour information; subtracting the distances at intervals of n degrees to obtain a gradient at the maximum angle of the gradient
Figure 100004_DEST_PATH_IMAGE002
Judging the position to be dislocation, wherein the corresponding gradient is dislocation amplitude, if the dislocation amplitude is smaller than a threshold value T, carrying out rough judgment, otherwise, carrying out rough judgmentAnd performing fine judgment. According to the invention, the fine dislocation and the fine dislocation amplitude are calculated, and then the fine dislocation is repaired according to the fine dislocation amplitude, so that the dislocation generated due to the displacement of the relative position of the catheter and the blood vessel wall during oct imaging can be accurately corrected and repaired.

Description

Oct image dislocation correction method
Technical Field
The invention relates to the technical field of oct image processing, in particular to an oct image dislocation correction method.
Background
Optical Coherence Tomography (OCT) is a non-invasive probing technique. OCT technology has been widely used for imaging the structure of a living body cross section of a biological tissue. By measuring the scattered light as a function of depth, OCT can provide high resolution, high sensitivity tissue structures.
When the OCT is applied to blood vessel imaging, near infrared light is emitted to the inner surface of a blood vessel by using a rotatable optical lens and an optical fiber, and reflected light waves are received by using an optical interferometer and imaged. Because the optical wave imaging is utilized, the OCT imaging resolution is high, the axial resolution can reach 10-20um, and the components and the microstructure on the plaque surface can be imaged. However, the near infrared light wave is not very permeable (about 1.0-2.5 mm), and blood cells, red thrombus, and plaque lipid core or plaque necrosis all affect OCT observation of vessel wall structure and estimation of plaque burden. Because of red light scattering by red blood cells, past OCT requires constant injection of contrast agent to wash blood away during imaging. Modern OCT systems partially reduce the interference of red blood cells and the like on imaging by techniques such as rapid rotational withdrawal, and complete imaging of a length of blood vessel in a few seconds.
The dislocation cause of OCT blood vessel imaging: in the process of imaging by emitting near infrared light to the inner surface of a blood vessel through the rapid rotation of the optical lens and the optical fiber on the catheter, the relative position of the catheter and the blood vessel wall is displaced due to the movement of the catheter, the contraction and the relaxation of the blood vessel and the like, and the displacement is compared with the initial position and is accumulated to generate obvious dislocation. The resulting misalignment, shown at the top left in fig. 1, affects the observation and judgment.
Disclosure of Invention
The invention aims to provide an oct image misalignment correction method aiming at the defects in the prior art.
In order to achieve the above object, the present invention provides an oct image misalignment correction method, which includes:
reading an oct image in a blood vessel, and acquiring lumen contour information of the blood vessel from the oct image;
using the oct image and the lumen contour information as the center (x) of the oct image0,y0) Polar coordinate conversion is carried out for the origin, and the distance from the lumen contour to the central point under each angle theta within the range of 0-360 DEG is obtained according to the lumen contour information
Figure DEST_PATH_IMAGE001
At a distance of n DEG each
Figure 650708DEST_PATH_IMAGE001
Subtracting to obtain a gradient
Figure 722569DEST_PATH_IMAGE002
At maximum angle of gradient
Figure DEST_PATH_IMAGE003
Judging the position of the dislocation, wherein,
Figure 306653DEST_PATH_IMAGE004
corresponding gradient
Figure DEST_PATH_IMAGE005
If the dislocation amplitude is smaller than the threshold value T, carrying out rough judgment, otherwise, carrying out fine judgment;
the rough judgment comprises the following steps: the angles are traversed from 0 to 360 in steps of 1 pixel at each angle θ, ranging from
Figure 313923DEST_PATH_IMAGE006
Traversing the distance r to the origin, calculating the spacing distance at each angle as
Figure DEST_PATH_IMAGE007
Gray scale difference of the points
Figure 308424DEST_PATH_IMAGE008
Distance from the point of maximum gradation gradient to the origin as gradation gradient
Figure DEST_PATH_IMAGE009
Namely:
Figure 59343DEST_PATH_IMAGE010
calculating the difference value of the distances from the point with the maximum gray gradient to the original point, and judging the angle with the maximum absolute value as the dislocation
Figure DEST_PATH_IMAGE011
Figure 557320DEST_PATH_IMAGE012
Wherein,
Figure DEST_PATH_IMAGE013
is composed of
Figure 266650DEST_PATH_IMAGE014
The distance from the point with the maximum gray scale gradient of the adjacent angles to the original point is judged as the dislocation
Figure DEST_PATH_IMAGE015
Then carrying out fine judgment;
the fine judgment comprises the following steps: the traversal range is
Figure 623813DEST_PATH_IMAGE016
2 DEG and/or
Figure DEST_PATH_IMAGE017
2 °, step size 0.1 °, at each angle
Figure 303056DEST_PATH_IMAGE018
Distance r from the lower traverse to the center point in the range of
Figure DEST_PATH_IMAGE019
Step size of 0.5 pixels, per angle
Figure 937431DEST_PATH_IMAGE018
Distance from point with maximum lower gray gradient to center point
Figure 614400DEST_PATH_IMAGE020
Namely:
Figure DEST_PATH_IMAGE021
calculating the difference value of the distances from the point with the maximum gray gradient to the origin, and judging the angle with the maximum absolute value as the fine dislocation
Figure 583493DEST_PATH_IMAGE022
Comprises the following steps:
Figure DEST_PATH_IMAGE023
the corresponding fine misalignment magnitudes are:
Figure 673284DEST_PATH_IMAGE024
and carrying out dislocation repair on the fine dislocation part, and specifically comprising the following steps: selecting direction by taking the angle of the fine dislocation as a starting point, traversing each pixel point (x, y) in a fan-shaped range which rotates by a certain angle p along the direction as an end point, and calculating to obtain the radius of the point
Figure DEST_PATH_IMAGE025
Angle of and
Figure 83537DEST_PATH_IMAGE026
wherein:
Figure 931407DEST_PATH_IMAGE027
Figure 387796DEST_PATH_IMAGE028
wherein,
Figure 284208DEST_PATH_IMAGE029
is a unit vector;
calculating the dislocation amplitude of the point
Figure 345705DEST_PATH_IMAGE030
If Δ r > 0, the angle is
Figure 161214DEST_PATH_IMAGE026
Radius of
Figure 100002_DEST_PATH_IMAGE031
And assigning a pixel value at the position to the point, not processing the catheter in the image within the range, if delta r is less than 0, indicating that the contour at the position needs to be expanded, and filling the area between the repaired lumen contour and the lumen contour before repair.
Further, the lumen contour information is obtained through manual segmentation of an interactive interface.
Further, the lumen contour information is obtained by automatic segmentation through a machine learning method.
Further, the direction selection comprises manually selecting a clockwise or counterclockwise direction by manually referring to the blood vessel image information of the upper and lower frames on the interactive interface.
Further, the direction selection includes calculating curvature along the lumen contour at each angle according to the lumen contour information, and dividing two types of angle sets according to the clustering information of the curvature, wherein the direction in which the set with the smaller number is located is the selected direction.
Further, the filling comprises calculating a pixel distribution histogram in a range between the outside of the catheter and the lumen contour, and randomly selecting one of a plurality of pixel values with the highest frequency for filling.
Further, the method further includes performing boundary smoothing processing on the image subjected to the dislocation restoration, and specifically includes: in a manner that
Figure 980266DEST_PATH_IMAGE032
And (4) traversing pixel points in the boundary range by taking +/-1 degrees as the boundary range, and carrying out median filtering by taking 5-5 pixels around the boundary point and eliminating the pixel points in the boundary range as masks.
Further, the angle θ includes 2 °.
Further, the n ° includes 5 °.
Further, the
Figure 742686DEST_PATH_IMAGE007
Including 9 pixels.
Has the advantages that: according to the invention, the fine dislocation and the fine dislocation amplitude are calculated, and then the fine dislocation is repaired according to the fine dislocation amplitude, so that the dislocation generated due to the displacement of the relative position of the catheter and the blood vessel wall during oct imaging can be accurately corrected and repaired.
Drawings
FIG. 1 is a schematic diagram of the conventional misalignment of oct images;
FIG. 2 is a flow chart of an oct image misalignment correction method according to an embodiment of the present invention;
FIG. 3 is a schematic illustration of the oct image after the restoration of the misalignment;
fig. 4 is a schematic diagram of the oct image after the boundary smoothing process.
Detailed Description
The present invention will be further illustrated with reference to the accompanying drawings and specific examples, which are carried out on the premise of the technical solution of the present invention, and it should be understood that these examples are only for illustrating the present invention and are not intended to limit the scope of the present invention.
As shown in fig. 2, an embodiment of the present invention provides a method for correcting an oct image misalignment, including:
and reading an oct image in the blood vessel, and segmenting the oct image to obtain lumen contour information of the blood vessel. The lumen contour information of the obtained blood vessel can be obtained by a certain blood vessel segmentation method. Specifically, an interactive interface can be created, and manual lumen segmentation is performed on the oct image in the blood vessel manually to acquire lumen contour information of the blood vessel. Or automatically segmenting the oct image based on a machine learning method to acquire lumen contour information of the blood vessel.
The oct image and the lumen contour information are taken as the center (x) of the oct image0,y0) Polar coordinate conversion is carried out for the origin, and the distance from the lumen contour to the central point under each angle theta within the range of 0-360 DEG is obtained according to the lumen contour information
Figure 721006DEST_PATH_IMAGE033
. Among them, the angle θ is preferably 2 °.
At a distance of n DEG each
Figure 910679DEST_PATH_IMAGE033
Subtracting to obtain a gradient
Figure 100002_DEST_PATH_IMAGE034
N ° is preferably 5 °, at the maximum angle of the gradient
Figure 748185DEST_PATH_IMAGE003
Judging the position of the dislocation, wherein,
Figure 720820DEST_PATH_IMAGE035
corresponding gradient
Figure 553647DEST_PATH_IMAGE005
If the dislocation amplitude is smaller than the threshold value T, rough judgment is carried out, and otherwise, fine judgment is carried out. The threshold T is typically 10 to 15 unit pixels and can be adjusted as desired.
The rough judgment of the embodiment of the invention comprises the following steps: the angles are traversed from 0 to 360 in steps of 1 pixel at each angle θ, ranging from
Figure 100002_DEST_PATH_IMAGE036
Traversing the distance r to the origin, calculating the spacing distance at each angle as
Figure 320746DEST_PATH_IMAGE007
Gray scale difference of the points
Figure 239023DEST_PATH_IMAGE037
As a result of the gradient of the gray scale,
Figure 343245DEST_PATH_IMAGE007
preferably 9 pixels, the distance from the point of maximum gray gradient to the origin
Figure 100002_DEST_PATH_IMAGE038
Namely:
Figure 905945DEST_PATH_IMAGE039
calculating the difference value of the distances from the point with the maximum gray gradient to the original point, and preliminarily judging the angle with the maximum absolute value of the difference value as a dislocation
Figure 100002_DEST_PATH_IMAGE040
Figure 234158DEST_PATH_IMAGE041
Wherein,
Figure 100002_DEST_PATH_IMAGE042
is composed of
Figure 240729DEST_PATH_IMAGE009
The distance from the point with the maximum gray gradient of the adjacent angles to the original point is judged as the dislocation
Figure 945380DEST_PATH_IMAGE015
And then carrying out fine judgment.
The fine judgment of the embodiment of the invention comprises the following steps: the traversal range is
Figure 690482DEST_PATH_IMAGE016
2 DEG and/or
Figure 330542DEST_PATH_IMAGE043
2 °, step size 0.1 °, at each angle
Figure 692253DEST_PATH_IMAGE018
Distance r from the lower traverse to the center point in the range of
Figure 100002_DEST_PATH_IMAGE044
Step size of 0.5 pixels, per angle
Figure 607120DEST_PATH_IMAGE018
Distance from point with maximum lower gray gradient to center point
Figure 878832DEST_PATH_IMAGE020
Namely:
Figure 752110DEST_PATH_IMAGE045
calculating the difference value of the distances from the point with the maximum gray gradient to the origin, and judging the angle with the maximum absolute value as the fine dislocation
Figure 866697DEST_PATH_IMAGE032
Comprises the following steps:
Figure 100002_DEST_PATH_IMAGE046
the corresponding fine misalignment magnitudes are:
Figure 585254DEST_PATH_IMAGE047
carry out dislocation restoration to meticulous dislocation department, specifically include: the direction selection is carried out by taking the angle of the fine dislocation as a starting point, and has two modes, specifically, the clockwise or anticlockwise direction can be manually selected by manually referring to the blood vessel image information of the upper frame and the lower frame on an interactive interface, or the clockwise or anticlockwise direction can be manually selected according to the lumen contourAnd calculating the curvature of the lumen contour under each angle by the information, dividing two types of angle sets according to the clustering information of the curvature, and taking the direction of the set with the smaller number as the selected direction. Traversing each pixel point (x, y) in a fan-shaped range which rotates for a certain angle p along the direction to be an end point, and calculating to obtain the radius of the point
Figure 39369DEST_PATH_IMAGE025
Angle of and
Figure 21232DEST_PATH_IMAGE026
wherein:
Figure 100002_DEST_PATH_IMAGE048
Figure 623114DEST_PATH_IMAGE028
wherein,
Figure 410942DEST_PATH_IMAGE029
is a unit vector and can be expressed as (0, 1).
Calculating the dislocation amplitude of the point
Figure 922826DEST_PATH_IMAGE049
If Δ r > 0, the angle is
Figure 137906DEST_PATH_IMAGE026
Radius of
Figure 100002_DEST_PATH_IMAGE050
And assigning a pixel value at the position to the point, not processing the catheter in the image within the range, if delta r is less than 0, indicating that the contour at the position needs to be expanded, and filling the area between the repaired lumen contour and the lumen contour before repair.
Filling the region between the lumen contour after repair and the lumen contour before repair comprises calculating a pixel distribution histogram in a range between the outside of the catheter and the lumen contour, and randomly selecting one of a plurality of pixel values with the highest frequency for filling. Can be expressed as:
Figure 227085DEST_PATH_IMAGE051
wherein H is a pixel distribution histogram in a range between the outside of the catheter and the lumen contour,
Figure 100002_DEST_PATH_IMAGE052
for a certain pixel in the pixel distribution histogram H,
Figure 225128DEST_PATH_IMAGE053
is the probability of the occurrence of the pixel,
Figure 100002_DEST_PATH_IMAGE054
for the b-th highest probability of occurrence, b is generally equal to 9, that is, one of the first 9 names with the highest probability of occurrence is randomly selected for filling.
As shown in fig. 3 and 4, a boundary after repairing may have an obvious boundary line, and boundary smoothing is required, so the embodiment of the present invention further includes performing boundary smoothing processing on the image after dislocation repairing, specifically including: in a manner that
Figure 119747DEST_PATH_IMAGE022
And (3) traversing pixel points in the boundary range by taking +/-1 degrees as the boundary range, taking 5-5 pixels around the boundary point and eliminating the pixel points in the boundary range as masks, carrying out median filtering, and then obtaining a finally repaired image.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that other parts not specifically described are within the prior art or common general knowledge to those of ordinary skill in the art. Without departing from the principle of the invention, several improvements and modifications can be made, and these improvements and modifications should also be construed as the scope of the invention.

Claims (10)

1. An oct image misalignment correction method, comprising:
reading an oct image in a blood vessel, and acquiring lumen contour information of the blood vessel from the oct image;
using the oct image and the lumen contour information as the center (x) of the oct image0,y0) Polar coordinate conversion is carried out for the origin, and the distance from the lumen contour to the central point under each angle theta within the range of 0-360 DEG is obtained according to the lumen contour information
Figure DEST_PATH_IMAGE002
At a distance of n DEG each
Figure 122176DEST_PATH_IMAGE002
Subtracting to obtain a gradient
Figure DEST_PATH_IMAGE004
At maximum angle of gradient
Figure DEST_PATH_IMAGE006
Judging the position of the dislocation, wherein,
Figure DEST_PATH_IMAGE008
corresponding gradient
Figure DEST_PATH_IMAGE010
If the dislocation amplitude is smaller than the threshold value T, carrying out rough judgment, otherwise, carrying out fine judgment;
the rough judgment comprises the following steps: the angles are traversed from 0 to 360 in steps of 1 pixel at each angle θ, ranging from
Figure DEST_PATH_IMAGE012
Traversing the distance r to the origin, calculating the spacing distance at each angle as
Figure DEST_PATH_IMAGE014
Gray scale difference of the points
Figure DEST_PATH_IMAGE016
Distance from the point of maximum gradation gradient to the origin as gradation gradient
Figure DEST_PATH_IMAGE018
Namely:
Figure DEST_PATH_IMAGE020
calculating the difference value of the distances from the point with the maximum gray gradient to the original point, and judging the angle with the maximum absolute value as the dislocation
Figure DEST_PATH_IMAGE022
Figure DEST_PATH_IMAGE024
Wherein,
Figure DEST_PATH_IMAGE026
is composed of
Figure 738097DEST_PATH_IMAGE018
The distance from the point with the maximum gray scale gradient of the adjacent angles to the original point is judged as the dislocation
Figure DEST_PATH_IMAGE027
Then carrying out fine judgment;
the fine judgment comprises the following steps: the traversal range is
Figure DEST_PATH_IMAGE028
2 DEG and/or
Figure DEST_PATH_IMAGE029
2 °, step size 0.1 °, at each angle
Figure DEST_PATH_IMAGE031
Distance r from the lower traverse to the center point in the range of
Figure DEST_PATH_IMAGE032
Step size of 0.5 pixels, per angle
Figure 17899DEST_PATH_IMAGE031
Distance from point with maximum lower gray gradient to center point
Figure DEST_PATH_IMAGE034
Namely:
Figure DEST_PATH_IMAGE036
calculating the difference value of the distances from the point with the maximum gray gradient to the origin, and judging the angle with the maximum absolute value as the fine dislocation
Figure DEST_PATH_IMAGE038
Comprises the following steps:
Figure DEST_PATH_IMAGE040
the corresponding fine misalignment magnitudes are:
Figure DEST_PATH_IMAGE042
and carrying out dislocation repair on the fine dislocation part, and specifically comprising the following steps: selecting direction by taking the angle of the fine dislocation as a starting point, traversing each pixel point (x, y) in a fan-shaped range which rotates by a certain angle p along the direction as an end point, and calculating to obtainRadius of the point
Figure DEST_PATH_IMAGE044
Angle of and
Figure DEST_PATH_IMAGE046
wherein:
Figure DEST_PATH_IMAGE048
Figure DEST_PATH_IMAGE050
wherein,
Figure DEST_PATH_IMAGE052
is a unit vector;
calculating the dislocation amplitude of the point
Figure DEST_PATH_IMAGE054
If Δ r > 0, the angle is
Figure 15375DEST_PATH_IMAGE046
Radius of
Figure DEST_PATH_IMAGE056
And assigning a pixel value at the position to the point, not processing the catheter in the image within the range, if delta r is less than 0, indicating that the contour at the position needs to be expanded, and filling the area between the repaired lumen contour and the lumen contour before repair.
2. The oct image-misalignment correction method according to claim 1, wherein the lumen contour information is obtained by manual segmentation via an interactive interface.
3. The oct image misalignment correction method according to claim 1, wherein the lumen contour information is obtained by automatic segmentation using a machine learning method.
4. The oct image-misalignment-correction method according to claim 1, wherein the direction selection includes manually selecting a clockwise or counterclockwise direction with reference to the vessel image information of the upper and lower frames on the interactive interface.
5. The oct image-misalignment correction method according to claim 1, wherein the direction selection includes calculating a curvature along the lumen contour at each angle according to the lumen contour information, and dividing two types of angle sets according to the curvature clustering information, wherein a direction in which a smaller number of types of sets are located is the selected direction.
6. The oct image-misalignment correction method according to claim 1, wherein the filling comprises calculating a pixel distribution histogram in a range between the outside of the catheter and the lumen contour, and randomly selecting one of a plurality of pixel values having the highest frequency for filling.
7. The oct image misalignment correction method according to claim 1, further comprising performing boundary smoothing on the image after misalignment restoration, specifically comprising: in a manner that
Figure DEST_PATH_IMAGE057
And (4) traversing pixel points in the boundary range by taking +/-1 degrees as the boundary range, and carrying out median filtering by taking 5-5 pixels around the boundary point and eliminating the pixel points in the boundary range as masks.
8. The oct image misalignment correction method of claim 1, wherein the angle θ comprises 2 °.
9. The oct image misalignment correction method of claim 1, wherein the n ° includes 5 °.
10. The oct image misalignment correction method of claim 1, wherein the oct image misalignment correction method is performed in a computer
Figure 78140DEST_PATH_IMAGE014
Including 9 pixels.
CN202110165934.0A 2021-02-07 2021-02-07 Oct image dislocation correction method Active CN112509095B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110165934.0A CN112509095B (en) 2021-02-07 2021-02-07 Oct image dislocation correction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110165934.0A CN112509095B (en) 2021-02-07 2021-02-07 Oct image dislocation correction method

Publications (2)

Publication Number Publication Date
CN112509095A true CN112509095A (en) 2021-03-16
CN112509095B CN112509095B (en) 2021-05-07

Family

ID=74953191

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110165934.0A Active CN112509095B (en) 2021-02-07 2021-02-07 Oct image dislocation correction method

Country Status (1)

Country Link
CN (1) CN112509095B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113349737A (en) * 2021-06-30 2021-09-07 深圳英美达医疗技术有限公司 Method for calibrating OCT image of intravascular dual-mode imaging system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101398935A (en) * 2007-09-29 2009-04-01 上海西门子医疗器械有限公司 Method for judging image dislocation
CN104237167A (en) * 2013-09-09 2014-12-24 深圳市斯尔顿科技有限公司 Correction method and system for distortion of scanning device during OCT sectional image scanning
US20150150447A1 (en) * 2012-06-15 2015-06-04 Oregon Health & Science University Non-invasive 3d imaging and measuring of anterior chamber angle of the eye
US20160360962A1 (en) * 2015-06-11 2016-12-15 Tomey Corporation Anterior Ocular Segment Optical Coherence Tomographic Imaging Device and Anterior Ocular Segment Optical Coherence Tomographic Imaging Method
CN110033496A (en) * 2019-03-27 2019-07-19 江苏理工学院 The motion artifact bearing calibration of time series three-dimensional retina SD-OCT image
CN111461961A (en) * 2020-03-27 2020-07-28 佛山科学技术学院 OCT blood vessel image dislocation correction method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101398935A (en) * 2007-09-29 2009-04-01 上海西门子医疗器械有限公司 Method for judging image dislocation
US20150150447A1 (en) * 2012-06-15 2015-06-04 Oregon Health & Science University Non-invasive 3d imaging and measuring of anterior chamber angle of the eye
CN104237167A (en) * 2013-09-09 2014-12-24 深圳市斯尔顿科技有限公司 Correction method and system for distortion of scanning device during OCT sectional image scanning
US20160360962A1 (en) * 2015-06-11 2016-12-15 Tomey Corporation Anterior Ocular Segment Optical Coherence Tomographic Imaging Device and Anterior Ocular Segment Optical Coherence Tomographic Imaging Method
CN110033496A (en) * 2019-03-27 2019-07-19 江苏理工学院 The motion artifact bearing calibration of time series three-dimensional retina SD-OCT image
CN111461961A (en) * 2020-03-27 2020-07-28 佛山科学技术学院 OCT blood vessel image dislocation correction method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
DAMIAN SIEDLECKI等: "Distortion Correction of OCT Images of the Crystalline Lens Gradient Index Approach", 《OPTOMETRY AND VISION SCIENCE》 *
VASILY MATKIVSKY等: "Determination and correction of aberrations in full field OCT using phase gradient autofocus by maximizing the likelihood function", 《JOURNAL OF BIOPHOTONICS》 *
张明蓉等: "OCT图像运动伪差校正算法的研究", 《中国医疗器械杂志》 *
高振玉等: "基于C扫描的视盘OCT图像运动伪差校正算法的研究", 《中国医疗器械杂志》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113349737A (en) * 2021-06-30 2021-09-07 深圳英美达医疗技术有限公司 Method for calibrating OCT image of intravascular dual-mode imaging system
CN113349737B (en) * 2021-06-30 2023-05-26 深圳英美达医疗技术有限公司 Calibration method of OCT (optical coherence tomography) image of intravascular dual-mode imaging system

Also Published As

Publication number Publication date
CN112509095B (en) 2021-05-07

Similar Documents

Publication Publication Date Title
US9955951B2 (en) Sensor coordinate calibration in an ultrasound system
CN102068281B (en) Processing method for space-occupying lesion ultrasonic images
US7787673B2 (en) Method and apparatus for airway detection and segmentation using 3D morphological operators
CN108537838B (en) Detection method for hip joint bony acetabulum angle
CN108765388B (en) Automatic segmentation method and system for esophageal endoscopic OCT (optical coherence tomography) image hierarchical structure
CN105488796A (en) Lung segmentation method
CN109829942A (en) A kind of automatic quantization method of eye fundus image retinal blood vessels caliber
Xu et al. Alignment of 3-D optical coherence tomography scans to correct eye movement using a particle filtering
US20210272291A1 (en) Method and computer program for segmentation of optical coherence tomography images of the retina
CN108257126A (en) The blood vessel detection and method for registering, equipment and application of three-dimensional retina OCT image
CN112509095B (en) Oct image dislocation correction method
US6748257B2 (en) Detection of ribcage boundary from digital chest image
Artaechevarria et al. Airway segmentation and analysis for the study of mouse models of lung disease using micro-CT
CN113012127A (en) Cardiothoracic ratio measuring method based on chest medical image
JP4091318B2 (en) X-ray CT system
CN111260673B (en) Visceral organ parenchyma segmentation method and device suitable for edge-breaking visceral organ radiography
CN110033496B (en) Motion artifact correction method for time sequence three-dimensional retina SD-OCT image
CN117274216A (en) Ultrasonic carotid plaque detection method and system based on level set segmentation
CN116277978A (en) Multimode bone joint digital 3D printing method
CN113012151B (en) OCT (optical coherence tomography) image correction method and system for SS-OCT operation navigation system
Nakamura et al. Automated segmentation and morphometric analysis of the human airway tree from multidetector CT images
Roy et al. MDL-IWS: multi-view deep learning with iterative watershed for pulmonary fissure segmentation
CN114494177A (en) IVOCT (in-vivo visual optical coherence tomography) branch blood vessel identification method by utilizing longitudinal section and withdrawal property
CN108269258B (en) Method and system for segmenting corneal structures from OCT corneal images
CN112085711A (en) Method for automatically tracking muscle pinnate angle by combining convolutional neural network and Kalman filtering

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
GR01 Patent grant
GR01 Patent grant