发明名称 修正PS-OCT的测量数据的程序及安装该程序的PS-OCT系统
摘要 将PS‑OCT中测量的数据非线性修正,提高PS‑OCT的定量解析能力,能够进行包含病变部位的病期在内的准确的定量诊断,可以作为用于计算机诊断的有效手段使用。即使在PS‑OCT(1)中的光程差包含误差而成为噪音,其分布既不为正态分布、也不在真值周围对称分布的情况下,也能使用通过用蒙特卡洛仿真解析噪音特性所获得的分布变换函数变换测量数据,由此去除系统性误差,推测被噪音掩盖的真值,将PS‑OCT(1)的图像修正得更清晰。
申请公布号 CN103814286B 申请公布日期 2016.10.05
申请号 CN201280034132.5 申请日期 2012.04.25
申请人 国立大学法人筑波大学 发明人 安野嘉晃;段炼;伊藤雅英
分类号 G01N21/17(2006.01)I 主分类号 G01N21/17(2006.01)I
代理机构 北京银龙知识产权代理有限公司 11243 代理人 曾贤伟;曹鑫
主权项 一种PS‑OCT系统,具备PS‑OCT和处理由该PS‑OCT获得的测量数据的装置,所述PS‑OCT系统的特征在于,所述装置作为以下单元发挥作用:预先将由分别关于0~π的多个相位差的分布组成的多个相位差的分布的组针对每个不同的有效信噪比生成多个组并存储在存储装置中的单元;确定将由PS‑OCT测量的相位差的分布数据变换为对称的相位差的分布数据的非线性变换函数的系数的单元;根据由PS‑OCT测量的相位差的分布数据估算有效信噪比,并从针对所述每个不同的有效信噪比生成的相位差的分布的多个组中选择与该估算的有效信噪比对应的相位差的分布的组的单元;在根据所述选择的相位差的分布的组确定的系数下,通过所述变换函数将由PS‑OCT测量的相位差的分布数据变换为对称的相位差的分布数据的单元;以及通过对所述对称的相位差的分布数据进行平均来计算真的相位差的单元,所述非线性变换函数为下面的式(2),确定所述变换函数的系数的单元的结构为,针对所述非线性变换函数,从0~π为止每隔π/m改变,生成下面的m+1的联立方程式即式(3),求解该联立方程式来制作变换函数的系数b<sub>i</sub>的表格,其中,式(2)中,δ<sub>E</sub>为具有对称分布的相位差,δ<sub>M</sub>为由PS‑OCT测量的相位差,b<sub>i</sub>为变换函数的系数,γ为有效信噪比,n为由PS‑OCT测量的相位差的个数,其中,i=0、1、2、3、......n,式(3)中的m是划分π的次数,<>的意思是表示由多次蒙特卡洛仿真得到的值的平均值,式(2)<maths num="0001"><math><![CDATA[<mrow><msub><mi>&delta;</mi><mi>E</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msub><mi>b</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><msubsup><mi>&delta;</mi><mi>M</mi><mi>i</mi></msubsup><mo>=</mo><msub><mi>b</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><msub><mi>&delta;</mi><mi>M</mi></msub><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><msubsup><mi>&delta;</mi><mi>M</mi><mn>2</mn></msubsup><mo>+</mo><mo>...</mo><mo>+</mo><msub><mi>b</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><msubsup><mi>&delta;</mi><mi>M</mi><mi>n</mi></msubsup></mrow>]]></math><img file="FDA0001001932520000011.GIF" wi="1581" he="198" /></maths>式(3)<maths num="0002"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mrow><mn>0</mn><mo>=</mo><msub><mi>b</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>&delta;</mi><mrow><mi>M</mi><mo>,</mo><mn>0</mn></mrow></msub><mo>></mo><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&lt;</mo><msubsup><mi>&delta;</mi><mrow><mi>M</mi><mo>,</mo><mn>0</mn></mrow><mn>2</mn></msubsup><mo>></mo><mo>+</mo><mo>...</mo><mo>+</mo><msub><mi>b</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&lt;</mo><msubsup><mi>&delta;</mi><mrow><mi>M</mi><mo>,</mo><mn>0</mn></mrow><mi>n</mi></msubsup><mo>></mo></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mi>&pi;</mi><mi>m</mi></mfrac><mo>=</mo><msub><mi>b</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>&delta;</mi><mrow><mi>M</mi><mo>,</mo><mi>&pi;</mi><mo>/</mo><mi>m</mi></mrow></msub><mo>></mo><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&lt;</mo><msubsup><mi>&delta;</mi><mrow><mi>M</mi><mo>,</mo><mi>&pi;</mi><mo>/</mo><mi>m</mi></mrow><mn>2</mn></msubsup><mo>></mo><mo>+</mo><mo>...</mo><mo>+</mo><msub><mi>b</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&lt;</mo><msubsup><mi>&delta;</mi><mrow><mi>M</mi><mo>,</mo><mi>&pi;</mi><mo>/</mo><mi>m</mi></mrow><mi>n</mi></msubsup><mo>></mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow><mi>&pi;</mi><mo>=</mo><msub><mi>b</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>&delta;</mi><mrow><mi>M</mi><mo>,</mo><mi>&pi;</mi></mrow></msub><mo>></mo><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&lt;</mo><msubsup><mi>&delta;</mi><mrow><mi>M</mi><mo>,</mo><mi>&pi;</mi></mrow><mn>2</mn></msubsup><mo>></mo><mo>+</mo><mo>...</mo><mo>+</mo><msub><mi>b</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>&lt;</mo><msubsup><mi>&delta;</mi><mrow><mi>M</mi><mo>,</mo><mi>&pi;</mi></mrow><mi>n</mi></msubsup><mo>></mo></mrow></mtd></mtr></mtable></mfenced><mo>.</mo></mrow>]]></math><img file="FDA0001001932520000021.GIF" wi="1644" he="360" /></maths>。
地址 日本茨城县