多波束测深数据处理：坐标转换、姿态改正与声线弯曲

简介《多波束水深测量数据处理》是一份面向海洋测绘、水深测量及多波束数据处理初学者与技术人员的专题文档围绕多波束面测深中的坐标系统、姿态改正与误差控制展开帮助读者理解从原始测深数据到成果归算的关键环节。压缩包仅含1个doc文件约152KB篇幅集中、便于检索适合作为坐标系定义、横摇纵摇改正、航向归算及系统偏差分析的学习底稿。文档结合测量船、多波束探头、姿态传感器、电罗经、水平坐标系与高斯测量坐标系等对象梳理坐标转换公式和声线处理思路并对传感器安装位置、系统偏差与测量误差的影响作出讨论。已有1125人学习可作为多波束水深测量数据处理入门梳理、公式查阅与误差排查的参考资料尤其适合需要补齐坐标转换和改正概念的读者。1. 一条测线跑完图上那块平地其实是礁石某次任务收工回看成果多波束全覆盖扫过的一片区域里单波束点和侧扫影像都报过的那个浅点在水深图上被抹得干干净净只剩一片顺滑的斜坡。设备没坏换能器没脏作业船速也压得住——问题出在从原始 Ping 到最终高斯坐标的那条归算链路上。多波束水深测量数据处理这件事难点从来不是会不会开机采集而是六个坐标系对不对得齐、姿态改正的符号有没有反、声线弯曲是不是还按常声速在算、系统偏差有没有真标过。这条链上一共就那么几个环节任何一个差一点点最后都糊在同一张图上看起来还特别像地形本来就长这样。下面按数据流的顺序把这条链路拆开走一遍。2. 六个坐标系先立住旋转矩阵错一个符号整片扇面就歪了多波束测深系统是一个多传感器综合系统测深仪、定位设备、测航向的电罗经、测纵摇横摇的姿态传感器、测上下起伏的涌浪滤波器缺一个都跑不起来。它们各自的安装位置、安装姿态都不一样最后要归到同一个水平面、同一个高斯坐标系里去。所以第一步不是写代码是把坐标系定义钉死。2.1 六套坐标的父子关系与各自职责工程上一般定义六套坐标系原点大多落在多波束探头发射中心轴向按右手系或按测绘习惯定义。这套约定一旦写进处理软件就不能中途改否则前后数据没法互相追溯。坐标系原点轴向约定主要用途测量船坐标系 XbYbZb探头发射中心向船艏看向左为 X 正沿艏艉线向船艏为 Y 正垂直向下为 Z 正描述所有传感器的安装位置探头坐标系 XcYcZc探头发射中心垂直扇面、指向船艏为 Y 正向左为 X 正中心波束方向为 Z 正波束向量起算姿态传感器坐标系 XdYdZd姿态传感器安装点与船坐标系对应轴同向平行横摇、纵摇观测电罗经坐标系 XeYeZe罗经基线与船坐标系对应轴同向平行航向观测水平坐标系 XaYaZa与平静时船坐标系重合XaYa 平面与当地水平面平行纵摇、横摇、航向倾斜改正的基准测量坐标系 XY高斯投影高斯平面直角坐标成果输出这张表的价值在于后面所有改正本质都是把某个坐标系下的向量转到另一个坐标系下。横摇改正和纵摇改正是船坐标系到水平坐标系航向归算是水平坐标系到高斯坐标系安装偏差改正是探头坐标系到水平坐标系。你只要在代码里把每一段转换写成独立函数出了问题就能逐段打日志看哪一段的旋转矩阵不是正交阵、哪一段的行列式不是 1。一个常见的错误是把 Z 轴方向搞反。测绘习惯里水深向下为正而姿态传感器和罗经往往按右手系、向上为正输出角度。这两个约定混在同一段代码里纵摇改正出来的水深会呈现越靠边越浅的假象而且量级不大很容易被当成声速剖面没选对。我一般会在处理流程最前面加一个断言静水状态下各传感器读数经转换后探头中心处的 Z 值应当等于静态吃水差超过 5 cm 就停下来查。2.2 横摇纵摇改正为什么工程上要用第二种写法姿态改正的数学表述很直接测量船坐标系 XbYbZb 绕水平面内某条直线 l 旋转一个角 ρ就能与水平坐标系 XaYaZa 重合。ρ 和等效旋转轴方位 κ 由横摇 α、纵摇 β 解出写成三次基本旋转的连乘形式就是一个 3×3 的旋转矩阵作用于波束向量。这套写法等效轴角推导漂亮但落到声线处理上不方便。原因是声线追踪要对每个波束逐个沿着声线积分中途每一步都要用当前入射角用等效轴角意味着每换一个波束就要重算一次轴和角中间过程还看不出是哪个姿态量在起作用。于是有了第二种写法先把水平坐标系绕 Za 轴旋转一个偏角 ε得到过渡坐标系 XaYaZa并且要求 Ya 轴正好落在 Yb 轴于 XaYa 平面上的投影方向再分别绕 Y即 Yb和 Xa 旋转两次旋转角 β、α 由 α、β 解出。三段旋转矩阵连乘形式上多了一个过渡坐标系好处是每一步都能单独验证第一段 Rz(ε) 只跟航向偏转有关静水直航时应当接近单位阵第二段 Ry(β) 只吃纵摇横摇为零时它单独起作用第三段 Rx(α) 只吃横摇纵摇为零时它单独起作用。偏转角 ε 有直接解算式不必反解矩阵。这一条在处理大批量 Ping 时很实用省掉一次矩阵求逆。实测下来两种写法在中低海况横摇 5° 以内下数值差异在毫米级取舍只看你的调试习惯要推导用第一种要排错用第二种。提示不管用哪种写法旋转顺序必须和设备内部定义一致。横摇在前还是纵摇在前结果差的是二阶小量但在大角度边缘波束上能放到厘米级验不出来就会一直怀疑声速。2.3 用 numpy 把姿态改正写成可单独测试的模块下面这段是姿态改正的最小闭环。核心思路是把三个基本旋转拆成三个函数组合顺序显式写出来中间不留隐式约定。import numpy as np def rx(a): 绕 X 轴旋转横摇方向 c, s np.cos(a), np.sin(a) return np.array([[1, 0, 0], [0, c, -s], [0, s, c]]) def ry(a): 绕 Y 轴旋转纵摇方向 c, s np.cos(a), np.sin(a) return np.array([[c, 0, s], [0, 1, 0], [-s, 0, c]]) def rz(a): 绕 Z 轴旋转航向方向 c, s np.cos(a), np.sin(a) return np.array([[c, -s, 0], [s, c, 0], [0, 0, 1]]) def boat_to_level(roll, pitch, orderroll_first): 把测量船坐标系向量转到水平坐标系。 roll : 横摇弧度右舷下沉为正须用静态试验确认设备符号 pitch : 纵摇弧度船艏抬起为正 order : 旋转顺序必须与设备内部定义一致 if order roll_first: return ry(-pitch) rx(-roll) # 先抵消横摇再抵消纵摇 return rx(-roll) ry(-pitch) def check_orthogonal(R, tol1e-9): 旋转矩阵自检正交且行列式为 1否则一定有符号或顺序写错了 assert np.allclose(R R.T, np.eye(3), atoltol), 矩阵不正交 assert abs(np.linalg.det(R) - 1.0) tol, 行列式不为 1可能是镜像 return True代码里的负号是抵消而不是施加姿态传感器给出的是船体相对水平面的倾角要把船体系下的向量搬到水平系旋转方向必须反过来。这一步符号搞反横摇改正会把倾斜放大一倍边缘波束的水深误差直接翻番。order参数看起来多余但它是最容易出问题的地方——同一批数据用两种顺序跑一遍边缘波束的差异如果超过 3 cm说明你的横摇角已经不小了得回头确认设备手册里的顺序定义。check_orthogonal这个断言别嫌麻烦。旋转矩阵一旦写错漏一个负号、把 sin 和 cos 写反、把转置当逆用矩阵仍然长得像个矩阵跑起来不报错只是地形慢慢变形。加上这个断言至少能保证错误不会以数据看起来还行的形式混进成果。验证办法很土但有效船靠泊静水记录一组横摇约 0°、纵摇约 0° 的姿态此时boat_to_level必须输出接近单位阵再人为把船压出一侧倾斜用压载或人员站位读出倾角喂进函数验证中心波束的方向余弦变化量和倾角的正弦值对得上。两关都过了姿态改正这一段基本可以封板。3. 航向归算与声线弯曲罗经读数和真实入射角之间隔着三层姿态改完波束向量从船体系转到了水平系接下来要解决两件事船头到底朝哪航向以及声线在水里到底拐了多少声线弯曲。这两件事都容易想当然也都在实测量级上吃过亏。3.1 光纤罗经和陀螺罗经测的不是同一个量两种罗经都靠地球自转工作但测出来的物理量完全不同处理方式也不同。光纤罗经基于光学环路中的相向传播光程差效应安装三个相互垂直的光纤陀螺分别测出地球自转角速度在这三个轴上的分量。地球自转角速度在当地可分解为水平分量和垂直分量两者都与当地纬度有关。船体水平静止时由三个分量就能直接解出航向。问题在于船在摇三个陀螺测到的是倾斜坐标系下的角速度分量必须先借姿态传感器测出的横摇、纵摇把这三个分量转到水平坐标系的两个水平方向上再解算真方位。顺序不能颠倒——先算航向再改正倾斜等于用错误的分量解了一个正确的方程。陀螺罗经走的是另一条路靠陀螺的定轴性和进动性让主轴自动找北。它的主罗经刻度盘随船体一起倾斜所以它测定的是倾斜面上的航向而不是水平面上的航向。如果设备本身没做倾斜修正就得自己按姿态角换算把陀螺主轴向量在倾斜坐标系下的分量经旋转矩阵变换到水平坐标系再解出真方位最后叠加上偏转角。这个改正的量级值得心里有数只考虑横摇、纵摇为零的情况下横摇 10° 时航向倾斜改正最大约 0.44°横摇 15° 时约 1°最大值出现在 45°、135°、225°、315° 四个航向上正横和正艏艉航向上为零。1° 的航向误差在 200 m 边缘波束上对应约 3.5 m 的点位偏移对 1:2000 比例尺的成图来说是不能接受的。所以罗经是进口的、精度 0.1°这种说法要打个问号——标称精度是水平状态下的倾斜状态下的残差得自己算。3.2 声线弯曲的二维分层积分与斯涅耳定律声线弯曲计算一般采用二维模式假定声速只在垂直深度方向变化水平方向均匀。这个假定在开阔水域成立得不错在河口锋面附近就会翻车那是后话。连续形式的表达是两条积分水深是声速乘入射角余弦沿时间的积分侧向中心距是声速乘入射角正弦沿时间的积分。实际处理一定离散化把声速剖面按深度分成若干层每层内声速取常数或线性变化逐层累加旅行时间、层厚和水平位移。关键在于层与层之间的入射角要按斯涅耳定律递推入射角正弦与声速之比在整条声线上是常数用换能器表面声速和初始入射角就能一路推下去。import numpy as np def ray_trace(svp, theta0, t_total, dt0.002, max_iter20000): 二维声线追踪声速只随深度变化。 svp : [(depth_m, sound_speed_mps), ...]按深度递增排列 theta0 : 换能器处的波束入射角弧度相对垂线 t_total : 单程旅行时间秒 dt : 时间步长秒2 ms 对应内插误差约 5 cm 量级 返回 : (D, X) 水深、侧向中心距单位米 D X t 0.0 theta theta0 v svp[0][1] p np.sin(theta) / v # 斯涅耳常数全程不变 for _ in range(max_iter): if t t_total: break # 当前深度落在哪一层取该层声速 depth D v _layer_velocity(svp, depth) sin_t min(1.0, p * v) # 防止数值越界 theta np.arcsin(sin_t) step v * dt D step * np.cos(theta) X step * np.sin(theta) t dt return D, X def _layer_velocity(svp, depth): 在声速剖面里做线性插值取当前深度的声速 depths np.array([s[0] for s in svp]) speeds np.array([s[1] for s in svp]) return float(np.interp(depth, depths, speeds))p sin(theta)/v是整段代码的骨架斯涅耳常数一旦算错后面全错。dt取 2 ms 不是随手写的按声速变化 15 m/s、时间步长 2 ms 估算步长中点处声线长度的变化小于 1.5 cm对垂直波束水深的影响在厘米以内。如果你的作业水域声速梯度特别大比如河口淡水楔这个步长要减半再验证一次收敛性。还有一个容易被忽略的细节theta0是波束入射角不是波束单元入射角。换能器各阵元的指向角是相对换能器垂直轴定义的船一摇这个角就不再等于波束中心线与水平面法线的夹角。必须先按姿态角把波束单元入射角换算成真正的波束入射角才能进声线追踪。换算过程涉及横摇、纵摇和波束方位角三者的组合边缘波束上差异能到好几度直接拿阵元角当入射角用边缘条带会整体上翘或下压。3.3 查表加线性内插把每 Ping 上千个波束算快一条测线每 Ping 少则几十、多则上千个波束逐波束跑积分在船载计算机上吃不消。常规做法是预先按垂线声速剖面对一组离散的入射角和旅行时间组合算好水深和中心偏距做成二维表实测时按每个波束的实际入射角和旅行时间在表里插值。插值一般用双线性先在时间方向插值再在角度方向插值或者反过来两种顺序在小步长下差异可以忽略。关键参数是两个步长入射角步长和旅行时间步长。误差估算有两个典型场景可以对照。边缘波束取斜距 200 m、入射角步长对应半角 1°、入射角 75°偏距误差约 0.03 m、水深误差约 0.01 m近中央波束取同样斜距和步长、入射角 1°偏距误差接近 0、水深误差约 0.03 m。再叠加声速在步长内变化带来的误差2 ms 内声速变 15 m/s声线长度变化小于 1.5 cm结论是入射角步长取 2°、时间步长取 2 ms 时线性内插对偏距和水深的影响不大于 5 cm。这个数字可以直接写进技术设计书的误差预算里。注意5 cm 是内插方法本身的误差不含声速剖面测错、姿态观测误差、安装偏差。别把它当成整体精度承诺。3.4 波束入射角算错边缘波束的水深就是假的波束入射角误差对成果的影响有解析式偏距误差正比于斜距乘入射角误差再乘余弦水深误差正比于斜距乘入射角误差再乘正弦。也就是说越靠边缘入射角越大同样的角度误差对水深的影响越大越靠中央对偏距的影响越大。斜距 200 m、入射角 60°、入射角误差 1° 时水深误差约 3 m。这个量级说明一件事入射角不是差不多就行的中间量它是决定水深的核心观测量之一。换算到水平坐标系时还要按波束方位把水深和偏距分解到 Xa、Ya 两个方向再叠加偏转角做一次平面旋转。这几步都是矩阵乘法写代码时建议每一步单独出一个中间数组并保留后处理发现条带异常时可以直接比对是哪一步的数值发散了。4. 系统偏差改正三条测线、两种水域把安装偏差逼出来坐标系定义得再漂亮传感器装上去总有偏差。罗经和姿态传感器的安装偏差还能靠静态测定校正多波束探头的安装偏差几乎没法静态测——你没法把探头拆下来放到水平台上量。所以只能靠用波束去打已知地形把偏差反演出来。4.1 横摇综合偏差为什么必须挑平坦水域往返测横摇综合偏差的测定要选平坦水域这个要求不是图省事。在平坦水域里航向偏差只影响水深点的平面位置不影响水深值纵摇偏差虽然影响水深但对同一 Ping 内各水深点是等幅度的改完之后水深还是平的。只有横摇偏差会让扇面左右两侧一高一低在平坦海底上表现为一个V形或倒 V形一眼就能看出来。测定方式是同测线往返测量。往返两趟的横摇偏差符号相反把两条测线的水深剖面叠在一起交叉点就是偏差为零的位置两侧的张开量直接对应偏差角。我一般会取重叠段内 200 个 Ping 做统计中误差超过 0.05° 就重测。顺序上有一条硬规矩横摇综合偏差必须先测先改。原因是它影响的是整个扇面的形状如果带着横摇偏差去测纵摇偏差纵摇的判读会被扇面倾斜污染反演出来的角度是错的。这条在实操里被违反的概率不低尤其是赶工期的时候。4.2 纵摇与艏向综合偏差斜坡和孤立目标怎么用纵摇综合偏差要选水深有明显变化的水域比如有一定坡度的斜坡航向垂直于斜坡走向。同样往返测量比较同一水深点在航向上的点位变化。为什么斜坡能凸显纵摇因为纵摇偏差会让波束在航向方向上提前或滞后打到底斜坡上同样的水深对应不同的平面位置往返一比较就出来了。平坦水域里这个效应被淹没了测不出来。这里有个顺序上的讲究纵摇偏差宜在艏向偏差之前测定。原因是艏向偏差会让探头扇面整体偏转一个角度如果没有纵摇和横摇影响同一船位处相邻两 Ping 是平行的测到的水深相同艏向偏差在平坦地形上几乎不可见。但一旦带着艏向偏差去测纵摇往返测线的同水深点根本对不上判读就废了。艏向综合偏差要选水深变化明显、最好有孤立目标点的水域用两条平行测线同向测量比较重叠带内同一目标点在航向上的位移。有沉船、暗礁这类目标点最理想因为平面位置可以精确到单个波束脚印。三种综合偏差的覆盖范围不一样这一点必须在报告里写清楚偏差类型包含的安装偏差测定水域要求测量方式横摇综合偏差探头 姿态传感器平坦水域同测线往返纵摇综合偏差探头 姿态传感器有明显坡度的斜坡同测线往返艏向综合偏差探头 电罗经水深变化明显、有孤立目标两平行测线同向还有一项容易被漏掉采用卫星定位的设备必须先做定位时延测定。船在动定位解算有延迟位置和水深点对不上表现出来就是整条测线沿航迹方向系统性偏移和艏向偏差的症状很像但改正手段完全不同。时延没测就去做艏向标定等于在一个错误的基础上再叠一层错误。4.3 涌浪改正的倾斜归算船头船尾的起伏不是一回事涌浪滤波器测的是它自己安装位置处的上下起伏但我们要的是探头发射中心处的起伏。船在浪里不是一个刚体平移而是绕重心转动船中后部起伏小船头起伏大。所以滤波器离探头越远直接拿观测值当探头起伏用的误差越大。正确处理是把滤波器的观测值先按安装偏差改正到探头坐标系再做倾斜转换到水平坐标系得到探头中心处的涌浪值。这一步需要滤波器的安装坐标相对探头原点和当时的横摇、纵摇。量级上设滤波器安装在探头后方 10 m 处横摇 10°、纵摇 5°姿态角本身有 1° 的系统偏差算下来涌浪改正值的差异能到 0.2 m。0.2 m 在浅水区是致命的在深水区会被当成噪声忽略掉然后某一天在浅滩上翻车。所以只要滤波器安装位置与探头的距离超过 5 m这一步就必须做不能因为平时都这么干就跳过。姿态传感器和罗经一体化的设备可以用综合偏差直接完成从探头坐标系到水平坐标系的转换姿态传感器和罗经充当中间过渡。分离式设备则要先静态测定姿态传感器与罗经坐标系之间的关系再对探头做动态标定。这里有个实际操作中的疑问点必须提前确认姿态传感器坐标系的 Y 轴到底指向哪、罗经坐标系是不是本身带倾斜。这两个问题不问清楚后面所有换算都建立在一个假设上。5. 五条踩坑记录从 1° 横摇到 1.66 米水深偏差下面这几条都是能在大批量数据里反复见到的症状按现象—原因—解决展开。前两条偏外业标定后三条偏内业归算。5.1 标定与安装阶段现象一平坦海底的成果图上出现一条贯穿全测区的V形折痕顶点在测线正下方。原因是横摇综合偏差没有改正或者改了但符号反了。平坦水域本来应该是一张平面横摇偏差让扇面左右两半各偏一个角度叠加起来就是 V 形。解决办法是回到平坦水域重做往返测用两条测线的交叉点确定零偏位置把改正量按左半扇加、右半扇减的方式重新施加。符号确认的土办法把改正量人为放大 10 倍看 V 形是变深还是变浅变浅说明符号反了。现象二往返两条测线的重叠带里水深明显分成两层差值几十厘米且沿整条测线稳定存在。大概率是纵摇综合偏差未改或者定位时延未测。两者的区别在于纵摇偏差造成的分层在往返方向上符号相反掉头回来后差值反向定位时延造成的偏移在同一航向上方向一致与航向无关。区分办法是把两条测线按航向分段统计偏移量随航向变化的是纵摇不随航向变化的是时延。解决办法分别对应重做斜坡往返标定和静态/动态时延测定。现象三涌浪改正加上去之后深水区数据反而抖得更厉害浅水区看不出问题。原因是涌浪滤波器安装位置离探头远观测值直接当成了探头处的起伏没有做倾斜归算。深水区波束斜距长同样的起伏误差在深度的投影上被放大浅水区斜距短误差被淹没在噪声里。解决办法是按滤波器安装坐标和姿态角做一次完整的归算把观测值搬到探头中心再施加改正。改完之后深水区的抖动应当明显收敛如果不收敛说明滤波器本身的观测在这一段就失真了风浪过大那要考虑限差作业而不是继续算。5.2 归算与后处理阶段现象四边缘波束整体上翘中央波束正常声速剖面换了几版都没改善。先别怀疑声速八成是波束入射角用错了。换能器各阵元的指向角是相对换能器垂直轴定义的船一摇这个角就不再是波束中心线与水平面法线的夹角。入射角取小了声线追踪算出来的水深偏深、偏距偏大扇面边缘就翘起来。解决是按姿态角和波束方位把阵元入射角换算成真实入射角再做声线追踪。验证办法取一段平直海底的数据边缘波束和中央波束的水深应当落在同一个统计分布里中误差不超过 0.1 m。现象五单波束和侧扫都报过的浅点多波束成果图上没有。这是最难受的一种。原因有两个方向一是系统探测能力的问题硬件层面就测不到二是测到了但处理不当。后者的典型场景是波束脚印内取了最快回波而不是最浅点——单波束习惯用最快回波代表水深多波束如果照搬这个逻辑在大入射角波束上会漏掉脚印内的高点浅点正好被抹平。另一个场景是粗差过滤把真实的浅点当噪声滤掉了尤其是孤立的单 Ping 浅点。解决办法是确认处理软件能测定波束脚印内的最浅点并反演其声线和入射角同时把粗差判别的阈值调到能保留单 Ping 异常、又能滤掉随机跳点的位置这一步只能靠人工复核加统计分布一起看。6. 验收怎么验一张误差预算表加一条往返自检线数据处理做完最怕的不是精度不够是不知道精度到底是多少。我的习惯是每条测线收工前跑一遍误差预算把各分项的贡献列清楚加起来和实测的往返较差对一对对不上就说明有一项被漏算了。误差分项典型量级主要控制手段横摇/纵摇观测误差0.05°0.1°设备标称精度 静态率定航向观测与倾斜改正残差0.1°0.4°光纤罗经倾斜改正、罗经基线校准声线弯曲与内插误差≤0.05 m入射角步长 2°、时间步长 2 ms声速剖面代表性误差0.10.3 m剖面站点密度、时段匹配综合偏差残差0.05°0.1°三条标定测线 重复观测涌浪与潮位改正残差0.050.2 m倾斜归算、潮位站同步定位时延残差0.10.5 m动态时延测定表里的数字是参考区间实际项目要用自己的标定结果替换。这张表真正的用处不是填完交差而是当实测往返较差超出预算时能快速定位是哪一项在冒头。再说一个我认为最值得坚持的验证动作用一段往返测线的重叠带做自检。不是看平均差值——平均差值会被系统性偏差吃掉看不出问题要看差值的分布形态。如果差值随入射角单调变化问题在声线弯曲或入射角换算如果差值随航向变化问题在纵摇或艏向偏差如果差值是常数偏移问题在时延或潮位如果差值在 Ping 之间跳变问题在姿态观测的时间同步或者粗差没滤干净。四种形态对应四类根因比逐项排查快得多。还有一个必须确认的功能处理软件要具备质量检查和水深恢复能力。所谓恢复是指原始观测值、改正量、中间结果都要能回溯。多波束水深测量是动态测量很多改正是在不确定的状态下施加的一旦发现某条测线的偏差改错了能不能把原始 Ping 重新拉出来、只重跑那一段改正决定了返工成本是一小时还是一星期。我见过不止一次因为中间结果没落盘一条测线从头重做到天亮。血泪经验是早些年做一条浅水测线横摇标定的符号反了图上 V 形特别规整规整到所有人都以为是真实地形直到和相邻测区拼接时对不上才被发现前后返工两天。从那以后我每次做完姿态改正都强制走一遍平坦水域自检 旋转矩阵正交性断言 往返重叠带分布形态这三关一关不过就不往下走。这套从坐标系到误差预算的完整拆解连同公式推导、参数表和可运行的追踪代码我整理成了一份可以直接照着走的资料需要的话可以拿去对着自己的处理流程逐段比对看哪一段的约定和你的软件不一致。希望帮到你。本文还有配套的精品资源点击获取