血管的三重建-讲解.docVIP

A题 血管的三维重建 断面可用于了解生物组织、器官等的形态。例如，将样本染色后切成厚约1? m的切片，在显微镜下观察该横断面的组织形态结构。如果用切片机连续不断地将样本切成数十、成百的平行切片， 可依次逐片观察。根据拍照并采样得到的平行切片数字图象，运用计算机可重建组织、器官等准确的三维形态。 假设某些血管可视为一类特殊的管道，该管道的表面是由球心沿着某一曲线（称为中轴线）的球滚动包络而成。例如圆柱就是这样一种管道，其中轴线为直线，由半径固定的球滚动包络形成。 现有某管道的相继100张平行切片图象，记录了管道与切片的交。图象文件名依次为0.bmp、1.bmp、…、 99.bmp，格式均为BMP，宽、高均为512个象素（pixel）。为简化起见，假设：管道中轴线与每张切片有且只有一个交点；球半径固定；切片间距以及图象象素的尺寸均为1。 取坐标系的Z轴垂直于切片，第1张切片为平面Z=0，第100张切片为平面Z=99。Z=z切片图象中象素的坐标依它们在文件中出现的前后次序为 （-256，-256，z），（-256，-255，z），…（-256，255，z）， （-255，-256，z），（-255，-255，z），…（-255，255，z）， …… （ 255，-256，z），（ 255，-255，z），…（255，255，z）。 试计算管道的中轴线与半径，给出具体的算法，并绘制中轴线在XY、YZ、ZX平面的投影图。 Z=1 Z=1 Z=0 Z=99 Z=98 Z=49 Z=50 一、思考步骤： 将图形迭加，画出三维图形，获得感性认识； 获得每一副图片的中心坐标（细化）； 对中心坐标进行3个平面的投影，结合拟合获得空间的三维曲线方程。 细化 怎么获得每一副图片的中心坐标？ 读取数据； 抽空、获得边界； 证明每一副图片中含有一个大圆------必有中心（写定理并证明）； 寻找中心（细化） 分步算法 读取数据： 读 x= imread(filename,fmt) x是什么？ 绘画 image(x) 抽空： 如果与一点相邻的8个点颜色为黑，则把该点的颜色边为白； 这样除一些特殊点外，基本抽空 寻找中心 第一副图片近似为圆，可以求出圆心以及半径 方法1：以其中一点为圆心，以已知半径画圆，圆上的点全部在图片内则该圆心即为所求； 方法2：从一点出发找另一点，满足1、两点的距离为直径；2、两点的切线平行（怎么判断平行？） 二、推广 以上方法可以求得中心坐标，但不准确，可以对起进行细化每一点细化为4点，再按上述操作。 垃圾运输问题 某城区有３６个垃圾集中点，每天都要从垃圾处理厂（第３７号节点）出发将垃圾运回。现有一种载重 6吨的运输车。每个垃圾点需要用１０分钟的时间装车，运输车平均速度为40公里／小时（夜里运输，不考虑塞车现象）；每台车每日平均工作 4小时。运输车重载运费1.8元/吨公里；运输车和装垃圾用的铲车空载费用0.4元/公里；并且假定街道方向均平行于坐标轴。请你给出满意的运输调度方案以及计算程序。问题： 1. 运输车应如何调度（需要投入多少台运输车，每台车的调度方案，运营费用） 2. 铲车应如何调度（需要多少台铲车，每台铲车的行走路线，运营费用） 3. 如果有载重量为4吨、6吨、8吨三种运输车，又如何调度？ ? 垃圾点地理坐标数据表 ? 序号 站点 编号 垃圾量T 坐标(km) 序号 站点 编号 垃圾量T 坐标(km) x y x y 1 1 1.50 3 2 20 15 1.40 19 9 2 2 1.50 1 5 21 32 1.20 22 5 3 3 0.55 5 4 22 22 1.80 21 0 4 4 1.20 4 7 23 23 1.40 27 9 5 6 0.85 0 8 24 24 1.60 15 19 6 5 1.30 3 11 25 25 1.60 15 14 7 7 1.20 7 9 26 26 1.00 20 17 8 8 2.30 9 6 27 27 2.00 21 13 9 9 1.40 10 2 28 28 1.00 24 20 10 10 1.50 14 0 29 29 2.10 25 16 11 11 1.10 17 3 30 30 1.20 28 18 12 12 2.70 14 6 31 31 1.90 5 12 13 13 1.80 12 9 32 21 1.30 17 16 14 14 1.80 10 12 33 33 1.60 25 7 15 20 0.60 7 14 34 34 1.20 9 20 16 16 1.50 2 16 35 35 1.50 9 15 17 17 0.80 6 18 36 36 1.30 30 12 18 18