CN112826443A - 一种眼底图像畸变的矫正方法以及共焦激光眼底成像仪 - Google Patents
一种眼底图像畸变的矫正方法以及共焦激光眼底成像仪 Download PDFInfo
- Publication number
- CN112826443A CN112826443A CN202011640430.1A CN202011640430A CN112826443A CN 112826443 A CN112826443 A CN 112826443A CN 202011640430 A CN202011640430 A CN 202011640430A CN 112826443 A CN112826443 A CN 112826443A
- Authority
- CN
- China
- Prior art keywords
- image
- pixel point
- fundus
- pixel
- original image
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 210000001525 retina Anatomy 0.000 claims abstract description 14
- 238000003384 imaging method Methods 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 6
- 201000010099 disease Diseases 0.000 description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 2
- 230000002207 retinal effect Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 208000003174 Brain Neoplasms Diseases 0.000 description 1
- 206010065559 Cerebral arteriosclerosis Diseases 0.000 description 1
- 206010008111 Cerebral haemorrhage Diseases 0.000 description 1
- 208000010412 Glaucoma Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 206010038933 Retinopathy of prematurity Diseases 0.000 description 1
- 206010064930 age-related macular degeneration Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 206010008118 cerebral infarction Diseases 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 210000003161 choroid Anatomy 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 201000005851 intracranial arteriosclerosis Diseases 0.000 description 1
- 208000017169 kidney disease Diseases 0.000 description 1
- 208000002780 macular degeneration Diseases 0.000 description 1
- 210000001328 optic nerve Anatomy 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000004256 retinal image Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/14—Arrangements specially adapted for eye photography
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Ophthalmology & Optometry (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明提供了一种眼底图像畸变的矫正方法以及共焦激光眼底成像仪,所述眼底图像畸变的矫正方法包括:S1:获取原图像中每一像素点对应的扫描光斑在视网膜上的光斑像的像高;S2:计算原图像中每一像素点对应的光斑像的累计像高差;S3:根据原图像中每一像素点对应的累计像高差,计算原图像中每一像素点在新图像中的坐标;S4:在新图像中,若相邻两个像素点之间的坐标距离大于1,则对所述相邻两个像素点之间进行双线性插值,以使得插值后的相邻两个像素点之间的坐标距离小于1。相较于现有技术,本发明眼底图像畸变矫正方法可以很好地修正眼底图像畸变,从而使得眼底图像可以准确地呈现眼底视网膜的大小,进而方便医护人员判断病灶的大小。
Description
技术领域
本发明涉及一种眼底图像畸变的矫正方法以及共焦激光眼底成像仪。
背景技术
眼底成像属于医学成像领域,用于获取人眼视网膜图像,以便医疗人员检查眼底病或者辅助医疗人员判断其它器官的病情。由于眼底的血管是人体唯一可以通过体表直接观察到的血管,医疗人员通过眼底相机可以检查眼底的视神经、视网膜、脉络膜以及屈光介质是否存在病变,同时还可以通过眼底相机协助对其它系统疾病进行诊断和病情判断,例如通过筛选视网膜照片检测脑梗塞、脑溢血、脑动脉硬化、脑肿瘤、糖尿病、肾病、高血压、早产儿视网膜病变、青光眼、老年化黄斑变性等。越早检测出这些疾病越有利于临床治疗,因此眼底照相机普遍用于临床筛查眼底疾病,成为不可或缺的医疗器械。
为了获取高清、大视场的眼底图像,通常采用共焦激光眼底成像仪来获取眼底图像。然而,通过共焦激光眼底成像仪获取的眼底图像会有部分畸变。产生图像畸变的原因包括:光学系统自身对图像造成的畸变、眼底视网膜曲面在成像时导致的图像畸变以及振镜扫描线速度的变化导致的图像畸变等等。其中,光学系统自身导致的图像畸变可以在安装前进行校准从而解决该问题;而眼底视网膜曲面导致的图像畸变,在设计光学系统时,已经考虑到该畸变了。因此这两者导致的图像畸变很小,可以忽略不计。然而,在成像系统中,由于振镜扫描的角速度不变,而视网膜却是曲面,从而使得入射光照在视网膜中心和边缘位置时,光程不相等,存在光程差。而该光程差导致振镜扫描时产生不同的线速,进而使得图像产生畸变。
鉴于上述问题,有必要提供一种新的眼底图像畸变的矫正方法,以解决上述问题。
发明内容
本发明的目的在于提供一种眼底图像畸变的矫正方法,该眼底图像畸变矫正方法可以很好地修正眼底图像畸变,从而使得眼底图像可以准确地呈现眼底视网膜的大小,进而方便医护人员判断病灶的大小。
为实现上述目的,本发明提供了一种眼底图像畸变的矫正方法,包括:S1:获取原图像中每一像素点对应的扫描光斑在视网膜上的光斑像的像高;S2:计算原图像中每一像素点对应的光斑像的累计像高差;S3:根据原图像中每一像素点对应的累计像高差,计算原图像中每一像素点在新图像中的坐标;S4:在新图像中,若相邻两个像素点之间的坐标距离大于1,则对所述相邻两个像素点之间进行双线性插值,以使得插值后的相邻两个像素点之间的坐标距离小于1。
作为本发明的进一步改进,所述像高为光斑像在横坐标上的像高。
作为本发明的进一步改进,所述原图像在横坐标方向上的像素点数为nRows,原图像中像素点P(i,j)在新图像中的对应像素点为P(x,y);其中,像素点P(x,y)的横坐标为Px,则
Px=δri*(nRows-1)/δr(nRows-1)。
作为本发明的进一步改进,原图像中像素点P(i+1,j+1)在新图像中的对应像素点为P(x+1,y+1);若|Px+1-Px|>1,则根据像素点P(i,j)的像素灰度值和像素点P(i+1,j+1)的像素灰度值,在所述像素点P(x,y)和像素点P(x+1,y+1)之间进行双线性插值,以使得新图像中在横坐标方向,相邻的两个像素点之间的横坐标之差不大于1。
作为本发明的进一步改进,将在所述像素点P(x,y)和像素点P(x+1,y+1)之间插入的像素点记为P(m,n);像素点P(x,y)的像素灰度值为H(x,y);像素点P(x+1,y+1)的像素灰度值为H(x+1,y+1);像素点P(m,n)的像素灰度值为H(m,n),横坐标为Pm则:
作为本发明的进一步改进,所述像高为光斑像在纵坐标上的像高。
作为本发明的进一步改进,所述原图像在纵坐标方向上的像素点数为nCols,原图像中像素点P(i,j)在新图像中的对应像素点为P(x,y);其中,像素点P(x,y)的纵坐标为Py,则
Py=δci*(nCols-1)/δc(nCols-1)。
作为本发明的进一步改进,原图像中像素点P(i+1,j+1)在新图像中的对应像素点为P(x+1,y+1);若|Py+1-Py|>1,则根据像素点P(i,j)的像素灰度值和像素点P(i+1,j+1)的像素灰度值,在所述像素点P(x,y)和像素点P(x+1,y+1)之间进行双线性插值,以使得新图像中在纵坐标方向,相邻的两个像素点之间的纵坐标之差不大于1。
作为本发明的进一步改进,将在所述像素点P(x,y)和像素点P(x+1,y+1)之间插入的像素点记为P(m,n);像素点P(x,y)的像素灰度值为H(x,y);像素点P(x+1,y+1)的像素灰度值为H(x+1,y+1);像素点P(m,n)的像素灰度值为H(m,n),纵坐标为Pn则:
本发明还揭示了一种共焦激光眼底成像仪,所述共焦激光眼底成像仪实施前述眼底图像畸变的矫正方法。
本发明的有益效果是:本发明眼底图像畸变矫正方法可以很好地修正眼底图像畸变,从而使得眼底图像可以准确地呈现眼底视网膜的大小,进而方便医护人员判断病灶的大小。
附图说明
图1是本发明眼底图像畸变的矫正方法的流程示意图。
图2是振镜转动角度与光斑像的坐标位置的关系图。
图3是原图像映射到新图像中的关系示意图。
图4是本发明共焦激光眼底成像仪的模块示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面结合附图和具体实施例对本发明进行详细描述。
请参阅图1所示,本发明揭示了一种眼底图像畸变的矫正方法,包括如下步骤:
S1:获取原图像中每一像素点对应的扫描光斑在视网膜上的光斑像的像高。
首先,根据共焦激光眼底成像系统以及眼底模型,获取振镜的转动角度以及扫描光斑在视网膜上的光斑像的坐标位置,如图2所示。然后,根据振镜的转动角度以及光斑像的坐标位置,获取原图像中每一像素点对应的光斑像的像高。其中,原图像是共焦激光成像系统获取的眼底图像。所述像高包括光斑像在横坐标上的像高以及光斑像在纵坐标上的像高。
S2:计算原图像中每一像素点对应的光斑像的累计像高差。
S3:根据原图像中每一像素点对应的累计像高差,计算原图像中每一像素点在新图像中的坐标。
所述新图像是指对所述原图像进行畸变矫正后的图像。所述原图像在横坐标方向上的像素点数目为nRows,在纵坐标上的像素点数目为nCols。原图像中像素点P(i,j)在新图像中的对应像素点为P(x,y)。像素点P(x,y)的横坐标为Px,纵坐标为Py。其中,
Px=δri*(nRows-1)/δr(nRows-1);
Py=δci*(nCols-1)/δc(nCols-1)。
S4:在新图像中,若相邻两个像素点之间的坐标距离大于1,则对所述相邻两个像素点之间进行双线性插值,以使得插值后的相邻两个像素点之间的坐标距离小于1。
请参阅图3所示,当修正原图像中的畸变从而得到新图像后,原图像中相邻的两个像素点P(i,j)和P(i+1,j)在新图像中的对应像素点为P(x,y)和P(x+1,y)。由于振镜扫描线速度的变化,从而导致原图像发生畸变,因此修正后的新图像中,像素点P(x,y)和P(x+1,y)则不一定相邻,即:像素点P(x,y)和P(x+1,y)在横坐标方向上的距离或者纵坐标方向上的距离可能大于1个像素。同理,原图像中相邻的两个像素点P(i,j)和P(i,j+1)在新图像中的对应像素点为P(x,y)和P(x,y+1),其中像素点为P(x,y)和P(x,y+1)也不一定相邻。
具体地,原图像中像素点P(i+1,j)在新图像中的对应像素点为P(x+1,y);若|Px+1-Px|>1,则根据像素点P(i,j)的像素灰度值和像素点P(i+1,j)的像素灰度值,在所述像素点P(x,y)和像素点P(x+1,y)之间进行双线性插值,以使得新图像中在横坐标方向,相邻的两个像素点之间的横坐标之差不大于1。原图像中像素点P(i,j+1)在新图像中的对应像素点为P(x,y+1);若|Py+1-Py|>1,则根据像素点P(i,j)的像素灰度值和像素点P(i,j+1)的像素灰度值,在所述像素点P(x,y)和像素点P(x,y+1)之间进行双线性插值,以使得新图像中在纵坐标方向,相邻的两个像素点之间的纵坐标之差不大于1。
进一步地,将在所述像素点P(x,y)和像素点P(x+1,y)之间插入的像素点记为P(m,y)。像素点P(x,y)的像素灰度值为H(x,y);像素点P(x+1,y)的像素灰度值为H(x+1,y);像素点P(m,y)的像素灰度值为H(m,y),横坐标为Pm,纵坐标为Pn。其中,
进一步地,将在所述像素点P(x,y)和像素点P(x,y+1)之间插入的像素点记为P(x,n);像素点P(x,y)的像素灰度值为H(x,y);像素点P(x,y+1)的像素灰度值为H(x,y+1);像素点P(x,n)的像素灰度值为H(x,n),纵坐标为Pn则:
如此设置,可以有效消除原图像中的畸变,从而使得新图像中的病灶的各部分与视网膜上的真实病灶之间的比例基本一致,从而有利于医护人员根据新图像判断病灶的大小,进而避免因图像失真而使得医护人员出现判断失误的问题。
相较于现有技术,本发明眼底图像畸变矫正方法可以很好地修正眼底图像畸变,从而使得眼底图像可以准确地呈现眼底视网膜的大小,进而方便医护人员判断病灶的大小。
请参阅图4所示,本发明还揭示了一种共焦激光眼底成像仪100,包括共焦激光成像模块10以及图像处理模块20。所述图像处理模块20根据所述眼底图像畸变矫正方法对所述共焦激光成像模块20获取的原图像进行图像处理,从而获得修正畸变后的新图像。
以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围。
Claims (12)
1.一种眼底图像畸变的矫正方法,其特征在于,包括:
S1:获取原图像中每一像素点对应的扫描光斑在视网膜上的光斑像的像高;
S2:计算原图像中每一像素点对应的光斑像的累计像高差;
S3:根据原图像中每一像素点对应的累计像高差,计算原图像中每一像素点在新图像中的坐标;
S4:在新图像中,若相邻两个像素点之间的坐标距离大于1,则对所述相邻两个像素点之间进行双线性插值,以使得插值后的相邻两个像素点之间的坐标距离小于1。
2.如权利要求1所述的眼底图像畸变的矫正方法,其特征在于:所述像高为光斑像在横坐标上的像高。
4.如权利要求3所述的眼底图像畸变的矫正方法,其特征在于:所述原图像在横坐标方向上的像素点数为nRows,原图像中像素点P(i,j)在新图像中的对应像素点为P(x,y);其中,像素点P(x,y)的横坐标为Px,则
Px=δri*(nRows-1)/δr(nRows-1)。
5.如权利要求4所述的眼底图像畸变的矫正方法,其特征在于:原图像中像素点P(i+1,j)在新图像中的对应像素点为P(x+1,y);若|Px+1-Px|>1,则根据像素点P(i,j)的像素灰度值和像素点P(i+1,j)的像素灰度值,在所述像素点P(x,y)和像素点P(x+1,y)之间进行双线性插值,以使得新图像中在横坐标方向,相邻的两个像素点之间的横坐标之差不大于1。
7.如权利要求1所述的眼底图像畸变的矫正方法,其特征在于:所述像高为光斑像在纵坐标上的像高。
9.如权利要求8所述的眼底图像畸变的矫正方法,其特征在于:所述原图像在纵坐标方向上的像素点数为nCols,原图像中像素点P(i,j)在新图像中的对应像素点为P(x,y);其中,像素点P(x,y)的纵坐标为Py,则
Py=δci*(nCols-1)/δc(nCols-1)。
10.如权利要求9所述的眼底图像畸变的矫正方法,其特征在于:原图像中像素点P(i,j+1)在新图像中的对应像素点为P(x,y+1);若|Py+1-Py|>1,则根据像素点P(i,j)的像素灰度值和像素点P(i,j+1)的像素灰度值,在所述像素点P(x,y)和像素点P(x,y+1)之间进行双线性插值,以使得新图像中在纵坐标方向,相邻的两个像素点之间的纵坐标之差不大于1。
12.一种共焦激光眼底成像仪,其特征在于:实施如权利要求1~11中任意一项所述的眼底图像畸变的矫正方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011640430.1A CN112826443B (zh) | 2020-12-31 | 2020-12-31 | 一种眼底图像畸变的矫正方法以及共焦激光眼底成像仪 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011640430.1A CN112826443B (zh) | 2020-12-31 | 2020-12-31 | 一种眼底图像畸变的矫正方法以及共焦激光眼底成像仪 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112826443A true CN112826443A (zh) | 2021-05-25 |
CN112826443B CN112826443B (zh) | 2022-02-22 |
Family
ID=75926947
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011640430.1A Active CN112826443B (zh) | 2020-12-31 | 2020-12-31 | 一种眼底图像畸变的矫正方法以及共焦激光眼底成像仪 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112826443B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07271961A (ja) * | 1994-03-30 | 1995-10-20 | Canon Inc | 画像処理装置 |
US20080204656A1 (en) * | 2007-02-22 | 2008-08-28 | Kowa Company Ltd. | Image processing apparatus |
CN103908222A (zh) * | 2013-01-08 | 2014-07-09 | 华晶科技股份有限公司 | 眼底图像检测装置及获取动物眼睛底部图像的方法 |
WO2020102734A1 (en) * | 2018-11-16 | 2020-05-22 | Arizona Board Of Regents On Behalf Of The University Of Arizona | Identification and control of myopic progression using distortion |
-
2020
- 2020-12-31 CN CN202011640430.1A patent/CN112826443B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07271961A (ja) * | 1994-03-30 | 1995-10-20 | Canon Inc | 画像処理装置 |
US20080204656A1 (en) * | 2007-02-22 | 2008-08-28 | Kowa Company Ltd. | Image processing apparatus |
CN103908222A (zh) * | 2013-01-08 | 2014-07-09 | 华晶科技股份有限公司 | 眼底图像检测装置及获取动物眼睛底部图像的方法 |
WO2020102734A1 (en) * | 2018-11-16 | 2020-05-22 | Arizona Board Of Regents On Behalf Of The University Of Arizona | Identification and control of myopic progression using distortion |
Non-Patent Citations (1)
Title |
---|
JIAPAN GUO ET AL.: "Automatic Determination of Vertical Cup-to-Disc Ratio in Retinal Fundus Images for Glaucoma Screening", 《IEEE ACCESS》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112826443B (zh) | 2022-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7322698B2 (en) | Method for judging changes in images of the eye or its component parts | |
EP1976424B1 (en) | A method of eye examination by optical coherence tomography | |
US8636364B2 (en) | Systems and methods for widefield mapping of the retina | |
EP2360641B1 (en) | Image processing method and image processing device | |
Seo et al. | Measurement of ocular torsion using digital fundus image | |
US8184149B2 (en) | Ophthalmic apparatus and method for increasing the resolution of aliased ophthalmic images | |
CN112826443B (zh) | 一种眼底图像畸变的矫正方法以及共焦激光眼底成像仪 | |
CN114092405A (zh) | 一种针对黄斑水肿oct图像的视网膜层自动分割方法 | |
US11642019B2 (en) | Ophthalmic device | |
CN112426333A (zh) | 一种偏盲患者立体视觉辅助电子设备及其控制方法 | |
Fu et al. | A GAN-based deep enhancer for quality enhancement of retinal images photographed by a handheld fundus camera | |
CN103919525A (zh) | 拼接图像的方法及摄像系统 | |
CN112669357B (zh) | 一种眼底图像合成方法以及眼底成像仪 | |
CN115690183A (zh) | 一种图像配准数据处理方法及系统 | |
TWI524876B (zh) | 拼接影像的方法及攝像系統 | |
CN113384236A (zh) | 一种眼部光学相干断层扫描的方法及其装置 | |
CN111652805A (zh) | 一种用于眼底图像拼接的图像预处理方法 | |
CN117495864B (zh) | 基于图像处理的影动方向计算系统及屈光度估计系统 | |
Xian et al. | Image registration based on camera calibration for dual-wavelength retinal oximetry | |
Bernardes et al. | Multimodal functional and morphological nonrigid image registration | |
Mujat et al. | Non-Rigid Registration for High-Resolution Retinal Imaging. Diagnostics 2023, 13, 2285 | |
CN118864384A (zh) | 基于oct图像的目标区域识别方法和系统 | |
CN110660035A (zh) | 一种显微系统图像畸变矫正方法 | |
Ho et al. | Dacryomeniscometry: an early attempt | |
CN117176868A (zh) | 高精度扫描成像系统 |
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 |