讨论建模a题！！！！！！！！

[programhunter]

现在得到半径30，如何快速找到中轴点？？？？？？？？

[SoftWare1999]

我用三段六次多项式拟合曲线，半径选择29.5
用Delphi做了个检测程序，误差不重合率：4%-9%感觉最大值大了一点。
不知道是数据有问题？？？

能否贴一下，找到的中轴点？

[starfish]

我的结果:
血管平均半径为:
29.77
每一切片上的轴线坐标为和对应的半径为：
-161.00 0.00 0.00 29.83
-161.00 0.00 1.00 29.83
-161.00 0.50 2.00 29.72
-161.00 0.50 3.00 29.72
-161.00 0.50 4.00 29.72
-161.00 1.00 5.00 29.61
-161.00 1.00 6.00 29.61
-161.00 1.50 7.00 29.72
-161.00 1.50 8.00 29.72
-161.00 2.00 9.00 29.61
-161.00 2.50 10.00 29.72
-161.00 3.00 11.00 29.61
-161.50 4.00 12.00 29.77
-161.50 4.50 13.00 29.84
-161.50 5.50 14.00 30.01
-161.50 6.00 15.00 30.10
-161.50 7.50 16.00 30.01
-161.50 8.50 17.00 30.21
-161.50 10.00 18.00 30.10
-161.50 11.50 19.00 30.01
-161.50 12.50 20.00 30.01
-161.50 13.50 21.00 30.01
-161.00 14.50 22.00 29.52
-161.00 19.00 23.00 29.07
-160.50 22.50 24.00 29.54
-160.00 24.50 25.00 29.52
-160.00 26.50 26.00 29.35
-159.50 29.50 27.00 29.54
-159.50 30.00 28.00 29.57
-159.00 33.50 29.00 29.35
-158.50 36.50 30.00 29.54
-158.00 39.00 31.00 29.43
-157.00 43.00 32.00 29.43
-156.00 47.00 33.00 29.43
-155.00 50.00 34.00 29.61
-154.00 53.50 35.00 29.72
-153.50 55.50 36.00 29.84
-152.00 59.50 37.00 29.72
-150.00 64.00 38.00 29.55
-148.00 69.00 39.00 29.83
-146.50 72.00 40.00 29.89
-144.00 77.00 41.00 29.73
-142.00 80.50 42.00 29.70
-140.00 84.00 43.00 29.83
-136.00 90.50 44.00 29.75
-134.50 92.50 45.00 29.77
-131.00 97.50 46.00 29.77
-126.50 103.50 47.00 29.84
-123.00 107.50 48.00 30.27
-119.00 112.00 49.00 30.53
-115.00 116.00 50.00 30.53
-110.50 120.00 51.00 29.70
-104.00 125.50 52.00 29.60
-95.50 132.00 53.00 29.77
-90.50 135.50 54.00 29.77
-85.00 139.00 55.00 29.83
-79.00 142.50 56.00 29.82
-72.00 146.00 57.00 30.00
-64.50 149.50 58.00 29.74
-57.50 152.50 59.00 29.74
-49.50 155.50 60.00 29.74
-43.00 157.50 61.00 29.62
-31.50 160.50 62.00 29.47
-26.00 161.50 63.00 29.60
-20.50 162.50 64.00 29.61
-10.50 164.00 65.00 29.95
-3.50 164.50 66.00 29.61
4.50 165.00 67.00 29.95
12.50 165.00 68.00 29.95
19.50 164.50 69.00 29.50
27.00 164.00 70.00 29.83
36.50 163.00 71.00 29.90
42.00 162.00 72.00 29.73
49.50 160.50 73.00 29.81
57.00 158.50 74.00 29.62
68.50 155.00 75.00 29.75
76.50 152.00 76.00 29.77
82.50 149.50 77.00 29.74
93.50 144.50 78.00 29.81
100.50 140.50 79.00 29.81
108.00 136.00 80.00 29.83
115.50 131.00 81.00 29.77
123.50 125.00 82.00 30.10
131.00 118.50 83.00 30.10
136.50 113.50 84.00 29.91
144.50 105.00 85.00 30.10
149.00 100.00 86.00 30.02
155.00 92.50 87.00 29.77
158.50 87.50 88.00 29.74
163.50 79.50 89.00 29.81
166.50 74.50 90.00 29.81
170.50 67.00 91.00 29.89
173.50 60.50 92.00 29.74
176.00 54.50 93.00 29.57
178.50 48.00 94.00 29.92
181.00 40.50 95.00 29.57
182.50 35.50 96.00 29.71
184.50 27.50 97.00 29.71
185.50 22.00 98.00 29.65
187.00 14.00 99.00 30.00

[starfish]

我们做出来的血管在三个面的投影图在
http://starfish.vip.sina.com/XY.bmp
http://starfish.vip.sina.com/YZ.bmp
http://starfish.vip.sina.com/XZ.bmp

[andrew80]

累死了，一觉睡了10个小时。呵呵。

我说一下我们的想法:
首先,可以采用穷举法,对于边界上的所有点任取两个验证是不是最大内切圆.就可以找到圆心.这种方法只是对每个切片单独考虑,会丢失很多信息.我们又采取了第二种方法,设法找到与某一切片的最大内切圆相对应的另外切片上的小圆圆心在该切片上的投影.然后对这些圆心坐标进行平均,这样就可以得到较准确的圆心坐标了.