Nurbs样条线的原理及实现例码

2024-03-01 08:08

本文主要是介绍Nurbs样条线的原理及实现例码,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

最近做项目要用到到nurbs样条线,于是在网上找了一些Nurbs样条线的资料进行了研究,参照例子用自己的方法实现了Nurbs样条线,使用了几组例子进行测试。语言环境是OSG库。

(1)nurbs样条线的原理简介:

公式:

其中是控制顶点(构成控制多边形),为权因子,为定义于非均匀控制矢量上的次B-样条基函数,如下:


若未经说明,一般假定>0;

令:

上式可写为:


在这里给出链接地址,感兴趣的朋友可以看看:http://give.zju.edu.cn/cgcourse/new/book/8.6.htm

(2)实现的代码:

<NURBS基函数>

double N_Base_New(int i, int k, int t, double u, std::vector<float> U)//NURBS basis function
{
    double N[20];//保存基函数的中间变量
    double saved, temp;
    N[0] = 1.0;
    for (int j = 1; j <= t; j++)
    {
        saved = 0.0;
        for (int r = 0; r < j; r++)//
        {
            if (N[r] == 0 && (U[i + r + 1] - U[i + 1 - j + r]) == 0)
            {
                temp = 0;
            }
            else
            {
                temp = N[r] / (U[i + r + 1] - U[i + 1 - j + r]);
            }
            N[r] = saved + (U[i + r + 1] - u)*temp;
            saved = (u - U[i + 1 - j + r])*temp;

        }
        N[j] = saved;
    }
    return N[k];
}

<生成线顶点>

osg::ref_ptr<osg::Vec3Array> NB_Spline_New(int t, osg::ref_ptr<osg::Vec3Array> CP, std::vector<float> U)//参数:度数,控制点,节点向量
{
    osg::ref_ptr<osg::Vec3Array> curvePs = new osg::Vec3Array;//保存生成的顶点
    int Num = 50;
    auto Unumm = U.size();
    std::vector<float> NewU;//[0,1]之间的节点向量U
    std::vector<int> CPIndexs;
    for (int i = 0; i < Unumm; i++)
    {
        float tempu = U[i] / U[Unumm - 1];
        NewU.push_back(tempu);
    }
    for (double u = 0; u <= Num; u++)
    {
        float deltau = float(1.0) / Num*u;
        osg::Vec3 tempP(0.0, 0.0, 0.0);
        int i, j;

        if (u == Num)
        {
            auto Uvecnum = NewU.size();
            for (j = 0; j < Uvecnum - 1; j++)
            {
                if (NewU[j] <= deltau && deltau <= (NewU[j + 1]))
                {
                    i = j;
                    break;
                }
            }
        }
        else
        {
            auto Uvecnum = NewU.size();
            for (j = 0; j < Uvecnum - 1; j++)
            {
                if (NewU[j] <= deltau && deltau < NewU[j + 1])
                {
                    i = j;
                    break;
                }
            }
        }
        std::vector<osg::Vec3> tempPs;
        for (int k = 0; k <= t; k++)
        {
            osg::Vec3 tempPP = CP->at(i + k - t)*N_Base_New(i, k, t, deltau, NewU);//i:节点向量的下标索引
            tempPs.push_back(tempPP);
            CPIndexs.push_back(i + k - t);
        }
        auto Numm = tempPs.size();
        for (int k = 0; k < Numm; k++)
        {
            tempP += tempPs[k];
        }
        curvePs->push_back(tempP);
    }
    return curvePs;
}

(3)效果展示:

测试数据:


效果如图(红色为nurbs曲线,经过两端顶点,绿色为控制多边形顶点线):




这篇关于Nurbs样条线的原理及实现例码的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/761743

相关文章

hdu1043(八数码问题,广搜 + hash(实现状态压缩) )

利用康拓展开将一个排列映射成一个自然数,然后就变成了普通的广搜题。 #include<iostream>#include<algorithm>#include<string>#include<stack>#include<queue>#include<map>#include<stdio.h>#include<stdlib.h>#include<ctype.h>#inclu

深入探索协同过滤:从原理到推荐模块案例

文章目录 前言一、协同过滤1. 基于用户的协同过滤(UserCF)2. 基于物品的协同过滤(ItemCF)3. 相似度计算方法 二、相似度计算方法1. 欧氏距离2. 皮尔逊相关系数3. 杰卡德相似系数4. 余弦相似度 三、推荐模块案例1.基于文章的协同过滤推荐功能2.基于用户的协同过滤推荐功能 前言     在信息过载的时代,推荐系统成为连接用户与内容的桥梁。本文聚焦于

【C++】_list常用方法解析及模拟实现

相信自己的力量,只要对自己始终保持信心,尽自己最大努力去完成任何事,就算事情最终结果是失败了,努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一:什么是list? •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二:list常用接口 •🌰1.默认成员函数 🔥构造函数(⭐) 🔥析构函数 •🌰2.list对象

【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯,先赞后看哦~🎈🎈 🏆 作者简介:景天科技苑 🏆《头衔》:大厂架构师,华为云开发者社区专家博主,阿里云开发者社区专家博主,CSDN全栈领域优质创作者,掘金优秀博主,51CTO博客专家等。 🏆《博客》:Python全栈,前后端开发,小程序开发,人工智能,js逆向,App逆向,网络系统安全,数据分析,Django,fastapi

hdu4407(容斥原理)

题意:给一串数字1,2,......n,两个操作:1、修改第k个数字,2、查询区间[l,r]中与n互质的数之和。 解题思路:咱一看,像线段树,但是如果用线段树做,那么每个区间一定要记录所有的素因子,这样会超内存。然后我就做不来了。后来看了题解,原来是用容斥原理来做的。还记得这道题目吗?求区间[1,r]中与p互质的数的个数,如果不会的话就先去做那题吧。现在这题是求区间[l,r]中与n互质的数的和

让树莓派智能语音助手实现定时提醒功能

最初的时候是想直接在rasa 的chatbot上实现,因为rasa本身是带有remindschedule模块的。不过经过一番折腾后,忽然发现,chatbot上实现的定时,语音助手不一定会有响应。因为,我目前语音助手的代码设置了长时间无应答会结束对话,这样一来,chatbot定时提醒的触发就不会被语音助手获悉。那怎么让语音助手也具有定时提醒功能呢? 我最后选择的方法是用threading.Time

Android实现任意版本设置默认的锁屏壁纸和桌面壁纸(两张壁纸可不一致)

客户有些需求需要设置默认壁纸和锁屏壁纸  在默认情况下 这两个壁纸是相同的  如果需要默认的锁屏壁纸和桌面壁纸不一样 需要额外修改 Android13实现 替换默认桌面壁纸: 将图片文件替换frameworks/base/core/res/res/drawable-nodpi/default_wallpaper.*  (注意不能是bmp格式) 替换默认锁屏壁纸: 将图片资源放入vendo

C#实战|大乐透选号器[6]:实现实时显示已选择的红蓝球数量

哈喽,你好啊,我是雷工。 关于大乐透选号器在前面已经记录了5篇笔记,这是第6篇; 接下来实现实时显示当前选中红球数量,蓝球数量; 以下为练习笔记。 01 效果演示 当选择和取消选择红球或蓝球时,在对应的位置显示实时已选择的红球、蓝球的数量; 02 标签名称 分别设置Label标签名称为:lblRedCount、lblBlueCount

Kubernetes PodSecurityPolicy:PSP能实现的5种主要安全策略

Kubernetes PodSecurityPolicy:PSP能实现的5种主要安全策略 1. 特权模式限制2. 宿主机资源隔离3. 用户和组管理4. 权限提升控制5. SELinux配置 💖The Begin💖点点关注,收藏不迷路💖 Kubernetes的PodSecurityPolicy(PSP)是一个关键的安全特性,它在Pod创建之前实施安全策略,确保P

工厂ERP管理系统实现源码(JAVA)

工厂进销存管理系统是一个集采购管理、仓库管理、生产管理和销售管理于一体的综合解决方案。该系统旨在帮助企业优化流程、提高效率、降低成本,并实时掌握各环节的运营状况。 在采购管理方面,系统能够处理采购订单、供应商管理和采购入库等流程,确保采购过程的透明和高效。仓库管理方面,实现库存的精准管理,包括入库、出库、盘点等操作,确保库存数据的准确性和实时性。 生产管理模块则涵盖了生产计划制定、物料需求计划、