CN112509095B - Oct image dislocation correction method - Google Patents

Oct image dislocation correction method Download PDF

Info

Publication number
CN112509095B
CN112509095B CN202110165934.0A CN202110165934A CN112509095B CN 112509095 B CN112509095 B CN 112509095B CN 202110165934 A CN202110165934 A CN 202110165934A CN 112509095 B CN112509095 B CN 112509095B
Authority
CN
China
Prior art keywords
oct image
point
angle
dislocation
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.)
Active
Application number
CN202110165934.0A
Other languages
Chinese (zh)
Other versions
CN112509095A (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
And judging the position to be misplaced, wherein the corresponding gradient is the misplaced amplitude, if the misplaced amplitude is smaller than a threshold value T, rough judgment is carried out, and if not, fine judgment is carried out. 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 100002_DEST_PATH_IMAGE001
At a distance of n DEG each
Figure 631491DEST_PATH_IMAGE001
Subtracting to obtain a gradient
Figure 100002_DEST_PATH_IMAGE002
At maximum angle of gradient
Figure 100002_DEST_PATH_IMAGE003
Judging the position of the dislocation, wherein,
Figure 100002_DEST_PATH_IMAGE004
corresponding gradient
Figure 100002_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 100002_DEST_PATH_IMAGE006
Traversing the distance r to the origin, calculating the spacing distance at each angle as
Figure 100002_DEST_PATH_IMAGE007
Gray scale difference of the points
Figure 100002_DEST_PATH_IMAGE008
Distance from the point of maximum gradation gradient to the origin as gradation gradient
Figure 100002_DEST_PATH_IMAGE009
Namely:
Figure 100002_DEST_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 100002_DEST_PATH_IMAGE011
Figure 100002_DEST_PATH_IMAGE012
Wherein the content of the first and second substances,
Figure 100002_DEST_PATH_IMAGE013
is composed of
Figure 100002_DEST_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 100002_DEST_PATH_IMAGE015
Then carrying out fine judgment;
the fine judgment comprises the following steps: the traversal range is
Figure 100002_DEST_PATH_IMAGE016
2 DEG and/or
Figure 100002_DEST_PATH_IMAGE017
2 °, step size 0.1 °, at each angle
Figure 100002_DEST_PATH_IMAGE018
Distance r from the lower traverse to the center point in the range of
Figure 100002_DEST_PATH_IMAGE019
Step size of 0.5 pixels, per angle
Figure 889735DEST_PATH_IMAGE018
Distance from point with maximum lower gray gradient to center point
Figure 100002_DEST_PATH_IMAGE020
Namely:
Figure 100002_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 100002_DEST_PATH_IMAGE022
Comprises the following steps:
Figure 100002_DEST_PATH_IMAGE023
the corresponding fine misalignment magnitudes are:
Figure 100002_DEST_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 100002_DEST_PATH_IMAGE025
Angle of and
Figure 100002_DEST_PATH_IMAGE026
wherein:
Figure 100002_DEST_PATH_IMAGE027
Figure 100002_DEST_PATH_IMAGE028
wherein the content of the first and second substances,
Figure 100002_DEST_PATH_IMAGE029
is a unit vector;
calculating the dislocation amplitude of the point
Figure 100002_DEST_PATH_IMAGE030
If Δ r > 0, the angle is
Figure 42630DEST_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 DEST_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 370230DEST_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 DEST_PATH_IMAGE033
. Among them, the angle θ is preferably 2 °.
At a distance of n DEG each
Figure 425911DEST_PATH_IMAGE033
Subtracting to obtain a gradient
Figure DEST_PATH_IMAGE034
N ° is preferably 5 °, at the maximum angle of the gradient
Figure 773716DEST_PATH_IMAGE003
Judging the position of the dislocation, wherein,
Figure DEST_PATH_IMAGE035
corresponding gradient
Figure 228837DEST_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 DEST_PATH_IMAGE036
Traversing the distance r to the origin, calculating the spacing distance at each angle as
Figure 424326DEST_PATH_IMAGE007
Gray scale difference of the points
Figure 512892DEST_PATH_IMAGE008
As a result of the gradient of the gray scale,
Figure 692201DEST_PATH_IMAGE007
preferably 9 pixels, the distance from the point of maximum gray gradient to the origin
Figure DEST_PATH_IMAGE037
Namely:
Figure DEST_PATH_IMAGE038
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 DEST_PATH_IMAGE039
Figure DEST_PATH_IMAGE040
Wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE041
is composed of
Figure 75646DEST_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 125642DEST_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 991967DEST_PATH_IMAGE016
2 DEG and/or
Figure DEST_PATH_IMAGE042
2 °, step size 0.1 °, at each angle
Figure 642260DEST_PATH_IMAGE018
Distance r from the lower traverse to the center point in the range of
Figure DEST_PATH_IMAGE043
Step size of 0.5 pixels, per angle
Figure 329198DEST_PATH_IMAGE018
Distance from point with maximum lower gray gradient to center point
Figure DEST_PATH_IMAGE044
Namely:
Figure DEST_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 482968DEST_PATH_IMAGE032
Comprises the following steps:
Figure DEST_PATH_IMAGE046
the corresponding fine misalignment magnitudes are:
Figure DEST_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, two modes are provided for the direction selection, 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 each angle is calculated according to the lumen contour informationAnd dividing two types of angle sets according to the clustering information of the curvatures, wherein the direction of the set with the smaller number is taken 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 644828DEST_PATH_IMAGE025
Angle of and
Figure 736412DEST_PATH_IMAGE026
wherein:
Figure DEST_PATH_IMAGE048
Figure 743551DEST_PATH_IMAGE028
wherein the content of the first and second substances,
Figure 627193DEST_PATH_IMAGE029
is a unit vector and can be expressed as (0, 1).
Calculating the dislocation amplitude of the point
Figure DEST_PATH_IMAGE049
If Δ r > 0, the angle is
Figure 431726DEST_PATH_IMAGE026
Radius of
Figure 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 DEST_PATH_IMAGE051
wherein H is a pixel distribution histogram in a range between the outside of the catheter and the lumen contour,
Figure DEST_PATH_IMAGE052
for a certain pixel in the pixel distribution histogram H,
Figure DEST_PATH_IMAGE053
is the probability of the occurrence of the pixel,
Figure 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 712403DEST_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_IMAGE001
At a distance of n DEG each
Figure 506329DEST_PATH_IMAGE001
Subtracting to obtain a gradient
Figure DEST_PATH_IMAGE002
At maximum angle of gradient
Figure DEST_PATH_IMAGE003
Judging the position of the dislocation, wherein,
Figure DEST_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 DEST_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 DEST_PATH_IMAGE008
Distance from the point of maximum gradation gradient to the origin as gradation gradient
Figure DEST_PATH_IMAGE009
Namely:
Figure DEST_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 DEST_PATH_IMAGE012
Wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE013
is composed of
Figure 709953DEST_PATH_IMAGE009
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_IMAGE014
Then carrying out fine judgment;
the fine judgment comprises the following steps: the traversal range is
Figure DEST_PATH_IMAGE015
2 DEG and/or
Figure DEST_PATH_IMAGE016
2 °, step size 0.1 °, at each angle
Figure DEST_PATH_IMAGE017
Distance r from the lower traverse to the center point in the range of
Figure DEST_PATH_IMAGE018
Step size of 0.5 pixels, per angle
Figure 443685DEST_PATH_IMAGE017
Distance from point with maximum lower gray gradient to center point
Figure DEST_PATH_IMAGE019
Namely:
Figure DEST_PATH_IMAGE020
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_IMAGE021
Comprises the following steps:
Figure DEST_PATH_IMAGE022
the corresponding fine misalignment magnitudes are:
Figure DEST_PATH_IMAGE023
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_IMAGE024
Angle of and
Figure DEST_PATH_IMAGE025
wherein:
Figure DEST_PATH_IMAGE026
Figure DEST_PATH_IMAGE027
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE028
is a unit vector;
calculating the dislocation amplitude of the point
Figure DEST_PATH_IMAGE029
If Δ r > 0, the angle is
Figure 730572DEST_PATH_IMAGE025
Radius of
Figure DEST_PATH_IMAGE030
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_IMAGE031
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 672508DEST_PATH_IMAGE007
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 CN112509095A (en) 2021-03-16
CN112509095B true 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)

Families Citing this family (1)

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

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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

Family Cites Families (1)

* 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

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
Determination and correction of aberrations in full field OCT using phase gradient autofocus by maximizing the likelihood function;Vasily Matkivsky等;《Journal of Biophotonics》;20200708 *
Distortion Correction of OCT Images of the Crystalline Lens Gradient Index Approach;Damian Siedlecki等;《Optometry and Vision Science》;20120531;第89卷(第5期) *
OCT图像运动伪差校正算法的研究;张明蓉等;《中国医疗器械杂志》;20151231;第39卷(第1期) *
基于C扫描的视盘OCT图像运动伪差校正算法的研究;高振玉等;《中国医疗器械杂志》;20161231;第40卷(第2期) *

Also Published As

Publication number Publication date
CN112509095A (en) 2021-03-16

Similar Documents

Publication Publication Date Title
US10561403B2 (en) Sensor coordinate calibration in an ultrasound system
US10154820B2 (en) Method and apparatus for navigating CT scan with a marker
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
CN108257126B (en) Blood vessel detection and registration method, equipment and application of three-dimensional retina OCT image
CN109829942A (en) A kind of automatic quantization method of eye fundus image retinal blood vessels caliber
CN107481251A (en) A kind of method that terminal bronchi tree is extracted from lung CT image
US20210272291A1 (en) Method and computer program for segmentation of optical coherence tomography images of the retina
CN112509095B (en) Oct image dislocation correction method
CN111553849B (en) Cone beam CT geometric artifact removing method and device based on local feature matching
US6748257B2 (en) Detection of ribcage boundary from digital chest image
CN104361554A (en) Intravascular ultrasound image based automatic adventitia detection method
Artaechevarria et al. Airway segmentation and analysis for the study of mouse models of lung disease using micro-CT
Kiraly et al. 3D human airway segmentation for virtual bronchoscopy
CN113012127A (en) Cardiothoracic ratio measuring method based on chest medical image
JP4091318B2 (en) X-ray CT system
CN116277978B (en) Multimode bone joint digital 3D printing method
CN110033496B (en) Motion artifact correction method for time sequence three-dimensional retina SD-OCT image
CN115969400A (en) Apparatus for measuring area of eyeball protrusion
CN113012151B (en) OCT (optical coherence tomography) image correction method and system for SS-OCT operation navigation system
CN111260673B (en) Visceral organ parenchyma segmentation method and device suitable for edge-breaking visceral organ radiography
CN114820678A (en) Automatic extraction method and system for inner contour of blood vessel wall based on OCT image
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