计算几何06_B样条曲面

最新推荐文章于 2024-06-22 09:53:00 发布
最新推荐文章于 2024-06-22 09:53:00 发布
阅读量3.4k 收藏 33
点赞数 7
分类专栏： 计算几何 文章标签： 计算几何
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_35312463/article/details/112427112
版权

在这里插入图片描述

一、B样条曲面定义

在这里插入图片描述

二、双三次均匀B样条曲面片

2.1 理论

已知曲面控制顶点Pij, i=0,1,2,3; j= 0,1,2,3,即m=n=3,参数(u, v) [0, 1)x[0, 1], p=q=3。依次用线段连接点Pij阵列中的相邻两点,所形成的空间网格称为控制网格。由m=n=3构成了4x4=16个顶点的控制网格,其相应的曲面称为双三次B样条曲面。
具体曲面定义表达式，如下：
在这里插入图片描述
在这里插入图片描述

2.2 代码实现

2.2.1 定义曲面类

#pragma once
#include"P3.h"

class CBicubicUniformBSplinePatch
{
public:
	CBicubicUniformBSplinePatch(void);
	~CBicubicUniformBSplinePatch(void);
	void ReadControlPoint(CP3 P[4][4]);//读入16个控制点
	void DrawCurvedPatch(CDC* pDC);//绘制双三次均匀B样条曲面片
	void DrawControlGrid(CDC* pDC);//绘制控制网格
private:
	void LeftMultiplyMatrix(double M[][4],CP3 P[][4]);//左乘顶点矩阵
	void RightMultiplyMatrix(CP3 P[][4],double M[][4]);//右乘顶点矩阵
	void TransposeMatrix(double M[][4]);//转置矩阵
	CP2 ObliqueProjection(CP3 Point3);//斜二测投影
public:
	CP3 P[4][4];//三维控制点
};

#include "StdAfx.h"
#include "BicubicUniformBSplinePatch.h"
#include<math.h>
#define ROUND(d) int(d+0.5)

CBicubicUniformBSplinePatch::CBicubicUniformBSplinePatch(void)
{
}

CBicubicUniformBSplinePatch::~CBicubicUniformBSplinePatch(void)
{
}

void CBicubicUniformBSplinePatch::ReadControlPoint(CP3 P[4][4])
{
   for (int i=0;i<4;i++)
	   for(int j=0;j<4;j++)
   			this->P[i][j]=P[i][j];
}
void CBicubicUniformBSplinePatch::DrawCurvedPatch(CDC* pDC)
 {
	CPen NewPen, *pOldPen;
	NewPen.CreatePen(PS_SOLID, 1, RGB(255, 0, 0));
	pOldPen = pDC->SelectObject(&NewPen);
	double M[4][4];//系数矩阵Mb
	M[0][0]=-1;M[0][1]=3; M[0][2]=-3;M[0][3]=1;
	M[1][0]=3; M[1][1]=-6;M[1][2]=3; M[1][3]=0;
	M[2][0]=-3;M[2][1]=0; M[2][2]=3; M[2][3]=0;
	M[3][0]=1; M[3][1]=4; M[3][2]=1; M[3][3]=0;
	CP3 P3[4][4];//曲线计算用控制点数组
	for(int i=0;i<4;i++)
		for(int j=0;j<4;j++)
			P3[i][j]=P[i][j];
	LeftMultiplyMatrix(M,P3);//数字矩阵左乘三维点矩阵
	TransposeMatrix(M);//计算转置矩阵
	RightMultiplyMatrix(P3,M);//数字矩阵右乘三维点矩阵
	double Step=0.1;//步长
	double u0,u1,u2,u3,v0,v1,v2,v3;//u，v参数的幂
	for(double u=0;u<=1;u+=Step)
		for(double v=0;v<=1;v+=Step)
		{
			u3=u*u*u;u2=u*u;u1=u;u0=1;v3=v*v*v;v2=v*v;v1=v;v0=1;
			CP3 pt=(u3*P3[0][0]+u2*P3[1][0]+u1*P3[2][0]+u0*P3[3][0])*v3
			      +(u3*P3[0][1]+u2*P3[1][1]+u1*P3[2][1]+u0*P3[3][1])*v2
			      +(u3*P3[0][2]+u2*P3[1][2]+u1*P3[2][2]+u0*P3[3][2])*v1
			      +(u3*P3[0][3]+u2*P3[1][3]+u1*P3[2][3]+u0*P3[3][3])*v0;
			pt=pt/36;
			CP2 Point2=ObliqueProjection(pt);//斜投影      
			if(v==0)
				pDC->MoveTo(ROUND(Point2.x),ROUND(Point2.y));
			else
				pDC->LineTo(ROUND(Point2.x),ROUND(Point2.y));
		}		  
	for(double v=0;v<=1;v+=Step)
		for(double u=0;u<=1;u+=Step)
		{
			u3=u*u*u;u2=u*u;u1=u;u0=1;v3=v*v*v;v2=v*v;v1=v;v0=1;
			CP3 pt=(u3*P3[0][0]+u2*P3[1][0]+u1*P3[2][0]+u0*P3[3][0])*v3
			      +(u3*P3[0][1]+u2*P3[1][1]+u1*P3[2][1]+u0*P3[3][1])*v2
			      +(u3*P3[0][2]+u2*P3[1][2]+u1*P3[2][2]+u0*P3[3][2])*v1
			      +(u3*P3[0][3]+u2*P3[1][3]+u1*P3[2][3]+u0*P3[3][3])*v0;
			pt=pt/36;
			CP2 Point2=ObliqueProjection(pt);//斜投影    
			if(0==u)
				pDC->MoveTo(ROUND(Point2.x),ROUND(Point2.y));
			else
				pDC->LineTo(ROUND(Point2.x),ROUND(Point2.y));
		}
	pDC->SelectObject(pOldPen);
	NewPen.DeleteObject();
}

void CBicubicUniformBSplinePatch::LeftMultiplyMatrix(double M[][4],CP3 P[][4])//左乘矩阵M*P
{
	CP3 T[4][4];//临时矩阵
	for(int i=0;i<4;i++)
		for(int j=0;j<4;j++)
			T[i][j]=M[i][0]*P[0][j]+M[i][1]*P[1][j]+M[i][2]*P[2][j]+M[i][3]*P[3][j];	
	for(int i=0;i<4;i++)
		for(int j=0;j<4;j++)
			P[i][j]=T[i][j];
}

void CBicubicUniformBSplinePatch::RightMultiplyMatrix(CP3 P[][4],double M[][4])//右乘矩阵P*M
{
	CP3 T[4][4];//临时矩阵
	for(int i=0;i<4;i++)
		for(int j=0;j<4;j++)
			T[i][j]=P[i][0]*M[0][j]+P[i][1]*M[1][j]+P[i][2]*M[2][j]+P[i][3]*M[3][j];
	for(int i=0;i<4;i++)
		for(int j=0;j<4;j++)
			P[i][j]=T[i][j];
}

void CBicubicUniformBSplinePatch::TransposeMatrix(double M[][4])//转置矩阵
{
	double T[4][4];//临时矩阵
	for(int i=0;i<4;i++)
		for(int j=0;j<4;j++)
			T[j][i]=M[i][j];
	for(int i=0;i<4;i++)
		for(int j=0;j<4;j++)
			M[i][j]=T[i][j];
}

CP2 CBicubicUniformBSplinePatch::ObliqueProjection(CP3 Point3)//斜二测投影
{
	CP2 Point2;
	Point2.x=Point3.x-Point3.z*sqrt(2.0)/2.0;
	Point2.y=Point3.y-Point3.z*sqrt(2.0)/2.0;
	return Point2;
}

void CBicubicUniformBSplinePatch::DrawControlGrid(CDC *pDC)//绘制控制网格
{
	CP2 P2[4][4];//二维控制点
	for(int i=0;i<4;i++)
		for(int j=0;j<4;j++)
			P2[i][j]=ObliqueProjection(P[i][j]);
	CPen NewPen,*pOldPen;
	NewPen.CreatePen(PS_SOLID,2,RGB(0,0,255));
	pOldPen=pDC->SelectObject(&NewPen);
	for(int i=0;i<4;i++)
	{
		pDC->MoveTo(ROUND(P2[i][0].x),ROUND(P2[i][0].y));
		for(int j=1;j<4;j++)
			pDC->LineTo(ROUND(P2[i][j].x),ROUND(P2[i][j].y));
	}
	for(int j=0;j<4;j++)
	{
		pDC->MoveTo(ROUND(P2[0][j].x),ROUND(P2[0][j].y));
		for(int i=1;i<4;i++)
			pDC->LineTo(ROUND(P2[i][j].x),ROUND(P2[i][j].y));
	}
	pDC->SelectObject(pOldPen);
	NewPen.DeleteObject();
}

2.2.2 读取16个控制点

CP3 P[4][4];//控制点

void CGeometricfiguretestView::ReadPoint()//读入控制多边形顶点
{
	P[0][0].x=20;  P[0][0].y=0;  P[0][0].z=200;
	P[0][1].x=0;   P[0][1].y=100;P[0][1].z=150;
	P[0][2].x=-130;P[0][2].y=100;P[0][2].z=50;
	P[0][3].x=-250;P[0][3].y=50; P[0][3].z=0;
	P[1][0].x=100; P[1][0].y=100;P[1][0].z=150;
	P[1][1].x= 30; P[1][1].y=100;P[1][1].z=100;
	P[1][2].x=-40; P[1][2].y=100;P[1][2].z=50;
	P[1][3].x=-110;P[1][3].y=100;P[1][3].z=0;
	P[2][0].x=280; P[2][0].y=90; P[2][0].z=140;
	P[2][1].x=110; P[2][1].y=120;P[2][1].z=80;
	P[2][2].x=0;   P[2][2].y=130;P[2][2].z=30;
	P[2][3].x=-100;P[2][3].y=150;P[2][3].z=-50;
	P[3][0].x=350; P[3][0].y=30; P[3][0].z=150;
	P[3][1].x=200; P[3][1].y=150;P[3][1].z=50;
	P[3][2].x=50;  P[3][2].y=200;P[3][2].z=0;
	P[3][3].x=0;   P[3][3].y=100;P[3][3].z =-70;
}

2.2.3 调用绘制

void CGeometricfiguretestView::OnDraw(CDC* pDC)
{
	CTestDoc* pDoc = GetDocument();
	ASSERT_VALID(pDoc);
	if (!pDoc)
		return;
	// TODO: add draw code for native data here
	CRect rect;//定义矩形
	GetClientRect(&rect);//获得客户区的大小
	pDC->SetMapMode(MM_ANISOTROPIC);//pDC自定义坐标系
	pDC->SetWindowExt(rect.Width(),rect.Height());//设置窗口范围
	pDC->SetViewportExt(rect.Width(),-rect.Height());//设置视区范围,x轴水平向右，y轴垂直向上
	pDC->SetViewportOrg(rect.Width()/2,rect.Height()/2);//客户区中心为原点
	rect.OffsetRect(-rect.Width()/2,-rect.Height()/2);
	DrawGraph(pDC);//绘制图形
}

void CGeometricfiguretestView::DrawGraph(CDC* pDC)//绘制图形
{
	CBicubicUniformBSplinePatch patch;
	patch.ReadControlPoint(P);
	patch.DrawCurvedPatch(pDC);
	patch.DrawControlGrid(pDC);
}

编译，运行。
在这里插入图片描述

三、非均匀B样条曲面

3.1 理论

在这里插入图片描述

在这里插入图片描述

3.2 代码实现

3.2.1 定义曲面类

#pragma once
#include"P3.h"

class CBSplineSurface
{
public:
	CBSplineSurface(void);
	~CBSplineSurface(void);
	void Initialize(CP3** ptr,int m,int p,int n,int q);//初始化
	void GetKnotVector(double* T,int nCount,int num,int order, BOOL bU);//获取节点矢量
	void DrawCurvedSurface(CDC* pDC);//绘制曲面
	double BasisFunctionValue(double u,int i,int k,double* T);//根据节点矢量T，计算基函数
	void DrawControlGrid(CDC* pDC);//绘制控制网格
	CP2 ObliqueProjection(CP3 Point3);//斜二测投影
public:
	int m,p;//u向的顶点数减1与次数    
	int n,q;//v向的顶点数减1与次数
	double **V;//u向节点矢量数组
	double **U;//v向节点矢量数组
	CP3 **P;//三维控制点
};

#include "StdAfx.h"
#include "BSPlineSurface.h"
#include<math.h>
#define ROUND(d) int(d+0.5)
#define PI 3.1415926

CBSplineSurface::CBSplineSurface(void)
{
	P = NULL;
	U = NULL;
	V = NULL;
}

CBSplineSurface::~CBSplineSurface(void)
{
	if(NULL != P)
	{
		for(int i=0;i<n+1;i++)
		{
			delete []P[i];
			P[i]=NULL;
		}
		delete []P;
		P=NULL;
	}
	if(NULL!=U)
	{
		for(int i=0;i<n+1;i++)
		{
			delete []U[i];
			U[i]=NULL;
		}
		delete []U;
		U=NULL;
	}
	if(NULL!=V)
	{
		for(int i=0;i<m+1;i++)
		{
			delete []V[i];
			V[i]=NULL;
		}
		delete []V;
		V=NULL;
	}
}

void CBSplineSurface::Initialize(CP3** ptr,int n,int q,int m,int p)//初始化
{
	P=new CP3*[n+1];//建立动态二维数组
	for(int i=0;i<n+1;i++)
		P[i]=new CP3[m+1];
	for(int i=0;i<n+1;i++)//二维数组初始化
		for(int j=0;j<m+1;j++)
			P[i][j]=ptr[i][j];
	this->m=m,this->p=p;
	this->n=n,this->q=q;
	U=new double*[n+1];//建立u向节点矢量动态数组
	for(int i=0;i<n+1;i++)
		U[i]=new double[m+p+2];	
	V=new double*[m+1];//建立v向节点矢量动态数组
	for(int i=0;i<m+1;i++)
		V[i]=new double[n+q+2];
}

void CBSplineSurface::DrawCurvedSurface(CDC* pDC)//绘制曲面
{
	for(int i=0;i<m+1;i++)
		GetKnotVector(V[i],i,n,q, FALSE);//获取V
	for(int i=0;i<n+1;i++)
		GetKnotVector(U[i],i,m,p, TRUE);//获取U
	CPen NewPen(PS_SOLID,1,RGB(255,0,0));//曲线颜色
	CPen* pOldPen=pDC->SelectObject(&NewPen);
	double Step=0.1;//步长
	for(double u=0.0;u<=1.0;u+=Step)//绘制v向曲线
	{
		for(double v=0.0;v<=1.0;v+=Step)
		{
			CP3 pt(0,0,0);
			for(int i=0;i<n+1;i++)
			{
				for(int j=0;j<m+1;j++)
				{
					double BValueU=BasisFunctionValue(u,j,p,U[i]);
					double BValueV=BasisFunctionValue(v,i,q,V[j]);
					pt=pt+P[i][j]*BValueU*BValueV;
				}
			}
			CP2 Point2=ObliqueProjection(pt);//斜投影
			if (v==0.0)
				pDC->MoveTo(ROUND(Point2.x),ROUND(Point2.y));
			else
				pDC->LineTo(ROUND(Point2.x),ROUND(Point2.y));
		}
	}
	for(double v=0.0;v<=1.0;v+=Step)//绘制u向曲线
	{
		for(double u=0.0;u<=1.0;u+=Step)
		{
			CP3 pt(0,0,0);
			for(int i=0;i<n+1;i++)
			{		
				for(int j=0;j<m+1;j++)
				{		
					double BValueU=BasisFunctionValue(u,j,p,U[i]);
					double BValueV=BasisFunctionValue(v,i,q,V[j]);
					pt=pt+P[i][j]*BValueU*BValueV;
				}
			}
			CP2 Point2=ObliqueProjection(pt);//斜投影
			if (u==0.0)
				pDC->MoveTo(ROUND(Point2.x),ROUND(Point2.y));
			else
				pDC->LineTo(ROUND(Point2.x),ROUND(Point2.y));
		}
	}
	pDC->SelectObject(pOldPen);
	NewPen.DeleteObject();
}

void CBSplineSurface::DrawControlGrid(CDC* pDC)//绘制控制网格
{
	CP2** P2=new CP2*[n+1];
	for(int i=0;i<n+1;i++)
		P2[i]=new CP2[m+1];
	for(int i=0;i<n+1;i++)
		for(int j=0;j<m+1;j++)
			P2[i][j]=ObliqueProjection(P[i][j]);
	CPen NewPen,*pOldPen;
	NewPen.CreatePen(PS_SOLID,2,RGB(0,0,255));
	pOldPen=pDC->SelectObject(&NewPen);
	for(int i=0;i<n+1;i++)
	{
		pDC->MoveTo(ROUND(P2[i][0].x),ROUND(P2[i][0].y));
		for(int j=1;j<m+1;j++)
			pDC->LineTo(ROUND(P2[i][j].x),ROUND(P2[i][j].y));
	}
	for(int j=0;j<m+1;j++)
	{
		pDC->MoveTo(ROUND(P2[0][j].x),ROUND(P2[0][j].y));
		for(int i=1;i<n+1;i++)
			pDC->LineTo(ROUND(P2[i][j].x),ROUND(P2[i][j].y));
	}
	pDC->SelectObject(pOldPen);
	NewPen.DeleteObject();
	if(NULL!=P2)
	{
		for(int i=0;i<n+1;i++)
		{
			delete []P2[i];
			P2[i] = NULL;
		}
		delete []P2;
		P2 = NULL;
	}
}

void CBSplineSurface::GetKnotVector(double* T,int nCount,int num,int order, BOOL bU)//Hartley-Judd算法获取节点矢量数组
{
	for(int i=0;i<=order;i++) //小于曲线次数order的节点值为0
		T[i]=0.0;
	for(int i=num+1;i<=num+order+1;i++)//大于控制点个数num的节点值为1
		T[i]=1.0;
	//计算num-order个内节点
	for(int i=order+1;i<=num;i++)  
	{
		double sum=0.0;
		for(int j=order+1;j<=i;j++)
		{ 
			double numerator=0.0;
			for(int loop=j-order;loop<=j-1;loop++)
			{
				if(bU)
					numerator+=(P[nCount][loop].x-P[nCount][loop-1].x)*(P[nCount][loop].x-P[nCount][loop-1].x)+(P[nCount][loop].y-P[nCount][loop-1].y)*(P[nCount][loop].y-P[nCount][loop-1].y);
				else
					numerator+=(P[loop][nCount].x-P[loop-1][nCount].x)*(P[loop][nCount].x-P[loop-1][nCount].x)+(P[loop][nCount].y-P[loop-1][nCount].y)*(P[loop][nCount].y-P[loop-1][nCount].y);
			}
			double denominator=0.0;//计算分母
			for(int loop1=order+1;loop1<=num+1;loop1++)
			{
				for(int loop2=loop1-order;loop2<=loop1-1;loop2++)
				{
					if(bU)
						denominator+=(P[nCount][loop2].x-P[nCount][loop2-1].x)*(P[nCount][loop2].x-P[nCount][loop2-1].x)+(P[nCount][loop2].y-P[nCount][loop2-1].y)*(P[nCount][loop2].y-P[nCount][loop2-1].y);
					else
						denominator+=(P[loop2][nCount].x-P[loop2-1][nCount].x)*(P[loop2][nCount].x-P[loop2-1][nCount].x)+(P[loop2][nCount].y-P[loop2-1][nCount].y)*(P[loop2][nCount].y-P[loop2-1][nCount].y);
				}
			} 
			sum+=numerator/denominator;			
		}
		T[i]=sum;
	}
}

double CBSplineSurface::BasisFunctionValue(double t,int i,int order,double* T)//计算B样条基函数
{
	double value1,value2,value;
	if(order==0)
	{
		if(t>=T[i]&&t<T[i+1])
			return 1.0;
		else
			return 0.0;
	}
	if(order>0)
	{
		if(t<T[i]||t>T[i+order+1])
			return 0.0;
		else
		{
			double coffcient1,coffcient2;//凸组合系数1，凸组合系数2
			double denominator=0.0;//分母
			denominator=T[i+order]-T[i];//递推公式第一项分母
			if(denominator==0.0)//约定0/0
				coffcient1=0.0;
			else
				coffcient1=(t-T[i])/denominator;
			denominator=T[i+order+1]-T[i+1];//递推公式第二项分母
			if(0.0==denominator)//约定0/0
				coffcient2=0.0;
			else
				coffcient2=(T[i+order+1]-t)/denominator;
			value1=coffcient1*BasisFunctionValue(t,i,order-1,T);//递推公式第一项的值
			value2=coffcient2*BasisFunctionValue(t,i+1,order-1,T);//递推公式第二项的值
			value=value1+value2;//基函数的值
		}
	}
	return value;
}


CP2 CBSplineSurface::ObliqueProjection(CP3 Point3)//斜二测投影
{
	CP2 Point2;
	Point2.x=Point3.x-Point3.z*sqrt(2.0)/2.0;
	Point2.y=Point3.y-Point3.z*sqrt(2.0)/2.0;
	return Point2;
}

3.2.2 初始化数据

CP3** P3;//二维控制网格
int m,n;//u向和v向的控制点数-1
int p,q;//u向和v向的次数
CBSplineSurface BSurface;//B样条曲面对象	

CGeometricfiguretestView::CGeometricfiguretestView()
{
	// TODO: add construction code here
	m=4, p=3;//u向顶点数5个，曲线次数为3
	n=3, q=2;//v向顶点数4个，曲线次数为2
	P3=new CP3*[n+1];//建立动态二维数组
	for(int i=0;i<n+1;i++)
		P3[i]=new CP3[m+1];
	P3[0][0]=CP3(20,0,200),   P3[0][1]=CP3(100,100,150),P3[0][2]=CP3(230,130,140),   P3[0][3]=CP3(350,100,150),P3[0][4]=CP3(450,30,150);
 	P3[1][0]=CP3(0,100,150),  P3[1][1]=CP3(30,100,120), P3[1][2]=CP3(110,120,80),    P3[1][3]=CP3(200,150,50), P3[1][4]=CP3(300,150,50);
 	P3[2][0]=CP3(-130,100,50),P3[2][1]=CP3(-40,120,50), P3[2][2]=CP3(30,150,30),     P3[2][3]=CP3(50,150,0),   P3[2][4]=CP3(150,200,0);
 	P3[3][0]=CP3(-250,50,0),  P3[3][1]=CP3(-150,100,0), P3[3][2]=CP3(-100,150,-50),  P3[3][3]=CP3(0,150,-70),  P3[3][4]=CP3(80,220,-70);
	BSurface.Initialize(P3,n,q,m,p);
}

3.2.3 调用绘制

void CGeometricfiguretestView::OnDraw(CDC* pDC)
{
	CTestDoc* pDoc = GetDocument();
	ASSERT_VALID(pDoc);
	if (!pDoc)
		return;
	// TODO: add draw code for native data here
	CRect rect;//定义矩形
	GetClientRect(&rect);//获得客户区的大小
	pDC->SetMapMode(MM_ANISOTROPIC);//pDC自定义坐标系
	pDC->SetWindowExt(rect.Width(),rect.Height());//设置窗口范围
	pDC->SetViewportExt(rect.Width(),-rect.Height());//设置视区范围,x轴水平向右，y轴垂直向上
	pDC->SetViewportOrg(rect.Width()/2,rect.Height()/2);//客户区中心为原点
	rect.OffsetRect(-rect.Width()/2,-rect.Height()/2);
	DrawGraph(pDC);//绘制图形
}

void CGeometricfiguretestView::DrawGraph(CDC* pDC)//绘制B样条曲面
{
	BSurface.DrawCurvedSurface(pDC);
	BSurface.DrawControlGrid(pDC);
}

编译，运行。

在这里插入图片描述

沉默的舞台剧
关注
  • 7
    点赞
  • 踩
  • 33
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
专栏目录
java nurbs几何库_NURBS曲线与曲面
weixin_29821699的博客
02-28 686
B样条方法在表示与设计自由型曲线曲面形状时显示了强大的威力，然而在表示与设计初等曲线曲面时时却遇到了麻烦。因为B样条曲线包括其特例的Bezier曲线都不能精确表示出抛物线外的二次曲线，B样条曲面包括其特例的Bezier曲面都不能精确表示出抛物面外的二次曲面，而只能给出近似表示。提出NURBS方法，即非均匀有理B样条方法主要是为了找到与描述自由型曲线曲面的B样条方法既相统一、又能精确表示二次曲线弧与...
matlab绘制B样条曲面代码
05-02
在MATLAB中，B样条（B-Spline）是一种非常...总之，MATLAB中的B样条曲面绘制涉及数学理论、数值计算和图形可视化等多个方面，通过理解和应用这些知识，可以创建出复杂而精确的几何模型，对于学习和研究具有很高的价值。
计算几何05_B样条曲线
qq_35312463的博客
01-16 4108
首先我们应该来了解一下为什么要有B样条曲线，主要原因是因为Bezier曲线的不足： 1.确定了控制多边形的顶点个数，也就确定了曲线的次数； 2.控制多边形与曲线的逼近程度较差，次数越高，逼进程度越差； 3.曲线不能局部修改，调整某一控制点将影响到整条曲线，原因是伯恩斯坦基函数在整个[0,1]区间内有支撑，所以曲线在区间内任何一点的值都将受到全部顶点的影响，调整任何控制点的位置，将会引起整条曲线的改变； 4.Bezier曲线的拼接比较复杂。 为了保留Bezier方法的优点，仍采用控制顶点定义曲线。改用B.
B样条曲线曲面--拟合技术
最新发布
x13262608581的博客
06-22 1015
几何
计算几何07_NURBS曲线与曲面
qq_35312463的博客
01-23 9977
B样条方法在表示与设计自由曲线曲面时展现出强大的威力，但在设计与表示初等曲线曲面时却遇到了麻烦。因为B样条曲线、及其特例的Bezier曲线都不能精确表示除抛物线以外的二次曲线弧，B样条曲面、及其特例的Bezier曲面都不能精确表示除抛物面以外的二次曲面，而只能给出近似表示。解决这一问题的途径是改进现有的B样条方法，在保留描述自由曲线曲面长处的同时，扩充其表示二次曲线弧与二次曲面的能力，这种方法就是有理B样条方法。 一、NURBS曲线的定义及几何性质 二、权因子对NURBS曲线形状的影响 三、NURBS曲线的
【XJTUSE计算机图形学】第三章 几何造型技术(3)——B样条曲线与曲面
weixin_47692652的博客
02-14 3627
【XJTUSE计算机图形学】第三章 几何造型技术(3)——B样条曲线与曲面
MFC绘制任意B样条曲面
ha0ha0xuexi的博客
06-19 950
MFC绘制任意B样条曲面 黑色为控制点形成的控制点网，红色为绘制的3*2次B样条曲面 参考《计算几何算法与实现》–孔令德 给定，u、v方向的控制点数与坐标及次数，根据Hartely-Judd法计算控制节点，根据de-Boor递推公式绘制曲面；完整的绘制类代码如下 #pragma once #include"P3.h" class CBSplineSurface { public: CBSplineSurface(void); ~CBSplineSurface(void); void Initia
B样条曲线与曲面相关知识点汇总
Qlz的博客
02-18 1550
本文整理自西安交通大学软件学院祝继华老师的计算机图形学课件，请勿转载 文章目录B样条曲线与曲面de Boor-Cox递推定义两个特殊观察性质节点数目计算局部性开曲线定义域变差缩减性造型灵活性 B样条曲线与曲面 de Boor-Cox递推定义 两个特殊观察 基函数Ni,k(t)N_{i,k}(t)Ni,k​(t)在[ti,ti+k)[t_i,t_{i+k})[ti​,ti+k​)上非零； 在任何一个节点区间[ti,ti+1)[t_i,t_{i+1})[ti​,ti+1​), 最多有k个(k-1)次基函.
【速览版】B样条曲线曲面定义及计算例题
Lindy_pass的博客
10-20 758
其中，称为控制点，是一组预先指定的坐标称为B样条的次数，是一个预先指定的常数称为B样条基函数，是关于u的分段函数，需要计算得到B样条基函数通过如下方式定义：先定义一个非递减（前数小于等于后数）的序列称为节点向量，其满足（即 控制点总数+B样条次数）定义了节点向量后，B样条基函数。
B样条曲线 及 B样条曲面
MFDY
08-01 5095
B样条曲线 在保留Bézie优点的同时，克服了他的弱点：采用样条——分段连续多项式 例如，现在有n+1个点，在每两点之间构造一条多项式，则有n个小区间。每个小区间构造一条三次多项式，变成了n段的三次多项式拼接在一起，段与段之间要C²连续，这就是三次样条 定义  P(t)=∑i=0nPiNi,k(t)P(t)=\sum_{i=0}^{n}P_{i}N_{i,k}(t)P(t)=∑i=0n​Pi​Ni,k​(t) 其中，Pi是控制点，Ni,k(t)N_{i,k}(t)Ni,k​(t)是调和函数，被称为k阶(k-
B样条曲线曲面拟合
younggung的博客
05-26 1901
B样条曲线曲面拟合(更新中...)
Bezier曲面和B样条曲面_B样条曲面_bezier曲面_
09-29
B样条曲线的一个关键特性是其分段性质，即相邻的曲线段可以在端点处平滑连接，这使得B样条曲面能够更好地适应复杂的几何形状。在OpenGL中，可以使用GLU库（OpenGL Utility Library）中的gluBuild2DMesh()或...
matlab绘制B样条曲面代码.rar_Bspline曲面MATLAB_B样条曲线matlab_B样条曲线曲面_am9i2_b样
07-14
在`matlab绘制B样条曲面代码`这个压缩包中，你将找到示例代码，这些代码将演示如何定义控制点、选择合适的参数化变量，并使用MATLAB的内置函数绘制B样条曲线和曲面。 总的来说，掌握B样条曲线和曲面的绘制技巧对...
Uniform-Rational-B-spline-fitting.rar_B曲面_B样条曲面_b spline_曲面插值mat
07-14
而在工程计算中，如结构分析、流体力学模拟等，B样条曲面可以用来近似复杂几何体的边界，帮助解决实际问题。 总的来说，这个压缩包提供的工具可以帮助用户进行数据驱动的三维建模，通过MATLAB的高效计算能力，实现...
B样条曲线曲面拟合精选代码汇总
06-19
- **Surface construction**：B样条曲面是由两组相互正交的B样条曲线通过逐行和逐列插值形成的。 在“B样条曲线曲面拟合汇总”这个压缩包中，可能包含了一些示例代码，这些代码可能涉及到上述的各种技术，如控制点...
计算几何02_三次样条曲线
qq_35312463的博客
01-07 8955
一、样条 样条（Spline）函数是由舍恩伯格于1946年提出的。样条是富有弹性的细木条或有机玻璃条，它的作用相当于“万能”曲线板。早期船舶、汽车、飞机放样时用铅压铁压住样条，使其通过一系列型值点，调整压铁达到设计要求后绘制其曲线，称为样条曲线。这样设计曲线的方法在20世纪六七十年代得到了广泛应用。 二、几何连续性 2.1 连续性条件 通常单一的曲线段或曲面片难以表达复杂的形状，必须将一些曲线段拼接成组合曲线，或将一些曲面片拼接成组合曲面，才能表达复杂的形状。为了保证在结合点处光滑过渡，需要满足连续性条件
计算几何03_三次参数样条曲线与Cardinal曲线
qq_35312463的博客
01-09 4457
一、三次参数样条曲线 三次样条曲线的唯一缺点就是缺乏几何不变形。即当型值点发生几何变换时不能保证参数递增。因此提出了以弦长为参数的三次参数样条曲线。 1.1 定义 已知n个型值点Pi(xi, yi), i = 1, 2,…, n且相邻型值点不重合；若 p(t) 满足下列条件： （1）型值点Pi在函数 p(t) 上。 （2）p(t) 在整个区间[P1, Pn]上二阶连续可导。 （3）在每个子区间[Pi, Pi+1], i = 1, 2,…, n-1上，分段函数 p(t) 都是参数t的三次多项式。则称函.
计算几何04_Bezier曲线与曲面
qq_35312463的博客
01-10 4099
一、Bezier曲线 1.1 Bezier曲线定义 给定n+1个控制点Pi（i＝0，1，2，…，n），则n次Bezier曲线定义为： 式中，t∈[0,1]，Pi（i＝0，1，2，…，n）是控制多边形的n+1个控制点。 Bi,n(t) 是Bernstein多项式，称为Bernstein基函数，其表达式为： 式中（i＝0，1，2，…，n）。 从式中可以看出：Bezier曲线是控制多边形的控制点关于Bernstein基函数的加权和。Bezier曲线的次数为n，则需要n+1个控制点来定义。在实际应用中，二次B.
代码实现由小四边形顶点计算B样条曲面得到3维几何造型
06-01
以下是一个使用Jzy3d库实现由小四边形顶点计算B样条曲面得到3维几何造型的示例代码： ```java import org.jzy3d.chart.Chart; import org.jzy3d.chart.controllers.mouse.camera.AWTMouseUtilities; import org.jzy3d.chart.factories.AWTChartComponentFactory; import org.jzy3d.chart.factories.ChartComponentFactory; import org.jzy3d.colors.Color; import org.jzy3d.maths.Coord3d; import org.jzy3d.plot3d.builder.Builder; import org.jzy3d.plot3d.builder.concrete.OrthonormalGrid; import org.jzy3d.plot3d.builder.delaunay.DelaunayTessellator; import org.jzy3d.plot3d.builder.delaunay.jdt.JdtTessellator; import org.jzy3d.plot3d.builder.mapper.Mapper; import org.jzy3d.plot3d.primitives.Shape; import org.jzy3d.plot3d.rendering.canvas.CanvasAWT; import org.jzy3d.plot3d.rendering.canvas.Quality; import java.awt.*; import java.util.ArrayList; import java.util.List; public class B_Spline_Surface { public static void main(String[] args) { // 定义小四边形顶点坐标和控制点权重 double[][] controlPoints = { {0, 0, 0, 1}, {1, 0, 0, 1}, {1, 1, 0, 1}, {0, 1, 0, 1}, {0.5, 0.5, 0.5, 1}, {0.5, 0.5, -0.5, 1}, }; // 定义B样条曲面的度数和节点矢量 int degree = 3; double[] knotVector = {0, 0, 0, 0, 0.5, 1, 1, 1, 1}; // 定义曲面的细分数 int steps = 20; // 计算曲面顶点坐标 double[][] vertices = calculateBSplineSurface(controlPoints, degree, knotVector, steps); // 将曲面顶点坐标转换为3维几何造型 List<Coord3d> points = new ArrayList<>(); for (int i = 0; i < vertices.length; i++) { points.add(new Coord3d(vertices[i][0], vertices[i][1], vertices[i][2])); } // 创建3维几何造型 Shape surface = Builder.buildDelaunay(points); surface.setColor(Color.RED); // 创建3维画布并添加几何造型 Chart chart = AWTChartComponentFactory.chart(Quality.Advanced, ChartComponentFactory.Toolkit.newt); chart.add(surface); // 创建AWT画布并设置尺寸 CanvasAWT canvas = (CanvasAWT) chart.getCanvas(); canvas.setPreferredSize(new Dimension(800, 600)); // 显示画布 Frame frame = new Frame("B-Spline Surface"); frame.setLayout(new BorderLayout()); frame.add(canvas, BorderLayout.CENTER); frame.pack(); frame.setVisible(true); // 添加鼠标控制器 AWTMouseUtilities.addMouseController(chart); } /** * 计算B样条曲面的顶点坐标 * * @param controlPoints 控制点的坐标和权重 * @param degree 曲面的度数 * @param knotVector 节点矢量 * @param steps 曲面的细分数 * @return 曲面顶点的坐标 */ private static double[][] calculateBSplineSurface(double[][] controlPoints, int degree, double[] knotVector, int steps) { // 创建B样条曲面的映射器 Mapper mapper = new Mapper() { @Override public double f(double x, double y) { return 0; } }; // 创建网格 OrthonormalGrid grid = new OrthonormalGrid(controlPoints[0][0], controlPoints[1][0], steps, controlPoints[0][1], controlPoints[3][1], steps); // 创建三角剖分器 DelaunayTessellator tessellator = new JdtTessellator(); // 计算B样条曲面的顶点坐标 Shape surface = Builder.buildOrthonormal(grid, mapper); surface.setColor(Color.RED); surface.setFaceDisplayed(true); surface.setWireframeDisplayed(false); surface.setWireframeColor(Color.BLACK); tessellator.build(surface); return surface.getVertices().toArray(new double[0][0]); } } ``` 在上述代码中，我们使用Jzy3d库创建了一个3维画布，并添加了一个B样条曲面几何造型。我们定义了小四边形顶点坐标和控制点权重，然后计算B样条曲面的顶点坐标，并将其转换为3维几何造型。最后，我们将几何造型添加到3维画布中，并显示到屏幕上。 需要注意的是，这里使用了Jzy3d库的DelaunayTessellator三角剖分器，用于将曲面网格化。Jzy3d库还提供了其他三角剖分器，如EarClippingTessellator和IncrementalTessellator，可以根据实际需求选择使用。
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值