C# 实现道路横断面自动生成与格式转换（最小二乘拟合 + 方向向量法）

摘要：本文介绍一款基于 C# WinForms 的横断面数据处理工具，重点剖析最小二乘直线拟合剔除离群点、方向向量判断左右侧、以及 CASS 与纬地格式转换的核心算法。

1. 引言

道路横断面测量中，外业采集的散点常存在偏差，且不同软件（CASS、纬地）的横断面数据格式迥异。本文工具实现：

根据中桩方向自动识别断面左右侧。
最小二乘拟合剔除超限点。
一键导出 CASS 绝对偏距格式或纬地增量格式。

2. 核心算法详解

2.1 方向向量与左右侧判断

复制代码

// 计算前进方向向量（中桩序列）
double dirE = next.E - prev.E;
double dirN = next.N - prev.N;
// 归一化
double len = Math.Sqrt(dirE*dirE + dirN*dirN);
dirE /= len; dirN /= len;

// 计算散点相对于中桩的偏距及左右
double dx = p.E - mid.E;
double dy = p.N - mid.N;
double cross = dirE * dy - dirN * dx;   // 叉积
double offset = (cross > 0) ? Math.Sqrt(dx*dx+dy*dy) : -Math.Sqrt(dx*dx+dy*dy);

原理：二维向量叉积 A×B = Ax*By - Ay*Bx，正值表示 B 在 A 的顺时针方向（右侧）。
优势：无需显式计算方位角，避免了角度象限判断的繁琐。

2.2 最小二乘直线拟合与离群点剔除

断面点理论上共线，实际存在测量误差。采用最小二乘法拟合直线 y = a*x + b：

复制代码

double sumX=0, sumY=0, sumXY=0, sumX2=0;
foreach(var p in points) {
    sumX += p.E; sumY += p.N;
    sumXY += p.E*p.N; sumX2 += p.E*p.E;
}
double a = (n*sumXY - sumX*sumY) / (n*sumX2 - sumX*sumX);
double b = (sumY - a*sumX) / n;

// 剔除距离大于阈值的点
foreach(var p in points) {
    double dist = Math.Abs(a*p.E - p.N + b) / Math.Sqrt(a*a + 1);
    if(dist <= threshold) validPoints.Add(p);
}

效果：有效去除粗差点，保证横断面地面线的准确性。

2.3 纬地增量格式转换

纬地格式要求平距和高差均为相对于前一点的增量，且左侧平距也为正数。

复制代码

// 左侧点按距中桩由近到远排序
var leftPoints = offsets.Where(o=>o.Offset<0).OrderBy(o=>Math.Abs(o.Offset));
double prevOffset = 0, prevH = midH;
foreach(var p in leftPoints) {
    double deltaDist = Math.Abs(p.Offset) - prevOffset; // 平距增量（正）
    double deltaH = p.H - prevH;
    // 写入 deltaDist, deltaH
    prevOffset = Math.Abs(p.Offset);
    prevH = p.H;
}

输出示例：

复制代码

K1+000
3 3.500 0.800 1.700 -0.500 0.900 0.200
2 2.000 0.600 1.500 -0.300

3. 可视化交互

双缓冲绘图避免闪烁。
鼠标滚轮缩放、中键平移、左键框选放大。
青色箭头指示路线方向，绿色虚线为拟合直线，黄色虚线为阈值范围，灰色叉号为离群点。

4. 结语

本文工具将复杂的断面处理流程自动化，算法简洁高效，适用于道路、铁路等线性工程的横断面数据处理。完整代码可联系作者获取。

关键词：横断面；最小二乘；C#；纬地格式；CASS；方向向量