计算机系统应用 ２０１ ０年第１ ９卷第６期 Ｄ Ｉ Ｃ Ｏ Ｍ医学图像显示算法改进与实现① 林红利１，２ 陈真诚１ 王伟胜２（中南大学信息物理工程学院湖南长沙４１ ００８３； ２．湖南大学计算机与通信学院湖南长沙４１ ００８２） 摘要： 为了解决医学图像处理过程中在需要频繁改变窗宽、窗位的情况下图像不能实时显示，显示过程存在 明显迟滞的问题，分析了ＤｌＣＯＭ医学图像的显示过程以及传统显示过程中可以改进的地方，提出了 使用统一ＬＵＴ（Ｌｏｏｋ Ｕｐ Ｔａｂｌｅ）的方法来提高显示效率的解决方法，并使用Ｂｒｅｓｅｎｈａｍ算法改进了 传统的ＬＵＴ的生成算法。实验结果表明该改进后的算法能提高医学影像的显示效率。 关键词： ＤｌＣＯＭ；医学图像显示；Ｂｒｅｓｅｎｈａｍ算法；查找表 Ｉｍｐｌｅｍｅｎｔａｔｉｏｎ ｏｆ ｔｈｅ Ｉｍｐｒｏｖｅｄ Ａｌｇｏｒｉｔｈｍ ｆｏｒ ＤＩＣＯＭ Ｍｅｄｉｃａｌ Ｉｍａｇｅ Ｄｉｓｐｌａｙ ＬＩＮ Ｈｏｎｇ—Ｌｉ ’ ，ＣＨＥＮ Ｚｈｅｎ—Ｃｈｅｎｇ ，ＷＡＮＧ Ｗｅｉ－Ｓｈｅｎｇ￣ （１．Ｓｃｈｏｏｌ ｏｆ ｌｎｆｏ－ｐｈ），ｓｉｃｓ ａｎｄ Ｇｅｏｍａｔｉｃｓ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｃｅｎｔｒａｌ Ｓｏｕｔｈ Ｕｎｉｖｅｒｓｉｔｙ，Ｃｈａｎｇｓｈａ ４ １ ００８３，Ｃｈｉｎａ； ２．Ｓｃｈｏｏｌ ｏｆＣｏｍｐｕｔｅｒ ａｎｄ Ｃｏｍｍｕｎｉｃａｔｉｏｎ，Ｈｕｎａｎ Ｕｎｉｖｅｒｓｉｔｙ，Ｃｈａｎｇｓｈａ ４１００８２，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｏ ｉｍｐｒｏｖｅ ｄｉｓｐｌａｙ ｓｐｅｅｄ ｏｆ ＤＩＣＯＭ ｍｅｄｉｃａｌ ｉｍａｇｅ ｗｉｔｈ ｃｏｎｔｉｎｕｏｕｓｌｙ ａｄｊｕｓｔｉｎｇ ｗｉｎｄｏｗ ｖａｌｕｅｓ，ｔｈｉｓ ｐａｐｅｒ ｄｅｓｃｒｉｂｅｓ ｔｈｅ ｔｒａｎｓｆｏｒｍａｔｉｏｎ ｏｆ ｍａｎｕｆａｃｔｕｒｅｒ ｄｅｐｅｎｄｅｎｔ ｐｉｘｅｌ ｖａｌｕｅｓ ｏｆ ＤＩＣＯＭ ｍｅｄｉｃａｌ ｉｍａｇｅ ｉｎｔｏ ｐｉｘｅｌ ｖａｌｕｅｓ ｔｈａｔ ａｒｅ ｍｅａｎｉｎｇｆｕｌ ｆｏｒ ｄｉｓｐｌａｙ，ａｎｄ ｐｒｏｐｏｓｅｓ ａ ｓｏｌｕｔｉｏｎ ｕｓｉｎｇ ＬＵＴ．Ｆｉｎａｌｌｙ，ｔｈｅ Ｂｒｅｓｅｎｈａｍ ａｌｇｏｒｉｔｈｍ ｉｓ ｐｒｏｐｏｓｅｄ ｔｏ ｉｍｐｒｏｖｅ ｔｈｅ ＬＵＴ ｇｅｎｅｒａｔｉｏｎ ａｌｇｏｒｉｔｈｍ．Ａ ｇＯｏｄ ｒｅｓｕｌｔ ｉｓ ｇｏｔ ａｎｄ ａｎ ａｆｉｆｒｍａｔｉｏｎ ｆｒｏｍ ｄｏｃｔｏｒｓ ｉｓ ｍａｄｅ ｉｎ ｐｒａｃｔｉｃｅ． Ｋｅｙｗｏｒｄｓ：ＤＩＣＯＭ；ｍｅｄｉｃａｌ ｉｍａｇｅ ｄｉｓｐｌａｙ；ｂｒｅｓｅｎｈａｍ ａｌｇｏｒｉｔｈｍ；ＬＵＴ（Ｌｏｏｋ Ｕｐ Ｔａｂｌｅ） １ 引言 实际医生操作的需要有着重要的实际意义。本文从讨 随着ｕｌ（超声成像）、ＣＴ（计算机断层成像）、ＭＲＩ 论ＤｌＣＯＭ格式的图像显示过程八手，分析了显示过 磁共振成像）、ＣＲ（计算机Ｘ线成像）、电子内窥镜、 程中的关键问题，提出了统一ＬＵＴ的解决方法，并对 Ｐ盯（正电子发射断层成像）和Ｍｌ（分子影像）等医学影 ＬＵＴ的生成算法进行了改进。 像设备不断涌现，利用计算机对医学影像设备采集到 的图像进行后处理的医学图像处理与分析技术也越来 ２ ＤＩＣＯＭ医学影像显示过程分析 越成为放射科医生和临床医生进行疾病诊断的重要辅 ＤＩＣＯＭ是美国放射学会和美国电器制造商协会 助手段。在医学图像的处理过程中医生经常需要频繁 组织指定的用于医学图像存储和传输的标准【１，２】。目前 地变换窗宽、窗位来对感兴趣的病灶部位进行细致的 几乎大部分的医学设备产生的医学图像都是采用这种 研究，但传统的医学图像的显示方法不能达到实时显 格式。作为一种特殊的图像编码格式，通用图像浏览 示图像的目的，图像的显示具有明显的迟滞现象，致 软件并不能直接显示，要显示ＤｌＣＯＭ格式的医学图 使医生很多的时间浪费在等待影像的显示中，工作效 像必须将原始图像数据经过一系列的转换才能得到可 率比较低，为此开发新的医学图像的显示算法来满足 直接在显示设备上显示的数据（称之为Ｐ—Ｖａｌｕｅｓ）。 ①基金项目：湖南省自然科学基金（３９４２０１ ０４，０７Ｊ６１ ３３） 收稿时间：２００９—０９—１ ８：收到修改稿时间：２００９－１１—２１ ２０４经验交流Ｅｘｐｅｒｉｅｎｃｅｓ Ｅｘｃｈａｎｇｅ ２０１ ０年第１ ９卷第６期 计算机系统应用 ＤＩＣＯＭ医学图像显示转换过程如图１所示，需要经 过Ｍｏｄａｌｉｔｙ ＬＵＴ、ＶＯＩ ＬＵＴ、Ｐｒｅｓｅｎｔａｔｉｏｎ ＬＵＴ三 斜率／截￣（Ｒｅｓｃａｌｅ／Ｉｎｔｅｒｃｅｐｔ）变换是一种线性 变换，使用的公式是：标准图像像素值一原始图像像 素值×斜率十截距，斜率和截距在ＤＩＣＯＭ文件中可以 读取。通过这一步转换后图像像素就从设备有关变成 设备无关了。 个转换过程，最终输出的Ｐ Ｖａｌｕｅｓ才是可以直接显示 的图像数据。 ＰＶａｌｕｅｓ ２．２ ＶＯＩ（Ｖａｌｕｅ ０ｆ Ｉｎｔｅｒｅｓｔ）ＬＵＴ转换（感兴趣区转换） 由于医学图像数据动态范围大（像素深度通常不低 于４０９６个灰度级），因此一般显示器很难提供如此高 圜Ｒｅｓｃａｌｅ／图圈 图１ ＤｌＣＯＭ图像像素转换过程 ２．１ Ｍｏｄａｌｉｔｙ ＬＵＴ转换（数据规范化转 通常不同生产厂商的设备很难保证在一种设备上 生成的图像和其他生产厂商的同类型设备上生成的图 像在度量上是一致的，为此就需要将不同设备厂家产 生的图像的原始数据转换到一个标准的度量空间， Ｍｏｄａｌｉｔｙ ＬＵＴ转换就是完成这个功能的。医疗设备的 生产厂商都会在自己的图像中采用ＤＩＣＯＭ标准规定 的格式说明如何将自己的数据转换为标准图像数据， ＤＩＣＯＭ中规定可以使用通过查找表（Ｌｏｏｋ Ｕｐ Ｔａｂｌｅ， 简称ＬＵ．ｒ）查找和通过斜率／截距（ＲｅｓｃａＩｅ／Ｉｎｔｅｒｃｅｐｔ） 转换两种方法中的一种。 查找表方法（Ｍｏｄａｌｉｔｙ ＬＵ＂１３是一种非线性变换算 法。一个查找表由很多数据项组成，每一个数据项为 相应原始数值转换后的数据值，同时还会提供～个被 转换为查找表第一个数据项的原始数据的数据值。在 进行数据转换时只要知道被转换为查找表第一项的原 始数据，比原始数据大１的数据则转换为查找表第二 项的值，比原始数据大２的数据转换后的数据则为查 找表的第三项的值，依次类推直到查找表的最后一项。 大于查找表能转换的最大的原始数据则转换为查找表 的最后一项的值，小于查找表能转换的最小的原始数 据则转换为查找表的第一项的值。用公式表示为： ｆ 印≤印一 Ｍ（Ｆｐ）：｛ 一‰砜 ＜印＜砜 ＋Ｅｃｏｕｎｔ 【‰一１ Ｆｐ 仇 ＋Ｅｃｏｕｎｔ 其中：Ｍ：规范化后的标准图像像素值，Ｆｐ：原始像 素值，Ｆｐｍｉｎ：被转换为查找表第一项的原始数据， Ｍ０：ＬＵＴ中的第一个数据项的值，Ｍｐ—ｐｍｉｎ：ＬＵＴ 中第Ｆｐ—ＦｐｍＩｎ个数据项的值，Ｅｃｏｕｎｔ：查找表的数 据项个数。 的动态范围一次显示整幅图像的全部信息细节，在图像 的处理中一般都是先选择一个操作者感兴趣的区域，然 后将该区域的图像信息映射到显示器能显示的整个数 据范围，这样就增加了该区域的图像信息的对比度。这 个过程ＤＩＣＯＭ标准中称之为感兴趣区（ｖａｌｕｅ ｏｆ ｉｎｔｅｒｅｓｔ简称ＶＯＩ）ＬＵＴ（Ｌｏｏｋ Ｕｐ Ｔａｂｌｅ）转换。临床医 生感兴趣的窗宽、窗位调节功能就是ＶＯｆ ＬＵＴ转换的 Ｇ 一种算法实现。ＶＯＩ ＬＵＴ转换可以使用设置窗宽、窗 ＩＩ ，●●●●●●●●●，【 ， ●● ●位的线性转换算法和通过查找查找表（ＬＵＴ）转换的非线 【 性算法两种算法中的一种且只能使用其中的一种，具体 翌 使用哪种算法在ＤＩＣＯＭ文件中有专门的标记来设置。 。 窗宽、窗位线性转换算法是根据预知的窗宽和窗 ２ 一位值获得需要显示的窗口大小（窗宽）和中心位置（窗 Ｃ 位），从而将窗口区域的图像数据线性地转换到显示器 的最大范围内，高于或低于窗口上、下限的图像数据 则分别设置为最高或最暗的显示值。用公式表示为： ＜Ｃ一 ／２ Ｃ一 ／２ Ｓ Ｃ＋ ／２ ＞Ｃ＋ ／２ 其中Ｖ：图像数据，Ｇ：显示器显示值，ｇｍ：显示器 的最大显示值，Ｗ：窗宽，Ｃ：窗位。在这里窗宽是指 需要显示图像的范围，调节窗宽主要影响对比度，窗 宽越大图像灰度层次多，组织对比度减少，细节显示 差，而窗位也称窗中心，表示显示区域的中心位置， 例如骨骼的窗位（Ｃ）、窗宽（Ｗ）分别为Ｃ＝４００， Ｗ＝２０００，那么就可以利用调窗处理将窗宽调节到骨 骼窗宽，窗位调节到骨骼的窗位，然后利用上面的公 式将图像Ｖ换算成显示器显示值，最终得到的结果是 只显示窗口范围内的图像，也就是骨骼。 ＶＯＩ ＬＵＴ是一种查找表算法，算法处理过程同 Ｍｏｄａｌｉｔｙ ＬＵＴ。 ２．３ Ｐｒｅｓｅｎｔａｔｉｏｎ ＬＵＴ转换 Ｐｒｅｓｅｎｔａｔｉｏｎ ＬＵＴ是对图像像素要做的最后一 Ｅｘｐｅｒｉｅｎｃｅｓ Ｅｘｃｈａｎｇｅ经验交流２０５ 计算机系统应用 ２０１ ０年第１ ９卷第６期 个变换，它用于特定图像的显示。这一模块转换完成 后的输出值为Ｐ—Ｖａｌｕｅｓ，Ｐ—Ｖａｌｕｅｓ是于任何显 以大大提高图像的调窗处理时间，达到实时调窗的性 能要求。在和医生沟通后我们发现这个想法是可行的， 因为在进行疾病诊断中医生一般会根据医学图像的类 示设备的特性曲线，与人的视觉反应近似相关的值， 可直接作为已经校正的软拷贝设备或硬拷贝的输入。 Ｐｒｅｓｅｎｔａｔｉｏｎ ＬＵＴ转换也有两种转换方法，一种是通 过Ｐｒｅｓｅｎｔａｔｉｏｎ ＬＵＴ进行转换的非线性转换方法， 一型和检查的身体部位不同去查看几种典型的窗宽、窗 位。例如在ＣＴ图像中医生频繁使用的是骨骼 （Ｃ＝４００，Ｗ＝２０００），胸腔（Ｃ＝５０，Ｗ＝３５０），肺部 （Ｃ＝－６００，Ｗ＝１ ５００），腹部（Ｃ＝４５，Ｗ＝２５０）等。而 Ｍｏｄａｌｉｔｙ ＬＵＴ考虑到图像操作时只使用一次，可以将 种是通过Ｐｒｅｓｅｎｔａｔｉｏｎ Ｓｈａｐｅ的转换方法，这两 种转换方法只能使用一种。Ｐｒｅｓｅｎｔａｔｉｏｎ ＬＵＴ转换的 过程基本同上面介绍的两种ＬＵＴ算法。Ｐｒｅｓｅｎｔａｔｉｏｎ 从原始像素规范化后的数据值保存在一个矩阵中，每 Ｓｈａｐｅ算法使用ＩＮＤＥＮＴＩＴＹ和ＬＩＮ ＯＤ两者中的一 次ＶＯＵ ＬＵＴ时只要直接使用这个矩阵中的值就可以 个。ＩＮＤＥＮＴＩＴＹ表示经过ＶＯＩ ＬＵＴ转换得来的数据 了，不需要进行算法改造了；Ｐｒｅｓｅｎｔａｔｉｏｎ ＬＵＴ的线 值即Ｐ—Ｖａｌｕｅｓ值，不需要转换，一般图像的软拷贝 性处理算法在软拷贝中只使用Ｐｒｅｓｅｎｔａｔｉｏｎ Ｓｈａｐｅ 使用这个值；ＬＩＮ ＯＤ是实现图像的硬拷贝时使用的， 中的ＩＮＤＥＮＴＩＴＹ值也不需要算法改造；所以只需要 在软件的开发中没有使用到就不讨论了Ｏ 改造ＶＯｌ ＬＵＴ的线性算法就可以了。 ． 分析Ｍｏｄａｌｉｔｙ ＬＵＴ、ＶＯＩ ＬＵＴ和Ｐｒｅｓｅｎｔａｔｉｏｎ ３算法改进 ＬＵＴ发现ＤＩＣＯＭ标准的ＬＵＴ是由多个数据项组成， 在图像后处理软件的开发过程中我们发现一幅 第一个数据项保存的是第一个被转换的原始数据在查 ＤＩＣＯＭ图像从调入，到显示需要的时间比较长，而且 找表的对应值，第二个数据项为第一个被转换的数据 医生在利用医学图像进行疾病诊断的过程中需要频繁 项＋１的原始数据项在查找表中的对应值，依次类推， 地变换图像的窗宽、窗位来对感兴趣的区域进行仔细 其数学表示方法为：Ｙ［Ｏ］＝Ｆ（ＸＯ），Ｙ［１】：Ｆ（×Ｏ＋１）， 的研究，每次调窗处理后图像的显示具有明显的迟滞 Ｙ［ｉ］＝Ｆ（ＸＯ＋ｉ），其中ｉ的取值范围为【０一ＬＵＴ数据项个 现象，达不到医生的临床诊断要求，为了能动态、实 数１，Ｘ０表示第一个被转换的原始数据。结合我们要 时地对图像进行调窗必须改进图像处理算法，缩短每 实现的算法，要将线性算法改造成ＬＵＴ算法首先必须 次图像调窗处理时间。 计算出第一个被映射的原始数据的值，被映射的原始 ３．１使用统一的ＬＵＴ算法 数据的个数（决定了要构造的ＬＵＴ中数据项的个数）， 通过分析图像显示过程，我们发现医学图像从原 然后是构造ＬＵＴ。考虑到经过ＶＯｌ ＬＵＴ转换后的灰度 始数据到Ｐ—Ｖａｌｕｅｓ要经过Ｍｏｄａｌｉｔｙ ＬＵＴ、Ｖ０Ｉ ＬＵＴ、 值最小为０，最大为２Ｓ５，灰度间隔为１，因此第一 Ｐｒｅｓｅｎｔａｔｉｏｎ ＬＵＴ三个处理步骤，而Ｍｏｄａｌｉｔｙ ＬＵＴ、 个被映射的原始数据的值为Ｘ０＝Ｃ－０．５一（Ｗ一１）／２， ＶＯＩ ＬＵＴ又有线性变换和非线性变换两种算法，其中 ＬＵＴ数据项的个数为Ｗ，其中Ｃ为窗位，Ｗ为窗宽。 非线性变换算法是采用查找表的方法，这种方法相比 ３．２ ＬＵＴ生成算法优化 线性变换算法来说的优点是变换效率高，而内存的占 最简单的ＬＵＴ生成算法是将取值范围在Ｃ—Ｗ／２ 用会大些；而线性变换算法的优点是内存占用少，但 －Ｃ＋Ｗ ／２之间的的整数（ＤＩＣＯＭ中图像图像为整数） 算法效率不高，例以５１ ２ Ｘ ５１ ２的一幅医学图像来说， 使用ＶＯＩ ＬＵＴ转换中的线性算法公式进行计算后保存 每一次ＶＯＩ ＬＵＴ变换需要进行５１ ２×５１ ２次的浮点 到ＬＵＴ中。分析发现这个线性算法需要进行大量的浮 数乘除法运算，对实时性要求比较高的医学图像处理 点数的乘法和除法运算，是比较耗时间的，如果能对 来说性能上是难以接受的，为此借鉴ＬＵＴ的非线性处 这个线性算法进行改进将能提高图像的处理效率。再 理方法考虑将线性处理算法改造成ＬＵＴ算法，同时找 看线性转换的公式： 出医生实际操作时频繁需要的调窗操作对应的窗宽、 ＜Ｃ一／ｋ＂／２ 窗位提取出来，在图像装入的时候创建其相应的ＬＵＴ， Ｃ一 ，２≤ ≤Ｃ ４－ ／２ 这样在医生进行调窗处理的时候不需要重新生成 ＞Ｃ＋ ／２ ＬＵＴ，直接从对应的ＬＵＴ表中查表就可以了，这样可 可表示为Ｇ＝ａｖ＋ｂ的形式，生成查找表的过程和 ２０６经验交流Ｅｘｐｅｒｉｅｎｃｅｓ Ｅｘｃｈａｎｇｅ ２０１ ０年第１ ９卷第６期 计算机系统应用 的值 画直线的过程是一样的，因此可以使用计算机图形学 中的直线扫描转换算法来对线性算法进行优化。直线 扫描转换算法可以将大量的浮点数乘法、除法运算转 换成加法和减法运算，提高运算效率。常用的直线扫 描转换算法有数值微分（Ｄｉｇｉｔａｌ ｄｉｆｆｅｒｅｎｔｉａｌ ａｎａｌｙ－ ｋ＝２５６／ｉ—Ｅｃｏｕｎｔ：１／２５６为显示器的能表示的灰度 值个数，即ＬＵＴ的最大值 ｅ＝一０．５： ｘ＝ｉ—ＦＶａｌｕｅ；／／第一个被转换的原始数据值 ｙ＝０：／／转换后的第～个数据值，ＬＵＴ的第一项数据值 ｆｏｒ（ｉｎｔ ｉ＝０；ｉ＜：ｉ—Ｅｃｏｕｎｔ：ｉ＋＋） ｚｅｒ，ＤＤＡ）算法Ｉ３，４１、中点画线算法【ｓ１、Ｂｒｅｓｅｎｈａｍ 算法【６，７１，其中以Ｂｒｅｓｅｎｈａｍ算法效率最高，因此考虑 使用Ｂｒｅｓｅｎｈａｍ算法来生成ＬＵＴ。 ３．３改造后算法 在系统设计时我们使用一个全局的链表来保存当前 ｛ Ｐ＝ｙ： ／／将数据项保存到向ＬＵＴ Ｐ＋＋／／指针前指 ｘ＝×＋１：／／被转换的原始数据值加１ ｅ＝ｅ＋ｋ：／／ｅ在这里表示对应的查找表数据项的值是 否需要加１，当ｅ＞０时表示查找表的数据项的值要加１ ｉｆ（ｅ＞０）／／如果误差ｅ＞０则查找表数据项值加１ ｛ｙ＋＋：ｅ＝ｅ一１：｝ ｝／／循环结束 ， 系统中设置了哪些窗宽、窗位的信息，链表中的每个结 点保存的数据类型是一个结构体，结构体的定义如下： ｔｙｐｅｄｅｆ ｓｔｒｕｃｔ ＷｉｎｄｏｗＰａｒ／／定义结构体 ｛ｉｎｔ ｉ＿ｃｏｕｎｔ；，／查找表数据项数 ｉｎｔ｜－ＦｉｒｓｔＶａｌｕｅ；／／映射到查找表第一项的原始数据的值 ｉｎｔ ＬＵＴＰ；，，指向查找表的指针 ） ３．４优化效果 通过优化后，效果比较明显，达到了医学图像实 时调窗的要求。另外算法还可以继续完善，下一步我 们计划借鉴计算机操作系统对内存的调度方法来进一 步优化ＬＵＴ的生成管理。 将Ｂｒｅｓｅｎｈａｍ算法应用于ＬＵＴ的生成的最终算法为： Ｖｏｉｄ ＣｒｅａｔｅＬｕｔＯｎｔ Ｃ。ｉｎｔ Ｗ）／／创建ＬＵＴ ｛／／Ｃ：窗位，Ｗ：窗宽 ｉｎｔ．＿Ｅｃｏｕｎｔ，．－ＦＶａｌｕｅ；／／ｉ—Ｅｃｏｕｎｔ：查找表数据项数 目，ｉ—ＦＶａｌｕｅ：映射到查找表第一项的原始数据 ｓｔｒｕｃｔ ＷｉｎｄｏｗＰａｒ ＷｉｎＰ；／／ＷｉｎＰ：查找表结构变量 ／／计算ＬＵＴ数据项的个数 ｉ—Ｅｃｏｕｎｔ＝Ｗ： ４结语 本文在分析ＤｌＣＯＭ医学图像显示过程的基础上， 使用ＬＵＴ来统一ＤｌＣＯＭ医学图像的处理过程，并对 ＬＵＴ的生成算法进行了改进，并给出了具体的代码实 现。实验证明改进效果明显，具有很好的实用价值。 参考，：献 ／／计算第一个被映射的原始数据的值 ｉ—ＦＶａｌｕｅ＝（ｉｎｔ）（Ｃ－０．５一（Ｗ一１）／２）？： ／／创建ＬＵＴ数据结构 ｉｎｔ ｐ＝ｎｅｗ ｉｎｔ［ｉ Ｅｃｏｕｎｔ］； ／／设置保存窗宽、窗位以及该窗宽、窗位下的ＬＵＴ的 首地址的信息到一个全局链表中，链表中的数组元素的类型 是一个结构体 ＷｉｎＰ．－－ｃｏｕｎｔ＝．－Ｅｃｏｕｎｔ； ＷｉｎＰ．．－ＦｉｒｓｔＶａｌｕｅ＝－＿ＦＶａｌｕｅ； ＷｉｎＰ．ＬＵＴＰ＝ｐ： ｌ田捷，包尚联，周明全．医学影像处理与分析．北京：电 子工业出版社，２００３．２０—４０． ２ ＡＣＲ／ＮＥＭＡ．Ｄｉｇｉｔａｌ Ｉｍａｇｉｎｇ ａｎｄ Ｃｏｍｍｕｎｉｃａｔｉｏｎ ｉｎ Ｍｅｄｉｃｉｎｅ，ｄｒａｆｔ Ｓｔａｎｄａｒｄ ３．０ ２００３，ｆｒｏｍ ｆＩｐ：／／邱． ｎｅｍａ．ｏｒｇ／ｍｅｄｉ２ｃａｌ／ｄｉｃｏｍ，２００８．７． ３ Ｈｅａｒｎ Ｄ，Ｂａｋｅｒ ＭＰ．Ｃｏｍｐｕｔｅｒ Ｇｒａｐｈｉｃｓ．２ｎｄ ｅｄ．， Ｂｅｒｉｎｇ：Ｔｓｉｎｇｈｕａ Ｕｎｉｖｅｒｓｉｔｙ Ｐｒｅｓｓ，１ ９９８．１ ２０—１ ２７． ４ Ｂｒｅｓｅｎｈａｍ ＪＥ．Ａｌｇｏｒｉｔｈｍ ｆｏｒ ｃｏｍｐｕｔｅｒ ｃｏｎｔｒｏｌ ｏｆ ａ ／／生成一个新的结点并加入到全局链表中，此处代码 省略： ／／循环执行Ｂｒｅｓｅｎｈａｍ算法 ／／在这里因为显示器能表示的灰度范围为０—２Ｓ５， 因此查找表的数据项的最小值为０，作为查找表的第一个数 据项，数据项的最大值为２ ５ ５为查找表的最后一个数据项。 ｆｌｏａｔ ｋ，ｅ： ｄｉｇｉｔａｌ ｐｌｏｔｔｅｒ．ＩＢＭ Ｓｙｓｔｅｍ Ｊｏｕｒｎａｌ，１９６５，４（１）：２５—３０． ５刘晶，李俊，孙涵．改进的Ｂｒｅｓｅｎｈａｍ直线生成算法．计 算机应用与软件，２００８，２５（１０）：１３３—１３７． ６贾银亮，张焕春，经亚枝．Ｂｒｅｓｅｎｈａｍ直线生成算法改 进．中国图形图像学报，２００８，１３（１）：２０—２４． ７孙岩，唐棣．并行的Ｂｒｅｓｅｎｈａｍ直线生成算法．计算机 工程与应用，２００１，３７（２１）：１３６—１４０． Ｅｘｐｅｒｉｅｎｃｅｓ Ｅｘｃｈａｎｇｅ经验交流２０７ ｉｎｔ Ｘ，ｙ；／／ｘ表示被转换的原始数据值，Ｙ表示转换后
