IPv6单栈模式是指仅使用IPv6不支持IPv4。向 IPv6单栈过渡是必然趋势。而IPv6单栈演进的可视化测量研究目前较少。设计IPv6网络可视化测量系统，既要考虑IPv6协议本身的标识、封装属性，又要考虑可视化测量系统自身在IPv6单栈环境下的运行。

       随着IPv6规模部署的推进，各类应用系统及设备对IPv6网络和协议支持需求的剧增成为IPv6单栈演进需要考虑的关键问题之一。支持IPv6的高精度卫星定位网络测量系统，作为信息物理系统^[1](Cyber-Physical Systems，CPS)的一部分，对各行各业开展更大范围、更深层次的IPv6协同创新，推动IPv6单栈大规模部署具有非常重要的意义。

      在此背景下，作为共性关键技术的IPv6单栈的演进和北斗系统^［2］的融合应用变得尤其重要。因此，IPv6可视化测量系统需要兼顾支持IPv6单栈和高精度卫星定位这两方面的设计问题，系统自身也能在IPv6单栈环境下稳定运行。

       北斗三号导航定位系统作为全球导航卫星系统^［3］（GlobalNavigation Satellite System，GNSS），发射卫星数量最多，覆盖120个国家和地区，轨道径向精度6.6厘米^［4］。北斗三号在高精度算法和高精度板卡制造方面，运用RTK(Real-time Kinematic，实时动态差分 ) 技术能将精度提升至厘米级。

       消息队列遥测传输^［5］(Message Queuing Telemetry Transport，MQTT)是ISO标准下基于发布/订阅范式的消息协议。它工作在TCP/IP协议族上，是为硬件性能低下的远程设备以及网络状况糟糕的情况而设计的发布 / 订阅型消息协议。
 

       MQTT 协议在受限的环境中，如卫星链路通信传感器、机器与机器通信和物联网中已广泛使用。通过MQTT协议，目前已经扩展出了数十个MQTT服务器端程序，可以通过PHP、Java、Python、C、C#等编程语言来向MQTT发送相关消息。

       本文针对移动终端运行关键数据，与全球卫星定位系统的经纬度坐标绑定，基于MQTT协议进行IPv6网络发布和订阅，为指挥控制中心全面了解和监测现场移动终端和被测对象的运行状态，提供可视化测量系统并进行验证和展示。

       全球卫星定位系统的基本原理是：测量出已知位置的卫星到用户接收机之间的距离，然后综合多颗卫星的数据，就可知道接收机的具体位置。
 

      要达到这一目的，卫星的位置可以根据星载时钟所记录的时间在卫星星历中查出。而用户到卫星的距离则通过记录卫星信号传播到用户所经历的时间，再将其乘以光速得到。由于大气层电离层的干扰，这一距离并不是用户与卫星之间的真实距离，而是伪距（PR）。
 

       当GPS卫星正常工作时，会不断地用1和0二进制码元组成的伪随机码（简称“伪码”）发射导航电文。然而，由于用户接收机使用的时钟与卫星星载时钟不可能总是同步，所以除了用户的三维坐标 x、y、z 外，还要引进一个Δt，即卫星与接收机之间的时间差作为未知数，然后用4个方程将这4个未知数解出来。因此，如果想知道接收机所处的位置，至少要能接收到4颗卫星的信号^［6］。
 

图1 系统组成结构

      接收机与RTCM服务器为同一台服务器，移动终端设备由UM482基站、USB转RS485模块、支持本地链路地址FE80::/10的通信模块组成，移动终端(传感器 + 嵌入式 OS)和MQTT服务器之间采用MQTT协议通信，订阅RTCMessage信息并发布DeviceMessage信息。

       IPv6协议单栈部署设计包括IPv6单栈网络设计和 IPv6单栈站点设计。

       IPv6单栈网络设计的重点为基于DHCPv6的地址生成系统，此系统支持基于用户信息和地理信息生成IPv6地址，系统组成原理如图2所示。

图2 单栈IPv6地址生成系统原理

      具体而言，DHCPv6服务器确定所述地理信息与DHCPv6服务器的地理信息的差值，并分别对该差值以及与所述用户信息相关的账号信息进行压缩处理。

      由此确定与所述差值相关的5个字节，以及与账号信息相关的3个字节。其中，5个字节包括与经度信息相关的2个字节，与纬度信息相关的2个字节以及与高度信息相关的1个字节。

      压缩处理后，确定与该DHCPv6服务器单元通信的用户端的IPv6地址的后64位，并判断该用户端的IPv6地址是否与数据库中存储的已有IPv6地址相同。

       若相同，则DHCPv6服务器单元随机生成一修正IPv6地址，并将其下发至该用户端；若不同，则保留该用户端的IPv6地址，并将其下发至该用户端。
 

       IPv6单栈站点设计的重点在于禁止IPv4协议的运行或在网卡上禁用IPv4地址，本文基于ubuntu20.04和Apache2.4.48实现。

       根据可视化测量系统的设计流程，通过MQTT服务器传输的消息结构体主要包括RTCMessage 和DeviceMessage两种，其中RTCM主要是增加长度字段转发。本文设计了移动终端设备到MQTT服务器的DeviceMessage数据结构，该结构体为56字节，如下所示，保留字段支持不同类型的时延地理敏感性指标参数填充。
 

RTCMessage消 息 发 布 订 阅 程 序 流 程 设 计 如 下：

图3 RTCM发布流程和MQTT订阅流程

DeviceMessage消息发布订阅程序流程设计如下：

       通过对单点定位和多点定位精度测量的对比可以发现，多点定位时经纬度位置偏移范围更小。接收的卫星定位经纬度信息可通过对保存的历史数据进行偏移分析来获取。在多点定位固定基站和移动站环境搭建完成后，可以部署 MQTT程序传输技术参数数据并进行应用展示。MQTT程序需安装开发库以支持RTCMessage消息和DeviceMessage的发布和订阅。
 

      DeviceMessage消息订阅程序根据移动终端发布的内容进行订阅，消息内容包括设备号、参数大小、定位模式、经纬度、设备状态等信息，其中经纬度信息是卫星多系统联合定位的结果，联合定位输出信息包含GNSS接收机的时间和定位相关数据。语句以GNGGA^［7］开头，表示GNSS定位信息，只有北斗参与时以BDGGA开头，有两个卫星系统及以上的卫星参与定位解算都以GNGGA形式输出。表1为GNGGA数据结构表。

表1 GNGGA数据结构

其中，GPS质量指示符qual字段表示定位解算的状态，主要有：

      GNGGA数据结构经度（DDDmm.mm）和纬度（DDmm.mm），其中，D是度，m是分，小数点前有4~5位，与传统的表示法不同。移动终端设备在发布之前转换成传统表示法，按小数点前两位截取，把m部分都除以60，放小数点后面。纬度范围是-90~+90，可能出现D在赤道上纬度接近于0的情况，仍然按小数点前两位，没有的位数补零，经纬度方向有正负，北纬为正，东经为正。

       Python编程语言通过格式化字符串（format strings）规定转化的方法和格式。根据 Python 类型与格式和大小的对应关系，上述Device消息56字节的结构体打包和解包的代码为：struct.pack('>IfcddBBBHHLLQQ',message.payload)struct.unpack('<ifcddbbbh HLLQQ',message.payload)，其中‘>’为字节流大端序，用于把数据发送到网络，‘<’为字节流小端序，用于从网络解包字节流。

       其中连接回调函数on_connect在客户端连接后被调用，依据返回值rc判断客户端是否连接成功。同时创建一个MQTT客户端，该客户端将连接到消息回调函数on_message，该函数将在客户端收到消息后被调用，在该函数中打印出订阅的topic名称以及接收到的消息内容。Device消息订阅程序文件命名为：Mqtt-sub-device-message.py。运行程序，整个系统正常工作，经纬度数据小数点后长度不同表示精度不同，精度不同分配的IPv6地址前缀不同。

       从显示的信息顺序可以看出，定位状态首先是0（不可用或无效），过渡到RTCM消息还没参与解算是1（单点定位），参数设备自身收到伪距差分信息后状态为2（伪距差分），然后根据RTCM消息和伪距差分消息结算为5（RTK 浮动解），一分钟后得到4（RTK 固定解），并保持稳定。固定解小数点后5位保持不变，表示定位精度已达到厘米级。

       参数的变化过程为：首先是未确定固定解时8NM 上下浮动的误差值，然后是7.82NM的设定目标值，中间校准过程分别为6.07NM、5.94NM和6.05NM，最后达到了7.82NM的设置目标值，参数校准过程完成。可预先设置阈值进行报警。

       定位状态收敛过程为：首先是状态1（单点定位）下的低精度范围，然后是状态5（RTK 浮动解）定位下的重精度范围，最后为状态4（RTK 固定解）定位下的高精度范围。

       Web服务器基于IPv6单栈技术建设，支持IPv6单 栈访问。系统通过MQTT协议读取不同精度的数据存入Web服务器，调用百度地图接口显示经纬度坐标，经过坐标转换生成html页面。通过IPv6地址访问http://[2402:e740::666]/site_ipv6/line.html站点，显示图4所示的效果，发现在100米长度的线路上，移动设备技术参数大小与位置变化的过程。

图4 IPv6移动终端订阅技术参数信息位置显示

      访问http://[2402:e740::666]/site_ipv6/bpos.html站点，显示图5所示的效果，因百度地图的精度为米或亚米级别，基本看不出精度是厘米级的点的区别。

图5 Web服务展示低精度参数信息位置显示

      通过高精度交互软件转换为离散轨迹图，显示出图6的效果，可以看出精度范围在6厘米内，满足测量场景中技术参数厘米级的要求。

图6 高精度交互软件订阅高精度卫星数据显示

       本文从IPv6创新应用向IPv6单栈演进的角度出发，设计并实现了一种支持IPv6的可视化测量系统。在系统中对经纬度信息、高度信息和用户身份信息如何映射到IPv6后64位（可变）进行了详细设计，提供了IPv6站点内的位置查询相关服务。尝试解决可视化测量系统在推进IPv6规模部署时面临的扩展性、实时性、移动性、管理性四个方面的重大挑战，提供一套综合解决方案。同时，在面向工业互联网环境下，尝试从信息物理系统的角度出发，为基于IPv6网络测量物理世界技术参数的安全性方面提供了参考，在IPv6创新应用领域中起到一定的应用示范作用。

参考文献（上下滑动查看全部内容）

本文来源：中国教育网络，作者：黄友俊（清华大学网络科学与网络空间研究院、赛尔网络有限公司），李想、徐绘凯、刘保君、张甲（清华大学网络科学与网络空间研究院）