CN100380406C - 放射治疗计划系统中获得三维解剖结构投影轮廓线的方法 - Google Patents
放射治疗计划系统中获得三维解剖结构投影轮廓线的方法 Download PDFInfo
- Publication number
- CN100380406C CN100380406C CNB200610022710XA CN200610022710A CN100380406C CN 100380406 C CN100380406 C CN 100380406C CN B200610022710X A CNB200610022710X A CN B200610022710XA CN 200610022710 A CN200610022710 A CN 200610022710A CN 100380406 C CN100380406 C CN 100380406C
- Authority
- CN
- China
- Prior art keywords
- point
- pixel
- projection
- tumour
- marking image
- 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
Images
Abstract
本发明公开了放射治疗计划系统中,制作肿瘤三维解剖结构的投影轮廓线方法,包括以下步骤:a)将用三角片结构表示的肿瘤三维解剖结构,投影到投影平面上建立标记图像;b)将投影平面上得到的每个三角片每条边的投影线段,映射到标记图像上画线,得到具有标记的像素点区域形成二值标记图像;c)在二值标记图像上对像素点按规定的方向进行扫描搜索,确定初始轮廓点后,进行各个轮廓点跟踪,以获取肿瘤三维解剖结构的投影轮廓线。本发明的优点是:用画线和轮廓跟踪方法代替了现有技术中对投影线段间的求交、排序、取舍等复杂运算,节省了大量计算时间,具有很高的适用性,为放射治疗计划系统的临床应用提供了便利条件。
Description
技术领域
本发明涉及一种图形处理方法,特别涉及一种在放射治疗计划系统中获取三维解剖结构的投影轮廓线的方法。
背景技术
目前,利用放射性射线杀死或抑制恶性肿瘤是目前治疗癌症患者的基本手段之一,专家估计大约有70%的肿瘤病人需要接受放射治疗。放射治疗的目标是最大限度地将射线集中到病变区内杀灭肿瘤细胞,而使周围正常组织和器官受到最大程度的保护。三维适形放射治疗是通过采用立体定位技术,在治疗机上附加特制形状的铅块或多叶准直器等技术实施非共面不规则野照射,使各照射野的形状在束眼观BEV上与靶区形状一致,使等剂量面在三维空间分布上紧扣靶区。在BEV上靶区形状即为靶区的投影轮廓线。
传统的射野设计的验证是在模拟机Simulator上进行,通过模拟机上得到的胶片,检验射野是否可行。这种模拟过程存在的缺陷是:整个模拟过程,包括机架、病床旋转角度放置,胶片制作等操作很耗时;x射线透视本身存在盲区,如肝、脾等,加上许多肿瘤位于人体内较深处,使得靶区和紧要器官在胶片上难以准确定位,从而影响最优放疗计划的获得。由于CT图像的灰度直接与人体组织密度对应,利用CT图像序列数据,通过复杂的投影算法,可以在计算机上得到传统X光机产生的X胶片图像,完成模拟机的验证过程,这被称为虚模拟,英文为Virtual Simulator。它的输出图像就是DRR,英文为Digitally Reconstructed Radiograph。虚模拟可以实现交互式的射野设计和验证。通常在DRR图像上需要同时显示紧要器官、肿瘤靶区的对应投影,这也需要计算其投影轮廓线。
在放射治疗计划系统中是以由三角片为单位构成的物体表面,来表达三维解剖结构的。将三维解剖结构向投影平面上投影,其投影轮廓线为包围投影区域的最大外轮廓。
人体器官结构复杂,构成器官表面的三角片数量通常很大,而且投影平面、投影方向都具有很大的任意性,因此采用图形学方法来计算投影轮廓线具有相当的复杂性。
发明内容
本发明的目的是:提供一种在放射治疗计划系统中获取肿瘤三维解剖结构投影轮廓线的方法,以实现照射野的适形设计和验证。
放射治疗计划系统中获取肿瘤三维解剖结构的投影轮廓线方法,包括以下步骤:
a)将用三角片结构表示的肿瘤三维解剖结构,投影到投影平面上建立标记图像;
b)将投影平面上得到的每个三个片每条边的投影线段,影射到标记图像上画线,得到具有标记的像素点区域形成二值标记图像;
c)在二值标记图像上对像素点按规定的方向进行扫描搜索,确定初始轮廓点后,进行各个轮廓点跟踪,以获取肿瘤三维解剖结构的投影轮廓线。
本发明的优点是:用画线和轮廓跟踪方法代替了现有技术中对投影线段间的求交、排序、取舍等复杂运算,因而节省了大量计算时间,具有很高的适用性,为放射治疗计划系统的临床应用提供了便利条件。
附图说明
图1是以三角片为单位构成的肿瘤三维解剖结构
图2是投影坐标系和标记图像;
图3是画投影线段的流程图;
图4是画线后的标记图像;
图5是方向链码示意图;
图6是轮廓跟踪示意图;
图7是投影轮廓线。
具体实施方式
下面结合附图,详细描述在放射治疗计划系统中,获取肿瘤三维解剖结构的投影轮廓线方法:
步骤一:将用三角片结构表示的肿瘤三维解剖结构,投影到投影平面上建立标记图像;
图1是以三角片为单位构成的肿瘤三维解剖结构,图2是投影坐标系和标记图像。在图2中,以投影平面为XOY平面,以投影方向即照射野的中心轴为Z轴,建立投影坐标系。将肿瘤三维解剖结构1的包围盒(所述包围盒即是包围肿瘤三维解剖结构的最小长方体)投影到投影坐标系的XOY平面,得到一个最大包围的矩形区域2,建立了与该矩形区域2大小相适应的标记图像,该图像原点(坐标为[0,0]的像素)在矩形区域的左上角3。
步骤二:将投影平面上的每个三个片每条边的投影线段,影射到标记图像上画线得到具有标记的像素点区域,形成二值标记图像;
图3是画投影线段的流程图。图中,按照该画投影线段的流程,将三角片的两个顶点投影到XOY平面,就可以得到投影线段的两个端点。然后将这两个端点映射到标记图像上,并在标记图像上画线,对其所经过的像素点做上标记,得到具体标记的像素点区域,形成二值标记图像。
具体方法是:
第一步ST1:计算所有三角片的顶点在投影坐标系下的坐标,将三角片个数n置于0;
第二步ST2:选取第n个三角片,将条边数m置于0;
第三步ST3:选取当前三角片的第m条边;
第四步ST4:计算该边在标记图像上投影线段的端点的映射坐标并画线;
第五步ST5:判断该三角片还有没有边未画线?如果有m加1,并返回第三步ST3步;如果没有,即m=3,进入第六步;
第六步ST6:判断还有没有三角片未画线?如果有n加1,并返回第二步ST2步;如果没有,即n等于所有三角片的个数,画投影线段结束。
图4是画线后的标记图像,它是一副标记过的像素点的像素值为1,其它像素点的像素值为0的二值图像。
步骤三:在二值标记图像上对像素点按规定的方向进行扫描搜索,确定初始轮廓点后进行各个轮廓点跟踪,以获取肿瘤三维解剖结构的投影轮廓线。
由于构成三维解剖结构的三角片在几何上具有拓扑一致性,在经过画线操作后的标记图像上,涂画标记过的像素所构成的区域,对应于肿瘤三维解剖结构的投影区域,其最外边缘就是投影轮廓线。
在本发明的实施方式中,给出了一个能适应于任意形状区域的轮廓跟踪算法。这种算法是基于方向链码的思想,图5给出了方向链码的示意图。从图中可以看出:对于任一个像素点P,其邻近的8个像素点,可以用方向链码的八个基本方向0、1、2、3、4、5、6、7表示。
设R是二值图像中一个连接区域,且区域上各像素点的灰度即像素值为1,R的边界点定义为R中这样的象素,它的8邻域中至少有一个点为不在R中的像素点,这些边界点的全体集合即构成区域R的轮廓线。
定义:direction为搜索过程中的当前方向。其值范围为0、1、2、3、4、5、6、7,分别对应于图5像素点P的方向链码的8个方向。
图6是轮廓跟踪示意图。跟踪方法如下:
第一步:搜索初始轮廓点;对二值标记图像自上而下,从左到右扫描每个象素点,直到找到一个像素点为轮廓点,将该点作为轮廓跟踪过程的初始点,并规定扫描搜索的方向;
第二步,扫描当前轮廓点的邻域点;以direction为扫描方向,按逆时针方向逐一测试当前轮廓点的邻域像素点,凡扫描得到像素值为0的邻域像素点都不是下一个轮廓点,只有扫描得到像素值为1的邻域像素点才是找到的轮廓点,重复这个跟踪过程可以找到一系列像素的轮廓点Pn;在本例中规定按逆时针行走方向搜索,接下来的扫描方向direction取像素点①的邻近链码5的方向;
第三步:测试得到的该目标区域的轮廓点是否为初始轮廓点①,如果是,该轮廓跟踪过程结束,跟踪得到的所有轮廓点构成轮廓线;如果该轮廓点不是初始轮廓点①,则以该轮廓点作为搜索下一个轮廓点的起始点,重复执行第二步。
以图6为例具体作法如下:
1)对图6自上而下、从左至右扫描每个像素点,直到找到一像素点①为轮廓点,将该点作为轮廓跟踪过程的起始点;
2)按方向链码图5所示,依图6中像素点①为起始,其扫描搜索的方向direction取该像素点①的邻近链码值5的方向,扫描了图6中的对应像素点②,因其像素值为0(未标记过),像素点②不是轮廓点,予以跳过;
3)接着direction按逆时针依该像素点①的邻近链码值6的方向,扫描图6中的对应像素点③,其像素值为1(标记过),所以是目标区域点,接受为轮廓点;
4)再将像素点③作为新的轮廓跟踪点,起始扫描方向direction取该像素点③的邻近链码值5的方向,扫描邻近像素点④,因其像素值为0(未标记过),所以该像素点④不是轮廓点;
5)接着direction按逆时针取该像素点③的邻近链码值6的方向,扫描下一个像素点⑤,其像素值为0,像素点⑤不是轮廓点;
6)direction再按逆时针取该像素点③的邻近链码值7的方向,扫描下一个像素点⑥,其像素值为1,该像索点⑥为轮廓点;
7)direction再按逆时针取该像素点⑥的邻近链码值5的方向,扫描下一个像素点⑦,其像素值为0,该像素点⑦不是轮廓点;
8)direction再按逆时针取该像素点⑥的邻近链码值6的方向,扫描下一个像素点⑧,其像素值为0,该像素点⑧不是轮廓点;
9)direction再按逆时针取该像素点⑥的邻近链码值7的方向,扫描下一个像素点⑨,其像素值为0,该像素点⑨不是轮廓点;
10)direction再按逆时针取像素点⑥的邻近链码值0的方向,扫描下一个像素点⑩,其像素值为1,该像素点⑩是轮廓点;
24)direction再按逆时针取像素点的邻近链码值4的方向,扫描下一个像素点,其像素值为1,该像素点是轮廓点;由于像素点与像素点①重合为一个像素点,因此由像素点①、像素点③、像素点⑥、像素点⑩、像素点、像素点、像索点各轮廓点,组成了标记图像上的肿瘤结构轮廓线。
从上例可以看出:在标记图像上对像素点进行轮廓跟踪扫描方法,即是对某一个进行搜索的像素点Pn,当确定它是当前轮廓点时,往下的跟踪搜索工作即是按规定的搜索方向,扫描该像素点Pn的周围邻近8个像素点,凡扫描得到像素值为0的像素点都不是下一个轮廓点,只有扫描像素值为l的像素点才是下一个轮廓点,接下去的轮廓跟踪扫描过程,是以新的轮廓点为扫描起点,将其周围邻近的8个像素点作为扫描对象逐一进行测试,通过其测试的像素值来确定下一个轮廓点,重复这个过程,直到搜索出的最后一个轮廓点与初始轮廓点相重合而得到轮廓线。图7是跟踪后得到的投影轮廓线,该投影轮廓线即为靶区的投影轮廓线。
Claims (1)
1.放射治疗计划系统中获得三维解剖结构投影轮廓线的方法,包括以下步骤:
步骤一,将用三角片结构表示的肿瘤三维解剖结构,投影到投影平面上建立标记图像;
步骤二,将投影平面上的每个三个片每条边的投影线段,影射到标记图像上画线得到具有标记的像素点区域,形成二值标记图像;
步骤三,在二值标记图像上对像素点按规定的方向进行扫描搜索,确定初始轮廓点后进行各个轮廓点跟踪,以获取肿瘤三维解剖结构的投影轮廓线;其中步骤一的具体作法是:
首先,制作以三角片为单位构成的肿瘤三维解剖结构;
接着,在投影平面上以投影方向即射野的中心轴为Z轴,建立投影坐标系;
然后,将肿瘤三维解剖结构的包围盒投影到投影坐标系的XOY平面上,得到一个最大包围的矩形区域,从而建立了与该矩形区域大小相适应的标记图像;
步骤二的具体作法是:
第一步ST1:计算所有三角片的顶点在投影坐标系下的坐标,将三角片个数n置于0;
第二步ST2:选取第n个三角片,将条边数m置于0;
第三步ST3:选取当前三角片的第m条边;
第四步ST4:计算该边在标记图像上投影线段的端点的映射坐标并画线;
第五步ST5:判断该三角片还有没有边未画线?如果有m加1,并返回第三步ST3步;如果没有,即m=3,进入第六步;
第六步ST6:判断还有没有三角片未画线?如果有n加1,并返回第二步ST2步;如果没有,即n等于所有三角片的个数,画投影线段结束;
步骤三的具体作法是:
第一步,搜索初始轮廓点;对二值标记图像自上而下,从左到右扫描每个象素点,直到找到一个像素点为轮廓点,将该点作为轮廓跟踪过程的初始轮廓点,并规定扫描搜索的方向;
第二步,扫描当前轮廓点的邻域点;以direction为扫描方向,按逆时针方向逐一测试当前轮廓点的邻域像素点,凡扫描得到像素值为0的邻域像素点都不是下一个轮廓点,只有扫描像素值为1的邻域像素点才是找到的轮廓点,重复这个跟踪过程可以找到一系列像素的轮廓点Pn;
第三步:测试得到的该目标区域的轮廓点Pn是否为初始轮廓点,如果是,该轮廓跟踪过程结束,跟踪得到的所有轮廓点构成轮廓线;如果该轮廓点Pn不是初始轮廓点,则以该轮廓点Pn作为搜索下一个轮廓点的起始点,重复执行第二步。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200610022710XA CN100380406C (zh) | 2006-12-29 | 2006-12-29 | 放射治疗计划系统中获得三维解剖结构投影轮廓线的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200610022710XA CN100380406C (zh) | 2006-12-29 | 2006-12-29 | 放射治疗计划系统中获得三维解剖结构投影轮廓线的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101000689A CN101000689A (zh) | 2007-07-18 |
CN100380406C true CN100380406C (zh) | 2008-04-09 |
Family
ID=38692647
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB200610022710XA Expired - Fee Related CN100380406C (zh) | 2006-12-29 | 2006-12-29 | 放射治疗计划系统中获得三维解剖结构投影轮廓线的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100380406C (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1946339A (zh) | 2004-02-20 | 2007-04-11 | 佛罗里达大学研究基金会公司 | 用于提供适形放射治疗同时对软组织进行成像的系统 |
JP4944007B2 (ja) * | 2007-12-21 | 2012-05-30 | 株式会社島精機製作所 | 輪郭抽出装置と輪郭抽出方法及び輪郭抽出プログラム |
US10561861B2 (en) | 2012-05-02 | 2020-02-18 | Viewray Technologies, Inc. | Videographic display of real-time medical treatment |
CN103530639B (zh) * | 2013-10-30 | 2017-10-31 | 湖南轻创科技有限公司 | 一种图像轮廓有序点集提取方法 |
CN105006002B (zh) * | 2015-08-31 | 2018-11-13 | 北京华拓金融服务外包有限公司 | 自动图像抠图方法和装置 |
EP3423153B1 (en) | 2016-03-02 | 2021-05-19 | ViewRay Technologies, Inc. | Particle therapy with magnetic resonance imaging |
EP3554635B1 (en) | 2016-12-13 | 2021-01-20 | ViewRay Technologies, Inc. | Radiation therapy systems |
CN116036499A (zh) | 2017-12-06 | 2023-05-02 | 优瑞技术公司 | 多模态放射疗法的优化 |
US11209509B2 (en) | 2018-05-16 | 2021-12-28 | Viewray Technologies, Inc. | Resistive electromagnet systems and methods |
CN112263788B (zh) * | 2020-11-02 | 2022-08-30 | 浙江省肿瘤医院 | 在放射治疗过程中形态变化的定量探测系统 |
CN116228634B (zh) * | 2022-12-07 | 2023-12-22 | 辉羲智能科技(上海)有限公司 | 用于图像检测的距离变换计算方法、应用、终端及介质 |
CN116652704B (zh) * | 2023-07-28 | 2023-10-31 | 中国人民解放军空军工程大学 | 一种自适应飞机结构外形的复合修理方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1250196A (zh) * | 1998-10-02 | 2000-04-12 | 国际商业机器公司 | 特征值抽取方法及设备,存储图象分析程序的存储介质 |
CN1804867A (zh) * | 2006-01-19 | 2006-07-19 | 上海交通大学 | 舌图像中舌体提取方法 |
-
2006
- 2006-12-29 CN CNB200610022710XA patent/CN100380406C/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1250196A (zh) * | 1998-10-02 | 2000-04-12 | 国际商业机器公司 | 特征值抽取方法及设备,存储图象分析程序的存储介质 |
CN1804867A (zh) * | 2006-01-19 | 2006-07-19 | 上海交通大学 | 舌图像中舌体提取方法 |
Non-Patent Citations (5)
Title |
---|
一种快速边缘跟踪与提取的新算法. 葛澎.微电子学与计算机,第22卷第8期. 2005 * |
三维放射治疗计划的坐标系统. 杨代伦,吴章文,罗正明.核技术,第27卷第2期. 2004 * |
图象边界的链码表示研究. 唐振军,张显全.微计算机信息,第23期. 2005 * |
基于标记信息的快速轮廓跟踪算法. 桑红石,傅勇,张天序,刘云生.华中科技大学学报(自然科学版),第33卷第9期. 2005 * |
面向RP骨骼CT图像的轮廓提取与精简. 任国成,王广春.山东大学学报(工学版),第34卷第3期. 2004 * |
Also Published As
Publication number | Publication date |
---|---|
CN101000689A (zh) | 2007-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100380406C (zh) | 放射治疗计划系统中获得三维解剖结构投影轮廓线的方法 | |
CN100496386C (zh) | 精确放射治疗计划系统 | |
US8818105B2 (en) | Image registration for image-guided surgery | |
Aird et al. | CT simulation for radiotherapy treatment planning | |
US7634304B2 (en) | Method and apparatus for lesion localization, definition and verification | |
EP2175931B1 (en) | Systems for compensating for changes in anatomy of radiotherapy patients | |
CN101600473B (zh) | 定量数据分析中的运动补偿及治疗装置 | |
US7453984B2 (en) | Real-time target confirmation for radiation therapy | |
CN103620646B (zh) | 用于生成图像配准图的系统、方法及计算机可读介质 | |
US20050047544A1 (en) | Apparatus and method for registering 2D radiographic images with images reconstructed from 3D scan data | |
EP3547257B1 (en) | Deformable image registration with awareness of computed tomography (ct) reconstruction area | |
CN105228527A (zh) | 使用从规划ct和轮廓生成的模板的用于igrt的肿瘤可见性的前瞻性评价 | |
Munbodh et al. | Automated 2D‐3D registration of a radiograph and a cone beam CT using line‐segment enhancement a | |
Liang et al. | A deep learning framework for prostate localization in cone beam CT‐guided radiotherapy | |
Mostafavi et al. | Detection and localization of radiotherapy targets by template matching | |
Arimura et al. | Computerized method for estimation of the location of a lung tumor on EPID cine images without implanted markers in stereotactic body radiotherapy | |
CN107480445A (zh) | 基于隐马尔可夫模型的肿瘤临床靶区侵犯概率计算方法 | |
Paskalev et al. | Daily target localization for prostate patients based on 3D image correlation | |
Van Der Put et al. | Contour propagation in MRI-guided radiotherapy treatment of cervical cancer: the accuracy of rigid, non-rigid and semi-automatic registrations | |
Jenkins et al. | Using a handheld stereo depth camera to overcome limited field‐of‐view in simulation imaging for radiation therapy treatment planning | |
Fox et al. | Performance evaluation of an automated image registration algorithm using an integrated kilovoltage imaging and guidance system | |
CN105816196A (zh) | 用于4dct成像的标记带及4dct成像方法 | |
Li et al. | A CT-based software tool for evaluating compensator quality in passively scattered proton therapy | |
Feng et al. | Invited Review—Image Registration in Veterinary Radiation Oncology: Indications, Implications, and Future Advances | |
Munbodh et al. | Automated 2D–3D registration of portal images and CT data using line‐segment enhancement |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080409 Termination date: 20171229 |