随笔 - 373  文章 - 1  评论 - 771  阅读 - 137万

数学图形(1.48)Cranioid curve头颅线

  这是一种形似乎头颅的曲线.这种曲线让我想起读研的时候,搞的医学图像三维可视化.那时的原始数据为脑部CT图像.而三维重建中有一种方式是面绘制,是将每一幅CT的颅骨轮廓提取出来,然后一层层地罗列在一起,生成一个3d的MESH.我的研究口味重吧.这里,我也会将一层层的头颅线转化为3D的MESH,也就是一个天灵盖的图形模型.

  头颅线的方程式为: (x2 + y2)2 - 2 y2 (y + x) + a2 (b - 1) x2 + (1 - a2) x2 = 0.

  相关软件参见:数学图形可视化工具,使用自己定义语法的脚本代码生成数学图形.该软件免费开源.QQ交流群: 367752815

 头颅线脚本代码:

复制代码
#http://www.2dcurves.com/trig/trigcr.html

vertices = 1000

t = from (0) to (2*PI)

a = 2
b = 0.5

r = sin(t) + a*sqrt(1 - b*cos(t)*cos(t))

x = r*cos(t)
y = r*sin(t)
复制代码

cranioid face A:

 

复制代码
#http://www.2dcurves.com/trig/trigcr.html

vertices = D1:360 D2:32

u = from (0) to (2*PI) D1
v = from (0) to 5 D2

b = 0.5

r = sin(u) + v*sqrt(1 - b*cos(u)*cos(u))

x = r*cos(u)
y = r*sin(u)
复制代码

 

 

cranioid face B:

复制代码
vertices = D1:360 D2:32

u = from (0) to (2*PI) D1
v = from (0) to 1 D2

a = 2

r = sin(u) + a*sqrt(1 - v*cos(u)*cos(u))

x = r*cos(u)
y = r*sin(u)
复制代码

MESH:

复制代码
vertices = D1:50 D2:50

u = from (0) to (2*PI) D1
v = from (0) to (PI/2) D2

a = 2
b = 0.5

r = sin(u) + a*sqrt(1 - b*cos(u)*cos(u))
r = r*cos(v)

x = r*cos(u)
y = r*sin(u) - 1.0*cos(v)
z = 0.8*a*sin(v)
复制代码

posted on   叶飞影  阅读(2281)  评论(1编辑  收藏  举报
编辑推荐:
· go语言实现终端里的倒计时
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
阅读排行:
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· 分享 3 个 .NET 开源的文件压缩处理库,助力快速实现文件压缩解压功能!
· Ollama——大语言模型本地部署的极速利器
· DeepSeek如何颠覆传统软件测试?测试工程师会被淘汰吗?
< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5

点击右上角即可分享
微信分享提示