background image

个叫有限元方法

(Finite Element Method),这种方法最开始就是用于处理实际材料力学中的

这种问题。但是

Finite Element Method 并非很简单的方式,需要解矩阵方程,运算量可不小,

自己编程也不是容易做的事情。实际上有的公司就是专门做有限元力学分析,对外出售有限
元计算包,所以可想而知

FEM 并不是容易的东西。所以大家看到别人的论文里提到 FEM,

就知道他们在做什么事情了。对于组织之间的相互应影响,侧滑,

FEM 建模等,有很多相

关论文,大家可以在

IEEE 或 google scholar 上找到一大堆,因为这样做很麻烦复杂,而发

论文往往是把简单问题复杂化,复杂问题更复杂化让大家看不懂,所以仅仅是研究生或
PHD 为了发论文而言,可以尝试做。在这一步,其实可以我们可以自己来设计计算方式,
虽然不是杨氏模量,只要是可以得到相对软硬信息,而且问题的关键是效率很高可以做到
实时处理就可以了。

第三步:错误压制和显示。第二步得到实时的软硬信息后,就可以做显示了。但还是面临问
题:错误!因为要做实际有意义的东西,前两步必须要能够实时处理,所以牺牲准确性是
必然的,到最后这步就需要压制这些错误,然后进行显示。软硬信息得到以后一般是把该信
息赋值,然后根据其值映射到伪彩显示。根据直接进行映射伪彩是根本不可取的,即使前两
步处理做得不错,直接映射

strain 结果的效果就四个字描述:惨不忍睹:).我在处理的时候,

在组织压缩或伸张量比较大的时候减少压制,在压缩或伸张量比较小的时候增大压制。这就
好像咱们要测一张纸的厚度,只测一张往往结果很不准,但如果测的是一万张纸的厚度,
然后除以一万,得到的结果就会要准确得多。这样处理以后再做一些简单的图像后处理,例
如控制弹性伪彩的总体亮度较均匀变化,以及

persistence 这样的基本处理,就可以覆盖在

原始

B-MODE 图像上进行显示了。

以上我概要的介绍了一下超声弹性成像的处理,在实际做工作时遇到的问题比我这里提到
的要多得多,我只是起个抛砖引玉的作用。希望有志于该技术国人自己商业化的朋友能够有
个大致的了解。对于有心实现可商用化产品的技术人员,我建议在

RF 上进行处理,我最开

始也是打算在

RF 上做,只是当时由于实验室的设备暂时不方便采集连续帧的 RF 数据,我

才转到在图像上做尝试。毕竟西门子和日立公司已经在

RF 上做出了商用化产品,我们就可

以确性在

RF 上做是肯定可以做出产品的。以前的 5,6 个月,我自己做了一些图像上的算法

尝试及编程验证,根据我的经验在图像上做,要同时获得实时的高效和可以接受范围内的
准确度面临的问题比较艰巨