简介Nurbs3.0.11是一套面向VS2010环境的开源NURBS非均匀有理B样条库适合需要精确建模曲线曲面的C开发者可应用于CAD、三维建模与工程计算场景。压缩包共253个文件以130个cpp源码、43个h头文件为主体辅以dsp/vcxproj工程文件、测试相关文件等整体仅828KB结构紧凑便于直接纳入项目。“TestWin”提供基于单文档界面的集成示例演示在GUI中显示、编辑与交互操作NURBS对象“Test”及其他用例覆盖创建、裁剪、求交等常见功能nurbs核心代码封装了控制点、权重、阶数、knot向量等关键算法实现。通过学习这些源码开发者能掌握NURBS数据结构的底层逻辑并快速完成库的配置与调用。该资源已被1049人下载学习特别适合希望深入理解NURBS原理并落地到VS2010工程的开发者。1. 直接能用还是黑匣子Nurbs 3.0.11 在 VS2010 下的真实定位拿到的这个 NURBS 3.0.11 源码包是 VS2010 时代留存下来的 CAD 底层库源代码很多人双击打开 .sln 后第一反应是“这玩意能编过吗”。它不是一个教学用的玩具工程而是一整套能参与产品级开发的 NURBS 几何内核覆盖曲线曲面的创建、求值、升阶、插值、求交以及 3DM 文件读写。对做 CAD/CAM 软件、数控编程后处理、模型转换工具的人来说这套源代码相当于把几何建模的轮子直接递到你手上省掉从论文公式到可运行代码之间的黑匣子距离。它适合三类人要往自研软件里嵌入 NURBS 数据结构的工程师、需要读写 3DM 文件的工具开发者以及想在老 Windows 工程里搭一套可靠几何算法的人。2. 把库拆开看核心类、数据组织与 3DM 文件协议的三个支点NURBS 相关的开源库不少但这个包的价值在于它的类设计足够“底”它不假设你用它做 CAD 还是做 CAM只给你数学对象和文件 IO。打开主头文件 opennurbs.h你会看到一串 ON_ 前缀的类——ON_NurbsCurve、ON_NurbsSurface、ON_Brep、ONX_Model。这些类之间的层次关系决定了你在工程里怎么组织代码。2.1 从控制点到节点向量ON_NurbsCurve 的数据结构在 NURBS 的世界里曲线不是一堆坐标点而是一组控制点、一条节点向量、一组权重。openNURBS 系代码里ON_NurbsCurve 这个类装的核心字段有下面这些m_dim曲线所在空间的维数三维曲线传 3二维曲线传 2m_is_rat是否有理曲线等于 1 时控制点带权重等于 0 则权重全部默认为 1m_order阶数等于次数加 1三次曲线的 order 4m_cv_count控制点个数m_knot节点向量长度为 order m_cv_count - 1很多人第一次写创建代码时把 order 当成 degree 传结果节点向量长度怎么都对不上这是这个库最典型的第一坑。节点向量的基本约定是首尾重复度等于 order形成 clamped 向量让曲线经过首尾控制点。整条向量长度必须严格等于 order m_cv_count - 1。比如 5 个控制点的三次曲线节点向量应该有 4 5 - 1 8 个数也就是两个方向各重复 4 次中间夹一个内部节点。ON_NurbsCurve* curve new ON_NurbsCurve(); int dim 3; int is_rat 0; // 0 表示非有理曲线控制点不带权重 int order 4; // 三次曲线order degree 1 int cv_count 5; // 控制点数量 // 节点向量clamped 均匀分布长度 order cv_count - 1 double knot[8] { 0.0, 0.0, 0.0, 0.0, 0.5, 1.0, 1.0, 1.0 }; // 控制点5 个三维点按 xyz 顺序平铺 double cv[15] { 0.0, 0.0, 0.0, 1.0, 2.0, 0.0, 2.0, -1.0, 0.0, 3.0, 1.0, 0.0, 4.0, 0.0, 0.0 }; BOOL ok curve-Create(dim, is_rat, order, cv_count, knot, cv); if (!ok) { // 参数不合法通常是节点向量不等长或控制点数组长度不对 delete curve; curve nullptr; }这段代码里的 Create 是这个库创建曲线的标准入口。注意前两个参数是 dim3三维和 is_rat0非有理紧跟着的是 order 而不是 degree这是最容易看走眼的地方。节点向量我写的是“0,0,0,0,0.5,1,1,1,1”对应 5 个控制点的三次曲线头四个 0 保证曲线从第一个控制点出发尾四个 1 保证曲线终止于最后一个控制点中间 0.5 是唯一的内部节点把参数区间切成两段使曲线在两个区间上分别由不同的基函数支撑。这段代码创建出来的曲线在节点 0.5 处是 C2 连续的因为内部节点没有重复。如果想让曲线在某处出现折角就把那个节点的重复度加到 3。控制点数组按 xyz 平铺5 个点就是 15 个 double。这个内存布局在整个库里是一致的不管是求值还是求交底层都在和这种平铺数组打交道。处理好这一步后续所有几何运算都建立在这个数据结构上。2.2 曲面与裁剪ON_NurbsSurface 的二维参数域曲面比曲线多一个维度。ON_NurbsSurface 里有两个方向的 order、cv_count、knot控制点从一维列表变成网格cv_count[0] x cv_count[1]。创建曲面时第一参数是维数第二参数是有理标志然后是 U 向 order、U 向 CV 数、V 向 order、V 向 CV 数接着是 U 向节点向量、V 向节点向量、控制点数组。控制点数组按 U 向优先平铺也就是先走完 U 向一行再换 V 向下一行。这里要补一个曲面里的关键概念裁剪。NURBS 曲面本身是一个矩形参数域但 Rhino 里的修剪面不是矩形网它由底层未裁剪曲面加一个二维参数域里的内外环组成。修剪信息不写在 ON_NurbsSurface 类里而是在 ON_Brep 的 trim 表里。ON_Brep 上挂着 ON_BrepLoop、ON_BrepTrim 这些对象它们负责记录“哪条曲线是边界、边界在 U/V 参数域的哪一侧”。对比项ON_NurbsCurveON_NurbsSurface阶数一个 m_orderU/V 两个 m_order节点向量一条 m_knotU/V 两条 m_knot控制点cv_count 个点cv_count[0] x cv_count[1] 点阵裁剪信息无无裁剪在 ON_Brep求值参数t(u, v)工程里做模型修复时必须明白这个边界想导出在 Rhino 里看到的那种带孔带边的面必须走 Brep 相关接口不能只拿一个 ON_NurbsSurface 硬拼。否则你会发现导出的 3DM 打开后是一个完整矩形面孔洞全没了。2.3 3DM 文件读写为什么这版库值得留一个副本3DM 是 Rhino 的原生格式也是这个库自带读写支持最好的格式之一。所谓“值得留副本”是因为新版本库读不进的老 3DM 文件换到这一版反而经常能进去。我遇到过几次这种情况老犀牛文件里有些历史遗留的对象属性新版 openNURBS 升级时把读取路径改掉了而这版还保留着旧路径。做历史数据转换工具时这个源码包是少数能两头兼顾的选择。3DM 读写围绕 ONX_Model 展开。ONX_Model 相当于一个内存里的模型容器往里塞 ON_Geometry 对象再整体写盘。写盘代码看起来很简单ONX_Model model; ONX_Model_Object obj; obj.m_object curve; // curve 是上面创建出来的对象 obj.m_bDeleteObject true; // true模型接管对象所有权 model.m_ModelObjectArray.Append(obj); FILE* fp ON_FileStream::Open(Ldemo.3dm, Lwb); if (fp) { model.Write(fp, 0, nullptr); // 第二参 0 表示当前库默认文件版本 ON_FileStream::Close(fp); }注意 m_bDeleteObject 这个字段。设为 true 后model 析构时会帮你 delete curve如果你后面又手动 delete就是双重释放。老版本库里这种内存所有权边界没有强制约定全靠调用方自觉。我在工程里的习惯是把它当成一种“所有权移交”一旦设了 true就不再碰原指针。ON_FileStream::Open 的宽字符版本在某些 VS2010 配置下要确认工程字符集是 Unicode否则文件打不开这个问题第 5 章会再展开。3. 在 VS2010 里把它编译成库工程选择、运行库开关与迁移到新 IDE 的注意项老代码到手第一件事不是写业务代码而是把工程跑通。NURBS 3.0.11 这套源码的构建不算难但有几个位置不注意会原地转圈。3.1 源码包里的目录组织与应该打开哪个工程文件这个版本源码包里一般能见到 opennurbs、zlib 两个核心目录外加按编译器版本区分的工程目录。VS2010 对应 vc10 目录下的 opennurbs.sln。不要看到有 vc12、vc14 就顺手打开VS2010 打不开高版本工程文件会直接报“需要较高版本的 Visual Studio”。nurbs311/ ├─ opennurbs/ # 核心源码opennurbs.h 在这里 ├─ zlib/ # 压缩依赖写 3DM 文件时用到 └─ vc10/ # VS2010 工程目录 ├─ opennurbs.sln ├─ opennurbs.vcxproj └─ opennurbs_dll.vcxproj用 VS2010 打开 opennurbs.sln 后解决方案里会有 opennurbs、opennurbs_dll、zlib 等多个项目。构建顺序有讲究先把 zlib 单独生成再生成 opennurbs。右键项目选择“仅用于项目 - 仅生成 zlib”通常会看到它生成出一个静态库然后回到 opennurbs 项目生成。如果你跳过 zlib 直接编 opennurbs链接时会报一堆无法解析的外部符号因为压缩库还没编出来。现在很多人电脑上已经没有 VS2010 了。我一般用 VS2019 或 VS2022 打开 vc10 工程装一个 v140 工具集就能编译更省事的做法是用 VSCode 配合 CMake 重写一个构建脚本把这批源码当成“参考实现”构建速度反而更快。CMakeLists 不用写复杂核心就是指定 opennurbs 源码目录、include 路径、zlib 路径三个事情。3.2 静态库与 DLL 的选择一个宏决定编译姿态openNURBS 系的 DLL 需要定义 OPENNURBS_DLL 宏否则默认按静态库编译。这个宏影响符号导入导出也影响类内变量的内存布局。宏不一致时最常见的结果是链接报 LNK2019一堆 ON_NurbsCurve 的成员函数无法解析。配置点静态库DLL宏定义不定义 OPENNURBS_DLL定义 OPENNURBS_DLL链接产物opennurbs.libopennurbs.dll opennurbs.lib运行库建议与调用方一致即可配合 /MD 使用适用场景独立 exe、内部工具多个插件共享几何代码选择标准我一般看产品形态独立工具用静态库链接最简单发布时少带一个 DLL插件化产品里多个模块都要用几何代码用 DLL 可以共享一份实现避免每个插件各自持有一份 ON_NurbsCurve 的静态拷贝导致内存堆不一致。但 DLL 方式要求所有插件和宿主程序的运行库一致否则跨模块传递对象时析构函数可能在错误的堆上执行。还有一个容易忽略的字符集问题老工程默认使用多字节字符集新建的调用方工程如果默认 Unicode你会发现 ON_String 和 std::string 互传时乱码或编译报错。在项目属性 - 常规 - 字符集里把调用方工程也切到多字节能省掉很多麻烦。3.3 编译过程中最容易翻车的三个地方下面三条是我在多个版本上反复遇到的按现象、原因、解决写清楚。现象 1编译时报 error C4996提示 fopen、strcpy 被标记为 deprecated。原因VS2010 之后编译器默认对这些 CRT 函数警告老代码里大量使用它们。解决在预处理定义里加 _CRT_SECURE_NO_WARNINGS一条宏解决不需要改源码。现象 2链接时报 LNK2019大量 ON_NurbsCurve::Create、ON_NurbsCurve::~ON_NurbsCurve 无法解析。原因库编译成静态库但调用方工程定义了 OPENNURBS_DLL或者反过来库编成 DLL 但调用方没定义。解决检查两边的预处理宏保证完全一致然后重新编译。现象 3程序运行时在对象析构处崩溃但编译链接都正常。原因调用方工程设置了内存对齐 /Zp1 或 /Zp2而库工程用默认对齐。类成员变量布局变了析构时按错误的偏移量释放内存。解决调用方工程关掉自定义对齐恢复默认 8 字节对齐。编译通过只是第一步。接下来要做的不是立刻写业务代码而是把你自己的调用工程先配好跑一个最小示例验证库的功能。4. 最小可用工程建一条三次 NURBS 曲线、求值、再导出 3DM任何一个库拿到手里我都会先写一个最小工程验证三件事能不能建对象、能不能算数据、能不能落盘。下面三节就是这套验证流程。4.1 创建一条精确半圆弧有理曲线的黄金用例验证曲线功能最好的用例是半圆弧因为它的结果可以直观检查而且只有有理曲线才能精确表示圆弧非有理曲线永远得到的是逼近弧线的多项式曲线。这里创建一个 3 个控制点的二维半圆弧ON_NurbsCurve* arc new ON_NurbsCurve(); int dim 3; int is_rat 1; // 有理曲线每个控制点带权重 int order 3; // 二次曲线即可精确表达圆弧 int cv_count 3; // 节点向量首尾各重复 order 次 double knot[6] { 0.0, 0.0, 0.0, 1.0, 1.0, 1.0 }; // 控制点xyz 权重 w按 xyzw 平铺 // 这是一个从 (1,0,0) 到 (0,1,0) 的四分之一圆弧 double cv[12] { 1.0, 0.0, 0.0, 1.0, // 起点权重 1 1.0, 1.0, 0.0, 0.7071067811865476, // 中间点权重 cos(45°) 0.0, 1.0, 0.0, 1.0 // 终点权重 1 }; BOOL ok arc-Create(dim, is_rat, order, cv_count, knot, cv); if (!ok) { delete arc; return; }中间控制点的权重是 0.7071067811865476也就是根号 2 除以 2。这是 NURBS 表示圆弧的标准公式这个权重值决定了曲线在中间点处“往里收”多少。如果你把权重写成 1得到的是一条二次多项式曲线形状会明显偏离圆弧往外鼓出去。权重写错是这个库里最典型的“看着差不多、实际差很多”的坑。节点向量只有三个 0 和三个 1表示整条曲线就在这一个节点区间上没有内部节点。验证方法很简单对曲线求值取 t0.5得到的点应该约等于 cos45° 和 sin45° 的位置也就是 (0.7071, 0.7071, 0)。如果偏离超过 1e-4权重或控制点大概率写错了。4.2 求值与导数给参数 t要点和切线工程里最刚性的需求是给一个参数 t算曲线上的点和切向。数控加工里刀位点就是沿着曲线按步长求值出来的切向决定走刀方向。这个库的求值接口是 Evaluatedouble t 0.5; ON_3dPoint pt; ON_3dVector tangent; // der_count1 表示算点 一阶导side1 表示从右侧趋近 BOOL ok arc-Evaluate(t, 1, 1, pt.x, tangent.x); if (ok) { // pt 是曲线上的点tangent 是参数方向的一阶导 double length tangent.Length(); // 加工用的单位切向需要自己归一化 ON_3dVector unit_tangent tangent / length; }这里的参数说明很重要。der_count 传 1表示同时计算位置和速度向量传 0 只算位置。side 参数是求导方向1 表示右极限-1 表示左极限。在节点重复度比较高的位置左右导数可能不同比如一条曲线的内部节点重复度为 3 时那个参数位置左右两侧的切线方向是断开的side 参数决定你拿到哪一侧的值。CAM 后处理里走刀方向必须连续这时要主动指定 side不能让库自己选。Evaluate 输出的 tangent 一般不是单位向量它的长度和参数化速度有关。做刀轨前先 normalize这是很多初学这个库的人容易漏掉的一步。4.3 导出 3DM让 Rhino 能直接打开你的结果建好对象、算完数据最后验证落盘。前面第 2 章写过写盘代码这里把读写完整串起来ONX_Model model; ONX_Model_Object obj; obj.m_object arc; obj.m_bDeleteObject true; // 模型接管所有权之后不能再 delete arc model.m_ModelObjectArray.Append(obj); FILE* fp ON_FileStream::Open(Larc.3dm, Lwb); if (fp) { model.Write(fp, 0, nullptr); ON_FileStream::Close(fp); } // 读取验证 ONX_Model in_model; fp ON_FileStream::Open(Larc.3dm, Lrb); BOOL read_ok in_model.Read(fp, 0, nullptr); ON_FileStream::Close(fp); if (read_ok) { for (int i 0; i in_model.m_ModelObjectArray.Count(); i) { const ONX_Model_Object obj_read in_model.m_ModelObjectArray[i]; const ON_NurbsCurve* c ON_NurbsCurve::Cast(obj_read.m_object); if (c) { // 拿到了读回来的曲线可以继续求值 } } }Cast 是这个库做安全类型转换的标准方式从 ON_Object 转到 ON_NurbsCurve转换失败返回 nullptr。读取老版本 3DM 文件时如果 Read 返回失败先别怀疑文件损坏——检查一下 ON_FileStream::Open 的字符参数到底用的是宽字符还是 UTF-8 窄字符这个库不同版本对文件名的编码处理不一样VS2010 工程默认多字节字符集下宽字符路径会静默失败。这是我在文件读写上踩过的最隐蔽的坑。5. 避坑与排查NURBS 库用稳的几个隐藏边界编译过了、例子跑通了不等于生产环境安全。下面四条都是我在实际工程里反复遇到的边界问题每条都按现象、原因、解决写清楚。5.1 节点向量重复度不足端点尖点是这么来的现象控制点和插值目标完全一致但求值出的曲线首末端点不经过控制点端部出现明显的尖角。原因一条 order4 的曲线首节点必须至少重复 4 次才能把首控制点“钉死”在端点。如果首节点只重复 3 次第一条曲线段的支撑区间还没收敛端点自然偏离控制点。解决打印节点向量前四个数和后四个数确认 clamp 是否正确。代码里用循环填节点向量时少写一个 0 是最隐蔽的我会写一个校验函数创建后立刻检查 m_knot[0] m_knot[order-1] 且尾部对称。5.2 对有理曲线直接做仿射变换形状悄悄变现象对半圆弧控制点做平移和缩放得到的曲线不再是圆弧而是一条不可预测的高次曲线。原因有理曲线的几何坐标是控制点除以权重后的结果。直接变换 xyz 忽略了齐次坐标里的 w 分量等于把曲线重新参数化了。解决调用 curve-Transform(ON_Xform)库内部会做齐次变换不要自己循环 control point 手工改坐标。这个库内置了 Transform 接口属于“后悔药”型功能能自动处理权重省去手写齐次运算的麻烦。5.3 内存所有权混乱谁接管谁释放必须让所有人知道现象程序退出时崩溃调用栈停在 ON_NurbsCurve 的析构函数里。原因ONX_Model_Object 的 m_bDeleteObject 设成了 true模型析构时把曲线对象删了但业务代码在后面又 delete 了一次同一指针导致双重释放。解决在代码注释第一行写明所有权归属。凡是把裸指针塞进 ONX_Model要么 m_bDeleteObject 设为 false 并自己管理要么设为 true 后绝对不再碰原指针。团队协作时这个字段的值应该写进代码评审清单。5.4 平台位数混用DLL 加载就崩现象编译好的工具在 32 位程序里运行正常放到 64 位目录下加载动态库立即崩溃。原因可执行程序是 x86DLL 是 x64。VS2010 工程默认活动平台是 Win32很多人后来在 x64 平台下也生成了库但调用方工程没切换两类文件混在一个目录里运行时按错误位数加载。解决打开“配置管理器”把调用方工程和库工程的活动解决方案平台保持一致x86 和 x64 的产物分别输出到单独目录不要混放。VS2010 版本里平台切换不显眼每次全量重建后要手工确认输出目录里的 DLL 位数。这四条坑有一个共同点它们都不是 NURBS 数学本身的问题而是数据结构边界和构建配置的问题。这也说明这套库的数学部分是可靠的出问题的大多在使用边界。6. 进阶用控制点修改代替重建做轻量级形状编辑6.1 SetCV局部修形最被低估的接口很多人在做形状调整时选择删掉整条曲线、改参数、重新 Create。其实对于局部微调直接改控制点就够。B 样条有一个局部支撑性质移动第 i 个控制点只影响周围 order 个节点区间内的形状其他部分完全不变。这个性质在工程里极其有用。int cv_index 2; ON_3dPoint newpt(2.5, 0.5, 0.0); double weight 1.0; // 只改第 2 个控制点的位置和权重 BOOL ok curve-SetCV(cv_index, newpt.x, newpt.y, newpt.z, weight);SetCV 之后不需要重建对象也不需要重新做参数化。这条曲线在修改点附近的形状会变化远处的控制点完全不受影响。相比删除重建这种方式速度快而且保留了原来的节点向量和参数化结构不会引入意外的参数分布变化。做交互式造型预览时我用它做实时控制点拖拽效果很直接。6.2 点云拟合的调用路径另一个常见需求是把测量点拟合成曲线。库里有插值方向的构造入口但不同小版本的函数名有出入我一般先在头文件里 grep 关键字 Interpolate 或 Fit找当前版本实际可用的函数签名。拟合时有两个参数最影响结果容差和参数化方式。容差设 1e-3 和设 1e-6出来的控制点数量差别很大刀轨密度也随之变化。参数化用弦长还是向心参数化决定曲线在急弯处的表现。这两项没有绝对标准必须结合你的数据采样密度试算。我做数控加工前处理时会针对每个工件单独调一遍容差再批量出刀路。那次在客户现场用老库读一批遗留 3DM内存所有权没理清楚程序崩得毫无脾气。从那以后我每次拿到一个新库版本第一件事就是先写一个“对象所有权表”再编译官方自带例子全部通过才碰业务代码。这个习惯帮我避开了很多不明不白的崩溃希望帮到你。本文还有配套的精品资源点击获取