读取CAD的DXF文件中的SPLINE 样条曲线数据 求绘制此曲线的算法与实现

[xgc2067598]

如题，大家帮帮忙啊

[buluori]

不，cad里面省略是因为我们在画spline的时候没有去编辑控制点权值的值，默认都是1的，所以生成DXF的时候cad会自动省略掉，一旦用户编辑过spline的权值，DXF里就会有权值这组数据了，所以你如果不想以后解析出来的曲线变形的话就在程序里加入进去

[buluori]

而且用的公式是不一样的 例子里用的公式是 Dj * Nj,k一般样条曲线， 而 NURBS用的是 Wj * Dj * Nj,k / Wj * Nj,k

[xgc2067598]

Dj * Nj,k一般样条曲线， 而 NURBS用的是 Wj * Dj * Nj,k / Wj * Nj,k 能解释下么 我对算法的公式不是很了解 而且公式体现在程序的那一块里

[xgc2067598]

朋友可否加QQ 或者MSN 来交流下 这里效率太慢


我的QQ 123504023
MSN [e-mail]
请注明 NURBS

[buluori]

通讯工具在公司用起来不方便，公司禁止，有事可以邮件留言

[buluori]

有没有邮箱，发段例子给你把，只要把里面的 哈特利-贾德方法 算法，改下就可以了， 加入 权值， 节点的计算用DXF读取的值替换掉

[xgc2067598]

好的 [e-mail]

[xgc2067598]

你发的那份代码 曲线逼近的算法 很不柔和，autocad里面 6个控制点 画出来的两个波段都是很柔和的

不过还是要谢谢你 不知道如何使得更加柔和 据我所致 autocad里面的样条曲线时 NURBS的算法貌似

[hn13723899236]

给我也发一份代码吧

[buluori]

对，曲线逼近那里 也要改下，他用的步进值好像是du为 0.1 的，太绝对了， 用2分法去切要好点

[buluori]

改成NURBS 要在 取点的 d布尔算法那函数里 加权值 的

CRPoint CRSpline::Dbordv(double u, int iDegree)
/*
输入参数:
k表示次数，l表示相异节点个数，u表示参数值，t[l]表示相异节点数组，r[l]表示重复度数组，
d[dim][kl]表示所涉及的控制顶点。
输出参数:
p表示B样条曲线上参数为u值的一点
*/
{
int k = Knot_Index(u) - iDegree; //相关节点的起始点索引号
if(k < 0)
return CRPoint(0,0);

CArray<double,double> dt0,dt1, dw;
dt0.SetSize(iDegree + 1);
dt1.SetSize(iDegree + 1);
dw.SetSize(iDegree + 1); //权值

for(int i = 0;i <= iDegree;i++)
{
dw[i] = m_pW[k + i];
dt0[i] = dw[i] * m_pCtrlPtList[k + i].x;
dt1[i] = dw[i] * m_pCtrlPtList[k + i].y;
}

for(int II = 1;II <= iDegree;II++)
{
for(int j = 0;j <= iDegree - II;j++)
{
double denom,alpha;
denom = Knot_Value(j+iDegree+1+k) - Knot_Value(j+II+k);

if(fabs(denom) <= 1e-5)
alpha = 0.0;
else
alpha = (u - Knot_Value(j+II+k))/denom;

dt0[j] = (1.0 - alpha) * dt0[j] + alpha * dt0[j+1];
dt1[j] = (1.0 - alpha) * dt1[j] + alpha * dt1[j+1];
dw[j] = (1.0 - alpha) * dw[j] + alpha * dw[j+1]; //权值
}
}

CRPoint p;
p.x = dt0.GetAt(0) / dw.GetAt(0); //权值
p.y = dt1.GetAt(0) / dw.GetAt(0);

return p;

}

[xgc2067598]

好的 我已经发你右键了

[xgc2067598]

非常感谢 buluori 的大力支持！ 结贴了 基本完成了

[LIUKUN1106]

麻烦你也帮我发一份代码吧 谢谢了 [e-mail]

[csuchen6]

楼主已完成？能否发份代码给我？多谢
[e-mail]

[Jorry721]

楼主能否发份代码给我？多谢
[e-mail]

[cuiyaer]

楼主能否发份代码。。谢谢！

[swmnet]

楼主可否也发发份代码。。谢谢！在下也在摸索。。。
[e-mail]

[zhaoyong1105]

小弟也在摸索，楼主能否发份代码，不胜感激
[e-mail]

[geruiguang]

能不能给我发一份啊，谢谢了，现在正在研究这个
[e-mail]