简介基于VB.NET实现非均匀B样条曲线NURBS的完整示例项目面向计算机图形学学习者与VB.NET开发者演示在.NET环境中构建2D/3D建模所需的曲线控制与求值功能。压缩包共59个文件体积仅213KB包含13个VB源码文件核心曲线类、节点插入、反算及2D绘图等模块、13个DLL、5个EXE可执行程序以及resources、xml等资源与工程文件sln/vbproj覆盖从算法实现到界面调用的完整链路。已有184人学习下载项目清晰地展示了控制点、权重、Knot向量与基函数如何协同作用并提供曲线求值、反算、节点插入等常用操作配套窗体界面可直接绘制2D曲线便于对照源码理解每一步计算。对于想深入NURBS算法或需要VB.NET图形编程参考的读者这份小巧的源码包是高效的学习与复用样本。1. 非均匀B样条节点向量一换曲线性格全变接手过图形算法相关工作的人大概率在某个版本里被B样条三个字卡住过。均匀B样条看着简单控制点一拖、次数一选、曲线就出来了可真要碰那些需要局部精细控制的设计场景——比如让曲线在某一段贴着控制多边形走、在另一段又保持光顺——均匀节点向量完全不够用。非均匀B样条的核心差异就一句话节点向量不再等距分布曲线在参数空间里有了重点关照区域。这份基于VB.NET的源码包解决的正是这个痛点它把节点向量构造、递推基函数、曲线求值这一整套流程封装成了可见、可调、可改的类代码适合正在做CAD/CAM小工具、或是需要把曲线算法搬进.NET桌面程序的人。打开工程你会发现它不只是一段能跑的Demo更是一个能让你看清非均匀三个字到底改了什么的教学型源码。2. 节点向量与基函数Cox-de Boor递推的VB.NET落地2.1 均匀与准均匀节点向量先跑通最简单的形态要理解非均匀得先看非均匀的反面。均匀节点向量就是参数区间内节点等距排列比如次数p3、控制点n17个节点向量就是[0, 1/10, 2/10, ..., 1]这样的等差数列。它的缺点是曲线在首尾段的行为不可控——真实工程里我们几乎总是希望曲线穿过首尾控制点这就引出了准均匀Clamped节点向量首尾各重复p1次节点值内部节点依然等距。这份VB.NET工程里构造准均匀节点向量的代码大致长这样Public Function BuildClampedKnotVector(n As Integer, p As Integer) As Double() Dim m As Integer n p 1 Dim knots(m) As Double Dim i As Integer For i 0 To p knots(i) 0.0 Next Dim interiorCount As Integer n - p For i 0 To interiorCount knots(i p 1) i / interiorCount Next For i m - p To m knots(i) 1.0 Next Return knots End Function注意看循环边界第一个循环把前p1个节点全部置0最后一个循环把末尾p1个节点全部置1内部节点从i/(n-p)生成——分母用n-p而不是n是因为首尾重复的节点不参与内部等距切分。这个函数返回的数组长度是m1等于np2个元素对应B样条曲线定义里的节点向量应有m1个节点。拿n6、p3举例节点向量是[0,0,0,0,0.25,0.5,0.75,1,1,1,1]首尾各重复4次。用这个Clamped向量跑出来的曲线首尾点必然落在第一个和最后一个控制点上这是做交互式曲线编辑器最基本的前提。2.2 非均匀节点分布核密度导向的切分策略真正意义上的非均匀节点向量是把上一步的等距内部节点换成由数据分布或设计要求决定的非等距值。常见做法有两种一种是根据控制点的弦长比例来分配节点值让参数区间的疏密和控制点折线的长度占比挂钩另一种是人工指定关键节点位置让曲线在某一段拥有更多的参数区间从而获得更灵活的形变能力。弦长参数化的核心思路先把所有相邻控制点之间的欧氏距离算出来并累加成总弦长然后按每个控制点累计弦长占总弦长的比例去决定对应节点值落在参数区间的哪个位置。代码实现如下Public Function BuildChordKnotVector(pts() As PointF, p As Integer) As Double() Dim n As Integer pts.Length - 1 Dim m As Integer n p 1 Dim knots(m) As Double Dim chordLen(n) As Double Dim total As Double 0.0 Dim i As Integer For i 0 To n - 1 chordLen(i) Distance(pts(i), pts(i 1)) total chordLen(i) Next For i 0 To p knots(i) 0.0 Next Dim acc As Double 0.0 For i 1 To n - p acc chordLen(i - 1) knots(i p) acc / total Next For i m - p To m knots(i) 1.0 Next Return knots End Function内部节点从索引p1开始赋值一直填到m-p-1为止。这里有个容易写错的地方累计弦长从第一条边开始累加但第一个内部节点对应的是第一个控制点之后的位置所以knots(ip) acc/total这条赋值语句里的偏移量p一定要与后面基函数计算里的节点下标保持一致否则基函数定义域与控制点对不上。弦长法做出来的曲线逼近效果比均匀节点向量好一大截尤其是控制点分布密度不均的时候——控制点密集的区域内部节点也密集曲线在这个区间会更有弹性。2.3 基函数递归计算下标范围与边界处理B样条基函数的递推是Cox-de Boor公式VB.NET里实现它需要格外注意递归深度和数组边界。第i个p次基函数N_i,p(u)依赖于u落在哪个节点区间[u_k, u_k1)里递推公式分两段N(i,0)(u) 1 若 u_i u u_{i1} N(i,p)(u) (u - u_i)/(u_{ip} - u_i) * N(i,p-1)(u) (u_{ip1} - u)/(u_{ip1} - u_{i1}) * N(i1,p-1)(u)代码落地时最实用的写法不是直接递归而是用循环从0次往p次推进Public Function BasisFunc(i As Integer, p As Integer, u As Double, knots() As Double) As Double Dim N(p) As Double Dim j As Integer, k As Integer For j 0 To p N(j) 0.0 Next For j 0 To p If u knots(i j) AndAlso u knots(i j 1) Then N(j) 1.0 Exit For End If Next For k 1 To p Dim saved As Double 0.0 For j 0 To p - k Dim left As Double u - knots(i j) Dim right As Double knots(i j k) - knots(i j) If right 0.00000001 Then Dim ratio As Double left / right N(j) ratio * N(j) (1.0 - ratio) * N(j 1) Else N(j) 0.0 End If Next Next Return N(0) End Function这里有个工程上容易出事的细节当节点区间宽度为0时即节点重复right会等于0直接除零导致程序崩溃。上面代码用了一个极小阈值0.00000001来做保护right小于这个值时直接跳过该段贡献——这是能稳定跑完所有非均匀节点组合的关键。另一个细节是参数u1.0时最后一个半开区间[u_k, u_k1)会把u排除在外所以最后一小段要做一个Clamp处理让u落在合法区间内否则基函数会计算出全0结果。实际工程里我习惯在调用BasisFunc之前统一做一次参数裁剪If u knots(knots.Length - 1) - 0.0000001 Then u knots(knots.Length - 1) ElseIf u knots(0) Then u knots(0) End If3. 从控制点到可绘制曲线工程结构与绘图流程3.1 工程结构与核心类一处计算、多处复用打开这份VB.NET源码工程建议先花十分钟看Solution Explorer里的文件结构不要一上来就翻Form1.vb。它的核心逻辑不在窗体里而是一个独立的BsplineCurve类。这种把算法与界面分离的结构正是做曲线工具时值得模仿的——算法类不依赖任何Form控件之后不管是把它搬进WPF、还是转成控制台测试都能直接复用。BsplineCurve类里的关键成员大致是这样的分工Public Class BsplineCurve Private _ctrlPts() As PointF 控制点数组 Private _knots() As Double 节点向量 Private _degree As Integer 曲线次数 Public Sub New(pts() As PointF, degree As Integer, knotType As KnotType) _ctrlPts pts.Clone() _degree degree Select Case knotType Case KnotType.Clamped _knots BuildClampedKnotVector(pts.Length - 1, degree) Case KnotType.Chord _knots BuildChordKnotVector(pts, degree) End Select End Sub Public Function Eval(u As Double) As PointF Dim result As PointF New PointF(0.0F, 0.0F) Dim n As Integer _ctrlPts.Length - 1 For i As Integer 0 To n Dim basis As Double BasisFunc(i, _degree, u, _knots) result.X basis * _ctrlPts(i).X result.Y basis * _ctrlPts(i).Y Next Return result End Function End ClassKnotType这个枚举是工程里一个很友好的设计——它把均匀、准均匀、弦长参数化这三种最常见的节点生成策略收敛到同一个入口切换节点类型只需要改一行构造参数。Eval函数是整条曲线的心脏它把所有控制点与基函数做加权和得到一个参数u对应的曲线坐标。理解这个函数后续想加权因子、想算导矢、想转NURBS都只是在这个求和式上做文章。3.2 曲线采样与绘制Graphics的缩放和平移曲线本身是连续的数学对象屏幕上只能画离散点。采样数量直接影响绘制效果太稀曲线显得折线化太密计算开销大。工程里默认的做法是按参数区间均匀取200个点对大多数桌面绘制场景这个密度足够。绘制流程如下Private Sub DrawCurve(g As Graphics, curve As BsplineCurve, scale As Single, offset As PointF) Dim step As Double 1.0 / 200 Dim prev As PointF ScalePoint(curve.Eval(0.0), scale, offset) Dim pts(199) As PointF Dim j As Integer 0 For u As Double step To 1.0 step * 0.5 Step step If u 1.0 Then u 1.0 Dim current As PointF ScalePoint(curve.Eval(u), scale, offset) pts(j) current prev current j 1 Next If j 1 Then g.DrawLines(Pens.Blue, pts.Take(j).ToArray()) End If End Function Private Function ScalePoint(src As PointF, scale As Single, offset As PointF) As PointF Return New PointF(src.X * scale offset.X, src.Y * scale offset.Y) End Function这里有一个绘制的坑直接把Eval出来的数学坐标丢给Graphics是不行的数学坐标原点在左下角或者设计数据里的任意位置而屏幕坐标原点在左上角、y轴向下。ScalePoint函数做的事就是把数学坐标线性映射到屏幕空间。工程里的offset通常取自PictureBox的尺寸一半scale则要根据控制点分布范围动态算——比如所有控制点X范围是0到1000那scale就取PictureBox宽度除以1000左右这样曲线才能完整可见不会飞出画布。4. 避坑指南节点重复度、边界与下标越界4.1 节点首尾重复不够曲线根本不经过首尾控制点现象用均匀节点向量构造p3、6个控制点的曲线期望首尾点落在控制点上结果曲线开头差了一大截看起来像被削掉了一块。原因均匀节点向量的首尾节点值不重复基函数在边界处不为1曲线端点与控制点没有约束关系。只有首尾节点各重复p1次端点基函数才会退化成只有第一个基函数为1其余全为0的情形。解决构造节点向量时首尾各重复p1次中间节点再按均匀或弦长分布。这是Clamped类型存在的根本价值。动手调试时可以先写个快速断言用代码输出u0.0处的曲线点坐标和第一个控制点坐标比对误差超过1e-6就说明节点向量没做Clamped。4.2 基函数递归里right等于零程序直接崩现象跑一组带重复节点的非均匀节点向量基函数计算抛出DivideByZeroException调用栈指向BasisFunc内部的那行除法。原因非均匀节点向量允许内部节点重复。当节点区间[u_i, u_ip]宽度为零时也就是两个节点值相等Cox-de Boor公式里的分数项分母恰好为0。这在均匀和Clamped情况下不会触发但一旦切换到真正的非均匀分布就必然碰到。解决在做除法之前判断分母绝对值是否大于一个极小阈值比如1e-8小于这个阈值时整项按0处理。这个阈值不能调太小否则浮点数精度不够时还是会偶尔触发异常也不能调太大否则曲线会失真。经验值是取1e-8到1e-10之间。4.3 参数u超出定义域基函数输出全零曲线点归零现象Eval函数传入u1.0结果返回一个(0,0)坐标点整条曲线端点塌掉。传入u-0.01程序不报错但曲线起点位置异常。原因B样条基函数的支撑区间是[u_i, u_{ip1}]半开区间定义导致u等于最大节点值1.0时最后一个基函数的值被排除所有基函数累加结果变成0。u小于起始节点值时对于Clamped类型前p1个节点都是0u-0.01落在任何区间之外基函数自然也全为0。解决在Eval入口处统一做参数裁剪u小于knots(0)时夹到knots(0)u大于knots(knots.Length-1)时夹到该最大值。只在绘制循环里写一次裁剪还不够凡是外部传入u的地方都要过一遍——否则后续加求导、加插值功能时这些脏参数还会继续作妖。4.4 采样点太稀曲线看起来像折线现象采样点数取50或80曲线整体形状对但过弯处有明显棱角和真实光滑曲线差距很大。原因B样条本身是光滑的但屏幕绘制是离散点连线。控制点密集、曲线曲率变化快的区域需要更多采样点来逼近真实曲线。均匀采样200个点在控制点不多时够用但控制点数量上到十几个后局部细节就会丢。解决采样策略从固定数量改成基于弦长误差的自适应采样——递归检查相邻两个采样点之间的弦与理论曲线中点的偏差偏差超过阈值时在这段中间补插一个采样点。工程里的简化方案是把采样数关联到控制点数量sampleCount max(200, (n1)*30)简单有效。4.5 浮点数累积误差多轮编辑后曲线偏移现象反复拖动控制点、切换节点类型、缩放平移之后曲线和最初形态相比出现了微小的整体偏移误差量级大约在1e-5左右。原因坐标数据经过多次浮点运算——控制点从屏幕坐标反算回数学坐标时做了除法曲线求值又做乘法和累加VB.NET的Single类型在累加求和时有效位数只有约7位误差在多次变换后累积。解决控制点内部存储统一用Double类型只在Graphics绘制时才转成PointF。Eval求和时优先用两两相加的树形求和代替顺序累加能把浮点误差再压掉一个量级。这一点在数据精度要求高的CAD二次开发场景尤其值得重视。5. 参数调节实战权因子、次数与参数化的配合5.1 权因子端点不变、内部可控非均匀B样条的非均匀主要通过节点向量实现但如果控制点带了一个权重w就升级成了NURBS——这是B样条的推广形式也是这份源码向实用工具迈进的关键一步。权因子的工程价值在于它能让曲线在完全不改动控制点位置的前提下向某个控制点靠近或远离。核心公式变化在求和式上C(u) Σ(N_i,p(u) * w_i * P_i) / Σ(N_i,p(u) * w_i)分母是全部加权基函数之和这个归一化步骤不能省。当所有w_i1时公式退化为普通B样条把某个w_i调大曲线在该控制点附近明显被吸过去调小则曲线远离该控制点。端点控制点的权因子建议永远保持为1否则曲线首尾会脱离控制点而且很难通过直觉预测偏移方向。5.2 次数与参数化的联动调试B样条次数p的选择直接影响曲线光滑性和局部支撑性。p1是折线p2是各段抛物线拼接p3是常用三次样条。次数调高会让曲线更光滑但计算量、节点区间对曲线的约束都会改变。最典型的联动现象把p从3改成5之后同一个控制点网络生成的曲线在端点处可能不再通过首尾控制点——因为回头检查节点向量首尾重复次数是p1如果构造时重复次数还是按原来p3硬编码的4那节点向量就不满足Clamped条件了。联调建议写在调试工具的界面上节点类型下拉框、次数输入框、权因子表格这三个控件放在同一行方便随时修改后刷新绘图。每条曲线绘制完成后在PictureBox里同时画出控制多边形——用灰色虚线连起来——这样能直观看到控制点没动但曲线变了的差异比只盯一个最终曲线判断参数效果要高效得多。6. 验证与进阶用解析解反查实现正确性拿到这份源码第一步不要把时间花在改界面皮肤上先做正确性验证。最省事的验证方式是利用一条已知解析解的曲线反查。任意取四个控制点例如(0,0)、(1,2)、(3,1)、(4,3)次数取2用Clamped节点向量生成曲线然后在u0.5处求值。手算虽然繁琐但你可以用另一个独立工具比如熟悉的数学库算出同一参数位置的理论值和程序输出对比误差在1e-4以内说明基函数和求和逻辑是对的。这个交叉验证是我每次拿到新的B样条代码都会强制过的一关。进阶方向上真正值得深挖的是插值。这份工程目前是逼近式B样条——曲线靠近但不一定通过所有控制点。工程实际中经常要求曲线精确通过给定型值点那就需要反求控制点给定数据点Q_k和参数值u_k建立线性方程组解出控制点P_i。方程组矩阵是带状的VB.NET里用高斯消元或者三元组存储的迭代法都能解不过注意节点向量这时要用弦长参数化重新生成直接套用均匀节点会导致方程组病态。延伸一层的做法是把权因子加进求值公式升级成NURBS曲线这样就能精确表示圆弧和二次曲线——这是很多CAD文件格式比如某跨平台图形交换格式里NURBS实体的基本要求。这份源码的类结构很适合做这个扩展Eval函数里增加一个w数组求值时分母分子分开算其余代码几乎不用动。最后想说一个习惯问题。早前我拆这类曲线源码时总喜欢先跑一遍看效果出问题了再回头翻代码结果花在找bug上的时间比读代码的时间还多。后来改成固定顺序先看节点向量构造函数、再用测试点验证求值结果、最后才上界面。从那以后我每次拿到新的B样条或NURBS代码都强制走一遍节点向量输出 → 已知点验证 → 边界参数测试这三步基本能过滤掉七八成实现错误。这样做不一定最快但最稳。希望帮到你。本文还有配套的精品资源点击获取