体积OCT图像数据处理制造技术

本申请涉及体积OCT图像数据处理。一种处理C扫描数据以生成校正数据的方法，该C扫描数据包括由OCT成像系统获取的对成像目标的B扫描序列，该校正数据用于补偿由OCT成像系统和成像目标的相对运动引起的在该序列中的B扫描之间的轴向位移，该相对运动在B扫描的获取期间使OCT成像系统和成像目标之间的距离变化，该校正数据是通过以下方式生成的：对于该序列中多对相邻B扫描中的每个相邻B扫描对，确定(S10)在相邻B扫描中的共同眼部特征的相应表示之间的轴移的相应指示符；以及根据所确定的指示符随着在该序列中相对应的相邻B扫描对的位置的变化，确定(S20)该变化的第一频率分量，其指示B扫描的获取期间的相对运动。其指示B扫描的获取期间的相对运动。其指示B扫描的获取期间的相对运动。

【技术实现步骤摘要】
体积OCT图像数据处理


[0001]本文的示例方面总体上涉及眼科光学相干断层扫描(OCT)成像系统的领域，并且更具体地，涉及一种用于处理由OCT成像系统生成的C扫描数据以生成校正数据的方法和装置，该校正数据用于补偿由OCT成像系统和成像目标的相对运动引起的获取的B扫描序列中的B扫描之间的轴向位移，该相对运动在OCT成像系统获取B扫描期间使OCT成像系统和成像目标之间的距离变化。
[0002]背景
[0003]对受试者眼睛的一部分(例如视网膜)的典型体积光学相干断层扫描(OCT)扫描(还称为C扫描)的获取可能需要大约1秒钟到5秒钟。在此期间，OCT成像系统在第一扫描方向上重复扫描OCT光束，以记录(二维)B扫描序列，每个B扫描包括一系列轴向扫描(A扫描)，该一系列轴向扫描是沿着第一扫描方向在视网膜表面的相应点处被记录的。B扫描序列中的B扫描通常在垂直于第一方向的方向上排列。在获取C扫描的期间，眼睛可能通常会由于受试者的不自主运动而轴向(沿着OCT光束的方向)移动。为了成功渲染C扫描图像，并确保基于C扫描图像执行的后续测量的准确性，可能有必要校正由本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种处理C扫描数据以生成校正数据的方法，所述C扫描数据包括由光学相干断层扫描OCT成像系统(30)获取的对成像目标(20)的B扫描(210，220)序列，所述校正数据用于补偿由所述OCT成像系统(30)和所述成像目标(20)的相对运动引起的在所述B扫描序列中的B扫描(210，220)之间的轴向位移(d)，所述相对运动在由所述OCT成像系统(30)获取所述B扫描(210，220)期间使所述OCT成像系统和所述成像目标之间的距离变化，所述方法包括通过以下步骤生成所述校正数据：对于所述序列中多对相邻B扫描中的每个相邻B扫描(210，220)对，确定(S10)所述相邻B扫描(210，220)中的共同眼部特征(230)的相应表示之间的轴移(d)的相应指示符(510)；和根据所确定的指示符(510)随着所述序列中的相对应的相邻B扫描(210，220)对的位置的变化，确定(S20)所述变化的第一频率分量，所述第一频率分量指示在由所述OCT成像系统(30)获取所述B扫描(210，220)期间所述OCT成像系统(30)和所述成像目标(20)的所述相对运动。2.根据权利要求1所述的方法，其中，对于每个所述相邻B扫描(210，220)对，所述轴移(d)的所述相应指示符(510)由以下步骤来确定：计算所述相邻B扫描(210，220)对之间的互相关，并且确定对应于所计算的互相关中的峰值的所述B扫描(210，220)之间的偏移作为所述指示符(510)，或者识别所述相邻B扫描(210，220)对的B扫描(210，220)中的共同眼部特征的相应位置，并且确定所识别的位置之间沿着所述B扫描(210，220)的轴的位移，所述B扫描的轴对应于所述成像系统(30)的轴向方向。3.根据权利要求1或权利要求2所述的方法，其中，所述成像目标(20)具有曲率，并且所述方法还包括：确定所述变化的第二频率分量，所述第二频率分量指示所述成像目标(20)的曲率。4.根据权利要求3所述的方法，其中：通过将m阶多项式拟合(S210)到所确定的指示符(510)的变化来确定所述第二频率分量，以及通过以下方式来确定所述第一频率分量：从所述确定的指示符(510)的变化中的所述指示符(510)中减去(S220)所述m阶多项式的值以生成所述指示符(510)的校正变化，并且将n阶多项式拟合到所述指示符(510)的所述校正变化，其中，m和n是整数，并且m小于n。5.根据任一前述权利要求所述的方法，还包括通过将基于所述第一频率分量的偏移应用于所述B扫描(210，220)序列中的B扫描(210，220)，使用所述校正数据来补偿由所述OCT成像系统(30)和所述成像目标(20)的所述相对运动引起的在所述B扫描(210，220)序列中的B扫描(210，220)之间的轴向位移(d)。6.根据权利要求4所述的方法，其中，通过进一步执行包括以下步骤的过程的至少两次迭代来确定所述第一频率分量：(i)计算(S240)多个残值，每个残值被计算为指示符(510)的变化中的指示符(510)和n阶多项式的对应值之间的差值；(ii)确定(S250)所述多个残值是否包括超出第一正阈值或低于第一负阈值的异常值；(iii)在确定所述多个残值包括所述异常值的情况下，从指示符(510)的所述变化中去
除(S260)对应于所述异常值的所述指示符(510)以生成指示符的更新变化，并且在确定所述残值不包括所述异常值的情况下，确定所述n阶多项式作为所述第一频率分量并结束所述过程；和(iv)将所述n阶多项式拟合(S270)到指示符(510)的所述更新变化，其中，在所述过程的第一次迭代中，所述多个残值中的每个残值被计算为所述指示符的所述校正变化中的指示符和拟合到所述指示符的所述校正变化的所述n阶多项式的对应值之间的差值，并且其中，在所述过程的剩余的一次或更多次迭代中的每次迭代中，所述多个残值中的每个残值被计算为在所述过程的前一次迭代中生成的指示符的所述更新变化中的指示符和拟合到所述过程的所述前一次迭代中生成的指示符的所述更新变化的所述n阶多项式的对应值之间的差值。7.根据权利要求6所述的方法，还包括：确定(S310)所述多个残值中具有大于第二正阈值或小于第二负阈值的量值的残值的数量，其中，所述第二正阈值小于所述第一正阈值，并且所述第二负阈值大于所述第一负阈值；在所确定的残值的数量小于第三阈值的情况下，通过将基于所述第一频率分量的偏移应用于所述B扫描(210，220)序列中的B扫描(210，220)来补偿(S30)所述B扫描(210，220)序列中的B扫描(210，220)之间的轴向位移(d)；和在所确定的残值的数量不...

【专利技术属性】
技术研发人员：普拉文，
申请(专利权)人：奥普托斯股份有限公司，
类型：发明
国别省市：