[发明专利]一种多功能天文经纬仪子午圈测纬的星径曲率改正方法

权利要求书

1.一种多功能天文经纬仪子午圈测纬的星径曲率改正方法，其特征在于，包括以下步骤：

步骤1、多功能天文经纬仪采用转轴观测、短曝光采集多幅恒星像；

步骤2、根据多功能天文经纬仪的特点，计算理想情况下的子午圈测纬的星径曲率改正；

步骤3、考虑仪器方位差影响，计算非理想状态下测纬的星径曲率改正公式，提高天文纬度测定值的精度；

所述步骤2中的根据多功能天文经纬仪的特点，计算理想情况下的子午圈测纬的星径曲率改正，具体为：

步骤2.1、转轴前后所有露光图像上的被测星在天球上的对应位置，全部对称于多功能天文经纬仪的子午圈排列，以观测者为天球的中心o，x轴指向西点，z轴指向观测点的天顶，x、y、z三轴满足右手坐标系法则；σσ为周日平行圈，H为子午面，B、F和C、E分别为对称于子午面的平行平面；周日平行圈σσ与平面E、F、H、B、C的交点分别为：σ_n、σ₁、σ₀、σ_-1、σ_-n；σ₀为星过子午圈时的位置，多功能天文经纬仪在转轴前后各进行n次观测，则会有n对平行于子午面的平行平面，它们与周日平行圈亦有对称于σ₀的n对交点；

星径曲率的基本改正公式为：

其中，k为星径曲率改正值，δ为恒星视赤纬；本多功能天文经纬仪采用CCD相机连续短曝光的方式采集多幅恒星图像，第i副图像的露光时刻为tⁱ，其星像相对于子午面的角距为：

σ₀σ_i=(tⁱ-T₀)*15″*cosδ， (2)

其中，T₀为恒星过子午圈的时刻，角距的单位为角秒；

利用转换关系cos²δ*tanδ＝sin2δ/2，得第i副图像上的星径曲率改正值：

kⁱ＝[(tⁱ-T₀)*15″]²*sin1″*sin2δ/4 (3)

假设转轴时间为T，转轴前后各采集n幅图像，每副图像露光时间为Δt，取转轴的平均时刻T₀为时间轴起算点，即星过子午圈的时刻，时间前进的方向为正向；

转轴前的平均露光时刻为则转轴后的平均露光时刻为按露光顺序排列，转轴前1～n幅恒星像的露光时刻t¹～tⁿ分别为：

则转轴前第i副图像的露光时刻为：

转轴后n+1～2n副恒星像的露光时刻tⁿ⁺¹～t²ⁿ分别为：

则转轴后第j副图像的露光时刻为

由于转轴观测能够消除零点误差，故而不单独计算转轴前第i副图像的星径曲率改正值，而是取转轴前第i副图像星像的时角和与其相对称的转轴后第(2n+1-i)副图像的恒星像的时角，计算出它们对应的星径曲率改正项的平均值：

为压缩随机误差的影响，该被测星总的星径曲率改正项应取转轴前后所有图像所对应的星径曲率改正项的平均值：

根据本仪器的观测流程，转轴时间T＝22s，转轴前后各采集的图像数n＝100副，每副图像露光时间为Δt＝45.44ms，将以上参数代入公式(7)中，得D＝79768，则星径曲率改正为：

公式(8)就是在理想情况下，即测时的各种误差综合影响为零的情况下，测纬的星径曲率改正公式；

所述步骤3中的考虑仪器方位差影响，计算非理想状态下测纬的星径曲率改正公式，具体为：

对于非理想情况下的星径曲率改正，假设仪器方位差为a，时秒，由于被测星的方位系数：z为天顶距，为测站的纬度采用值；则对于每颗被测星，转轴前的平均露光时刻应为：

转轴后的平均露光时刻应为：

将仪器参数T＝22s，n＝100副，Δt＝45.44ms代入公式(4)、(5)、(7)、(9)、(10)中，得：

D＝(198.7392+15″*a*A)²+(-198.7392+15″*a*A)²+774(11)

由公式(11)计算出中间参数D，代入公式(7)中，便可得到方位差影响下的子午圈测纬星径曲率改正项；在对测纬进行平差时，加入这一改正项。

2.根据权利要求1所述的星径曲率改正方法，其特征在于，所述步骤1中的多功能天文经纬仪采用转轴观测、短曝光采集多幅恒星像，具体为：当一颗被测星星像通过视场中垂线之前开始连续露光，用主光路终端CCD相机先得到50幅人造星像，再得到100幅恒星像，随后将仪器上盘及其以上部分，绕方位轴旋转180°，同时将镜筒随着高度轴旋转至预置天顶距2倍的角度，再次指向被测星开始连续露光，也得到100幅恒星像和50副人造星像。