CN104236520A - 一种整层大气湍流倾斜等晕角的测量方法 - Google Patents

Abstract

本发明公开了一种倾斜等晕角的测量方法，涉及光波大气传输领域，在接收双孔满足一定条件下，双孔接收到的平面波相对光强起伏协方差近似是湍流二阶矩的函数，通过测量相对光强起伏协方差来得出湍流二阶矩，再通过湍流二阶矩计算倾斜等晕角。通过对恒星的观测，利用本方法可以不分昼夜连续测量整层大气倾斜等晕角。本发明方法与广泛使用的恒星闪烁测量等晕角有异曲同工之处，都是通测量湍流矩来间接测量所需量；然而本发明方法的理论误差远小于等晕角测量的恒星闪烁法。本发明方法简单易行，测量精度高，在自适应光学、高分辨率、大视场成像等领域具有重要意义。

Description

一种整层大气湍流倾斜等晕角的测量方法

技术领域

本发明属于光波大气传输领域，涉及大气光学湍流参数的测量，具体为一种整层大气倾斜等晕角的测量方法。

背景技术

光波在大气中传输时，由于大气光学湍流(大气折射率的随机起伏)的影响，光波的波面会发生畸变、强度会产生起伏。倾斜非等晕性指两光源的倾斜抖动的差别。倾斜等晕角，也称像运动等晕角，表征倾斜非等晕性的大小，在这个角度范围内大气湍流造成的倾斜抖动(像运动)可认为基本相同，是自适应光学、激光大气传输、湍流成像中的一个重要参量。例如，装有激光导引星自适应光学系统的地基望远镜的长曝光角分辨率受大气湍流引起的倾斜项误差校正精度影响，而倾斜项校正精度与倾斜参考星与目标的角间距θ以及倾斜等晕角θ_ΤΑ有关；有了θ_ΤΑ，我们就可以评估望远镜的分辨性能。

Sasiela和Shelton在1993年(Sasiela and Shelton，“Transverse spectral filtering and Mellintransform methods applied to the effect of outer scale on tilt and tilt anisoplantism”.J.Opt.Soc.Am.A10，646-660，1993)给出了倾斜等晕角的表达式：

${\mathrm{\θ}}_{\mathrm{TA}}=\frac{{0.195\mathrm{\λD}}^{1/6}}{{\mathrm{\μ}}_2^{1/2}}$ Equation Chapter(Next)Section 1Equation Section 2(1)

由倾斜非等晕性引起的均方波前相位误差可表示为

${\mathrm{\σ}}_{\mathrm{TA}}^2={\left(\frac{\mathrm{\θ}}{{\mathrm{\θ}}_{\mathrm{TA}}}\right)}^2---\left(2\right)$

对于倾斜等晕角，简便易用且能昼夜连续测量的测量方法还未见报道。1995年A.Sivaramakrishnan等人利用2.4米望远镜观测星团内几颗恒星抖动的相关性来测量θ_ΤΑ(Sivaramakrishnan A，Weymann R J，Beletic J W.Measurements of the angular correlation ofstellar centroid motion[J].The Astronomical Journal，1995,110:430.)，1998年A.Ghedina等人利用0.2米的望远镜观测月亮边缘角间距几个角秒至几个角分的倾斜相关性来测量θ_ΤΑ(Ghedina A，Ragazzoni R，Baruffolo A.Isokinetic patch measurements on the edge of theMoon[J].Astronomy and Astrophysics Supplement Series，1998，130:561-566.)。这类方法都是直接测量倾斜相关，最后将倾斜相关曲线匹配经验公式来得出θ_ΤΑ，具有一定的误差。另外，如果用星团作为光源，需要超大口径望远镜，这对于常规观测是不现实的；如果用月亮边缘作为光源，观测时间将大受限制，因为高仰角的月亮并不是时常出现在观测天空的。

由倾斜等晕角的表达式(1)可知，可通过测量湍流二阶矩μ₂来间接测量θ_ΤΑ。众所周知，在孔径平滑效应下，大口径接收的恒星闪烁具有近似u²的路径权重，即可通过测量大口径的恒星闪烁来间接测量倾斜等晕角。然而，该方法实际运用比较困难，首先，需要大口径望远镜(口径越大，路径权重越趋近于u²)；其次，由于孔径平滑效应，大口径接收的恒星闪烁指数十分小，难以探测(测量易受噪声影响)。

发明内容

本发明的目的是提供一种简便易行且能昼夜连续监测的整层大气倾斜等晕角测量方法。

为了达到上述目的，本发明所采用的技术方案为：

1、一种整层大气倾斜等晕角的测量方法，利用安装在望远镜系统前端的两个圆形光瞳接收恒星光强，并计算相对光强起伏协方差C_R，通过选定合适的大气湍流层高度h₀、望远镜观测的天顶角φ确定相对光强起伏协方差C_R的湍流积分路径权重W_C(u)及其最大值，进而计算大气倾斜等晕角，其特征在于，包括以下步骤：

①、选定两个直径为D的接收光瞳安置于望远镜前端接收恒星的光强，利用望远镜后端的光电探测器测量两个光瞳分别接收到的恒星光强值I(x₁)和I(x₂)，并计算I(x₁)和I(x₂)的相对光强起伏协方差C_R：

$C_R=\frac{C_I\left(d\right)}{{<I>}^2}=\frac{<[I\left(x_1\right)-<I>][I\left(x_2\right)-<I>]>}{{<I>}^2}---\left(3\right)$

其中括号<>代表统计平均，x₁、x₂分别为两个接收光瞳的中心位置，d为两接收光瞳中心间距d＝|x₁-x₂|，<I>为恒星星光在各光瞳上的光强统计均值；因为光瞳大小一样，作为无限大平面波的星光在各光瞳上的光强统计均值应该是一样的，设为<I>。实时测量各光瞳上的瞬时光强值I(x₁)与I(x₂)，最后在一个统计时间内计算C_R。

②、将C_R转换为对数振幅起伏协方差C_χ：

$C_{\mathrm{\&chi;}}=\frac{\mathrm{In}(C_R+1)}{4}---\left(4\right)$

③、实验选定大气湍流层高度h₀、望远镜观测的天顶角φ，进而计算C_χ的湍流积分路径权重W_C(u)及其最大值A：

$A=\mathrm{Max}\left[W_C\left(u\right)\right],W_C\left(u\right)={\&Integral;}_0^{\&infin;}{\mathrm{\&kappa;}}^{-\frac{8}{3}}{\left[\frac{{2J}_1(\mathrm{\&kappa;}/2)}{\mathrm{\&kappa;}/2}\right]}^2{J}_0\left(\mathrm{\&kappa;}\frac{d}{D}\right){\sin }^2\left(\frac{{\mathrm{\&kappa;}}^{2}u}{4\mathrm{\&pi;}{F}_N}\right)\mathrm{d\&kappa;}---\left(5\right)$

其中，J_n为第一类n阶贝塞尔函数，F_N为菲涅尔数D²/λ₀L，λ₀为探测波长，L＝h₀secφ为所取湍流传播距离，其中h₀为所取湍流大气层高度，φ为天顶角；

④、将实验选定大气湍流层高度h₀、望远镜观测的天顶角φ，光瞳直径D以及以上计算所得的C_χ、A代入下式计算等晕角θ_ΤΑ

${\mathrm{\&theta;}}_{\mathrm{TA}}=\frac{1.40{\mathrm{\&lambda;}}_c{D_c}^{1/6}{D}^{5/6}{A}^{1/2}}{L{\mathrm{\&lambda;}}_0{C_{\mathrm{\&chi;}}}^{1/2}}=\frac{1.40{\mathrm{\&lambda;}}_c{D_c}^{1/6}{F}_N{A}^{1/2}}{D^{7/6}{C_{\mathrm{\&chi;}}}^{1/2}}---\left(\text{6}\right)$

其中，λ_c、D_c为与倾斜等晕角对应的波长、孔径直径，F_N为D²/λ₀ h₀secφ。

所述大气湍流层高度h₀、望远镜观测的天顶角φ的选定决定菲涅尔数F_N，进而决定大气湍流归一化积分路径权重W_C(u)/A，从而影响W_C(u)/A与倾斜等晕角θ_ΤΑ的路径权重u²的匹配程度，当W_C(u)/A与倾斜等晕角θ_ΤΑ的路径权重u²近似匹配时，可以将相应的C_χ、A代入公式计算等晕角θ_ΤΑ。

当所述的0.5<F_N<2.5时，一定存在合适的d/D使得W_C(u)/A近似匹配倾斜等晕角θ_ΤΑ的路径权重u²；当1≤F_N≤2时，W_C(u)/A近似匹配倾斜等晕角θ_ΤΑ的路径权重u²匹配误差的均方根RMSE较小，且对应的d/D均在0.53左右。

所述光电探测器为如CCD、CMOS面成像器件或者PMT、PD单元探测器件。

本发明的原理及依据是：两光瞳接收到的平面波相对光强起伏协方差与整层大气倾斜等晕角一样也是湍流路径的积分量，在两光瞳的直径D与间距d满足一定条件下，平面波相对光强起伏协方差的归一化路径权重W_C(u)/A能近似匹配倾斜等晕角的路径权重u²，从而可以通过测量相对光强起伏协方差来间接测量倾斜等晕角。相对光强起伏协方差C_χ(d)的归一化路径权重W_C(u)/A与大气层倾斜等晕角路径权重u²的匹配，即：

W_C(u)/A≈u²   (7)

其中，A＝Max[W_C(u)]。

对于平面波在弱起伏条件下传输，在Kolmogorov湍流谱模型与Rytov微扰近似假设下，孔径上对数振幅起伏协方差C_χ(d)可表示为

$C_{\mathrm{\&chi;}}\left(d\right)=1.303k^2{\mathrm{LD}}^{5/3}{\&Integral;}_0^1\mathrm{du}{C}_n^2\left(\mathrm{uL}\right)W_C\left(u\right)---\left(8\right)$

其中，WC(u)为C_χ的湍流积分路径权重：

$W_C\left(u\right)={\&Integral;}_0^{\&infin;}{\mathrm{\&kappa;}}^{-\frac{8}{3}}{\left[\frac{{2J}_1(\mathrm{\&kappa;}/2)}{\mathrm{\&kappa;}/2}\right]}^2{J}_0\left(\mathrm{\&kappa;}\frac{d}{D}\right){\sin }^2\left(\frac{{\mathrm{\&kappa;}}^{2}u}{4\mathrm{\&pi;}{F}_N}\right)\mathrm{d\&kappa;}---\left(9\right)$

由式(7)可知光强起伏协方差的归一化湍流积分路径权重W_C(u)/A的线型受F_N、d/D影响，对每一个不同的F_N，都对应存在一个d/D值使得W_C(u)/A与u²匹配的最好，在这里以匹配均方根误差RMSE来评价匹配效果，其中RMSE的计算方式为：

$\mathrm{RMSE}=\sqrt{\frac{\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{(W_C\left(u_i\right)-{u_i}^2)}^2}{N}}---\left(10\right)$

其中，N为数值计算中选取的路径等分段数，本次计算中取N＝1000。

本发明的优点为：

本发明所提出的一种整层大气倾斜等晕角的测量方法，利用双孔光瞳接收测量恒星闪烁协相关来间接测量倾斜等晕角，开辟了倾斜等晕角测量的新途径，具有如下优点：

(1)、只需一颗恒星作为光源，可不分昼夜连续监测。克服了直接测量倾斜抖动(像运动)相关需要星团或月亮边缘等非合作光源的问题。

(2)、克服了利用大孔径闪烁测量倾斜等晕角的困难，无需大口径接收望远镜，且易于探测。利用两个小口径上光强起伏的协相关可以替代超大口径闪烁的路径权重。

(3)、方法理论误差小。在F_N>1.1时，测量误差小于0.3％。与广泛使用的等晕角测量方法—恒星闪烁法相比，倾斜等晕角的路径权重匹配误差(RMSE)小约一个量级；与直接测量倾斜相关的方法相比，没有引入经验公式带来的误差。

本发明方法简单实用，测量精度高，在自适应光学、高分辨率成像、大视场成像等领域具有重要意义。

附图说明

图1为三个不同高空风速参数W(对应高空湍流强度)的HV模型下，各个高度湍流层对倾斜等晕角的贡献曲线图。

图2为不同菲涅尔数对应的最佳d/D，以及在这些组合(F_N，d/D)下路径权重匹配的均方根误差RMSE曲线图。

图3为光强起伏协方差测倾斜等晕角与恒星闪烁测等晕角两种路径权重匹配的RMSE曲线图。

图4为d/D＝0.53时不同菲涅尔数下的理论测量偏差Dev曲线图。

图5为菲涅尔数分布为0.78、1.11时对应的闪烁协方差的归一化路径权重函数曲线图。

图6为接收光瞳示意图。

具体实施方式

下面结合附图及具体实施方式详细介绍本发明。

一种整层大气倾斜等晕角的测量方法，利用安装在望远镜系统前端的两个圆形光瞳接收恒星光强，并计算相对光强起伏协方差C_R，通过选定合适的大气湍流层高度h₀、望远镜观测的天顶角φ确定相对光强起伏协方差C_R的湍流积分路径权重W_C(u)及其最大值，进而计算大气倾斜等晕角，包括以下步骤：

①、选定两个直径为D的接收光瞳安置于望远镜前端接收恒星的光强，利用望远镜后端的光电探测器测量两个光瞳分别接收到的恒星光强值I(x₁)和I(x₂)，并计算I(x₁)和I(x₂)的相对光强起伏协方差C_R：

$C_R=\frac{C_I\left(d\right)}{{<I>}^2}=\frac{<[I\left(x_1\right)-<I>][I\left(x_2\right)-<I>]>}{{<I>}^2}---\left(11\right)$

其中括号<>代表统计平均，x₁、x₂分别为两个接收光瞳的中心位置，d为两接收光瞳中心间距d＝|x₁-x₂|，<I>为恒星星光在各光瞳上的光强统计均值；因为光瞳大小一样，作为无限大平面波的星光在各光瞳上的光强统计均值应该是一样的，设为<I>。实时测量各光瞳上的瞬时光强值I(x₁)与I(x₂)，最后在一个统计时间内计算C_R。

②、将C_R转换为对数振幅起伏协方差C_χ：

$C_{\mathrm{\&chi;}}=\frac{\mathrm{In}(C_R+1)}{4}---\left(12\right)$

③、实验选定大气湍流层高度h₀、望远镜观测的天顶角φ，进而计算C_χ的湍流积分路径权重W_C(u)及其最大值A：

$A=\mathrm{Max}\left[W_C\left(u\right)\right],W_C\left(u\right)={\&Integral;}_0^{\&infin;}{\mathrm{\&kappa;}}^{-\frac{8}{3}}{\left[\frac{{2J}_1(\mathrm{\&kappa;}/2)}{\mathrm{\&kappa;}/2}\right]}^2{J}_0\left(\mathrm{\&kappa;}\frac{d}{D}\right){\sin }^2\left(\frac{{\mathrm{\&kappa;}}^{2}u}{4\mathrm{\&pi;}{F}_N}\right)\mathrm{d\&kappa;}---\left(13\right)$

其中，J_n为第一类n阶贝塞尔函数，F_N为菲涅尔数D²/λ₀L，λ₀为探测波长，L＝h₀secφ为所取湍流传播距离，其中h₀为所取湍流大气层高度，φ为天顶角；

④、将实验选定参数以及以上计算所得的C_χ、A代入下式计算等晕角θ_ΤΑ

${\mathrm{\&theta;}}_{\mathrm{TA}}=\frac{1.40{\mathrm{\&lambda;}}_c{D_c}^{1/6}{D}^{5/6}{A}^{1/2}}{L{\mathrm{\&lambda;}}_0{C_{\mathrm{\&chi;}}}^{1/2}}=\frac{1.40{\mathrm{\&lambda;}}_c{D_c}^{1/6}{F}_N{A}^{1/2}}{D^{7/6}{C_{\mathrm{\&chi;}}}^{1/2}}---\left(\text{14}\right)$

其中，λ_c、D_c为与倾斜等晕角对应的波长、孔径直径。

所述大气湍流层高度h₀、望远镜观测的天顶角φ的选定决定菲涅尔数F_N，进而决定大气湍流归一化积分路径权重W_C(u)/A，进而影响W_C(u)/A与倾斜等晕角θ_ΤΑ的路径权重u²的匹配程度，当W_C(u)/A与倾斜等晕角θ_ΤΑ的路径权重u²近似匹配时，可以将相应的C_χ、A代入公式计算等晕角θ_ΤΑ。

当所述的0.5<F_N<2.5时，一定存在合适的d/D使得W_C(u)/A近似匹配倾斜等晕角θ_ΤΑ的路径权重u²；当1≤F_N≤2时，W_C(u)/A近似匹配倾斜等晕角θ_ΤΑ的路径权重u²匹配误差的均方根RMSE较小，且对应的d/D均在0.53左右。

所述光电探测器为如CCD、CMOS面成像器件或者PMT、PD单元探测器件。

本发明选用口径356mm的Schmidt-Cassegrain望远镜(型号：Meada 14"f/10 LX200GPS)作为接收主体，具体步骤如下：

(1)、确定用于精确匹配的最大湍流层高度h₀。由于在高层大气密度随高度指数下降，高度20km以上的大气层十分稀薄，因而湍流强度非常弱。湍流模型SLC Day和SLC Night认为高于地面20.5km处光学湍流强度为零，H-V模型则通常认为高于地面30km处光学湍流强度为零。图1计算了不同高空湍流强度的H-V模型下各个高度湍流层对倾斜等晕角贡献的大小，也即对二阶湍流矩的贡献大小，具体计算方式为：

$R=\frac{C_n^2\left(h\right)h^2\mathrm{dh}}{{\&Integral;}_0^H{C}_n^2\left(h\right)h^2\mathrm{dh}}$    Equation Section 4(1)

H＝30km，dh取H/100，湍流分布采用H-V模型，固定近地面湍流强度，改变高空风速参数大小来改变高空湍流强度。从图1中可看出，20km以上的湍流层对倾斜等晕角贡献甚微，路径权重在这段路径的匹配精度对最后的θ_ΤΑ值影响不大。我们希望光瞳孔径D小点以减小接收望远镜的尺寸，考虑到菲涅尔数F_N＝D²cosφ/λ₀h₀，h₀越小，D就越小，为此我们选择h₀＝20km。即匹配0—20km路径段的路径权重。

(2)、选定接收光瞳大小D及间距d。为了为满足Rytov近似条件，通常要求天顶角φ小于45度；等效探测波长λ₀一般选在较长波段的650nm附近以减弱白天天空背景辐射。图2给出了不同菲涅尔数对应的最佳d/D，以及在这些组合(F_N，d/D)下路径权重匹配误差的均方根值RMSE。由图2可看出，当1≤F_N＝D²cosφ/λ₀h₀≤2时，匹配误差的均方根RMSE都比较小，且最佳比值d/D变化缓慢，在0.53左右。当φ＝45度时，D>13.6cm才能使得F_N>1，而对于口径356mm的Schmidt-Cassegrain望远镜，光瞳最大直径D约为12cm。因此D取最大值12cm，此时当φ从45度变化到0度时，0.78≤F_N≤1.11，属于可接收范围，经计算，此时测量理论误差不超过1％，测量偏差变化范围见图4。菲涅尔数分布为0.78、1.11时对应的闪烁协方差的归一化路径权重函数见图5，从图中5可看出，菲涅尔数为1.11的路径权重函数曲线与u²几乎完全重合，菲涅尔数为0.78的路径权重函数与u²也近似匹配。

实际观测中，两个光瞳间距与大小一般是固定的，即d/D固定，而观测恒星的天顶角φ是变化的(变化范围：0-45度)，即F_N＝D²/(λ₀h₀secφ)在一定范围内变化。

为此，我们数值计算了不同菲涅尔数对应的最佳d/D，以及在这些组合(F_N，d/D)下路径权重匹配的RMSE，如图2。RMSE最小值为0.117％，对应的最佳组合为(F_N＝1.60，d/D＝0.52)；当1≤F_N≤2，匹配误差的均方根RMSE都比较小，且最佳比值d/D变化缓慢，在0.53左右。

为了对RMSE的大小有个感性的认识，我们计算了利用平面波闪烁测量等晕角的闪烁权重匹配u^5/3的均方根误差RMSE，在最佳匹配时其RMSE约为1％，两种路径权重匹配的RMSE比较见图3。从图3可看出，协方差权重匹配的RMSE均小于用于测量等晕角的闪烁权重匹配的RMSE(小约一个量级)；利用闪烁测量等晕角一般要求RMSE小于1.5％(图中黑横线)，当F_N>0.5时，协方差权重匹配的RMSE总是小于1.3％(F_N＝0.5，RMSE＝1.3％)。因而可认为，当F_N>0.5时，可以利用恒星闪烁的协方差测量倾斜等晕角θ_ΤΑ，且F_N>0.7时测量理论误差均小于利用恒星闪烁测量等晕角在权重最佳匹配时的理论误差。

图4计算了d/D＝0.53时不同菲涅尔数下的理论测量偏差Dev，其计算方式为

$\mathrm{Dev}=\frac{{\mathrm{\&theta;}}_{\mathrm{TA}}-{\mathrm{\&theta;}}_{\mathrm{True}}}{{\mathrm{\&theta;}}_{\mathrm{True}}}---\left(2\right)$

θ_ΤΑ为协方差通过(4)式计算得出，θ_Τrue为通过定义式对不同湍流模型积分得出，计算中选取适用于不同地区的典型湍流模型：Middle East，GreenWood，Clear I，HAP，SLC Night，H-V5/7。计算得出，当F_N>1.1，测量偏差在±1％以内，且不同模型下的偏差值变化不大，在0.3％内；当F_N<1.1，测量值总是偏小，偏差是负值且随着菲涅尔数的减小而较快减小，不同模型下的偏差值较为发散，在3％内变化。若对偏差曲线减去各模型偏差平均得到的中值，可得测量误差Er为，当F_N>1.1，Er<0.3％，当F_N<1.1，随着菲涅尔数的减小，误差越来越大，在F_N＝0.5时将达到约3％。这个误差对光学湍流参数的测量来说，是十分微小的。

(3)、利用CCD探测I(x₁)与I(x₂)，计算两光瞳上的相对光强起伏协方差C_R。如图6，将两个光瞳化分为三个面积S₁、S₀、S₂，要想测量瞬时I(x₁)与I(x₂)，可以在S₁、S₂上放置楔镜，将三个面积的光分开聚焦到CCD靶面，形成三个光斑点，将S₁、S₀的光强相加得到I(x₁)，将S₂、S₀的光强相加得到I(x₂)。

(4)、将C_R转换为对数振幅起伏协方差C_χ：

$C_{\mathrm{\&chi;}}=\frac{\mathrm{In}(C_R+1)}{4}$    Equation Section 3(1)

(5)、数值计算C_χ的湍流积分路径权重W_C(u)的最大值A：

$A=\mathrm{Max}\left[W_C\left(u\right)\right],W_C\left(u\right)={\&Integral;}_0^{\&infin;}{\mathrm{\&kappa;}}^{-\frac{8}{3}}{\left[\frac{{2J}_1(\mathrm{\&kappa;}/2)}{\mathrm{\&kappa;}/2}\right]}^2{J}_0\left(\mathrm{\&kappa;}\frac{d}{D}\right){\sin }^2\left(\frac{{\mathrm{\&kappa;}}^{2}u}{4\mathrm{\&pi;}{F}_N}\right)\mathrm{d\&kappa;}---\left(2\right)$

其中，J_n为第一类n阶贝塞尔函数，L＝h₀secφ，为所取湍流传播距离，其中h₀为所取湍流大气层高度，φ为天顶角。

(6)代入相关参数计算倾斜等晕角θ_ΤΑ：

${\mathrm{\&theta;}}_{\mathrm{TA}}=\frac{1.40{\mathrm{\&lambda;}}_c{D_c}^{1/6}{D}^{5/6}{A}^{1/2}}{L{\mathrm{\&lambda;}}_0{C_{\mathrm{\&chi;}}}^{1/2}}---\left(3\right)$

其中，λ_c、D_c为与倾斜等晕角对应的波长、孔径直径；C_χ、Α从步骤4、5得出；L＝h₀secφ，h₀＝20km，φ取所用恒星的天顶角值(可从相关测量软件得出)；D＝12cm，λ₀根据所用恒星、滤光片、CCD探测器等的光谱特性来确定。

本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (4)

1.一种整层大气倾斜等晕角的测量方法，利用安装在望远镜系统前端的两个圆形光瞳接收恒星光强，并计算相对光强起伏协方差C_R，通过选定合适的大气湍流层高度h₀、望远镜观测的天顶角φ确定相对光强起伏协方差C_R的湍流积分路径权重W_C(u)及其最大值A，进而计算大气倾斜等晕角，其特征在于，包括以下步骤：

①、实验选定大气湍流层高度h₀、望远镜观测的天顶角φ，并选定两个直径为D的接收光瞳安置于望远镜前端接收恒星的光强，利用望远镜后端的光电探测器测量两个光瞳分别接收到的恒星光强值I(x₁)和I(x₂)，并计算I(x₁)和I(x₂)的相对光强起伏协方差C_R：

$C_R=\frac{C_I\left(d\right)}{{<I>}^2}=\frac{<[I\left(x_1\right)-<I>][I\left(x_2\right)-<I>]>}{{<I>}^2}---\left(1\right)$

其中括号<>代表统计平均，x₁、x₂分别为两个接收光瞳的中心位置，d为两接收光瞳中心间距d＝|x₁-x₂|，<I>为恒星星光在各光瞳上的光强统计均值；

②、将C_R转换为对数振幅起伏协方差C_χ：

$C_{\mathrm{\&chi;}}=\frac{\mathrm{In}(C_R+1)}{4}---\left(2\right)$

③、进而计算C_χ的湍流积分路径权重W_C(u)及其最大值A：

$W_C\left(u\right)={\&Integral;}_0^{\&infin;}{\mathrm{\&kappa;}}^{-\frac{8}{3}}{\left[\frac{{2J}_1(\mathrm{\&kappa;}/2)}{\mathrm{\&kappa;}/2}\right]}^2{J}_0\left(\mathrm{\&kappa;}\frac{d}{D}\right){\sin }^2\left(\frac{{\mathrm{\&kappa;}}^{2}u}{4\mathrm{\&pi;}{F}_N}\right)\mathrm{d\&kappa;},$

A＝Max[W_C(u)]   (3)

其中，J_n为第一类n阶贝塞尔函数，F_N为菲涅尔数D²/λ₀L，λ₀为探测波长，L＝h₀secφ为所取湍流传播距离，其中h₀为所取湍流大气层高度，φ为天顶角，

W_C(u)/A≈u²；

④、将实验选定大气湍流层高度h₀、望远镜观测的天顶角φ，光瞳直径D以及以上计算所得的C_χ、A代入下式计算等晕角θ_ΤΑ

${\mathrm{\&theta;}}_{\mathrm{TA}}=\frac{1.40{\mathrm{\&lambda;}}_c{D_c}^{1/6}{D}^{5/6}{A}^{1/2}}{L{\mathrm{\&lambda;}}_0{C_{\mathrm{\&chi;}}}^{1/2}}=\frac{1.40{\mathrm{\&lambda;}}_c{D_c}^{1/6}{F}_N{A}^{1/2}}{D^{7/6}{C_{\mathrm{\&chi;}}}^{1/2}}---\left(\text{4}\right)$

其中，λ_c、D_c为与倾斜等晕角对应的波长、孔径直径，F_N为D²/λ₀ h₀secφ。

2.根据权利要求1所述的一种整层大气倾斜等晕角的测量方法，其特征在于：所述大气湍流层高度h₀、望远镜观测的天顶角φ的选定决定菲涅尔数F_N，进而决定大气湍流归一化积分路径权重W_C(u)/A，从而影响W_C(u)/A与倾斜等晕角θ_ΤΑ的路径权重u²的匹配程度，当W_C(u)/A与倾斜等晕角θ_ΤΑ的路径权重u²近似匹配时，可以将相应的C_χ、A代入公式计算倾斜等晕角θ_ΤΑ。

3.根据权利要求1所述的一种整层大气倾斜等晕角的测量方法，其特征在于：当所述的0.5<F_N<2.5时，一定存在合适的d/D使得大气湍流归一化积分路径权重W_C(u)/A近似匹配倾斜等晕角θ_ΤΑ的路径权重u²；当1≤F_N≤2时，大气湍流归一化积分路径权重W_C(u)/A与倾斜等晕角θ_ΤΑ的路径权重u²匹配误差的均方根RMSE较小，且对应的d/D均在0.53左右。

4.根据权利要求1所述的一种整层大气倾斜等晕角的测量方法，其特征在于：所述光电探测器为面成像器件或者单元探测器件。