发明名称 用于检测非球面的高精度位相型计算全息图的编码方法
摘要 用于检测非球面的高精度位相型计算全息图的编码方法,涉及计算全息技术领域,消除了计算全息图在全息编码过程中产生的编码误差,阻断该误差从源头逐级传递和放大,提高计算全息图再现基准波面的精度。包括:对连续非球面数学模型进行抽样离散化,再以2π为模对其进行周期量化获得位相图,通过放大取整将其转换为灰阶图,生成位相型计算全息图;提取全息图每个量化周期的最小边界值,建立线性方程组,求最小边界值偏差和每个量化周期的调整系数;通过非球面量化周期边界的定义,提取每个量化周期数据;用最小边界值偏差和调整系数,调整对应量化周期数据;累加每个量化周期数据,生成校正后位相型计算全息图。本发明适用于计算全息法检测非球面。
申请公布号 CN103791855A 申请公布日期 2014.05.14
申请号 CN201410062280.9 申请日期 2014.02.24
申请人 哈尔滨理工大学 发明人 张洪鑫;燕宇;乔玉晶;司俊山;马薇
分类号 G01B11/24(2006.01)I;G03H1/08(2006.01)I 主分类号 G01B11/24(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张宏威
主权项 1.用于检测非球面的高精度位相型计算全息图的编码方法,其特征是:它由以下步骤实现:步骤一、根据被测旋转对称非球面的二次曲面偏心率e,确定旋转对称非球面的类型;然后根据公式:<maths num="0001"><![CDATA[<math><mrow><mi>W</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msup><mi>cS</mi><mn>2</mn></msup><mrow><mn>1</mn><mo>+</mo><msqrt><mn>1</mn><mo>-</mo><mrow><mo>(</mo><mi>L</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><msup><mi>c</mi><mn>2</mn></msup><msup><mi>S</mi><mn>2</mn></msup></msqrt></mrow></mfrac><mo>+</mo><msub><mi>A</mi><mn>1</mn></msub><msup><mi>S</mi><mn>4</mn></msup><mo>+</mo><msub><mi>A</mi><mn>2</mn></msub><msup><mi>S</mi><mn>6</mn></msup><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msub><mi>A</mi><mi>n</mi></msub><msup><mi>S</mi><mrow><mn>2</mn><mi>n</mi></mrow></msup></mrow></math>]]></maths>生成被测旋转对称非球面模型W(x,y),单位为微米(um);式中:S<sup>2</sup>=x<sup>2</sup>+y<sup>2</sup>;(x,y)为被测旋转对称非球面上任意一点在底面投影的坐标;<img file="FDA0000468970870000012.GIF" wi="196" he="160" />R<sub>0</sub>为被测旋转对称非球面的顶点曲率半径;A<sub>1</sub>、A<sub>2</sub>、LA<sub>n</sub>为高次非球面各高次项系数;L=-e<sup>2</sup>;步骤二、对被测旋转对称非球面模型W(x,y)采用位相型计算全息编码;具体为:首先对被测旋转对称非球面模型W(x,y)进行抽样离散化处理,获得W(i,j);<maths num="0002"><![CDATA[<math><mrow><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><mi>W</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>comb</mi><mrow><mo>(</mo><mfrac><mi>x</mi><mi>N</mi></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>comb</mi><mfrac><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mi>N</mi></mfrac><mo>&CenterDot;</mo><mi>rect</mi><mrow><mo>(</mo><mfrac><mi>x</mi><mi>D</mi></mfrac><mo>,</mo><mfrac><mi>y</mi><mi>D</mi></mfrac><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>式中:i,j=1,2,3LN,N为采样数,D为被测旋转对称非球面的口径;然后以2π为模对W(i,j)进行周期量化,获得位相图,通过放大取整后,将位相图转换为灰阶图,灰阶分布范围为0-255,生成位相型计算全息图W<sub>CGH</sub>(i,j);<maths num="0003"><![CDATA[<math><mrow><msub><mi>W</mi><mi>CGH</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>round</mi><mo>[</mo><mn>255</mn><mo>&CenterDot;</mo><mfrac><mn>1</mn><mi>&lambda;h</mi></mfrac><mo>&CenterDot;</mo><mi>mod</mi><msub><mrow><mo>[</mo><mi>W</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>h</mi><mo>]</mo></mrow><mrow><mn>2</mn><mi>&pi;</mi></mrow></msub><mo>]</mo><mrow><mo>(</mo><mi>h</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>&lambda;</mi></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>式中:λ为波长,单位为微米(um);步骤三、提取步骤二获得的位相型计算全息图W<sub>CGH</sub>(i,j)上包含对称中心点的一维数组B<sub>i</sub>,即:B<sub>1</sub>、B<sub>2</sub>……B<sub>N</sub>,从一维数组B<sub>i</sub>中提取每个量化周期的最小边界值b<sub>k</sub>,建立线性方程组:z<sub>k</sub>=b<sub>k</sub>-b<sub>0</sub>并求解最小边界值偏差z<sub>k</sub>;式中:b<sub>0</sub>为期望值,k为量化周期数,k=1,2,3…C;C是最大量化周期数;步骤四、建立线性方程组:<maths num="0004"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>k</mi></msub><mo>=</mo><mfrac><msub><mi>a</mi><mn>0</mn></msub><mrow><msub><mi>a</mi><mi>k</mi></msub><mo>-</mo><msub><mi>z</mi><mi>k</mi></msub></mrow></mfrac><mo>-</mo><msub><mi>b</mi><mn>0</mn></msub></mrow></math>]]></maths>求取:每个量化周期的调整系数β<sub>k</sub>;式中:a<sub>k</sub>为各周期最大边界值,a<sub>0</sub>为期望值;步骤五、通过非球面量化周期的边界定义式:<img file="FDA0000468970870000022.GIF" wi="1423" he="230" />提取位相型计算全息图W<sub>CGH</sub>(i,j)每个量化周期数据W<sub>k</sub>(i,j);式中:λ为波长,单位为微米(um);步骤六、用一维数组调整系数β<sub>k</sub>通过公式:R<sub>k</sub>(i,j)=(1+β<sub>k</sub>)[W<sub>k</sub>(i,j)-z<sub>k</sub>(i,j)]调整对应量化周期数据W<sub>k</sub>(i,j),得到调整后量化周期数据R<sub>k</sub>(i,j);其中:最小边界值偏差矩阵z<sub>k</sub>(i,j)由一维数组z<sub>k</sub>构造:<img file="FDA0000468970870000023.GIF" wi="1134" he="198" />再累加调整后每个量化周期数据:<maths num="0005"><![CDATA[<math><mrow><msub><mi>W</mi><mrow><mi>CGH</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>round</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>C</mi></munderover><msub><mi>R</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths>生成校正后位相型计算全息图W<sub>CGH1</sub>(i,j);步骤七、判断步骤六生成的校正后位相型计算全息图W<sub>CGH1</sub>(i,j)的最小边界值偏差是否满足|z<sub>k</sub>|&lt;z<sub>min</sub>,z<sub>min</sub>为最小边界值偏差期望值,如果判断结果为否,则返回执行步骤三;如果判断结果为是,则完成高精度位相型计算全息图的编码。
地址 150080 黑龙江省哈尔滨市南岗区学府路52号