CN106202928B - 一种用于计算组织切割后新生成模型顶点纹理坐标的方法 - Google Patents
一种用于计算组织切割后新生成模型顶点纹理坐标的方法 Download PDFInfo
- Publication number
- CN106202928B CN106202928B CN201610541938.3A CN201610541938A CN106202928B CN 106202928 B CN106202928 B CN 106202928B CN 201610541938 A CN201610541938 A CN 201610541938A CN 106202928 B CN106202928 B CN 106202928B
- Authority
- CN
- China
- Prior art keywords
- newly
- vertex
- generated
- model vertices
- relationship
- 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
- 238000005520 cutting process Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000008520 organization Effects 0.000 title claims description 8
- 239000013598 vector Substances 0.000 claims description 17
- 238000013507 mapping Methods 0.000 abstract description 8
- 230000000694 effects Effects 0.000 abstract description 2
- 210000000056 organ Anatomy 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 210000002837 heart atrium Anatomy 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Data Mining & Analysis (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Pathology (AREA)
- General Physics & Mathematics (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Image Generation (AREA)
Abstract
一种用于组织切割后计算新生成顶点纹理坐标方法,其特征在于:包括如下步骤:记录所有包含切割路径点的三角面片的索引号;新生成的模型顶点是否在某一个三角面片内,遍历判断新生成的模型顶点是否在记录的三角形面片内并且记录相应三角面片的索引号;计算新生成的模型顶点与所在三角形顶点的关系,如果模型顶点在三角形内,计算出新生成顶点的关系,如果不在所记录的三角形面片内部,则利用点在平面的投影点,计算出与三角形顶点的关系;计算新生成的模型顶点的纹理坐标与相应三角形面片顶点的纹理坐标之间关系,即可计算出模型顶点的纹理坐标。本发明采用了计算新生成模型顶点的纹理坐标与相应三角形面片的顶点纹理坐标的关系的方法,可以不用根据某种纹理映射关系求取新的模型顶点纹理坐标,从而可以忽略中间的纹理映射关系,同时保证了纹理坐标的准确性,纹理无失真的效果。
Description
技术领域
本发明涉及计算机图形学等领域,尤其是虚拟手术中生物组织器官纹理映射领域。
背景技术
在虚拟手术场景中,要建立真实的视觉反馈系统,除了要正确逼真地进行器官模型建模外,还要通过人体器官的各种纹理细节来丰富其表现特征。切割是外科手术操作中必不可少的一部分。切割就会使模型三角形网络发生变化,就会沿着切割路径去除一部分面片和生成新的顶点,组建新的三角形网络。这个切割过程就会造成新的生成的顶点没有纹理坐标,造成切割周边的部分没有纹理,失去了真实感。
在虚拟手术中的模型的纹理坐标是通过Maya软件经过几种不同的纹理映射方法和经过手动调整生成的,因此无法通过某一个纹理映射关系计算出新生成模型顶点的纹理坐标。
发明内容
为了克服先有的虚拟手术无法通过某一个纹理映射关系计算出新生成模型顶点的纹理坐标的不足,本发明提供了一种用于计算组织切割后新生成模型顶点纹理坐标的方法,在虚拟手术切割操作后切口附近的纹理不失真、不扭曲、不重叠,实现更好的视觉反馈,对切割中新生成的模型顶点,快速、基本复现切口周边原有纹理。
本发明所采用的技术方案如下所示:
一种用于计算组织切割后新生成模型顶点纹理坐标的方法,所述的计算新生成顶点纹理坐标的方法包括如下步骤:
(1)记录所有包含切割路径点集合R,并搜索包含切割路径的三角形面片的索引号,记录三角形面片索引号集合L;
(2)判断新生成的模型顶点是否在某一个三角形面片内:遍历判断新生成的模型顶点是否在集合L中的某一个三角形面片内,并且记录相应三角形面片的索引号I;
(3)计算新生成的模型顶点P与所在三角形三个顶点的关系,如果新生成模型顶点在三角形内,计算出与新生成顶点的关系,如果不在所记录的三角形面片内部,则利用点在平面的投影点,计算出与三角形顶点的关系,过程如下:
由步骤(2)可知,如果新生成的模型顶点在记录集合L中的某一个的三角形面片内,则记录满足条件的相应的m和k的值,则会有P=A+m*(C-A)+k*(B-A);如果新生成的模型顶点不在三角形面片索引集合L中,则P在所有索引集合L所组成平面的外部,当P位于平面的外部时,根据P在距离最近三角形面片上的投影点P′与所在三角形顶点的关系,间接得到P与所在三角形顶点的关系;
(4)计算新生成的模型顶点P的纹理坐标与相应三角形面片顶点{vi0,vi1,vi2}的纹理坐标之间关系,即可计算出新生成模型顶点的纹理坐标(u,v)。
进一步所述步骤(3)中,求投影P′的过程如下:
步骤3.1设集合L组成平面的法向量s:{a,b,c},集合L中每个面片的法向量为:si:{ai,bi,ci},则投影面的法向量为:
因此可以设每个三角形面片si可以采用三个顶点{vi0,vi1,vi2},顶点分别设为vi0(xvi0,yvi0,zvi0),vi1(xvi1,yvi1,zvi1),vi2(xvi2,yvi2,zvi2)且顶点顺序按逆时针排列,则平面的法向量si为:
步骤3.2由步骤3.1得到平面的法向量s:{a,b,c},接着下来,需要遍历找到P点沿着s负方向与L集合中某一个三角形面片相交,求相交过程如下:
经过P点s方向的射线方程:
F=P+ts
每个三角形面片的平面方程为:
si·(F-vi0)=0
联立两个方程式得到:
若t>=0,则射线与平面相交,且交点为P+ts,则P′=P+ts,再根据步骤(2)计算P′在所在三角形面片顶点的关系。则可得到新生成顶点P与所在投影三角形面片顶点的关系。
再进一步,所述步骤(4)中,计算出新生成模型顶点的纹理坐标(u,v)的过程如下:假设ΔABC顶点的纹理坐标Auv=(u0,v0),Buv=(u1,v1),Cuv=(u2,v2)新生成模型顶点P的纹理坐标为Puv(u,v),由步骤(3)所得新生成顶点和所在三角形面片的关系式:P=A+m*(C-A)+k*(B-A),则有如下关系式:
Puv=Auv+m*(Cuv-Auv)+k*(Buv-Auv)
因此得到新生成模型顶点P的纹理坐标为Puv(u,v)。
本发明的有益效果体现在:本发明采用了计算新生成模型顶点的纹理坐标与相应三角形面片的顶点纹理坐标的关系的方法,可以不用根据某种纹理映射关系求取新的模型顶点纹理坐标,从而可以忽略中间的纹理映射关系,同时保证了纹理坐标的准确性,纹理无失真的效果。
附图说明
图1是三角形网络选取的切割路径的示意图。
图2是切割后新生成的模型顶点的示意图。
图3是三角形内部点的关系矢量图。
具体实施方式
下面对本发明做进一步的说明。
参照图1~图3,一种用于计算组织切割后新生成模型顶点纹理坐标的方法,所述的计算新生成模型顶点纹理坐标的方法包括以下步骤:
(1)记录所有包含切割路径点集合R,并搜索包含切割路径的三角形面片的索引号,记录三角形面片索引号集合L;
如说明书附图,图1所示三角形网络选取的切割路径,设切割路径顶点p0p1p2…pn- 1pn,如说明书附图,图2所示切割后新生成的模型顶点,设切割后新生成的模型顶点集合R包含模型顶点q0q1q2…qn-1qn,记录包含切割路径顶点的所有三角形面片的索引号,设索引集合为L={t0t1…tn-1tn};
(2)判断新生成的模型顶点是否在某一个三角形面片内:遍历判断新生成的模型顶点是否在集合L中的某一个三角形面片内,并且记录相应三角形面片的索引号I,具体过程如下:
如说明书附图,图3所示三角形内部点的关系矢量图,分别设三角形ΔABC顶点A(x0,y0,z0),B(x1,y1,z1),C(x2,y2,z2)和平面内任意点P(x,y,z),对于平面内任意点P(x,y,z)都可以有如下方程表示:
P=A+m*(C-A)+k*(B-A)
那么如果P位于ΔABC内部,m和k必须满足三个条件m>=0,k>=0和m+k<=1,因此我们需要继续计算m,k的值。令 则有公式:
因为一个方程,两个未知数,无法解出m和k,等式两边分别点乘得到等式:
解得方程可以得到:
根据求解的结果,判断是否满足P位于ΔABC内部的条件,就可以得到P是否在三角形ΔABC内部;
(3)计算新生成的模型顶点与所在三角形顶点的关系,如果模型顶点在三角形内,计算出新生成顶点的关系,如果不在所记录的三角形面片内部,则利用点在平面的投影点,计算出与三角形顶点的关系,计算过程如下:
由步骤(2)可知,如果新生成的模型顶点在记录集合L中的某一个的三角形面片内,则记录满足条件的相应的m和k的值,则会有P=A+m*(C-A)+k*(B-A);如果新生成的模型顶点不在三角形面片索引集合L中,则P在所有索引集合L所组成平面的外部,当P位于平面的外部时,可以根据P在距离最近三角面片上的投影点P′与所在三角形顶点的关系,间接得到P与所在三角形顶点的关系。因此先求投影点P′,详细步骤如下:
步骤3.1设集合L组成平面的法向量s:{a,b,c},集合L中每个面片的法向量为:si:{ai,bi,ci},则投影面的法向量为:
因此可以设每个三角形面片si可以采用三个顶点{vi0,vi1,vi2},顶点分别设为vi0(xvi0,yvi0,zvi0),vi1(xvi1,yvi1,zvi1),vi2(xvi2,yvi2,zvi2)且顶点顺序按逆时针排列,则平面的法向量si为:
步骤3.2由步骤3.1得到平面的法向量s:{a,b,c},接着下来,需要遍历找到P点沿着s方向与L集合中某一个三角形面片相交,求相交过程如下:
经过P点s方向的射线方程:
F=P+ts
每个三角形面片的平面方程为:
si·(F-vi0)=0
联立两个方程式得到:
若t>=0,则射线与平面相交,且交点为P+ts,则P′=P+ts,再根据步骤(2)计算P′在所在三角形面片顶点的关系。则可得到新生成顶点P与所在投影三角形面片顶点的关系;
(4)计算新生成的模型顶点的纹理坐标与相应三角形面片顶点的纹理坐标之间关系,具体实现过程如下:
假设ΔABC顶点的纹理坐标Auv=(u0,v0),Buv=(u1,v1),Cuv=(u2,v2)新生成模型顶点P的纹理坐标为Puv(u,v),由步骤(3)所得新生成顶点和所在三角形面片的关系式:P=A+m*(C-A)+k*(B-A),则有如下关系式:
Puv=Auv+m*(Cuv-Auv)+k*(Buv-Auv)因此得到新生成模型顶点P的纹理坐标为Puv(u,v)。
Claims (3)
1.一种用于计算组织切割后新生成顶点纹理坐标的方法,其特征在于:所述方法包括如下步骤:
(1)记录所有包含切割路径点集合R,并搜索包含切割路径的三角面片的索引号,记录三角面片索引号集合L;
(2)新生成的模型顶点是否在某一个三角形面片内:遍历判断新生成的模型顶点是否在集合L中的某一个三角形面片内,并且记录相应三角形面片的索引号I;
(3)计算新生成的模型顶点与所在三角形顶点的关系,如果模型顶点在三角形内,计算出与新生成顶点的关系,如果不在所记录的三角形面片内部,则利用点在平面的投影点,计算出与三角形顶点的关系过程如下:
由步骤(2)可知,如果新生成的模型顶点在记录集合L中的某一个的三角形面片内,则记录满足条件的相应的m和k的值,则会有P=A+m*(C-A)+k*(B-A);如果新生成的模型顶点不在三角形面片索引集合L中,则P在所有索引集合L所组成平面的外部,当P位于平面的外部时,根据P在距离最近三角形面片上的投影点P′与所在三角形顶点的关系,间接得到P与所在三角形顶点的关系;
(4)计算新生成的模型顶点的纹理坐标与相应三角形面片顶点的纹理坐标之间关系,即可计算出新生成模型顶点的纹理坐标(u,v)。
2.如权利要求1所述的一种用于计算组织切割后新生成顶点纹理坐标的方法,其特征在于:所述步骤(3)中,求投影点P′的过程如下:
步骤3.1设集合L组成平面的法向量s:{a,b,c},集合L中每个面片的法向量为:si:{ai,bi,ci},则投影面的法向量为:
因此可以设每个三角形面片si采用三个顶点{vi0,vi1,vi2},顶点分别设为vi0(xvi0,yvi0,zvi0),vi1(xvi1,yvi1,zvi1),vi2(xvi2,yvi2,zvi2)且顶点顺序按逆时针排列,则平面的法向量si为:
步骤3.2由步骤3.1得到平面的法向量s:{a,b,c},接着下来,需要遍历找到P点沿着s方向与L集合中某一个三角形面片相交,求相交过程如下:
经过P点s方向的射线方程:
F=P+ts
每个三角形面片的平面方程为:
si·(F-vi0)=0
联立两个方程式得到:
若t>=0,则射线与平面相交,且交点为P+ts,则P′=P+ts,再根据步骤(2)计算P′在所在三角形面片顶点的关系,则得到新生成顶点P与所在投影三角形面片顶点的关系。
3.如权利要求1所述的一种用于计算组织切割后新生成顶点纹理坐标的方法,其特征在于:所述步骤(4)中,所述计算新生成的模型顶点的纹理坐标与相应三角形面片顶点的纹理坐标之间的关系,方法如下:
假设ΔABC顶点的纹理坐标Auv=(u0,v0),Buv=(u1,v1),Cuv=(u2,v2)新生成模型顶点P的纹理坐标为Puv(u,v),由步骤(3)所得新生成顶点和所在三角形面片的关系式:P=A+m*(C-A)+k*(B-A),则有如下关系式:
Puv=Auv+m*(Cuv-Auv)+k*(Buv-Auv)
因此得到新生成模型顶点P的纹理坐标为Puv(u,v)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610541938.3A CN106202928B (zh) | 2016-07-05 | 2016-07-05 | 一种用于计算组织切割后新生成模型顶点纹理坐标的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610541938.3A CN106202928B (zh) | 2016-07-05 | 2016-07-05 | 一种用于计算组织切割后新生成模型顶点纹理坐标的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106202928A CN106202928A (zh) | 2016-12-07 |
CN106202928B true CN106202928B (zh) | 2019-02-26 |
Family
ID=57476213
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610541938.3A Active CN106202928B (zh) | 2016-07-05 | 2016-07-05 | 一种用于计算组织切割后新生成模型顶点纹理坐标的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106202928B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106846489B (zh) * | 2017-01-18 | 2019-04-30 | 南京大学 | 一种基于vtk处理obj文件的方法 |
CN108876919A (zh) * | 2017-05-12 | 2018-11-23 | 北京中天易观信息技术有限公司 | 一种对3dtin不规则模型切割构网并计算体积的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103699716A (zh) * | 2013-12-01 | 2014-04-02 | 北京航空航天大学 | 一种个性化三维医学图像驱动的器官虚拟显示方法 |
EP2826023A1 (en) * | 2012-06-27 | 2015-01-21 | Landmark Graphics Corporation | Systems and methods for creating a three-dimensional texture atlas |
CN104318056A (zh) * | 2014-09-24 | 2015-01-28 | 北京航空航天大学 | 基于位置动力学的软组织变形和切割模拟方法 |
-
2016
- 2016-07-05 CN CN201610541938.3A patent/CN106202928B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2826023A1 (en) * | 2012-06-27 | 2015-01-21 | Landmark Graphics Corporation | Systems and methods for creating a three-dimensional texture atlas |
CN103699716A (zh) * | 2013-12-01 | 2014-04-02 | 北京航空航天大学 | 一种个性化三维医学图像驱动的器官虚拟显示方法 |
CN104318056A (zh) * | 2014-09-24 | 2015-01-28 | 北京航空航天大学 | 基于位置动力学的软组织变形和切割模拟方法 |
Non-Patent Citations (3)
Title |
---|
《基于四面体网格的软组织位置动力学切割仿真算法》;白隽瑄 等;;《北京航空航天大学学报》;20150731;第41卷(第7期);第1343-1352页; |
《虚拟手术中软组织形变与切割技术研究》;朱玲;《中国博士学位论文全文数据库医药卫生科技辑》;20130215(第2期);第E080-10页; |
《虚拟手术训练系统中软组织切割模型的研究进展》;程强强 等;;《中国生物医学工程学报》;20150831;第34卷(第4期);第464-474页; |
Also Published As
Publication number | Publication date |
---|---|
CN106202928A (zh) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104794758B (zh) | 一种三维图像的裁剪方法 | |
EP2391988B1 (en) | Visual target tracking | |
CN108335353A (zh) | 动态场景的三维重建方法、装置和系统、服务器、介质 | |
CN106600665B (zh) | 一种相机动画路径生成方法和装置 | |
CN108885690A (zh) | 用于产生头相关传递函数滤波器的布置 | |
JP7442494B2 (ja) | 光学式捕捉によるパーソナライズされたhrtf | |
CN103679783B (zh) | 基于几何变形的三维动画角色模型皮肤变形方法 | |
CN106202928B (zh) | 一种用于计算组织切割后新生成模型顶点纹理坐标的方法 | |
CN104867106B (zh) | 一种深度图超分辨率方法 | |
Schönstein et al. | HRTF selection for binaural synthesis from a database using morphological parameters | |
CN106875462A (zh) | 一种基于元球模型和混合驱动方法的实时数字器官切割方法 | |
CN110349152A (zh) | 人脸图像质量检测方法及装置 | |
CN109741382A (zh) | 一种基于Kinect V2的实时三维重建方法与系统 | |
CN107958446A (zh) | 信息处理设备和信息处理方法 | |
CN106251281B (zh) | 一种基于形状插值的图像渐变方法 | |
CN106803094A (zh) | 基于多特征融合的三维模型形状相似性分析方法 | |
WO2018207282A1 (ja) | 対象物認識方法、装置、システム、プログラム | |
US9165393B1 (en) | Measuring stereoscopic quality in a three-dimensional computer-generated scene | |
WO2023116430A1 (zh) | 视频与城市信息模型三维场景融合方法、系统及存储介质 | |
CN106651835A (zh) | 基于熵的双视点无参考立体图像质量客观评价方法 | |
US20060160637A1 (en) | Method for simulating performance of golf club head | |
CN112101247A (zh) | 一种人脸姿态估计方法、装置、设备及存储介质 | |
CN107292865B (zh) | 一种基于二维图像处理的立体显示方法 | |
TWI536316B (zh) | 用於產生三維景象之裝置及由電腦執行之產生三維景象之方法 | |
CN107688720A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |