质点跟踪在气象矢量场分析中的应用

       这篇文章写于2009年，在电脑里放这么久，快发霉了，拿出来晒晒。希望对做矢量场可视化方面研究的同学有帮助。

       1 引言

       由于大气探测手段不断提高，人们获得气象数据的途径越来越多，气象数据也更加庞大、复杂，运用常规的分析方法很难发现数据的内在规律，对于既有大小又有方向的矢量数据，分析起来难度更高。科学计算可视化技术可将大量的数据转换为图像在计算机屏幕上显示出来，从而使冗繁而枯燥的气象数据变得形象直观，使研究人员方便、快捷地获得气象数据所包含的物理信息并找出其中隐藏的规律，为天气现象的准确预报以及气候研究提供帮助。

       本文首先介绍了质点跟踪技术的基本原理，随后从跟踪精度与效率两方面综合考虑，选择了合适的数值积分方法，并采用了根据矢量场方向变化情况动态调整跟踪步距的策略，进一步提高了跟踪轨迹的精度，最后采用该技术实现了三维不稳定风场数据的可视化。

       2 质点跟踪技术

       质点跟踪技术来源于流场的物理实验，在流场中将烟雾、染料等注入流场，观察其运动轨迹的分布，以了解流场的内部结构，该技术被广泛地用于流场实验。一般情况下，质点被认为是一点状发光的运动质点，质点的运动轨迹是一条曲线，由该质点在矢量场中的一系列连续位置组成。质点运动方程为：

                                                                               

其中P(t)是在t时刻质点的位置，V(P(t))是在矢量场中P(t)处的矢量值，即速度场。

       2.1 跟踪过程
       质点跟踪的过程就是确定质点在计算空间中一系列具体位置的过程，从初始位置P(x,y,z)开始，具体步骤为：
       （1）计算质点位置，获取质点所在的体单元；
       （2）在体单元内，通过插值求出该质点的矢量值；
       （3）解微分方程，计算出质点下一位置；
       （4）保存新的质点位置，并转到步骤（1）。
       跟踪过程中，如果质点运动到达计算空间的边界，应立刻停止跟踪，这是正常的终止条件。此外，由于在特定时刻任何空间点上的速度值唯一，两个质点不存在碰撞的情况，所以当两个质点相遇时，跟踪将停止。为了提高整个系统的实时性，每个质点都被赋予一定的跟踪时间步数，当跟踪迭代次数超过一定值时，将视质点超过生命周期，应停止该质点的跟踪过程。

       2.2 跟踪轨迹
       跟踪轨迹由质点在跟踪过程中获得的一系列离散位置组成，质点的位置可通过质点运动方程积分获得，对公式（1），其解为：

                                                           

       其中△t为时间片的一个逼近。在跟踪过程中，需要通过数值积分来求解微分方程。数值积分求解微分方程的方法很多，最常用的有欧拉法、改进的欧拉法以及四阶Runge-Kutta法等。由于欧拉法计算比较粗糙，质点轨迹误差较大，而Runge-Kutta法计算的质点轨迹精度最高，但计算开销却很大。因此本文采用改进的欧拉法，又称中点法，即计算速度快，同时还可以达到二阶精度。改进的欧拉法计算过程如下：首先用欧拉法求出预测位置 P'，

                                                                               
然后再利用P'(t + △t)得出t +△t时刻的速度值V(P'(t +△t), t + △t)，最后计算P(t +△t)的准确位置为：

                                                            
即利用t时刻P位置的速度值和在t + △t时刻位置P'的速度值求出平均速度，再计算准确的位置。

       2.3 跟踪步距
       在跟踪过程中△t的选择比较重要，△t取得太小，计算开销太大；取得太大，则会带来比较大的误差，降低质点跟踪的精度。同时在整个跟踪过程中如果时间步距固定为某个特定的值，那么在矢量场方向变化剧烈的情况下，质点有可能偏离正常的轨迹而造成跟踪的误差。因此本文采用一种自适应的步距调整方法，具体方法为首先计算某一质点跟踪路径上相邻两点之间的矢量方向夹角，然后通过夹角的大小来估计该区域矢量场方向的变化情况以相应调整跟踪步距。由于在跟踪过程中，相邻两点的连线与两点中第一个点的速度方向一致，因此在实际应用中，只需计算相邻两线间的夹角即可，如图1所示：

                                                                   

       对于所有质点，其初始步距可由下列公式确定：

                                                                                                 

其中，|V|为质点初始位置上的速度大小，Volume为质点初始时刻所在体单元的体积。通过该公式，可以确保质点初始跟踪的距离与其所在体单元的大小相匹配。在动态调整步距过程中，需不断计算轨迹线片段之间的夹角，一般地，当夹角大于15度时，表示速度场方向变化较大，则应将步距△t减半；当夹角小于3度时，表示速度场方向变化较小，则应将步距△t加倍。

       3 应用实例
       在气象领域中，风场是一个最为典型的矢量场。传统的气象可视化软件中，风场数据集的可视化通常只能以二维图像的形式进行显示。对于由纬向风、经向风和垂直风组成的三维风场，这些软件即不能动态地反映风场变化，也不能表达大气的垂直运动，因此本文将质点跟踪技术应用于三维风场数据的可视化，来帮助研究人员更好的掌握天气现象和规律。

                                                                               

       图2为基于质点跟踪的风场可视化流程图，在Visual C++ 6.0/OpenGL环境下编程实现了基于质点跟踪的三维不稳定风场可视化。图3为风场可视化效果图，实验数据格式为144x73x12。从跟踪结果看，用该方法可视化三维风场具有比较逼真的效果，并且绘制速度较快，适合在气象可视化软件中应用。

                               

       4 结束语
       气象领域中可视化技术的发展有着广阔前景,可视化技术的运用将成为未来气象预报、数据分析中的重要环节。本文详细介绍了质点跟踪技术，并应用该技术对三维不稳定风场数据进行了可视化，实验证明相比传统绘图软件，可视化结果更加形象、有立体感，并能准确、实时的反映三维不稳定风场的运动变化规律，具有一定的应用意义。