CN103606192A - 一种基于三维虚拟地球的风场可视化展示方法 - Google Patents

Abstract

本发明涉及一种基于三维虚拟地球的风场可视化展示方法，包括以下步骤：1.气象自动站数据接入：服务端以气象局提供的气象自动站数据作为风场数据的数据来源，以符合NetCDF规范的数据格式，定时导入气象自动站数据；2.风场数据格点化：获取的所述气象自动站数据为离散点数据，服务端将气象自动站离散点数据根据风场覆盖的可视化区域插值为风场格点数据；3.用户交互响应：客户端根据用户交互操作从服务端下载指定时间与空间的风场格点数据；4.风场可视化：客户端根据下载的风场格点数据驱动粒子系统，通过控制粒子在三维空间中的运动实现实时风场可视化。该方法不仅展示效果好，而且运行高效，实时性强。

Description

一种基于三维虚拟地球的风场可视化展示方法

技术领域

本发明涉及风场可视化展示技术领域，特别是一种使用格点化风场数据驱动粒子系统模拟风场、通过烟雾运动形式展示风场可视化、在三维虚拟地球环境下真实、快速模拟多层次（宏观、中观、局部）风场的方法。

背景技术

随着科学技术尤其是计算机技术的快速发展，传统可视化技术已经不能满足社会需求的快速提高，计算可视化成为了将科学计算数据转化为人们容易理解的图形或图像形式的主要手段，而且向着实时、交互、多维度、虚拟现实及网络应用等方向深入发展。

地理信息系统（GIS）作为地理与环境信息可视化的最好载体，已经越来越被社会认知与接受，尤其是近年来兴起的三维数字地球技术，集合传统GIS功能与三维可视化技术，在许多领域使用日渐频繁，尤其是在气象应用上，常用于大气环境监控、灾害预报预警等，可以将海量气象原始数据经过分析、转化、整理、统计等手段，得到人们需要的结果，但由于大气环境的多变性、复杂性，使用怎样的方法在三维数字地球中对其进行动态展示，一直是研究的热点方向。目前常用的风场可视化展示方法大体包含以下两类：

标量场可视化技术，如将二维网格分布的风场数据抽取出等值线，再对等值线进行取样插值平滑的等值线法。使用等值线法时通常搭配系列色板将各等值区域填充以便于使用者观看；

矢量场可视化技术，例如使用简单符号如箭头标注风向，箭头长短表示风速的标注法；使用一组带有表示气流方向箭头的曲线构成流线图的矢量线法；将风场中的风矢量数据赋予粒子系统，控制其中的粒子运动，从而动态展示风场的粒子法；基于纹理的矢量场可视化算法，例如线积分卷积中的点噪声法、纹理卷积算法。

但是，现有风场可视化技术普遍为二维可视化，基本方法如箭标图、风场分布图、流线图、纹理图等使用抽象元素，如箭头、颜色、线条等对风场的数据进行描述，不存在真实感也无法给予普通观察者直观感受，而且这些方法只能简单反映风场的瞬时状态，无法动态展示其发展与演变过程。粒子法可以生成较为细致的风场效果，可以进行动态展示，但由于其前期需要根据风场数据追踪粒子轨迹，后期需要对计算结果进行实时动画渲染，涉及到大量数学计算过程，同时随着粒子数量增加，计算量将迅速增长，将给系统软、硬件带来很大压力。而且真实风场变化频繁且涉及范围较大，使用粒子法对采集到的原始数据直接进行可视化处理，所需传输的数据也将占用大量网络带宽。因此传统的粒子法难以满足在三维数字地球系统中动态展示大面积实时风场的要求。

发明内容

本发明的目的在于提供一种基于三维虚拟地球的风场可视化展示方法，该方法不仅展示效果好，而且运行高效，实时性强。

为实现上述目的，本发明的技术方案是：一种基于三维虚拟地球的风场可视化展示方法，包括以下步骤：

步骤一、气象自动站数据接入：服务端以气象局提供的气象自动站数据作为风场数据的数据来源，以符合NetCDF规范的数据格式，定时导入气象自动站数据；

步骤二、风场数据格点化：获取的所述气象自动站数据为离散点数据，服务端将气象自动站离散点数据根据风场覆盖的可视化区域插值为风场格点数据；

步骤三、用户交互响应：客户端根据用户交互操作从服务端下载指定时间与空间的风场格点数据；

步骤四、风场可视化：客户端根据下载的风场格点数据驱动粒子系统，通过控制粒子在三维空间中的运动实现实时风场可视化。

进一步的，在步骤一中，服务端定时通过数据导入接口连接气象自动站数据收集服务器，导入气象自动站数据，导入时间间隔与气象自动站记录更新间隔相适应；数据导入后按气象自动站记录时间戳划分批次。

进一步的，在步骤二中，数据导入完成后触发风场格点数据生成，将获取的气象自动站离散点数据插值为覆盖所有可视化区域的风场格点数据；该方法与三维数字地球系统结合，风场格点在WGS-84坐标系下，以经纬度为单位进行划分，风场格点编号的确定方法为：将可视化区域以经纬度为单位，按四叉树规则分为若干等级，每个四叉树分块以n×n网格划分格点；设在任意分块内，经度最小值为Lon₀，最大值为Lon₁，纬度最小值为Lat₀，最大值为Lat₁，格点经度间隔为Lon_i，格点纬度间隔为Lat_i，经度方向格点数为LonNum，纬度方向格点数为LatNum，则Lon_i=(Lon₁－Lon₀)/LonNum，Lat_i=(Lat₁－Lat₀)/LatNum，对于在风场格点数据范围内的任意位置Pos_X(lon_X,lat_X)，其经向编号LonNo=[(lon_X-Lon₀)/Lon_i]，纬向编号LatNo=[(lat_X-Lat₀)/Lat_i]，式中符号[]表示取整；将所述风场格点数据存储于服务端，当客户端需要时通过时间与风场格点编号获取任意风场格点数据。

进一步的，在步骤三中，客户端根据数据时间与区域向服务端请求数据，数据时间由用户使用时间轴设置；数据区域由用户当前浏览区域确定，具体步骤为：首先由当前相机高度确定对应的四叉树数据加载等级，然后使用相机视景体与四叉树相交计算获取视野范围对应风场数据分块与格点行列范围，最后根据行列编号向服务端请求数据；客户端与服务端交互使用WCF数据服务，可根据网络环境需要选择连接方式。

进一步的，在步骤四中，风场可视化包括数据更新、视图更新和数据显示；

1）数据更新过程包括格点数据更新与粒子更新；数据更新过程中使用的时间包括：数据更新时间t，对应数据采集时间，用于控制数据加载；粒子更新时间t_i，即粒子系统更新粒子所使用的时间，用于控制粒子的运动；当t变化时，格点数据更新，当t_i变化时，粒子风场系统更新；

1.1）格点数据更新：客户端浏览不同时间点风场时，根据浏览时间加载前后两个时间点的风场格点数据；设前一个时间点为t₀，后一个时间点为t₁，系统当前时间为t，设当前浏览的格点在格网中的编号为(LonNo,LatNo)，该格点在两个时间点上的风速分别为Speed₀、Speed₁，则任意t满足t₀≤t≤t₁，t时刻速度Speed_T=(1-k)Speed₀+kSpeed₁，其中k=t/(t₁-t₀)；当系统没有新时间点数据时，速度Speed_T=Speed₀；

1.2）粒子更新：粒子更新过程包括粒子创建、粒子位置更新和边界检测；

粒子创建：根据粒子数量上限、限制边界、限制高度等预设参数与当前粒子分布等条件生成新的粒子，使粒子覆盖整个区域；粒子创建使用RenderToTexture技术，具体过程为：对应风场可视化区域定义一幅空白纹理；在随机渲染粒子的同时把粒子渲染到空白纹理的像素点上；在同一像素点重叠的粒子越多，像素灰度值越高；通过设置像素灰度阈值，仅在低于阈值像素对应空间区域生成粒子；

粒子位置更新：使用粒子系统展现风场可视化，粒子运动速度与方向由所在格点的风速与风向决定，粒子位置更新间隔由粒子更新时间t_i决定，当风场数据停止更新时，粒子将根据t_i继续更新；与风场格点相同，粒子系统采用WGS-84坐标系标注粒子坐标，但粒子位置变化使用世界坐标系计算：设粒子的初始WGS-84坐标为GeoPos₀(Lon₀,Lat₀,Alt₀)，地球半径为R_E，则粒子的初始世界坐标WorldPos₀(x₀,y₀,z₀)为：

x₀=(R_E+Alt₀)cosLat₀cosLon₀

y₀=(R_E+Alt₀)cosLat₀sinLon₀

z₀=(R_E+Alt₀)sinLat₀

设该粒子所在格点风向角为θ，风速为Speed_x，粒子更新时间为t_i，则新坐标WorldPos₁(x₁,y₁,z₁)为：

x₁=x₀+t_iSpeed_xsin(θ+180°)

y₁=y₀+t_iSpeed_xcos(θ+180°)

z₁=z₀

粒子的新WGS-84坐标GeoPos₁(Lon₁,Lat₁,Alt₁)为：

Lon₁=atan(y₁/x₁)

${\mathrm{Lat}}_1=a\sin (z_1/\sqrt{{x_1}^2+{y_1}^2+{z_1}^2})$

${\mathrm{Alt}}_1=\sqrt{{x_1}^2+{y_1}^2+{z_1}^2}-R_E$

粒子更新过程采用OpenCL算法实现；

边界检测：当粒子超过有效数据区域时，对粒子进行越界判断以移除无效粒子；边界检测使用RenderToTexture技术和3D绘制中的模版测试技术，具体过程为：首先定义一张空白纹理，将可视化区域绘制到模版缓冲区上；然后将粒子通过模版测试绘制到纹理上，使用粒子的ARGB值作为粒子的ID；最后移除所有绘制到纹理上的粒子；

2）视图更新过程，即用户在场景中观察位置变化的过程；风场可视化使用同一格点风场数据驱动的双层粒子实现，即覆盖视野范围的全局风场层与覆盖镜头的局部风场层；全局风场层根据视野范围更新数据以更新视图，保证相机位置变化时，能够看到连续的风场；局部风场层作用于相机周围小范围区域，以相机位置为中心，环绕相机并跟随相机同步移动，同时更新区域内视图，在镜头前构建更加逼真的效果；

3）数据显示过程，即粒子渲染过程；通过对同一纹理进行随机变换尺寸、透明度等处理，将其复用至所有粒子，粒子渲染使用GeometryInstance方法，将粒子顶点和纹理信息写入通道0，将粒子位置大小等信息写入通道1，使用不同的通道写入渲染设备。

相较于现有技术，本发明的有益效果是：

1、运行高效。客户端使用相机视景体与四叉树相交计算视野范围对应风场格点数据，减少了与服务端数据传输对网络的压力。使用格点化风场数据直接驱动粒子系统，不需要对单个粒子追踪轨迹，在粒子创建和边界检测中使用RenderToTexture渲染技术取代几何计算，从各方面降低了所需计算量，提高总体运行效率。

2、实时性强。服务端实时获取最新风场数据后将其格点化并存储，客户端可以直接获取到最新风场格点数据，客户端格点数据更新与粒子更新独立运行，可以保证动态可视化过程平滑切换至新风场，不出现视觉中断、卡顿或重置。

3、动态交互性强。由于服务端按时间分批存储了风场格点数据，客户端可自由查看任意时间的真实风场可视化效果，或直接按时间顺序播放风场变化情况。

4、真实性强。系统使用真实采集的风场数据进行可视化，最终可视化风场速度与方向与实际风场数据保持一致。使用了双层粒子系统，保证近景与远景的可视化效果，进一步提升了视觉真实性。

下面结合附图及具体实施例对本发明作进一步的详细说明。

附图说明

图1是本发明方法的工作流程图。

图2是本发明方法中风场可视化的模块组成图。

具体实施方式

本发明提供一种基于三维虚拟地球的风场可视化展示方法，如图1所示，包括以下步骤：

步骤一、气象自动站数据接入：服务端以气象局提供的气象自动站风速、风向数据作为风场数据的数据来源，以符合NetCDF规范的数据格式，定时导入气象自动站数据，以保证风场可视化的实时性；

步骤二、风场数据格点化：获取的所述气象自动站数据为离散点数据，服务端将气象自动站离散点数据根据风场覆盖的可视化区域插值为风场格点数据，插值方法采用距离反比法；同一时刻的气象自动站数据生成一组覆盖有效数据区域的格点数据；格点数据存储于服务端；

步骤三、用户交互响应：客户端根据用户交互操作从服务端下载指定时间与空间的风场格点数据；

步骤四、风场可视化：客户端根据下载的风场格点数据驱动粒子系统，通过控制粒子在三维空间中的运动实现实时风场可视化。

在步骤一中，服务端定时通过数据导入接口连接气象自动站数据收集服务器，导入气象自动站数据，导入时间间隔与气象自动站记录更新间隔相适应；数据导入后按气象自动站记录时间戳划分批次。

在步骤二中，数据导入完成后触发风场格点数据生成，将获取的气象自动站离散点数据插值为覆盖所有可视化区域的风场格点数据；该方法与三维数字地球系统结合，风场格点在WGS-84坐标系下，以经纬度为单位进行划分，风场格点编号的确定方法为：将可视化区域以经纬度为单位，按四叉树规则分为若干等级（具体分级数视可视化区域大小而定），每个四叉树分块以10×10网格划分格点；设在任意分块内，经度最小值为Lon₀，最大值为Lon₁，纬度最小值为Lat₀，最大值为Lat₁，格点经度间隔为Lon_i，格点纬度间隔为Lat_i，经度方向格点数为LonNum，纬度方向格点数为LatNum，则Lon_i=(Lon₁－Lon₀)/LonNum，Lat_i=(Lat₁－Lat₀)/LatNum，对于在风场格点数据范围内的任意位置Pos_X(lon_X,lat_X)，其经向编号LonNo=[(lon_X-Lon₀)/Lon_i]，纬向编号LatNo=[(lat_X-Lat₀)/Lat_i]，式中符号[]表示取整；将所述风场格点数据存储于服务端，当客户端需要时通过时间与风场格点编号获取任意风场格点数据。

在步骤三中，客户端根据数据时间与区域向服务端请求数据，数据时间由用户使用时间轴自由设置；数据区域由用户当前浏览区域确定，具体步骤为：首先由当前相机高度确定对应的四叉树数据加载等级，然后使用相机视景体与四叉树相交计算获取视野范围对应风场数据分块与格点行列范围，最后根据行列编号向服务端请求数据；客户端与服务端交互使用WCF数据服务，可根据网络环境需要灵活选择HTTP、TCP等多种连接方式。

如图2所示，在步骤四中，风场可视化包括数据更新、视图更新和数据显示；上述过程独立执行，无时间先后顺序；

1）数据更新过程包括格点数据更新与粒子更新；数据更新过程中使用的时间包括：数据更新时间t，对应数据采集时间，用于控制数据加载；粒子更新时间t_i，即粒子系统更新粒子所使用的时间，用于控制粒子的运动；当t变化时，格点数据更新，当t_i变化时，粒子风场系统更新；

1.1）格点数据更新：客户端浏览不同时间点风场时，根据浏览时间加载前后两个时间点的风场格点数据；设前一个时间点为t₀，后一个时间点为t₁，系统当前时间为t，设当前浏览的格点在格网中的编号为(LonNo,LatNo)，该格点在两个时间点上的风速分别为Speed₀、Speed₁，则任意t满足t₀≤t≤t₁，t时刻速度Speed_T=(1-k)Speed₀+kSpeed₁，其中k=t/(t₁-t₀)；当系统没有新时间点数据时，速度Speed_T=Speed₀；

1.2）粒子更新：粒子更新过程包括粒子创建、粒子位置更新和边界检测；

粒子创建：根据粒子数量上限、限制边界、限制高度等预设参数与当前粒子分布等条件生成新的粒子，使粒子覆盖整个区域；为了合理利用系统资源并提高显示效果，粒子不应在同一区域内过度重叠；由于粒子数量较大，使用常规算法在避免粒子重叠方面性能消耗较大，所以本发明中粒子创建使用RenderToTexture（渲染到纹理）技术以提高计算速度，具体过程为：对应风场可视化区域定义一幅空白纹理；在随机渲染粒子的同时把粒子渲染到空白纹理的像素点上；在同一像素点重叠的粒子越多，像素灰度值越高；通过设置像素灰度阈值，仅在低于阈值像素对应空间区域生成粒子；该方法主要优点在于，将RenderToTexture渲染技术作为计算方法使用，规避了几何计算过程，执行效率高，过程简单；

粒子位置更新：使用粒子系统展现风场可视化，粒子运动速度与方向由所在格点的风速与风向决定，粒子位置更新间隔由粒子更新时间t_i决定，当风场数据停止更新时，粒子将根据t_i继续更新，不会影响浏览效果；与风场格点相同，在三维数字地球系统中粒子系统采用WGS-84坐标系标注粒子坐标，但粒子位置变化使用世界坐标系（三维直角坐标系）计算：设粒子的初始WGS-84坐标为GeoPos₀(Lon₀,Lat₀,Alt₀)，地球半径为R_E，则粒子的初始世界坐标WorldPos₀(x₀,y₀,z₀)为：

x₀=(R_E+Alt₀)cosLat₀cosLon₀

y₀=(R_E+Alt₀)cosLat₀sinLon₀

z₀=(R_E+Alt₀)sinLat₀

设该粒子所在格点风向角为θ，风速为Speed_x，粒子更新时间为t_i，则新坐标WorldPos₁(x₁,y₁,z₁)为：

x₁=x₀+t_iSpeed_xsin(θ+180°)

y₁=y₀+t_iSpeed_xcos(θ+180°)

z₁=z₀

粒子的新WGS-84坐标GeoPos₁(Lon₁,Lat₁,Alt₁)为：

Lon₁=atan(y₁/x₁)

${\mathrm{Lat}}_1=a\sin (z_1/\sqrt{{x_1}^2+{y_1}^2+{z_1}^2})$

${\mathrm{Alt}}_1=\sqrt{{x_1}^2+{y_1}^2+{z_1}^2}-R_E$

由于需要对大量粒子进行实时更新，经测试，如果单纯使用CPU计算，在粒子数量超过10000时，将无法保证粒子更新的实时性，因此粒子更新过程采用OpenCL算法实现；该算法提供异构计算能力，可以利用GPU并发处理单元进行计算，使计算时间降低为原来的几十分之一，同时减小CPU占用；OpenCL技术支持当前各类主流集成、独立显卡，对硬件要求很低；

边界检测：当粒子超过有效数据区域时，对粒子进行越界判断以移除无效粒子；当使用区域边界较为复杂，例如行政区划分界等，使用常规算法将会导致大量的几何计算，所以边界检测使用RenderToTexture技术和3D绘制中的模版测试技术，具体过程为：首先定义一张空白纹理，将可视化区域绘制到模版缓冲区上；然后将粒子通过模版测试绘制到纹理上，使用粒子的ARGB值作为粒子的ID；最后移除所有绘制到纹理上的粒子；该方法与粒子创建方法相似，使用3D渲染技术规避复杂的几何计算过程；

2）视图更新过程，即用户在场景中观察位置变化的过程；本发明中，为了保证可视化效果，同时降低资源占用，风场可视化使用同一格点风场数据驱动的双层粒子实现，即覆盖视野范围的全局风场层与覆盖镜头的局部风场层；全局风场层根据视野范围更新数据以更新视图，保证相机位置变化时，能够看到连续的风场；局部风场层作用于相机周围小范围区域，以相机位置为中心，环绕相机并跟随相机同步移动，同时更新区域内视图，在镜头前构建更加逼真的效果；

3）数据显示过程，即粒子渲染过程；由于粒子数量较大，本发明中通过对同一纹理进行随机变换尺寸、透明度等处理，将其复用至所有粒子，最大化纹理利用率、减少视觉重复感，系统占用较低；粒子渲染使用GeometryInstance方法，将粒子顶点和纹理信息写入通道0，将粒子位置大小等信息写入通道1，使用不同的通道写入渲染设备；本方法适用于渲染大量近似目标，可极大的提高渲染效率。

本发明提出了将离散的数据进行插值并格点化，使用格点化的数据驱动粒子系统，不需要分别追踪单个粒子的轨迹。风场格点数据更新与粒子系统更新相互独立，采用分时控制异步模式，可以实现动态可视化过程平滑切换至新风场场景，不会出现视觉中断、卡顿或重置。

本发明还通过烟雾运动将不可见的风的运动形象化，增强了真实感；在三维虚拟地球环境中，创建双层粒子系统，保证了近景与远景的可视化效果，进一步提升了视觉真实性；用户交互时，系统根据用户的相机镜头及视野范围的变化，动态地在宏观、中观、局部风场效果之间连贯切换，带来连贯的视觉效果。

以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。

Claims (5)

1.一种基于三维虚拟地球的风场可视化展示方法，其特征在于，包括以下步骤：

步骤一、气象自动站数据接入：服务端以气象局提供的气象自动站数据作为风场数据的数据来源，以符合NetCDF规范的数据格式，定时导入气象自动站数据；

步骤二、风场数据格点化：获取的所述气象自动站数据为离散点数据，服务端将气象自动站离散点数据根据风场覆盖的可视化区域插值为风场格点数据；

步骤三、用户交互响应：客户端根据用户交互操作从服务端下载指定时间与空间的风场格点数据；

步骤四、风场可视化：客户端根据下载的风场格点数据驱动粒子系统，通过控制粒子在三维空间中的运动实现实时风场可视化。

2.根据权利要求1所述的一种基于三维虚拟地球的风场可视化展示方法，其特征在于，在步骤一中，服务端定时通过数据导入接口连接气象自动站数据收集服务器，导入气象自动站数据，导入时间间隔与气象自动站记录更新间隔相适应；数据导入后按气象自动站记录时间戳划分批次。

3.根据权利要求1所述的一种基于三维虚拟地球的风场可视化展示方法，其特征在于，在步骤二中，数据导入完成后触发风场格点数据生成，将获取的气象自动站离散点数据插值为覆盖所有可视化区域的风场格点数据；该方法与三维数字地球系统结合，风场格点在WGS-84坐标系下，以经纬度为单位进行划分，风场格点编号的确定方法为：将可视化区域以经纬度为单位，按四叉树规则分为若干等级，每个四叉树分块以n×n网格划分格点；设在任意分块内，经度最小值为Lon₀，最大值为Lon₁，纬度最小值为Lat₀，最大值为Lat₁，格点经度间隔为Lon_i，格点纬度间隔为Lat_i，经度方向格点数为LonNum，纬度方向格点数为LatNum，则Lon_i=(Lon₁－Lon₀)/LonNum，Lat_i=(Lat₁－Lat₀)/LatNum，对于在风场格点数据范围内的任意位置Pos_X(lon_X,lat_X)，其经向编号LonNo=[(lon_X-Lon₀)/Lon_i]，纬向编号LatNo=[(lat_X-Lat₀)/Lat_i]，式中符号[]表示取整；将所述风场格点数据存储于服务端，当客户端需要时通过时间与风场格点编号获取任意风场格点数据。

4.根据权利要求1所述的一种基于三维虚拟地球的风场可视化展示方法，其特征在于，在步骤三中，客户端根据数据时间与区域向服务端请求数据，数据时间由用户使用时间轴设置；数据区域由用户当前浏览区域确定，具体步骤为：首先由当前相机高度确定对应的四叉树数据加载等级，然后使用相机视景体与四叉树相交计算获取视野范围对应风场数据分块与格点行列范围，最后根据行列编号向服务端请求数据；客户端与服务端交互使用WCF数据服务，可根据网络环境需要选择连接方式。

5.根据权利要求1所述的一种基于三维虚拟地球的风场可视化展示方法，其特征在于，在步骤四中，风场可视化包括数据更新、视图更新和数据显示；

1）数据更新过程包括格点数据更新与粒子更新；数据更新过程中使用的时间包括：数据更新时间t，对应数据采集时间，用于控制数据加载；粒子更新时间t_i，即粒子系统更新粒子所使用的时间，用于控制粒子的运动；当t变化时，格点数据更新，当t_i变化时，粒子风场系统更新；

1.1）格点数据更新：客户端浏览不同时间点风场时，根据浏览时间加载前后两个时间点的风场格点数据；设前一个时间点为t₀，后一个时间点为t₁，系统当前时间为t，设当前浏览的格点在格网中的编号为(LonNo,LatNo)，该格点在两个时间点上的风速分别为Speed₀、Speed₁，则任意t满足t₀≤t≤t₁，t时刻速度Speed_T=(1-k)Speed₀+kSpeed₁，其中k=t/(t₁-t₀)；当系统没有新时间点数据时，速度Speed_T=Speed₀；

1.2）粒子更新：粒子更新过程包括粒子创建、粒子位置更新和边界检测；

粒子创建：根据粒子数量上限、限制边界、限制高度等预设参数与当前粒子分布等条件生成新的粒子，使粒子覆盖整个区域；粒子创建使用RenderToTexture技术，具体过程为：对应风场可视化区域定义一幅空白纹理；在随机渲染粒子的同时把粒子渲染到空白纹理的像素点上；在同一像素点重叠的粒子越多，像素灰度值越高；通过设置像素灰度阈值，仅在低于阈值像素对应空间区域生成粒子；

粒子位置更新：使用粒子系统展现风场可视化，粒子运动速度与方向由所在格点的风速与风向决定，粒子位置更新间隔由粒子更新时间t_i决定，当风场数据停止更新时，粒子将根据t_i继续更新；与风场格点相同，粒子系统采用WGS-84坐标系标注粒子坐标，但粒子位置变化使用世界坐标系计算：设粒子的初始WGS-84坐标为GeoPos₀(Lon₀,Lat₀,Alt₀)，地球半径为R_E，则粒子的初始世界坐标WorldPos₀(x₀,y₀,z₀)为：

x₀=(R_E+Alt₀)cosLat₀cosLon₀

y₀=(R_E+Alt₀)cosLat₀sinLon₀

z₀=(R_E+Alt₀)sinLat₀

设该粒子所在格点风向角为θ，风速为Speed_x，粒子更新时间为t_i，则新坐标WorldPos₁(x₁,y₁,z₁)为：

x₁=x₀+t_iSpeed_xsin(θ+180°)

y₁=y₀+t_iSpeed_xcos(θ+180°)

z₁=z₀

粒子的新WGS-84坐标GeoPos₁(Lon₁,Lat₁,Alt₁)为：

Lon₁=atan(y₁/x₁)

${\mathrm{Lat}}_1=a\sin (z_1/\sqrt{{x_1}^2+{y_1}^2+{z_1}^2})$

${\mathrm{Alt}}_1=\sqrt{{x_1}^2+{y_1}^2+{z_1}^2}-R_E$

粒子更新过程采用OpenCL算法实现；

边界检测：当粒子超过有效数据区域时，对粒子进行越界判断以移除无效粒子；边界检测使用RenderToTexture技术和3D绘制中的模版测试技术，具体过程为：首先定义一张空白纹理，将可视化区域绘制到模版缓冲区上；然后将粒子通过模版测试绘制到纹理上，使用粒子的ARGB值作为粒子的ID；最后移除所有绘制到纹理上的粒子；

2）视图更新过程，即用户在场景中观察位置变化的过程；风场可视化使用同一格点风场数据驱动的双层粒子实现，即覆盖视野范围的全局风场层与覆盖镜头的局部风场层；全局风场层根据视野范围更新数据以更新视图，保证相机位置变化时，能够看到连续的风场；局部风场层作用于相机周围小范围区域，以相机位置为中心，环绕相机并跟随相机同步移动，同时更新区域内视图，在镜头前构建更加逼真的效果；

3）数据显示过程，即粒子渲染过程；通过对同一纹理进行随机变换尺寸、透明度等处理，将其复用至所有粒子，粒子渲染使用GeometryInstance方法，将粒子顶点和纹理信息写入通道0，将粒子位置、大小等信息写入通道1，使用不同的通道写入渲染设备。

Images