CN103295259B - 一种自适应采样的最小梯度夹角预积分光照方法 - Google Patents
一种自适应采样的最小梯度夹角预积分光照方法 Download PDFInfo
- Publication number
- CN103295259B CN103295259B CN201310213788.XA CN201310213788A CN103295259B CN 103295259 B CN103295259 B CN 103295259B CN 201310213788 A CN201310213788 A CN 201310213788A CN 103295259 B CN103295259 B CN 103295259B
- Authority
- CN
- China
- Prior art keywords
- sampling
- integration
- formula
- section
- lambda
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000010354 integration Effects 0.000 title claims abstract description 48
- 238000005286 illumination Methods 0.000 title claims abstract description 30
- 238000005070 sampling Methods 0.000 claims abstract description 89
- 230000008569 process Effects 0.000 claims abstract description 19
- 230000008859 change Effects 0.000 claims abstract description 12
- 239000012634 fragment Substances 0.000 claims abstract description 12
- 230000000694 effects Effects 0.000 claims abstract description 8
- 230000008676 import Effects 0.000 claims abstract description 6
- 238000009825 accumulation Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012546 transfer Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000012800 visualization Methods 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003760 hair shine Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Landscapes
- Image Generation (AREA)
Abstract
一种自适应采样的最小梯度夹角预积分光照方法,包括以下步骤:1)、导入医学体数据;2)、首先建立多张不同步长的预积分表,表的数量可以根据实际应用的精度需求来决定,在采样过程中,根据标量场的变化趋势进行自适应采样,在标量值变化剧烈的区域,遇到极值,可以采用较小的采样步长;而在标量值变化相对平缓的区域,采样步长可以相应增大;3)、通过插值重构的方式对预积分片段内的梯度进行重构,并计算预积分采样段中与光照方向夹角最小的梯度,并将该梯度用于预积分光照渲染中,从而加强了局部细节的光照效果;4)、根据自适应采样策略得到的不同步长的预积分片段,通过查找不同步长的预积分表获得相应的颜色值和不透明度,然后进行体绘制积分。
Description
技术领域
本发明涉及一种自适应采样的最小梯度夹角预积分光照方法
背景技术
医学可视化是科学可视化中最为活跃的研究领域之一,在临床医学上有着重要的实际价值,标志着一个国家的医学水平的先进程度,是一个国家综合实力的一种体现。随着核磁共振成像技术(MRI),计算机断层投影技术(CT)等医学设备的快速发展,通过医学影像设备已经可以得到人体内部某一组织比较清晰的二维数字图像序列,这些图像序列包含了人体内部各个器官的三维信息。然而医生只能通过这些二维序列图像,凭借自己的临床经验对病人的病情进行诊断,如此做法可能会造成诊断治疗的偏差。为了提高诊断的准确性和科学性,旨在提供人体组织三维空间信息以及三维观察手段的医学可视化技术应运而生。医学可视化就是运用计算机图形学和图像处理技术,将一系列的二维数字图像构造出人体组织的三维图形,从而在屏幕上形象逼真地显示人体组织内部的复杂结构,医护人员通过旋转、缩放、平移以及图像分割的操作,可以详细的了解病灶组织的空间信息,帮助医生做出准确的诊断和制定正确的手术方案,提高医学治疗的准确性和科学性。
医学可视化技术中一种重要的方法是直接体绘制技术,相比面绘制技术,直接体绘制技术可以清晰的显示人体内部的组织信息,使得绘制结果具有极强的透视效果,具有较高的保真性。近年来,随着医学成像技术的进步和仪器设备精度的提高,直接体绘制技术面临的医学体数据量越来越大,尤其在增加了光照渲染之后,采样过程的计算开销过大,梯度和光照强度值的计算都会增加计算的复杂度,因此,在实际应用中很难获得良好的交互性。
另一方面,在实际采样过程中,由于采样率过低造成采样不足,导致体数据中的某些细节特征的丢失。预积分方法能够在一定程度上解决上述的问题,该方法能够有效的克服奈奎斯特采样定律的限制,并且对于非线性变化的传输函数也有良好的绘制效果。然而,为了获得高质量的绘制结果,使用预积分方法仍然需要较小的采样步长,无法达到良好的交互性能。
综合上述,如何获得高质量的体绘制结果,同时保证用户交互的实时性,成为目前直接体绘制中一个亟待解决的重要问题。
发明内容
为了在高质量的绘制结果和良好的交互性之间得到较好的统一,本发明提供一种交互性好,准确性高,可靠性强的自适应采样的最小梯度夹角预积分光照算法。
本发明解决其技术问题所采用的技术方案是:
一种自适应采样的最小梯度夹角预积分光照方法,所述的自适应采样的预积分光照方法包括以下步骤:
1)、导入医学体数据;
2)、首先建立多张不同步长的预积分表,表的数量可以根据实际应用的精度需求来决定,在采样过程中,根据标量场的变化趋势进行自适应采样,在标量值变化剧烈的区域,遇到极值(极小值和极大值),可以采用较小的采样步长;而在标量值变化相对平缓的区域,采样步长可以相应增大。
3)、通过插值重构的方式对预积分片段内的梯度进行重构,并计算预积分采样段中与光照方向夹角最小的梯度,并将该梯度用于预积分光照渲染中,从而加强了局部细节的光照效果。4)、根据自适应采样策略得到的不同步长的预积分片段,通过查找不同步长的预积分表获得相应的颜色值和不透明度,然后进行体绘制积分。
作为优选的一种方案:所述的步骤2)包括以下步骤:
Step1根据实际应用的需求,确定多张不同采样步长的预积分表,这些表的步长按指数级变化,并且满足如下的约束条件:
Stepn=2-n*Stepn-1(1)
其中Step0表示第n级别的预积分表的采样步长,公式(1)初始条件为Step0=L(L为采样过程中的步长),并且满足n=0,1,2,3...
Step2构建不同采样步长的预积分表,积分公式如下:
其中sf和sb分别表示采样段Slab的入点Pa和出点Pb的标量值,L表示采样段Slab的长度;
Step3计算采样光线s(x)上di处的方向导数Dv(g(x)),我们根据公式(3)来计算方向导数。
其中梯度值可以通过常规的梯度计算方法得到,即中心作差法,向前作差法等。采样点的梯度值可以预先存放在3D纹理中,或是在绘制时实时计算得到。v是光线的采样方向。
Step4根据采样段[di,di+1]的入点和出点的方向导数来判断该采样段中是否含有标量极值,从而决定是否需要对该段采样进一步细分,当入点di和出点di+1的方向导数乘积小于0时,即满足公式(4),采样段[di,di+1]需要进一步细分。
Dv(g(di))Dv(g(di+1))<0(4)
其中,Dv(g(di))和Dv(g(di+1))分别是采样光线s(x)上di和di+1处的方向导数。
Step5若公式(4)成立,则对该采样段进行自适应的划分。假定方向导数在较小的范围内是线性变化的,即在一个采样段内可以认为是线性变化的,则在方向导数为0的位置dx取得标量值的极值。dx的位置可以通过公式(5)计算得到:
Step6判定dx的大小,若dx小于1/2,则极值点位于采样段[di,di+1]的前半段,否则位于采样段[di,di+1]的后半段;若位于前(后)半段,则计算前(后)半段两个端点的方向导数,然后计算前(后)半段内方向导数为0的位置dx;
Step7重复step6的过程,直到采样段内最小划分的采样步长小于或等于预先设定的最小步长。
作为优选的另一种方案:所述的步骤3)还包括以下步骤:
Step1本文采用Catmull-Rom曲线来重构位于采样段Slab内的梯度值,插值公式如下:
N(t)=N0*B0(t)+N1*B1(t)+N2*B2(t)+N3*B3(t)(6)
B0(t)=0.5(-t3+2t2-t)
B1(t)=0.5(3t3-5t2+2)
B2(t)=0.5(-3t3+4t2+t)
B3(t)=0.5(t3-t2)
0≤t≤1
其中N0,N1,N2,N3分别是一条光线上连续的四个采样点各自的梯度值,并且已经归一化处理。
公式(6)经过整理之后,得到公式(7)
N(t)=0.5*(-N0+3.0*N1-3.0*N2+N3)*t3
+0.5*(2.0*N0-5.0*N1+4.0*N2-N3)*t2(7)
+0.5*(-N0+N2)*t
+0.5*(2.0*N1)
Step2最小梯度夹角预积分算法的关键在于求L·N(t)和H·N(t)的最大值,即L、H与N的夹角最小。根据公式(7)得到公式(9)
L·N(t)=0.5*(-L·N0+3.0*L·N1-3.0*L·N2+L·N3)*t3
+0.5*(2.0*L·N0-5.0*L·N1+4.0*L·N2-L·N3)*t2(8)
+0.5*(-L·N0+L·N2)*t
+0.5*(2.0*L·N1)
0≤t≤1
令a=0.5*(-L·N0+3.0*L·N1-3.0*L·N2+L·N3)
b=0.5*(2.0*L·N0-5.0*L·N1+4.0*L·N2-L·N3)
c=0.5*(-L·N0+L·N2)
d=0.5*(2.0*L·N1),并且a,b,c,d都是常数
上述问题转为求F(t)=at3+bt2+ct+d在区间0≤t≤1的最大值F(t)max。
H·N(t)求最大值的过程同理可得。
Step3在得到F(t)max的函数之后,接下来计算体绘制的光照渲染,在引入Blinn-Phong局部光照模型之后,体绘制的积分公式如下:
其中Ka,Kd,Ks分别代表的是环境光ambient,漫反射光diffuse,镜面高光specular的系数,nδ表示光亮规模系数,cs表示镜面高光的颜色,L表示光源到采样点方向的单位向量,H表示L与视线方向夹角的一半,表示采样点的梯度值。
作为优选的另一种方案:所述的步骤4)还包括以下步骤:
Step1根据步骤2)中的自适应采样策略,原始采样段Slab的长度为l,假设被细分成两个采样小段Slab1和Slab2,长度分别为l1和l2,通过查找预积分表可以得到Slab1和Slab2颜色值分别为和,以及得到不透明度分别为和。
Step2通过将和进行融合,可以得到采样段Slab的颜色值,如下:
Step3通过将和进行融合,可以得到采样段Slab的颜色值如下:
Step4将和代入到原始的体绘制积分公式中,得到如下公式:
其中Ci+1和Ai+1分别表示积累颜色值和积累不透明度值。
Step5指定积累不透明度阈值A_acc,当Ai+1≥A_acc时,结束当前的光线的采样。
本发明的技术构思为:导入需要绘制的医学体数据;首先建立多张不同步长的预积分表,表的数量可以根据实际应用的精度需求来决定,在采样过程中,根据标量场的变化趋势进行自适应采样,在标量值变化剧烈的区域,遇到极值(极小值和极大值),可以采用较小的采样步长;而在标量值变化相对平缓的区域,采样步长可以相应增大。通过插值重构的方式对预积分片段内的梯度进行重构,并计算预积分采样段中与光照方向夹角最小的梯度,并将该梯度用于预积分光照渲染中,从而加强了局部细节的光照效果。根据自适应采样策略得到的不同步长的预积分片段,通过查找不同步长的预积分表获得相应的颜色值和不透明度,然后进行体绘制积分。
本发明的有益效果是:实现一种自适应采样的最小梯度夹角预积分光照方法,可以在临床应用中将高质量的绘制结果和良好的交互性之间得到较好的统一,同时这个过程不需要或是尽可能少的要求用户去干预。
附图说明
图1是自适应采样的最小梯度夹角预积分光照方法的系统总体结构图。
图2是自适应采样的预积分光照方法的程序简要流程。
图3是自适应采样过程的简要流程图。
图4是计算最小梯度夹角的流程图。
图5是自适应采样的预积分体绘制积分的简要流程图。
具体实施方式
下面结合附图对本发明作进一步描述:
参照图1~图5,一种自适应采样的最小梯度夹角预积分光照方法,所述自适应采样的预积分光照算法包括以下步骤:
1)、导入医学体数据;
2)、首先建立多张不同步长的预积分表,表的数量可以根据实际应用的精度需求来决定,在采样过程中,根据标量场的变化趋势进行自适应采样,在标量值变化剧烈的区域,遇到极值(极小值和极大值),可以采用较小的采样步长;而在标量值变化相对平缓的区域,采样步长可以相应增大。
3)、通过插值重构的方式对预积分片段内的梯度进行重构,并计算预积分采样段中与光照方向夹角最小的梯度,并将该梯度用于预积分光照渲染中,从而加强了局部细节的光照效果。4)、根据自适应采样策略得到的不同步长的预积分片段,通过查找不同步长的预积分表获得相应的颜色值和不透明度,然后进行体绘制积分。
图1中,自适应采样的最小梯度夹角预积分光照方法系统总体结构图是本方法的总体方案,数据导入以及处理部分使用C++语言实现,自适应采样及光照渲染部分利用GLSL语言编写,整个程序在VS2008平台上实现。
图2中,自适应采样的最小梯度夹角预积分光照方法的步骤是:首先导入医学体数据,然后根据标量场的变化趋势进行自适应采样,将预积分采样段进一步细分,通过插值重构的方式对预积分片段内的梯度进行重构,计算得到预积分片段中与光照方向夹角最小的梯度,最后查找不同部长的预积分表获得相应的颜色值和不透明度,接着进行体绘制积分。
在上述的自适应采样的预积分光照方法的步骤中,如何确定预积分采样段内极值点的位置,对后续进一步细分采样步长起着至关重要的作用,也即直接影响着绘制结果的质量,是该方法的核心步骤。
图3中描述了如何对预积分采样段进行自适应采样的划分,首先计算采样光线s(x)上di处的方向导数Dv(g(x))
其中梯度值可以通过常规的梯度计算方法得到,即中心作差法,向前作差法等。v是光线的采样方向。
根据采样段[di,di+1]的入点和出点的方向导数来判断该采样段中是否含有标量极值,从而决定是否需要对该段采样进一步细分,当入点di和出点di+1的方向导数乘积小于0时,采样段[di,di+1]需要进一步细分。
Dv(g(di))Dv(g(di+1))<0
其中,Dv(g(di))和Dv(g(di+1))分别是采样光线s(x)上di和di+1处的方向导数。
假定方向导数在较小的范围内是线性变化的,即在一个采样段内可以认为是线性变化的,则在方向导数为0的位置dx取得标量值的极值。dx的位置可以通过如下公式计算得到:
判定dx的大小,若dx小于1/2,则极值点位于采样段[di,di+1]的前半段,否则位于采样段[di,di+1]的后半段;若位于前(后)半段,则计算前(后)半段两个端点的方向导数,然后计算前(后)半段内方向导数为0的位置dx,直到采样段内最小划分的采样步长小于或等于预先设定的最小步长。
图4描述了计算最小梯度夹角的过程,首先采用Catmull-Rom曲线来重构位于采样段Slab内的梯度值:
N(t)=N0*B0(t)+N1*B1(t)+N2*B2(t)+N3*B3(t)
B0(t)=0.5(-t3+2t2-t)
B1(t)=0.5(3t3-5t2+2)
B2(t)=0.5(-3t3+4t2+t)
B3(t)=0.5(t3-t2)
0≤t≤1
其中N0,N1,N2,N3分别是一条光线上连续的四个采样点各自的梯度值,并且已经归一化处理。
整理之后,得到如下公式:
N(t)=0.5*(-N0+3.0*N1-3.0*N2+N3)*t3
+0.5*(2.0*N0-5.0*N1+4.0*N2-N3)*t2
+0.5*(-N0+N2)*t
+0.5*(2.0*N1)
接下来计算L·N(t)和H·N(t)的最大值,即L、H与N的夹角最小。
L·N(t)=0.5*(-L·N0+3.0*L·N1-3.0*L·N2+L·N3)*t3
+0.5*(2.0*L·N0-5.0*L·N1+4.0*L·N2-L·N3)*t2
+0.5*(-L·N0+L·N2)*t
+0.5*(2.0*L·N1)
0≤t≤1
令a=0.5*(-L·N0+3.0*L·N1-3.0*L·N2+L·N3)
b=0.5*(2.0*L·N0-5.0*L·N1+4.0*L·N2-L·N3)
c=0.5*(-L·N0+L·N2)
d=0.5*(2.0*L·N1),并且a,b,c,d都是常数
上述问题转为求F(t)=at3+bt2+ct+d在区间0≤t≤1的最大值F(t)max。
H·N(t)求最大值的过程同理可得。
在得到F(t)max的函数之后,接下来计算体绘制的光照渲染,在引入Blinn-Phong局部光照模型之后,体绘制的积分公式如下:
其中Ka,Kd,Ks分别代表的是环境光ambient,漫反射光diffuse,镜面高光specular的系数,nδ表示光亮规模系数,cs表示镜面高光的颜色,L表示光源到采样点方向的单位向量,H表示L与视线方向夹角的一半,表示采样点的梯度值。
图5描述了自适应采样的预积分体绘制积分的过程。首先根据自适应采样策略,原始采样段Slab的长度为l,假设被细分成两个采样小段Slab1和Slab2,长度分别为l1和l2,通过查找预积分表可以得到Slab1和Slab2颜色值分别为和,以及得到不透明度分别为和。
通过将和进行融合,可以得到采样段Slab的颜色值如下:
通过将和进行融合,可以得到采样段Slab的颜色值如下:
将和代入到原始的体绘制积分公式中,得到如下公式:
其中Ci+1和Ai+1分别表示积累颜色值和积累不透明度值。
最后根据设定的积累不透明度阈值A_acc,当Ai+1≥A_acc时,结束当前的光线的采样。
Claims (1)
1.一种自适应采样的最小梯度夹角预积分光照方法,所述的自适应采样的预积分光照方法包括以下步骤:
1)、导入医学体数据;
2)、首先建立多张不同步长的预积分表,表的数量可以根据实际应用的精度需求来决定,在采样过程中,根据标量场的变化趋势进行自适应采样,在标量值变化剧烈的区域,遇到极值(极小值和极大值),可以采用较小的采样步长;而在标量值变化相对平缓的区域,采样步长可以相应增大;所述的步骤2)包括以下步骤:
Step21根据实际应用的需求,确定多张不同采样步长的预积分表,这些表的步长按指数级变化,并且满足如下的约束条件:
Stepn=2-n*Stepn-1(1)
其中Stepn表示第n级别的预积分表的采样步长,公式(1)初始条件为Step0=L(L为采样过程中的步长),并且满足n=0,1,2,3...
Step22构建不同采样步长的预积分表,积分公式如下:
其中sf和sb分别表示采样段Slab的入点Pa和出点Pb的标量值,L表示采样段Slab的长度;
Step23计算采样光线s(x)上di处的方向导数Dv(g(x)),我们根据公式(3)来计算方向导数;
其中梯度值可以通过常规的梯度计算方法得到,即中心作差法,向前作差法等;采样点的梯度值可以预先存放在3D纹理中,或是在绘制时实时计算得到;v是光线的采样方向;
Step24根据采样段[di,di+1]的入点和出点的方向导数来判断该采样段中是否含有标量极值,从而决定是否需要对该段采样进一步细分,当入点di和出点di+1的方向导数乘积小于0时,即满足公式(4),采样段[di,di+1]需要进一步细分;
Dv(g(di))Dv(g(di+1))<0(4)
其中,Dv(g(di))和Dv(g(di+1))分别是采样光线s(x)上di和di+1处的方向导数;
Step25若公式(4)成立,则对该采样段进行自适应的划分;假定方向导数在较小的范围内是线性变化的,即在一个采样段内可以认为是线性变化的,则在方向导数为0的位置dx取得标量值的极值;dx的位置可以通过公式(5)计算得到:
Step26判定dx的大小,若dx小于1/2,则极值点位于采样段[di,di+1]的前半段,否则位于采样段[di,di+1]的后半段;若位于前(后)半段,则计算前(后)半段两个端点的方向导数,然后计算前(后)半段内方向导数为0的位置dx;
Step27重复step6的过程,直到采样段内最小划分的采样步长小于或等于预先设定的最小步长;
3)、通过插值重构的方式对预积分片段内的梯度进行重构,并计算预积分采样段中与光照方向夹角最小的梯度,并将该梯度用于预积分光照渲染中,从而加强了局部细节的光照效果;步骤3)还包括以下步骤:
Step31本文采用Catmull-Rom曲线来重构位于采样段Slab内的梯度值,插值公式如下:
N(t)=N0*B0(t)+N1*B1(t)+N2*B2(t)+N3*B3(t)(6)
B0(t)=0.5(-t3+2t2-t)
B1(t)=0.5(3t3-5t2+2)
B2(t)=0.5(-3t3+4t2+t)
B3(t)=0.5(t3-t2)
0≤t≤1
其中N0,N1,N2,N3分别是一条光线上连续的四个采样点各自的梯度值,并且已经归一化处理;
公式(6)经过整理之后,得到公式(7)
N(t)=0.5*(-N0+3.0*N1-3.0*N2+N3)*t3
+0.5*(2.0*N0-5.0*N1+4.0*N2-N3)*t2
(7)
+0.5*(-N0+N2)*t
+0.5*(2.0*N1)
Step32最小梯度夹角预积分算法的关键在于求L·N(t)和H·N(t)的最大值,即L、H与N的夹角最小;根据公式(7)得到公式(8)
L·N(t)=0.5*(-L·N0+3.0*L·N1-3.0*L·N2+L·N3)*t3
+0.5*(2.0*L·N0-5.0*L·N1+4.0*L·N2-L·N3)*t2
(8)
+0.5*(-L·N0+L·N2)*t
+0.5*(2.0*L·N1)
0≤t≤1
令a=0.5*(-L·N0+3.0*L·N1-3.0*L·N2+L·N3)
b=0.5*(2.0*L·N0-5.0*L·N1+4.0*L·N2-L·N3)
c=0.5*(-L·N0+L·N2)
d=0.5*(2.0*L·N1),并且a,b,c,d都是常数
上述问题转为求F(t)=at3+bt2+ct+d在区间0≤t≤1的最大值F(t)max;
H·N(t)求最大值的过程同理可得;
Step33在得到F(t)max的函数之后,接下来计算体绘制的光照渲染,在引入Blinn-Phong局部光照模型之后,体绘制的积分公式如下:
其中Ka,Kd,Ks分别代表的是环境光ambient,漫反射光diffuse,镜面高光specular的系数,nδ表示光亮规模系数,cs表示镜面高光的颜色,L表示光源到采样点方向的单位向量,H表示L与视线方向夹角的一半,表示采样点的梯度值;
4)、根据自适应采样策略得到的不同步长的预积分片段,通过查找不同步长的预积分表获得相应的颜色值和不透明度,然后进行体绘制积分;步骤4)还包括以下步骤:
Step41根据步骤2)中的自适应采样策略,原始采样段Slab的长度为l,假设被细分成两个采样小段Slab1和Slab2,长度分别为l1和l2,通过查找预积分表可以得到Slab1和Slab2颜色值分别为和以及得到不透明度分别为和
Step42通过将和进行融合,可以得到采样段Slab的颜色值如下:
Step43通过将和进行融合,可以得到采样段Slab的颜色值如下:
Step44将和代入到原始的体绘制积分公式中,得到如下公式:
其中Ci+1和Ai+1分别表示积累颜色值和积累不透明度值;
Step45指定积累不透明度阈值A_acc,当Ai+1≥A_acc时,结束当前的光线的采样。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310213788.XA CN103295259B (zh) | 2013-05-31 | 2013-05-31 | 一种自适应采样的最小梯度夹角预积分光照方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310213788.XA CN103295259B (zh) | 2013-05-31 | 2013-05-31 | 一种自适应采样的最小梯度夹角预积分光照方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103295259A CN103295259A (zh) | 2013-09-11 |
CN103295259B true CN103295259B (zh) | 2016-05-18 |
Family
ID=49096059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310213788.XA Active CN103295259B (zh) | 2013-05-31 | 2013-05-31 | 一种自适应采样的最小梯度夹角预积分光照方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103295259B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108460825A (zh) * | 2018-03-15 | 2018-08-28 | 中辰远瞻(北京)照明设计有限公司 | 一种夜景照明方案效果图制作方法 |
CN110728744B (zh) * | 2018-07-16 | 2023-09-19 | 海信视像科技股份有限公司 | 一种体绘制方法、装置及智能设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1818974A (zh) * | 2006-03-08 | 2006-08-16 | 杭州电子科技大学 | 一种多模态医学体数据三维可视化方法 |
CN102147932A (zh) * | 2011-03-30 | 2011-08-10 | 北京航空航天大学 | 一种基于可移动欧拉网格的模型驱动的烟雾模拟方法 |
-
2013
- 2013-05-31 CN CN201310213788.XA patent/CN103295259B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1818974A (zh) * | 2006-03-08 | 2006-08-16 | 杭州电子科技大学 | 一种多模态医学体数据三维可视化方法 |
CN102147932A (zh) * | 2011-03-30 | 2011-08-10 | 北京航空航天大学 | 一种基于可移动欧拉网格的模型驱动的烟雾模拟方法 |
Non-Patent Citations (2)
Title |
---|
Phong Normal Interpolation Revisited;Overveld C V et al;《ACM Transactions on Graphics》;19971031;第16卷(第4期);397-419 * |
自适应多预积分表体绘制算法;郑迟 等;《中国图形学进展2012:第十七届全国计算机辅助设计与图形学学术》;20120701;第233页第3段,第2.1节第3段,第2.2节第1段,图1 * |
Also Published As
Publication number | Publication date |
---|---|
CN103295259A (zh) | 2013-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Usher et al. | A virtual reality visualization tool for neuron tracing | |
US11395715B2 (en) | Methods and systems for generating and using 3D images in surgical settings | |
Ritter et al. | Real-time illustration of vascular structures | |
CN106659455A (zh) | 用于患者肺部的多张ct扫描片的统一坐标系 | |
WO2022109903A1 (zh) | 三维血管合成方法、系统及冠状动脉分析系统和存储介质 | |
CN102074039A (zh) | 一种体绘制裁剪面绘制方法 | |
CN106570859A (zh) | 一种牙齿图像分析系统 | |
Bouaoud et al. | DIVA, a 3D virtual reality platform, improves undergraduate craniofacial trauma education | |
CN103295259B (zh) | 一种自适应采样的最小梯度夹角预积分光照方法 | |
Li et al. | Design and evaluation of personalized percutaneous coronary intervention surgery simulation system | |
CN111105476B (zh) | 一种基于Marching Cubes的CT图像三维重建方法 | |
CN103345774A (zh) | 一种三维多尺度矢量化的建模方法 | |
CN112967786A (zh) | 基于多模影像与交互设备的解剖导航的构建方法及系统 | |
Bai et al. | A framework for computing angle of progression from transperineal ultrasound images for evaluating fetal head descent using a novel double branch network | |
Tian et al. | Efficient tooth gingival margin line reconstruction via adversarial learning | |
CN108510506A (zh) | 一种管状结构图像分割方法 | |
Chen et al. | The research and practice of medical image enhancement and 3D reconstruction system | |
US8675941B2 (en) | Method for visualization of airway anatomical changes | |
CN103218813A (zh) | 一种基于边界元的dot/xct双模式成像的图像重建算法 | |
CN104299254A (zh) | 一种基于cuda的移动式c形臂三维可视化方法 | |
CN107170009A (zh) | 一种基于医学图像的护目镜基弯数据测量方法 | |
CN101996419B (zh) | 超声图像三维重建方法与系统 | |
Bohak et al. | Neck veins: an interactive 3D visualization of head veins | |
Yan et al. | Virtual laryngoscopy: a noninvasive tool for the assessment of laryngeal tumor extent | |
CN109242964A (zh) | 三维医学模型的处理方法和装置 |
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 |