天文选址大气光学参数测量仪的数值模拟方法
【技术领域】
[0001]本发明涉及大气光学领域，具体涉及天文选址大气光学参数测量仪的数值模拟方法。采用该方法可以搭建虚拟实验平台，一方面可以仿真天文选址大气光学参数测量仪的实际工作过程，另一方面可以数值模拟多种天文选址大气光学湍流测量仪间的性能对比，为其实际对比试验提供参考。
【背景技术】
[0002]天文选址是国际天文界发展的重要领域之一，其理论研究和实测技术在近十几年中得到快速发展。天文选址目的:利用小型望远镜获得天文观测成果，实际检验台址条件，避免直接投入大型天文设备的风险，同时提升国际显示度和社会影响力。
[0003]天文光学湍流是天文选址中需要评估、测量的重要参数之一。目前主要有两种方法直接测量天文光学湍流。一种利用光波波前相位起伏，包括波前法向光程差(piston)、波前倾斜(tip-tilt)和波前离焦(defocus)等起伏。另一种利用光波波前闪耀。两种方法的观测目标包括恒星、月亮和太阳。
[0004]目前天文选址光学参数测量仪主要有:差分图像运动测量仪(DifferentialImage Mot1n Monitor, DIMM,基于星光波前倾斜起伏)和广义视宁度测量仪(GeneralizedSeeing Monitor, GSM,基于星光波前倾斜起伏)、干涉视宁度测量仪(InterferentialSeeing Monitor, ISM,基于星光波前法向光程差起伏)、快速离焦测量仪(Fast Defocus,FADE,基于星光波前离焦起伏)、多孔径闪耀传感器(Mult1-Aperture Scintillat1nSensor, MASS, 基于星光波前闪耀)和单星闪耀仪(Scintillat1n Detect1n andRanging, SCIDAR,基于星光波前闪耀)、外尺度廓线仪(Monitor of Outer Scale profile,M0SP，基于月亮边缘的光波波前倾斜起伏)和月亮边缘廓线仪(Profiler of Moon Limb,PML,基于月売边缘的光波波如倾斜起伏)等。
[0005]数值模拟技术是考察天文选址大气光学湍流测量仪性能的重要方法。国际天文选址工作组都用数值模拟技术仿真天文选址大气光学湍流测量仪。例如，托洛洛山美洲际天文台(Cerro Tololo Inter-American Observatory, CT10)天文选址工作组利用数值模拟技术研究DIMM、MASS和FADE的性能。南极冰穹C (Dome C)天文选址工作组利用数值模拟技术研究MOSP和PML的性能。但国际和国内天文选址工作组都没有提供统一数值模拟平台，用于上述的天文选址大气光学湍流测量仪性能的模拟。本发明提供虚拟实验平台，可实现数值模拟多种天文选址大气光学湍流测量仪间的性能对比，为实际对比试验提供参考。
[0006]具体地说，现有技术存在下述问题:
(1)现有技术没有提供统一数值模拟平台用于多种天文选址大气光学湍流测量仪的性能模拟；
(2)现有技术无法实现数值模拟更多天文选址大气光学湍流测量仪间的性能对比；
(3)现有技术模拟速度一直是瓶颈。
【发明内容】
[0007]针对现有技术中存在的上述问题，本发明的目的是提供一种天文选址大气光学参数测量仪的数值模拟方法。该方法能够解决:(1)提供统一数值模拟平台用于多种天文选址大气光学湍流测量仪的性能模拟；(2)实现数值模拟更多天文选址大气光学湍流测量仪间的性能对比；(3)提高模拟速度，解决模拟速度的瓶颈问题。
[0008]本发明中所提及的一种天文选址大气光学参数测量仪的数值模拟方法，该方法设有以下模块:配置参数模块、光波在大气光学湍流中传播模型、光学仪器物理模型、探测器模型、数据处理模块和运行方案模块；其特征在于如下步骤:
⑴.在配置参数模块中，配置模拟输入条件；
⑵.光学仪器跟踪各自目标，其目标包括:恒星、太阳或月亮；每个目标具有对应的唯一标识号；光学仪器除了自身属性外，附有目标和探测器属性，并具有对应的唯一标识号;
(3).配置大气光学湍流:配置大气光学湍流层高度分布，湍流权重高度分布，风速和风向高度分布，及相位屏；
⑷.在运行方案模块中，编写命令流控制文件；
(5).启动运行，并顺序执行如下三个模块:
(a).执行大气光学模块:依据光学仪器属性和大气光学湍流属性，采用相位屏技术和菲涅尔衍射方法，模拟目标光波在大气光学瑞流中传播，并到达相应的光学仪器入瞳处；
(b).执行光学仪器模块:依据光学仪器属性，采用夫琅和费衍射或菲涅尔衍射方法，模拟光学仪器功能；
(C).执行探测器模块:依据光学仪器属性，采用光波能量粒子化方法，模拟探测器功倉泛;
(6).执行数据处理模块:依据光学仪器属性，采用相应的数据处理方法，以友好的数据格式，输出模拟结果，同时也输出性能对比结果；
(7).对上述步骤调用的模块，设计加速算法；
(8).在计算机CPU上实现主要控制逻辑算法,在图形处理器(GraphicsProcessingUnit, GPU)上实现并行度高、数据吞吐量大的算法。
[0009]更优化和更具体地说，各步骤的操作方法是:
步骤⑴.在配置参数模块中，配置模拟输入条件:所述的多台光学仪器可以重复，但其位置或高度必需不同；
步骤⑵.软件为各光学仪器指定对应的唯一标识号。光学仪器跟踪各自目标，其目标可以相同。光学仪器配有探测器。即光学仪器除了自身属性外，附有目标和探测器属性，并具有对应的唯一标识号。另一方面，可以配置大气光学湍流层高度分布，湍流权重高度分布，风速和风向高度分布，其中湍流权重高度分布可以人工配置或由大气光学湍流分布模型自动生成。软件为大气光学瑞流指定对应的唯一标识号。大气光学瑞流配有相位屏模拟方法，即附有相位屏属性。
[0010]步骤⑶.在运行方案模块中，编写命令流控制文件。检查配置参数的一致性，检查整个流程，初始化观测目标、大气湍流和监测仪器。
[0011]步骤⑷.启动运行，并顺序执行如下三个模块: (a).执行大气光学模块。依据光学仪器属性和大气光学湍流属性，采用相位屏技术和菲涅尔衍射方法，模拟目标光波在大气光学湍流中传播，并到达相应的光学仪器入瞳处。
[0012](b).执行光学仪器模块。依据光学仪器属性，采用夫琅和费衍射或菲涅尔衍射方法，模拟光学仪器功能。
[0013](C).执行探测器模块。依据光学仪器属性，采用光波能量粒子化方法，模拟探测器功能。
[0014]步骤(5).执行数据处理模块。依据光学仪器属性，采用相应的数据处理方法，以友好的数据格式，输出模拟结果，同时也输出性能对比结果。
[0015]步骤(6).对上述步骤调用的模块，设计加速算法；
步骤(7).在计算机CPU上实现主要控制逻辑算法，在GPU上实现并行度高、数据吞吐量大的算法。
[0016]以上方案中，对所述的各模块的进一步说明是:
(I)配置参数模块:提供数字模拟方法的输入参数，包括目标天体信息、大气光学湍流信息、光学仪器信息、探测器信息、相位屏信息、以及其它必需的参数信息。
[0017](2)光波在大气中传播模型:提供相位屏模拟方法、相位屏到光波波前投影的模拟方法，光波在湍流层间传播的模拟方法。
[0018](3)光学仪器物理模型:提供天文选址大气光学参数测量仪的光学仪器物理模型，包括模拟薄透镜成像的夫琅和费衍射方法、模拟光场在自由空间传播的菲涅尔衍射方法、模拟瞳孔配置的瞳孔分割和瞳孔函数的方法、像平面插值配比CCD的方法。
[0019](4)探测器模型:提供能量粒子化方法、依据探测器类型添加噪声。
[0020](5)数据处理模块:提供天文选址大气光学参数测量仪的数据处理方法，以友好的数据格式，输出模拟结果，同时输出性能对比结果。
[0021](6)运行方案模块:编写命令流控制上述步骤的执行，检查配置参数的一致性，实现不同仪器在不同条件下监测的实时仿真。
[0022]本发明的方法是针对天文选址大气光学参数测量仪。这些测量仪包括DIMM、GSM、ISM、FADE、MASS、SCIDAR、MOSP和PML等，并且其光学参数也模拟输入条件之一。
[0023]本发明的优点:(I)具有扩展性，以上各模块均可增加新的内容；(2)可以搭建虚拟实验平台，实现数值模拟天文选址大气光学湍流测量仪间的性能对比；(3)设计加速算法，实现湍流监测仪器的实时仿真，主要控制逻辑算法在计算机CPU上实现，而并行度高、数据吞吐量大的算法在GPU上实现，例如快速变化的大范围大气湍流相位屏模拟等。
【附图说明】
[0024]图1是本发明的实施结构框图。
[0025]图2是本发明的实施结构框图中各模块的具体组成。
[0026]图3是本发明的软件流程图。
【具体实施方式】
[0027]下面结合附图及实施例对本发明作进一步说明。如图1所示，本发明主要由配置参数模块、运行方案模块、大气光学湍流遥感模块和数据处理模块的功能实现，其中所述的大气光学湍流遥感模块包括大气光学模块、光学仪器模块和探测器模块。上述功能模块的具体组成，如图2所示如下:配置参数