基于bSpline插值的直纹面和旋转曲面

最新推荐文章于 2024-05-17 11:58:07 发布
最新推荐文章于 2024-05-17 11:58:07 发布
阅读量1.2k 收藏 2
点赞数 1
分类专栏： OpenGL学习 文章标签： vector float delete null vb class
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ITvsET/article/details/5924858
版权
该程序通过bSpline插值算法创建旋转曲面和直纹面，包括生成点、法线和面片。旋转曲面类（Revsurf）与直纹面类（Rulesurf）继承自bspline类，并实现相应的几何形状。程序通过视图界面进行交互，左侧视图更新参数，右侧视图实时绘制。视图间的通信由消息传递完成，当参数改变时，触发重绘以更新图形。
摘要由CSDN通过智能技术生成

此程序主要分为两部分，一部分是抽象旋转曲面和直纹面的类，另一部分是视图界面的类组织，二者通过在视图界面定义数据成员建立联系。

 

第一部分：

旋转曲面类：

 

class Revsurf:public bspline,public Mesh                         //旋转曲面

{

public:

       Revsurf(Point3* pt1,int num,int uu,int vv,int mm);             //构造函数

       ~Revsurf();

public:

       void writeMesh(char *fname);                             //将生成的点、法线、面片写入文件

       void makefaces();                                       //生成面片及索引

       void makenormals();                                     //生成法线及索引

       int ctlpt_num;                                          //控制点的个数

       int v;                                                  //经度方向

       int u;                                                  //纬度方向

       Point3 *ctlpt;                                           //保存控制点

       void makepoints();                                       //生成点及索引

};

 

Revsurf::Revsurf(Point3 *pt1,int num,int uu,int vv,int mm)                                          //旋转曲面

:ctlpt(pt1),ctlpt_num(num),u(uu),v(vv),bspline(mm,num)

{

       numVerts=u*v;

       numNorms=u*v;

       numFaces=(u-1)*v;

      

    pt=new Point3[numVerts];

       norm=new Vector3[numNorms];

       face=new Face[numFaces];

}

 

Revsurf::~Revsurf()

{

}

 

/******************************旋转曲面生成顶点*********************************/

void Revsurf::makepoints()

{

       int i=0,j=0,k=0;

       float t=0,x=0,y=0,z=0;

       Point3 *p;                                        //记录两条u线的点

       p=new Point3[u];

 

       //细节处理1e-5

       for(t=knot[0]+1e-5;i<u;t+=knot[ctlpt_num+m-1]/(float)(u-1)-1e-5)                           //采集两条u线的点

       {

              p[i].set(0,0,0);

              for(k=0;k<ctlpt_num;k++)

              {

                     p[i].x+=ctlpt[k].x*getbsplineValue(k,m,t,knot);

                     p[i].y+=ctlpt[k].y*getbsplineValue(k,m,t,knot);

                     p[i].z+=ctlpt[k].z*getbsplineValue(k,m,t,knot);

              }

              i++;

       }

 

       glBegin(GL_LINE_STRIP);

       for(i=0;i<u;i++)

              glVertex3f(p[i].x,p[i].y,p[i].z);

       glEnd();

 

       for(j=0;j<v;j++)                                           //生成网格点i*v+j                                  

              for(i=0;i<u;i++)

              {

                     x=p[i].x*cos(j/(float)v*2*PI);

                     y=p[i].y;

                     z=p[i].x*sin(j/(float)v*2*PI);

                     pt[j*u+i].set(x,y,z);

                     pt[j*u+i].set(j*u+i);

              }

              delete[] p;

              p=NULL;

}

 

/******************************旋转曲面生成法线*********************************/

void Revsurf::makenormals()

{

       Vector3 va(0,0,0),vb(0,0,0);

       int i=0,j=0;

 

       for(;j<v-1;j++)                                        

       {

              for(i=0;i<u-1;i++)

              {

                     va.set(pt[j*u+i+1].x-pt[j*u+i].x,pt[j*u+i+1].y-pt[j*u+i].y,pt[j*u+i+1].z-pt[j*u+i].z);

                     vb.set(pt[(j+1)*u+i].x-pt[j*u+i].x,pt[(j+1)*u+i].y-pt[j*u+i].y,pt[(j+1)*u+i].z-pt[j*u+i].z);

 

            va.normalize();

                     vb.normalize();

 

                     norm[j*u+i].set(va.cross(vb));

            norm[j*u+i].set(j*u+i);

              }

 

              va.set(pt[j*u+i].x-pt[j*u+i-1].x,pt[j*u+i].y-pt[j*u+i-1].y,pt[j*u+i].z-pt[j*u+i-1].z);

              vb.set(pt[(j+1)*u+i].x-pt[j*u+i].x,pt[(j+1)*u+i].y-pt[j*u+i].y,pt[(j+1)*u+i].z-pt[j*u+i].z);

 

              va.normalize();

              vb.normalize();

 

              norm[j*u+i].set(va.cross(vb));

              norm[j*u+i].set(j*u+i);

       }

 

       for(i=0;i<u-1;i++)

       {

              va.set(pt[j*u+i+1].x-pt[j*u+i].x,pt[j*u+i+1].y-pt[j*u+i].y,pt[j*u+i+1].z-pt[j*u+i].z);

              vb.set(pt[0*u+i].x-pt[j*u+i].x,pt[0*u+i].y-pt[j*u+i].y,pt[0*u+i].z-pt[j*u+i].z);

             

              va.normalize();

              vb.normalize();

             

              norm[j*u+i].set(va.cross(vb));

              norm[j*u+i].set(j*u+i);

       }

 

       va.set(pt[j*u+i].x-pt[j*u+i-1].x,pt[j*u+i].y-pt[j*u+i-1].y,pt[j*u+i].z-pt[j*u+i-1].z);

       vb.set(pt[0*u+i].x-pt[j*u+i].x,pt[0*u+i].y-pt[j*u+i].y,pt[0*u+i].z-pt[j*u+i].z);

      

       va.normalize();

       vb.normalize();

      

       norm[j*u+i].set(va.cross(vb));

       norm[j*u+i].set(j*u+i);

}

 

/******************************旋转曲面生成面片*********************************/

void Revsurf::makefaces()

{

       int i=0,j=0;

      

       for(;j<v-1;j++)

       {

              for(i=0;i<u-1;i++)

              {

                     face[j*(u-1)+i].nVerts=4;

                     face[j*(u-1)+i].vert=new VertexID[4];

                    

                     face[j*(u-1)+i].vert[0].normIndex=j*u+i;

                     face[j*(u-1)+i].vert[1].normIndex=j*u+i+1;

            face[j*(u-1)+i].vert[2].normIndex=(j+1)*u+i+1;

                     face[j*(u-1)+i].vert[3].normIndex=(j+1)*u+i;

                    

            face[j*(u-1)+i].vert[0].vertIndex=j*u+i;

                     face[j*(u-1)+i].vert[1].vertIndex=j*u+i+1;

            face[j*(u-1)+i].vert[2].vertIndex=(j+1)*u+i+1;

                     face[j*(u-1)+i].vert[3].vertIndex=(j+1)*u+i;

              }

       }

      

       for(i=0;i<u-1;i++)

       {

              face[j*(u-1)+i].nVerts=4;

              face[j*(u-1)+i].vert=new VertexID[4];

             

              face[j*(u-1)+i].vert[0].normIndex=j*u+i;

              face[j*(u-1)+i].vert[1].normIndex=j*u+i+1;

              face[j*(u-1)+i].vert[2].normIndex=0*u+i+1;

              face[j*(u-1)+i].vert[3].normIndex=0*u+i;

             

              face[j*(u-1)+i].vert[0].vertIndex=j*u+i;

              face[j*(u-1)+i].vert[1].vertIndex=j*u+i+1;

              face[j*(u-1)+i].vert[2].vertIndex=0*u+i+1;

              face[j*(u-1)+i].vert[3].vertIndex=0*u+i;

       }

}

 

void Revsurf::writeMesh(char *fname)

{

       Mesh::writeMesh(face,norm,pt,numFaces,numNorms,numVerts,fname);

}

 

直纹面类：

 

class Rulesurf:public bspline,public Mesh                            //直纹面

{

public:

       Rulesurf(Point3* pt1,Point3* pt2,int num,int uu,int vv,int mm);       //构造函数

       ~Rulesurf();                                                //析构函数

public:

       void writeMesh(char *fname);                                  //将生成点、法线、面片 写入文件

最低0.47元/天 解锁文章
ITvsET
关注
专栏目录
B样条插值（B-spline interpolation）
04-12
该文档讲述了B样条插值在图像处理中的应用，内容翔实，可供参考。
Bspline
weixin_40662209的博客
06-29 1506
样条函数 我们插值的时候，常见的高次多项式插值会遇到龙格现象： 故一般插值避免超过7次的高次插值； 所以有分段线性插值，分段线性插值光滑性差，即节点处可能不连续； 分段三次埃尔米特插值虽然一阶光滑，但插值条件过多； 三次样条插值函数就是分段三次多项式，在每个节点处有二阶连续导数，在实际中被广泛使用。 b样条 其中bj,n−1(t)b_{j,n-1}(t)bj,n−1​(t)就是其系数，m+1...
bspline_B样条插值_
09-30
一种新的插值拟合点的样条曲线，称为B样条曲线 不同于三次样条曲线
B-spline_B样条插值程序
04-09
很多人都想要的B样条插值程序,VC++界面写的，一看就会。。
b-spline:B样条插值
05-01
b样条 B样条插值 使用对任意维度的控制点进行B样条插值。 插值器可以采用可选的权重向量，如果需要，可以使所得曲线成为非均匀有理B样条曲线（NURBS）曲线。 结向量也是可选的，并且当未提供时，将在内部生成未夹紧的均匀结向量。 安装 $ npm install b-spline 例子 钳结矢量 var bspline = require ( 'b-spline' ) ; var points = [ [ - 1.0 , 0.0 ] , [ - 0.5 , 0.5 ] , [ 0.5 , - 0.5 ] , [ 1.0 , 0.0 ] ] ; var degree = 2 ; // As we don't provide a knot vector, one will be generated // internally and have the followin
曲线和曲面Bspline拟合
04-25
曲线和曲面Bspline拟合代码，基于python numpy开发
BSpline插值算法
01-19
BSpline插值算法具有良好的平滑性和局部逼近性能，因此在图像处理和计算机图形学中得到广泛应用。 BSpline插值算法的基本思想是通过在每个数据点附近定义一个局部的多项式函数来逼近数据。这些多项式函数称为基函数...
ScatteredTransform:这是3D Slicer扩展，可通过使用多级BSpline插值算法从散布点处定义的位移场创建BSpline变换
04-17
【ScatteredTransform】是一个专为3D Slicer设计的扩展模块，它的主要功能是利用多级BSpline插值算法从散乱分布的点上的位移场数据生成BSpline变换。这个工具对于处理非均匀分布的点云数据尤其有用，比如在医学图像...
matlab经典全集（内含插值源带码）B-spline插值代码
05-26
matlab经典全集 （内含插值源带码） B-spline插值代码
计算机图形学-Bezier与Bspline曲线，曲面
01-05
MIT计算机图形学作业 C++实现,Bezier曲线，Bspline曲线
三次样条插值代码
06-12
#include <iostream> #include <cmath> #include <fstream> using namespace std; class scyt { float *x,*y,*d,*h,*u,*q,*a,*b,*c,*l,*r,*o,*M; int m; float y0,y3; public: scyt(); void qiudao(); void zgf(); void qiujie(); ~scyt(); }; void main() { scyt hello; hello.qiudao(); hello.zgf(); hello.qiujie(); } scyt::scyt() { ifstream fin("三次样条插值.txt"); for(float j;fin>>j;) { m=int(j); break; } x=new float[m]; y=new float[m]; d=new float[m]; h=new float[m-1]; u=new float[m-2]; q=new float[m-2]; a=new float[m-1]; b=new float[m]; c=new float[m-1]; l=new float[m]; r=new float[m-1];//此处的r为追赶法中的u； o=new float[m];//此处o为追赶法中的y M=new float[m];//此处M为追赶法中的x； int jishu=0; for(j;fin>>j;) { if(jishu<=m-1) x[jishu]=j; if(jishu>m-1&&jishu;<2*m) { y[jishu-m]=j; } if(jishu==2*m) { y0=j; } if(jishu==2*m+1) { y3=j; } jishu++; } fin.close(); } void scyt::qiudao() { for(int i=0;i<m-1;i++) { h[i]=x[i+1]-x[i]; } for(i=0;i<m-2;i++) { u[i]=h[i] / (h[i] + h[i+1]); } for(i=0;i<m-2;i++) { q[i]=1-u[i]; } d[0]=6/h[0]*((y[1]-y[0])/h[0]-y0); for(i=1;i<m-1;i++) { d[i]=6/(h[i-1]+h[i])*((y[i+1]-y[i])/h[i]-((y[i]-y[i-1])/h[i-1])); } d[m-1]=6/h[m-2]*(y3-(y[m-1]-y[m-2])/h[m-2]); } void scyt::zgf() { u[m-2]=1; for(int i=0;i<m;i++) { b[i]=2; } c[0]=1; for(i=1;i<m-1;i++) { c[i]=q[i-1]; } //........................................ l[0]=b[0]; for(i=0;i<m-1;i++) { r[i]=c[i] / l[i]; l[i+1]=b[i+1] - (u[i] * r[i]); } o[0]=d[0] / l[0]; for(i=1;i<m;i++) { o[i]=(d[i]-u[i-1]*o[i-1]) / l[i]; } M[m-1]=o[m-1]; for(i=m-2;i>=0;i--) { M[i]=o[i]-r[i] * M[i+1]; } cout<<m<<"个点的导数值分别是："<<endl; for(i=0;i<m;i++) { cout<<"M"<<i+1<<"="; cout<<M[i]<<endl; } //M的值求出。。。。。。追赶法调用完毕 } void scyt::qiujie() { float S; for(;;) { float f; cout<<"请输入待求x的值(输入1000)时退出："; cin>>f; if(f==1000) break; for(int i=0;i<m;i++) { if(f>x[i]&&f<x[i+1]) { S=pow((x[i+1]-f),3)*M[i]/(6*h[i]) + pow(f-x[i],3)*M[i+1]/(6*h[i]) + (x[i+1]-f)*(y[i]-h[i]*h[i]*M[i]/6)/h[i] + (f-x[i])*(y[i+1]-h[i]*h[i]*M[i+1]/6)/h[i]; cout<<"S["<<f<<"]="<<S<<endl; } } } } scyt::~scyt() { delete []x,y,d,h,u,q,a,b,c,l,r,o,M; }
B-spline Surfaces
Zlp19970106的博客
08-21 1110
文章目录1. B样条曲面的构建方法2. 性质3. De Boor's Algorithm 1. B样条曲面的构建方法 2. 性质 3. De Boor’s Algorithm
数值计算 --- B样条函数（B-spline）
热门推荐
松下J27录放机
09-27 1万+
这篇文章给出了1至3阶B样条曲线的定义，并配有详细的图示说明。最后给出了相应的MATLAB code。
Bezier, B-Spline NURB曲线曲面
huangzhipeng的专栏
08-09 3455
地址： svn co https://curve-surface.svn.sourceforge.net/svnroot/curve-surfacecurve-surface 参考书籍：《Focus on Curves and Surfaces》 Bezier曲线 1. 令控制点个数为N (t + (1.f - t))^(N-1) 基函数系数 B-Spli...
KAN保姆级解析(1)--Bspline曲线、Bspline插值、Bspline拟合概念辨析
最新发布
罟霖同学的站点
05-17 1189
这一系列文章主要是为了解决KAN论文中对于Bspline介绍语焉不详的问题，为没有图形学、动态路径规划等领域基础知识的人提供一些BSpline的相关只知识。网络上现有的关于Bspline的介绍往往是从图形学角度出发，会略过许多的简单概念辨析，默认读者有相关图形学的基础，且以曲线生成为导向。而本文则尽量从新手的角度出发，解释KAN源码中所用到的Bspline，因此尽量简单明确，一切条件使用最基础、最理想的情况。
B-spline曲面的法线贴图
GAIA LABORATORY FOR GRAPHICS
06-15 8842
[教学视频]计算机图形学基础  在线学习教程       *原创文章，转载请注明出处*     B-spline曲面的法线贴图 Normal Mapping on B-spline Surface   在openGL Cg系列教程中，有一篇文章详细的介绍了法线贴图的原理和实现。这种法线贴图的实现是基于平面的，也就是最常看到的在四边形上的法线贴图。如果想在一些几何体上贴图，比如to...
B-Spline曲线拟合
hanmingjunv5的博客
05-15 1万+
B-Spline曲线拟合 – Python实现 1.曲线定义 定义曲线为p阶样条曲线 给定n+1个控制点 P0,P1,...Pn 节点向量 U = {u0,u1....um}，且m = n+p+1 B-Spline曲线定义如下： N(i,p)为样条基函数 P(i)为控制顶点 参数u为参数节点，一般取0 ≤ u ≤ 1 上式为样条基函数的递推公式 2.De Boor 算法 参考文章： https://pages.mtu.edu/~shene/COURSES/cs3621/NOTES/spli
冷却塔直纹面的matlab,直纹面 - calculus的日志 - 网易博客
weixin_39978101的博客
03-19 69
this.p={ m:2,b:2,loftPermalink:'',id:'fks_085074083084085075081084094095086083080064085095087069',blogTitle:'直纹面',blogAbstract:'所有图形见：',blogTag:'单叶双曲面,高等数学,直纹面,曲面,空间解析几何,微积分',blogUrl:'blog/static/251...
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值