CN101617342B - 用于动态信息的图形表示的方法和系统 - Google Patents
用于动态信息的图形表示的方法和系统 Download PDFInfo
- Publication number
- CN101617342B CN101617342B CN2007800498369A CN200780049836A CN101617342B CN 101617342 B CN101617342 B CN 101617342B CN 2007800498369 A CN2007800498369 A CN 2007800498369A CN 200780049836 A CN200780049836 A CN 200780049836A CN 101617342 B CN101617342 B CN 101617342B
- Authority
- CN
- China
- Prior art keywords
- value
- interpolation
- multidate information
- absolute value
- carried out
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 85
- 238000001914 filtration Methods 0.000 claims description 22
- 239000008280 blood Substances 0.000 abstract description 5
- 210000004369 blood Anatomy 0.000 abstract description 5
- 241001465754 Metazoa Species 0.000 abstract description 2
- 239000007788 liquid Substances 0.000 abstract description 2
- 238000002604 ultrasonography Methods 0.000 description 13
- 230000000694 effects Effects 0.000 description 10
- 239000000523 sample Substances 0.000 description 10
- 230000017531 blood circulation Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 8
- 230000014509 gene expression Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 238000013507 mapping Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 239000003086 colorant Substances 0.000 description 4
- 238000012800 visualization Methods 0.000 description 4
- 210000002216 heart Anatomy 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000002203 pretreatment Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 239000003643 water by type Substances 0.000 description 2
- 241001062009 Indigofera Species 0.000 description 1
- 240000001439 Opuntia Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000004217 heart function Effects 0.000 description 1
- 210000003709 heart valve Anatomy 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 210000001644 umbilical artery Anatomy 0.000 description 1
- 210000003606 umbilical vein Anatomy 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Generation (AREA)
- Digital Computer Display Output (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种用于对来自动态2D或者3D对象的动态信息进行图形表示的方法和系统。具体地,本发明涉及一种对具有不同符号的值的动态信息进行插值的新颖方法。本发明的所述方法对于动态信息的图形表示特别有用,特别是对于用彩色多普勒方法已经获得的流动数据特别有用。本发明的所述方法可以用于表示液体(例如血液)的流动,并且因此对于人类和动物的诊断来说是有用的。
Description
技术领域
本发明涉及一种用于对来自动态2D或者3D对象的动态信息进行图形表示的方法和系统。具体地,本发明涉及一种对具有不同符号的值的动态信息进行插值的新颖方法。
背景技术
本发明的所述方法对于动态信息的图形表示特别有用,具体地,对于用彩色多普勒方法已经获得的流动数据(flow data)特别有用。本发明的所述方法可以用于表示液体(例如血液)的流动并且因此对于人类和动物诊断来说是有用的。
尽管原则上可以表示任意动态信息,但是本发明意义上的动态信息具体意味着流动数据,比如从多普勒超声波扫描测量法得到的彩色流动数据。可以使用多普勒测量法来增强超声波扫描,多普勒测量法使用了多普勒效应来估计组织(比如血液)是向探针移动还是远离探针移动,并估计它的相对速度。通过计算特定样本体的频率偏移,例如流过心脏瓣膜的血流(a jet of blood),可以确定它的速度和方向并且将该速度和方向视觉化。这在心血管研究和其它器官(例如肝等等)的研究中特别有用。
可以在连续波系统和脉冲波系统中使用多普勒超声波扫描。尽管如此,连续波多普勒超声波不能确定波束(beam)中速度的特定位置,并且因此不能用于产生彩色流动图像。另一方面,脉冲波超声波方法允许对流动地点的深度或者范围进行测量。对于彩色流动图像来说,使用脉冲波超声波来生成用于多普勒超声波扫描图的数据。
为了显示,将动态信息方便地映射为不同颜色。颜色的色调代表了速度的方向并且亮或者暗代表了速度的量级。通常地,使用两种不同或者对比的颜色(比如红和蓝)来代表两个相反的流动方向(正或者负)。一般地,红颜色代表向探针流动并且蓝色代表远离探针流动。当在两个具有不同符号的值之间执行插值运算时,这样的颜色插值可能产生速度颜色图之外的颜色,导致了误导的彩色流动信息以及它们相应的图形,比如黑或者白的边界线。
为了将表示为采样体数据(sampled volume data)的彩色流动数据(colour flow data)进行视觉化,采用了特别的措施来重建所谓的“混叠效应(aliasing effects)”,该效应发生在高负-高正流动区域之间的边界上。普通的三线性插值方法不适用于这样的边界区域,这是由于它们不充分强调边界区域。当以给定采样频率或者脉冲重复频率传送脉冲时,能够明确测量到的最大多普勒频率是脉冲重复频率的一半。如果被测量的血液速度和波束/流动角度结合起来给出了大于脉冲重复频率一半的多普勒频率值,则在多普勒信号中发生模糊(ambiguity)。该模糊被认为是混叠(aliasing)。在现有技术中,为了避免混叠,已经进行了增加脉冲重复频率的尝试。尽管如此,增加脉冲重复频率导致的问题是不能标识低速度。因此,避免混叠效应的其它方式是使用不同的插值方法。
因此,现有技术提出了最小弧插值(minimum arc interpolation)。在由Sumanaweera等人的US 2006/0094963 A1中描述了这样的插值方法和系统。该文档教导了最小弧插值的精确建模,该插值导致为了插值而实现相当复杂的符号计算。在图形硬件上很难有效率的实施这样的计算。
在US 2006/0094963 A1中披露了另一种方法,该方法将速度分为实数/复数表示,分别对它们进行插值并且稍后使用预计算的查找表(LUT)将它们结合起来。当对LUT采样时,由于有限的精度和插值误差,该被采样的查找表可以导致数值假象(numerical artefacts)。
特别是对于三角函数(arctan2),普通字节值不提供足够的精度并且将产生假象。使用较高的精度仅在最新的GPU(图像处理单元)上才可能并且将导致严重的性能下降。
现有技术方法的另一个不足之处是使用最小弧插值本身。输入值中的小变化将导致速度的激烈变化。当沿着弧的插值方向突然改变时就是这种情况。使用基于合成(composition)的体-渲染(volume-rendering)(将穿过体的许多切片从后至前或者从前至后的结合起来),在方向上的这些突然改变可以导致不能忽视的假象(“灰化效应”(grizzling-effect))。
发明内容
因此本发明的目标是提供一种用于图像渲染的方法,该方法至少在某种程度上避免了现有技术的缺点。
本目标由一种用于对来自动态二维或者三维对象的动态信息进行图形表示的方法来解决,该方法包括:
(a)从所述动态对象提供动态信息,
(b)通过下述步骤将所述动态信息编码:将所述动态信息分为绝对值和方向值,并且将所述绝对值和所述方向值映射到一个或更多字节的分别的位,
(c)分别对所述绝对值和所述方向值进行插值,以给出两个分别的值,
(d)将在步骤(c)中获得的两个分别的值结合起来,以及
(e)表示来自步骤(d)的结合的值。
本发明还涉及一种系统,所述系统包括提供动态信息的超声波系统,以及图像处理子系统,该图像处理子系统包括图形处理单元(GPU)和显示设备,其中对所述图像处理子系统编程以执行本发明的方法。本发明还涉及一种用于对动态信息进行图形表示的计算机程序产品,该产品存储在计算机可用介质或本发明的系统上,并且所述计算机程序产品包括能够执行本发明的方法的软件代码部分。
本发明可以用于将任意类型的动态信息进行视觉化,具体地,速度信息,例如流动数据,特别地是从多普勒超声波扫描得到的彩色流动数据。
本发明的一个主要应用是对彩色流动数据的表示,该彩色流动数据是从血液流动的分析中获得的。人类心脏的血液流动的精确表示对于心脏功能的检查是极其有用的。除了心脏之外,也可以使用多普勒超声波扫描以及利用本发明的插值方法将其它器官中(比如静脉和动脉、脐带等)的血液流动进行图形化表示。
尽管如此,本发明的插值方法不仅限于多普勒超声波扫描,而是可以用于可能出现上述模糊和混叠效应的任意数据集。该数据可以通过例如恰当的超声波探针来获得,或者来自离线的源。本发明特别适用于硬件加速的渲染(rendering)系统(GPU图形处理单元),但是它也适用于运行在标准CPU上的基于软件的重建和渲染。
本发明的方法产生了令人视觉愉快的图像,其图形对于用户更友善,并且可以在GPU上非常有效的实施该方法。此外,本发明具有考虑到紧凑数据表示(compact data representation)的优点,可以在存储器消耗严格的系统中使用该紧凑数据表示。
具体地,本发明避免了在现有技术中使用的复杂并且消耗存储器的最小弧插值,并且允许对绝对速度值使用简单的三线性插值作为替代。通过对绝对速度值和方向值进行分别插值,使得给出两个分别的值成为可能,之后将所述两个分别的值结合起来以给出有符号(正或者负)的速度值。
可以用多种方式来实现本发明的方法。例如,在一个特定实施例中,在单一字节中对动态信息进行编码,这样具有低存储器消耗的优点。在本实施例中,将指示速度方向的符号映射至八个位的最低有效位(通常最低有效位是最低位),并且将绝对速度值映射至剩余的较高位。因此,编码是简单的位操作,不需要查找表。
发明人惊讶地发现,将符号在最低位进行编码并且使用其它较高的七个位用于所述速度值是有利的。如果对该表示进行线性插值,例如
y=t*V1+(1-t)*V2
其中t是插值因子并且V1和V2代表两个字节,V1和V2中的最低位(方向或者符号)不具有显著的影响。由于较高的七个位代表绝对值,线性插值适用于该值。对于该表示来说(LUT),必须确定方向(正或者负)。由于线性插值对方向位进行“掺水冲淡(waters down)”,需要第二次插值。为此,使用近邻插值,这意味着在V1和V2以及插值因子t的情况中,取决于t是更接近于V1还是V2,使用或者V1或者V2。如果t小于0.5,则使用V1,如果t等于或者大于0.5,则使用V2。作为与线性插值的对比,选择或者V1或者V2,而不选择在它们之间的值。这意味着从该值,可以提取出方向位(最低位)。将两个分别插值的值结合为最后的结果。较高七个位(绝对值)来自线性插值的V1/V2,符号来自最近插值的V1/V2的最低位。
一般地,如果根据现有技术在一个字节中对颜色进行编码,当使用3D纹理映射(texture mapping)对该值进行访问时,使用简单三线性插值是不利的,这是由于在负和正方向之间的插值将导致零相交(zero crossing)。在图2(A)所示的图形表示中描述了该零相交。
本发明通过下述方式来避免该问题:将给定动态信息值的符号值编码为一个位并且其它七个位用于映射绝对速度值,并且访问该数据两次:第一次对绝对速度使用三线性插值并且第二次对方向值进行插值。取决于实施本发明的系统,可以同时或者顺序地执行分别的插值步骤。通过近邻滤波来对方向值进行插值。如上所述,由于如果V1和V2的符号不同则仅能使用V1和V2中的一个的符号,所以所述符号值必须通过近邻滤波来得到。
最后,将两个值结合为一个带符号值。为此,可以使用直接方法(例如掩码屏蔽(mask out))或者可以使用查找表。在直接位操作的情况中,两个值的结合是容易的。在GPU实施的情况中,本步骤较难。尽管如此,可以在此使用查找表方法,即具有预先计算好的所有可能组合的二维表,该表允许省却位操作。
也可以使用备选组合模式,该模式使用了预先计算的二维查找表。当对颜色进行编码使得最低位具有符号(即方向)并且较高的七个位对绝对值进行编码时,可以使用图形硬件以取得平滑的彩色过渡。优选地,对于最大负和正方向来说,abs(dir)应当映射至0xFe(==254),其可以包含一个方向的反方向。然后GPU取得插值后的颜色值的较高七个位,同时取决于最近滤波的纹理的最低位来选择符号(方向)。
优选地,使用两个纹理级以从三维彩色数据体中提取标量体(scalar volume)。这些级具有与输入相同的3D-体,但是级0使用线性滤波(三线性插值),而级1使用最近滤波(选择与给定查找位置最接近的三维像素(voxel))。
然后着色器(shader)(GPU片段程序)使用来自级0的插值后的值以通过对最低位进行掩码屏蔽来确定绝对速度(无符号)。该插值是可工作的,因为最低位(包含符号)对于插值结果没有做出多大贡献,所以可以在插值期间保持它。
使用来自级1的通过使用最近滤波获得的值(因此没有被插值),可以提取出恰当的符号(方向)值。
最后,将两个值结合为带符号的插值后的速度值,可以将该速度值写入帧缓冲器或使用彩色查找表(LUT)进行其它处理。该LUT通常将正方向映射至一种颜色并且将负方向映射至第二种颜色(通常使用红色代表正并且蓝色代表负),因此使用新速度值可以对简单一维LUT进行采样。当使用alpha混合合成来渲染体时,这是优选模式。
利用所述方法,当渲染彩色部分时进行最大彩色速度投影是可能的,这是因为当将新速度值写入帧缓冲器时需要使用最大值运算(maximum opeartion)。该运算使用如上所述的相同效应:由于最低位(包含符号)对插值结果不做出过多贡献,最大值运算符的使用将选择较高七个位中所表示的最大值。
在最终的渲染步骤中,当已经结合了全部切片时,通过使用一维查找表,对所产生的速度进行映射,以显示颜色值(例如红和蓝)。
由于符号和绝对值的结合可能是在GPU上不易实施的操作(由于GPU通常对定点或浮点数据操作,所以不存在对于位运算的直接支持),本发明还提出了另一个结合方法。将首先计算符号值替换为
output=sign?(absVal|0x1):absVal;
然后使用该值访问一维彩色映射LUT,以将速度映射到屏幕颜色,可以修改本发明的方法以使用二维彩色LUT,在该方法中使用两个纹理级的结果(绝对值和符号值)作为进入2维LUT的索引,在预处理步骤中生成该2维LUT。
在本发明的另一个优选实施例中,由两个字节而不是单一字节对动态信息进行编码。这使得存储器消耗加倍,所以可以将所述体存储为双信道体纹理,例如在Open GL或者DirectX中所谓的亮度-alpha格式。对于第一优选实施例来说,此处也执行向绝对值和符号值的简单分离。这避免了现有技术中描述的复杂的实数/复数转换。将一个字节赋予带符号的速度值并且添加第二字节,在第二字节中仅存储原始速度值的绝对值。在存储器中最好将这两个值交错存储。在基于GPU的实施的情况中,双分量无符号字节格式是优选的,比如亮度-ALPHA(Luminance-ALPHA)。也可以交错附加信道,比如紊乱(turbulence),b模式(b-mode),等等。
然后仅需要使用单一的三线性插值查找以从该数据集中提取两个插值后的值。新的第二信道给出了插值后的绝对速度v_abs(如在第一实施例中),但是由于不需要在该字节中携带符号位,维持了完整的八位精度。
基于第一字节来选择恰当的符号(原始负dir/+dir速度值)。执行简单的检查:插值后的值是具有正还是负方向,并且将该符号应用至v_abs的值,即插值后的绝对值。图5示出了该处理的示意表示。
这样,具有较高速度的样本点在方向上具有更多的影响。
本运算是非常简单的比较运算,并且可以在几乎所有GPU上实施。如果使用带符号的速度来对一维彩色映射LUT进行索引,还可以应用在第一实施例中提出的二维LUT机制。这样,二维索引将由插值后的绝对速度和插值后的带符号速度构成。
由于所有东西都由插值后的值而不是由插值和非插值的纹理像素(texel)的混合(比如由根据现有技术US 2006/0094963 A1所述的方法所生成的)所隔开,所以上述第二实施例导致了较好的图像质量。修改后的符号/值映射是平滑的,这意味着速度的较小改变决不会导致大的颜色改变。这在图2(C2)(相比于使用第一实施例的图2(C1))中示出。
根据第一和第二实施例所述的方法都使用了可以在GPU本身上实现的相当简单的编码机制。由于大多数GPU是基于统一存储器模型的,可以将着色器应用于速度转换并且执行所谓的纹理位块传输(blit)。这样,通过将工作负荷移至GPU,减少了宝贵的CPU功率的消耗。
对于将数据从一个字节表示加倍至两个字节表示的本发明的第二实施例的插值方法来说,由于仅传输原始数据,本方法还将减少从系统存储器至GPU存储器的带宽。
对于两个方法来说,将速度数据编码为包装的符号值或者两个字节的符号值是非常简单的运算,并且因此可以高性能地实施。与现有技术相比,在转换期间不需要查找表或者复杂的三角运算。可以用简单的处理器指令(比较和位变换)来实施根据本发明所述的转换。该转换的变体将在GPU上进行转换。
另外,还可以使用提出的插值方法重建来自声学表示(原始数据)的笛卡尔数据集。在该重建期间,需要针对速度的插值方案。
再一次地,还可以在此处使用上述的两个实施例。使用单一字节表示的第一实施例的主要优点是,将彩色流动数据进行重排列和/或修改,使得可以利用图形卡的三线性插值单元来处理彩色流动插值的性能最关键的部分,该三线性插值单元取得八个邻近的三维像素并且将它们结合为一个标量值。使用该三线性插值后的值,能够使用非常简单的着色器以提取流动值以及方向,或者能够将该运算结合到彩色LUT中。当已经存在例如可以使用的快速三线性查找引擎时,所提出的方法也适用于纯粹基于CPU的方法。另一个优点是,该方法具有非常低的存储器消耗,这对于低端系统或具有有限数量存储器的系统(比如移动设备或者膝上型计算机中使用的GPU等等)来说是有利的。
根据使用符号值映射的第二实施例所述的方法的优点是:插值总体上是平滑的,这是由于不存在例如流动的小改变会导致视觉化的大改变这种情况。相比于现有技术,不存在由于使用三角映射的分离导致的数值问题和由于具有高数值精度和大小的查找表导致的性能问题。此外,可以在GPU上有效地实施本发明的方法,无需附加的查找表,这是由于带符号的速度的结合是非常简单的。
根据本发明的另一个实施例,还可以使用所提出的插值方法对彩色流动数据进行滤波(低通、减少噪声...)。利用彩色流动数据(带符号的动态信息)到可以使用标准线性插值进行插值的数据格式的转换,能够对原始或者重建数据应用滤波方法。首先将原始彩色流动数据转换为符号/值格式,并取决于已经选择了一字节符号/值表示还是两字节符号/值表示,在一个或者两个信道上运行标准滤波器实现。在滤波之后,可以将结果转换回原始的彩色流动格式,或者以一字节或两字节格式的符号/值表示进行进一步的处理。
一个特殊变体适用于两字节符号/值表示:对于大多数应用来说,仅对包含绝对速度的信道进行滤波就足够了,不对符号信道进行滤波。这将导致细微不同的图像,仅对一个信道而不是两个信道进行滤波的好处是例如更少的计算时间。
利用本发明的该实施例,可以使用原始形式的彩色流动数据,应用噪声减小滤波器(例如基于卷积的高斯平滑滤波器),将数据转换回原始形式,并且通过使用最小弧技术对数据进行重建或者插值。对于彩色流动超声波来说,去掉体数据中的高频是重要的,从而为了实现良好的图像质量,适合的滤波步骤是必须的。利用所提出的符号/值(两信道)表示,可以使用所有已知的标准滤波技术。该滤波可以用于三维数据、二维数据或者一维数据(例如原始超声波束)。
在本发明所有的实施例中提出的方法特别适用于实时超声波成像。可能的实现可以包括超声波系统,该超声波系统包括超声波探针,以获取速度数据(彩色流动)。备选地,可以使用用于稍后查看而存储的所获取的数据(没有探针的离线查看系统等等)。此外,存在图像处理子系统,包括图形处理单元(GPU)和显示设备。备选地,可以采用所有处理由CPU执行的图像处理子系统。
在本例中,给出了作为三维数据集的速度数据。然后在CPU上或者在GPU上将该数据转换为两种表示中的一种。在CPU处理的情况中,将转换的三维数据集传输至GPU。为了将结果可视化,GPU使用该输入数据提取(渲染)一个或者更多切片。由于该步骤需要对任意空间位置处的输入数据重采样,因此能够使用两种插值方法中的一种来计算带符号的速度值。对于三维数据来说,GPU或者CPU取得八个相邻三维像素数据样本,并执行插值。然后使用该插值结果作为显示彩色查找表的输入,或者对该插值结果进一步处理以执行例如最大速度投影。应当注意到,尽管本发明特别适用于三维数据,它也适用于二维图像,例如在仅要求单一切片的情况下。
附图说明
现在将在下面的图和例子中更详细地描述本发明。
图1(A)至(C)示出了在三种插值方法之间的差异。图1(A)示出了普通线性插值,图1B示出了最小弧插值,并且图1(C)示出了根据本发明所述的符号/值映射。图2示出了代表了四种不同插值方法的四幅不同的图。图2(A)示出了在两个值A和B之间的普通三线性插值。图2(B)示出了最小弧插值。图2(C1)和图2(C2)示出了根据本发明的第一和第二实施例所述的两个插值图。
图3示出了使用位运算将符号和绝对值相结合的流程图。
图4示出了使用二维查找表将符号和绝对值相结合的流程图,在预处理步骤中生成该查找表。
图5示出了使用来自双信道3D数据集的插值绝对速度和插值符号速度将符号和绝对值相结合的流程图。
图6示出了使用本发明的插值方法由多普勒表示的人类心脏中血液流动区域的示意图。
具体实施方式
在图1中,三个图像(A)至(C)示出了在不同插值方法中的差异。在全部三个情况中,在A和B之间的区域是待插值的。两个值A和B是具有相反方向(例如相反符号)的高值。相反方向由不同颜色所指示,红色代表正方向以及蓝色代表负方向。阴影强度(shadingintensity)指示了绝对速度值。对于使用标准纹理滤波(texture filtering)的图1(A)所示的普通插值来说,插值导致了零相交区域。A的正值减少直到它为零并且然后再次升高至B的值。缺点是当对这样的结果进行图形化表示时,零相交线将导致图形的颜色方案之外的颜色,即白线或者黑线,并且这将使得图像对于用户来说不易理解。
对于图1(B)中所示的最小弧插值来说,当将A和B按照圆形方式排列时,插值方向是沿着从A至B的最小距离。所以A的值升高直到它达到最大程度,然后调换符号并且然后从最高的蓝值向B的值移动。
最后,图1(C)示出了本发明的符号/值映射插值。在左侧,在灰色渐变长方形中示出了两个值A和B,在该长方形中将A和B放在与符号无关的相应强度上。在右侧的长方形中,根据方向(正或者负)放置值A和B,该方向是用红色或者用蓝色示出。分别执行插值。作为与最小弧插值的对比,任一值的强度永远不会低于或者高于A或者B的实际值。由近邻滤波(nearest neighbour filtering)来确定方向(符号)改变,所以转变点永远在|AB|区域的中间。
图2示出了数学图表,该图表描述了如图1所示的不同插值方法。图2(A)和图2(B)示出了普通插值和最小弧插值,然而图2(C1)和图2(C2)示出了根据本发明所述的插值,图2(C1)示出了根据第一实施例所述的插值并且图2(C2)示出了根据第二实施例所述的插值。在图2(C1)所示的第一实施例中,由近邻滤波来确定符号,即从一个方向到另一个方向的交叉发生在绝对值A和B之间的中途。因此,交叉永远是t等于0.5。在图2(C2)中,对符号进行线性插值。在本情况中,高于值B的值A因而具有对于符号的较大的影响。因此,正和负之间的交叉发生在A的绝对值变得小于B的绝对值的点处。
图3示出了使用位运算将符号和绝对值相结合的流程图。在步骤10,获得三维速度数据,然后通过分离每一个数据点的绝对值和符号值并且在一个位上对符号值进行编码并且使用其它七个位对绝对值进行编码,将三维速度数据编码为一个字节。随后,执行如步骤12和14所示的两个运算。在步骤12,使用三线性插值对绝对值进行插值。在步骤14,使用近邻插值对符号值进行插值。获得两个值,即绝对速度和符号值,在长方形16中以两个分别的长方形示出这两个值。然后将这两个值结合起来以给出结合的插值输出值,然后使用该输出值进行合成(在步骤18)。
图4示出了具有基本相同的过程的流程图,不同之处在于使用结合的2D查找表两个单独的插值后的值结合起来以获得将绝对值和符号值。
图5示出了代表本发明的方法的第二实施例的流程图,在该实施例中将三维速度数据在两个字节上进行编码,一个字节代表符号值并且一个字节代表绝对值。在长方形20中示出了这点。在长方形22中,执行与符号值无关的三线性插值。仅在步骤26中考虑符号值。如果符号是正,那么从三线性插值获得的输出值是正值,而如果符号是负,则插值后的值具有负的符号。步骤28与图3和4中的步骤18相对应。
图6示出了由多普勒超声波扫描所生成的彩色流动图像的示意图。图6A示出了图像的示意表示,其中已经使用了现有技术的最小弧插值(比如由Sumanaweera等人的US 2006/0094063 A1中披露的)对值进行插值。不同渐变区域表示了不同的颜色。黑条代表了心脏瓣膜并且阴影线区域代表了通过该心脏瓣膜的血液流动。如图6A所示,在血流的中心处存在区域,在该区域中在大体为黄色的区域中出现了青色区域的孤岛(isolated island)。黄色和青色之间的交越是该符号的反向。因此,标记为黄色和红色的血液流动与标记为青色或者蓝色的血液流动方向相反。在图6B中示出的放大的中心区域示出了在青色孤岛的周围存在双阴影线的区域。这代表出现所谓的灰化效应(grizzling-effect)的区域。这意味着由于插值误差,存在彼此相邻的不同颜色的非常小的细长区域,对于裸眼来说该细长区域看起来像是线(例如黄色和青色)。这意味着对于流动的方向来说存在模糊。当使用根据本发明所述的插值方法时,现在在图6C中示出的相同放大的区域中没有灰化效应。相反,很好地限定出黄色和青色区域。
Claims (14)
1.一种用于对来自动态二维或三维对象的动态信息进行图形表示的方法,该方法包括:
(a)从所述动态对象提供动态信息,
(b)通过下述步骤对所述动态信息进行编码:将所述动态信息分为绝对值和方向值,并且将所述绝对值和所述方向值映射到一个或更多字节的分别的位,
(c)分别对所述绝对值和所述方向值进行插值,以给出两个分别的值,
(d)将在步骤(c)中获得的所述两个分别的值相结合,
(e)对来自步骤(d)的相结合的值进行表示。
2.根据权利要求1所述的方法,其中,在步骤(d)中使用直接方法或者查找表将所述两个分别的值相结合。
3.根据权利要求1或者2所述的方法,其中,在步骤(c)中针对所述绝对值使用三线性插值。
4.根据权利要求1或2所述的方法,其中,在步骤(c)中针对所述方向值的插值使用近邻滤波。
5.根据权利要求1或2所述的方法,其中,在单一字节中对所述动态信息进行编码,使用一个位来映射所述方向值。
6.根据权利要求5所述的方法,其中,将所述方向值映射到最低有效位。
7.根据权利要求1或2所述的方法,其中,将所述绝对值与所述方向值映射到两个字节,针对每一个值使用一个字节。
8.根据权利要求2所述的方法,其中,在步骤(d)中使用的查找表将正方向值映射为第一颜色并将负方向值映射为第二颜色。
9.根据权利要求1或2所述的方法,其中,所述动态信息源自三维对象。
10.根据权利要求1或2所述的方法,其中,所述动态信息是速度信息。
11.根据权利要求1或2所述的方法,其中,使用图形处理单元(GPU)来执行所述方法。
12.根据权利要求1或2所述的方法,其中,把来自步骤(d)的相结合的插值后的值表示为极坐标图,或表示为二维或三维图像。
13.根据权利要求1或2所述的方法,其中,通过将所述动态信息转化为能够使用标准线性插值进行插值的数据格式并向原始或者重建数据应用滤波方法,对所述动态信息进行滤波。
14.一种用于对来自动态二维或三维对象的动态信息进行图形表示的设备,该设备包括:
用于从所述动态对象提供动态信息的装置;
用于通过将所述动态信息分为绝对值和方向值、并且将所述绝对值和所述方向值映射到一个或更多字节的分别的位而对所述动态信息进行编码的装置;
用于分别对所述绝对值和所述方向值进行插值以给出两个分别的值的装置;
用于将获得的所述两个分别的值相结合的装置;以及
用于对相结合的值进行表示的装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/EP2007/000336 WO2008086820A1 (en) | 2007-01-16 | 2007-01-16 | A method and a system for graphic representation of dynamic information |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101617342A CN101617342A (zh) | 2009-12-30 |
CN101617342B true CN101617342B (zh) | 2012-06-13 |
Family
ID=37909752
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007800498369A Expired - Fee Related CN101617342B (zh) | 2007-01-16 | 2007-01-16 | 用于动态信息的图形表示的方法和系统 |
Country Status (7)
Country | Link |
---|---|
US (1) | US8427477B2 (zh) |
EP (1) | EP1982308B1 (zh) |
JP (1) | JP4958978B2 (zh) |
CN (1) | CN101617342B (zh) |
AT (1) | ATE503237T1 (zh) |
DE (1) | DE602007013356D1 (zh) |
WO (1) | WO2008086820A1 (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8717364B2 (en) | 2008-11-15 | 2014-05-06 | New BIS Safe Luxco S.a.r.l | Data visualization methods |
EP2315046B1 (en) | 2009-10-20 | 2012-11-07 | TomTec Imaging Systems GmbH | Method and device for pre-processing Doppler ultrasound data |
EP2719336B1 (en) * | 2012-10-12 | 2021-03-17 | Samsung Medison Co., Ltd. | Method for displaying ultrasound image using doppler data and ultrasound medical apparatus thereto |
CN103209307B (zh) * | 2013-04-18 | 2016-05-25 | 清华大学 | 编码重聚焦计算摄像方法及装置 |
JP7152958B2 (ja) * | 2019-01-10 | 2022-10-13 | 富士フイルムヘルスケア株式会社 | 超音波撮像装置及び画像処理方法 |
JP7377016B2 (ja) * | 2019-07-23 | 2023-11-09 | フクダ電子株式会社 | 超音波画像生成装置およびその制御方法 |
CN110751957B (zh) * | 2019-09-25 | 2020-10-27 | 电子科技大学 | 一种使用堆叠多尺度模块的语音增强方法 |
CN113989395A (zh) * | 2021-10-29 | 2022-01-28 | 北京字节跳动网络技术有限公司 | 一种绒毛效果生成方法、系统、设备、介质及产品 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5402533A (en) * | 1993-04-22 | 1995-03-28 | Apple Computer, Inc. | Method and apparatus for approximating a signed value between two endpoint values in a three-dimensional image rendering device |
CN1340791A (zh) * | 2000-08-29 | 2002-03-20 | 朗迅科技公司 | 对三维的图像重建执行线性内插的方法和装置 |
CN1720005A (zh) * | 2002-12-04 | 2006-01-11 | 皇家飞利浦电子股份有限公司 | 高帧频三维超声成像仪 |
CN1722178A (zh) * | 2004-06-30 | 2006-01-18 | 通用电气公司 | 图像的三维重建方法和系统 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2939594B2 (ja) * | 1990-03-20 | 1999-08-25 | ジーイー横河メディカルシステム株式会社 | カラーフロー表示装置 |
US5628321A (en) * | 1995-12-18 | 1997-05-13 | Diasonics Ultrasound, Inc. | Processing velocity information in an ultrasonic system |
US5724974A (en) * | 1996-03-22 | 1998-03-10 | Acuson Corporation | Two-dimensional ultrasound display system |
US5735797A (en) * | 1996-12-30 | 1998-04-07 | General Electric Company | Method and apparatus for combining topographic flow power imagery with a B-mode anatomical imagery |
US7678048B1 (en) * | 1999-09-14 | 2010-03-16 | Siemens Medical Solutions Usa, Inc. | Medical diagnostic ultrasound system and method |
JP2002113001A (ja) * | 2000-10-10 | 2002-04-16 | Toshiba Corp | 超音波診断装置および超音波信号処理方法 |
FI109311B (fi) * | 2001-03-16 | 2002-06-28 | Nokia Corp | Menetelmä informaatioelementin reunan määrittämiseksi, järjestelmä ja elektroniikkalaite |
US7245746B2 (en) * | 2001-06-12 | 2007-07-17 | Ge Medical Systems Global Technology Company, Llc | Ultrasound color characteristic mapping |
US7044913B2 (en) * | 2001-06-15 | 2006-05-16 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus |
US7206457B2 (en) * | 2001-11-27 | 2007-04-17 | Samsung Electronics Co., Ltd. | Method and apparatus for encoding and decoding key value data of coordinate interpolator |
US20040246252A1 (en) * | 2002-01-14 | 2004-12-09 | Morrow Jesse James | Method and apparatus for visualizing data |
US7136518B2 (en) * | 2003-04-18 | 2006-11-14 | Medispectra, Inc. | Methods and apparatus for displaying diagnostic data |
US7744538B2 (en) * | 2004-11-01 | 2010-06-29 | Siemens Medical Solutions Usa, Inc. | Minimum arc velocity interpolation for three-dimensional ultrasound imaging |
US7840623B2 (en) * | 2005-09-26 | 2010-11-23 | Dai Nippon Printing Co., Ltd. | Interpolator and designing method thereof |
-
2007
- 2007-01-16 US US11/912,406 patent/US8427477B2/en active Active
- 2007-01-16 AT AT07702793T patent/ATE503237T1/de active
- 2007-01-16 WO PCT/EP2007/000336 patent/WO2008086820A1/en active Application Filing
- 2007-01-16 EP EP07702793A patent/EP1982308B1/en not_active Not-in-force
- 2007-01-16 CN CN2007800498369A patent/CN101617342B/zh not_active Expired - Fee Related
- 2007-01-16 DE DE602007013356T patent/DE602007013356D1/de active Active
- 2007-01-16 JP JP2009545068A patent/JP4958978B2/ja not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5402533A (en) * | 1993-04-22 | 1995-03-28 | Apple Computer, Inc. | Method and apparatus for approximating a signed value between two endpoint values in a three-dimensional image rendering device |
CN1340791A (zh) * | 2000-08-29 | 2002-03-20 | 朗迅科技公司 | 对三维的图像重建执行线性内插的方法和装置 |
CN1720005A (zh) * | 2002-12-04 | 2006-01-11 | 皇家飞利浦电子股份有限公司 | 高帧频三维超声成像仪 |
CN1722178A (zh) * | 2004-06-30 | 2006-01-18 | 通用电气公司 | 图像的三维重建方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
US20090315888A1 (en) | 2009-12-24 |
ATE503237T1 (de) | 2011-04-15 |
EP1982308B1 (en) | 2011-03-23 |
CN101617342A (zh) | 2009-12-30 |
DE602007013356D1 (de) | 2011-05-05 |
US8427477B2 (en) | 2013-04-23 |
EP1982308A1 (en) | 2008-10-22 |
JP2010515495A (ja) | 2010-05-13 |
WO2008086820A1 (en) | 2008-07-24 |
JP4958978B2 (ja) | 2012-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101617342B (zh) | 用于动态信息的图形表示的方法和系统 | |
CN101357077B (zh) | 医学诊断成像中单模式数据的多体积绘制 | |
Rezk-Salama et al. | Interactive exploration of volume line integral convolution based on 3D-texture mapping | |
US7355597B2 (en) | Method, apparatus and computer program product for the interactive rendering of multivalued volume data with layered complementary values | |
CN110811686B (zh) | 超声流体成像方法及超声流体成像系统 | |
CN101542535B (zh) | 用于在2d图像中描绘3d图像的数据组的方法和设备 | |
Everitt | Interactive order-independent transparency | |
CN101689306B (zh) | 有效的二维及三维图形处理 | |
US6421057B1 (en) | Configurable volume rendering pipeline | |
Díaz et al. | Real-time ambient occlusion and halos with summed area tables | |
EP2124197A1 (en) | Image processing apparatus and computer program product | |
US20060256111A1 (en) | System and method for fast 3-dimensional data fusion | |
JP2000149064A (ja) | カットプレ―ン変更バ―ジョンを表示するためのボリュ―ムデ―タ集合を実時間で変更する装置及び方法 | |
EP2058672A1 (en) | Ultrasound imaging system including a graphic processing unit | |
CN107945176B (zh) | 一种彩色ivoct成像方法 | |
WO2021165165A1 (en) | Rendering three dimensional overlays on two dimensional images | |
Drouin et al. | PRISM: An open source framework for the interactive design of GPU volume rendering shaders | |
Razdan et al. | Volume visualization of multicolor laser confocal microscope data | |
Tatarchuk et al. | Advanced interactive medical visualization on the GPU | |
Nikolov et al. | Real time 3D visualization of ultrasonic data using a standard PC | |
US20030122821A1 (en) | Overlapping triangle meshes | |
Rossl et al. | Visualization of volume data with quadratic super splines | |
Holub et al. | Enabling real-time volume rendering of functional magnetic resonance imaging on an iOS device | |
CN109754869A (zh) | 着色的超声图像对应的着色描述符的呈现方法和系统 | |
Cohen | Focus and context for volume visualization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120613 |
|
CF01 | Termination of patent right due to non-payment of annual fee |