CN106708981A - 基于WebGL的MPR三维重建方法 - Google Patents
基于WebGL的MPR三维重建方法 Download PDFInfo
- Publication number
- CN106708981A CN106708981A CN201611121799.5A CN201611121799A CN106708981A CN 106708981 A CN106708981 A CN 106708981A CN 201611121799 A CN201611121799 A CN 201611121799A CN 106708981 A CN106708981 A CN 106708981A
- Authority
- CN
- China
- Prior art keywords
- image
- tangent plane
- dimensional
- webgl
- plane picture
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/95—Retrieval from the web
- G06F16/957—Browsing optimisation, e.g. caching or content distillation
- G06F16/9574—Browsing optimisation, e.g. caching or content distillation of access to content, e.g. by caching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种基于WebGL的MPR三维重建方法,包括:步骤一、通过canvas将各轴位图像绘制到一个指定区域,并缓存成WebGL的texture对象;步骤二、根据所需要显示的切面图像在由多个轴位图像所构建出的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标;再根据各点的三维坐标确定该点所在的轴位图像,进而得到该轴位图像所在的texture对象及该texture对象所在区域,使用GLSL的texture2D函数从该texture对象中取出该点的像素值;步骤三、将所需要显示的切面图像以二维图像的形式显示出来。本发明首次真正在浏览器中实现了MPR的三维重建,计算速度快,响应速度快。
Description
技术领域
本发明涉及图像处理技术,尤其涉及一种基于WebGL的MPR三维重建方法。
背景技术
多平面重组(multi-planar reconstruction,MPR),指多层螺旋CT后处理技术之一,在横断面扫描的基础上,对某些标线指定的组织进行不同方位的重组,从而得到冠状、矢状、斜位、曲面等方位的二维图像。
一般MPR都是由桌面应用程序或者APP实现,其中,MPR的桌面应用程序有用ITK等库实现的。但目前还没有在Web上实现MPR的。其原因在于,浏览器的容量较小,当需要在浏览器中查看大量网页时,将会导致大量图像缓存很慢,严重影响对大量图像的获取和查看速率,也更加无法实现从大量图像中生成所需要的切面图像。
发明内容
针对上述技术问题,本发明提供了一种基于WebGL的MPR三维重建方法,其首次真正在浏览器中实现了MPR的三维重建,计算速度快,响应速度快,对大数据的序列的兼容性好。
本发明提供的技术方案为:
一种基于WebGL的MPR三维重建方法,包括:
步骤一、提供多个轴位图像,通过canvas将各轴位图像绘制到一个指定区域,并缓存成WebGL的texture对象,其中,每个轴位图像上的各点对应有三维坐标和像素值;
步骤二、根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标;再根据所需要显示的切面图像上的各点的三维坐标确定该点所在的轴位图像,进而得到该轴位图像所在的texture对象及该texture对象所在区域,使用GLSL的texture2D函数从该texture对象中取出该点的像素值;
步骤三、根据所需要显示的切面图像上的各点的像素值,将所需要显示的切面图像以二维图像的形式显示出来。
优选的是,所述的基于WebGL的MPR三维重建方法中,所述步骤一中,通过canvas将各轴位图像绘制到一个指定区域,并缓存成WebGL的texture对象,其具体过程为:
构建一个canvas对象,并从该canvas对象取得WebGL context;从WebGL Conext中取出当前设备支持的texture个数的最大值和texture大小的最大值;根据texture大小的最大值生成Canvas对象的最大值,并将需要缓存的轴位图像使用canvas的drawimage方法绘制到指定的区域;当该canvas绘制满后,将该canvas转换成WebGL的texture对象。
优选的是,所述的基于WebGL的MPR三维重建方法中,所述步骤二中,所需要显示的切面图像为冠状位图像或矢状位图像。
优选的是,所述的基于WebGL的MPR三维重建方法中,所述步骤二中,当所需要显示的切面图像为斜位图像时,根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标,其具体过程为:
先根据与所需要显示的切面图像成一定角度的冠状位图像或矢状位图像在由多个轴位图像所构建的三维对象中的位置,确定该冠状位图像或矢状位图像上的各点的三维坐标;再根据该角度,通过三维坐标变换计算出所需要显示的切面图像上的各点的三维坐标。
优选的是,所述的基于WebGL的MPR三维重建方法中,所述步骤二中,当所需要显示的切面图像为曲面图像时,根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标,其具体过程为:
根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置以及该切面图像的弧度,确定所需要显示的切面图像上的各点的三维坐标。
本发明所述的基于WebGL的MPR三维重建方法具有以下有益效果:
(1)首次实现真正在浏览器中实现了MPR的三维重建;
(2)本发明中,大量的计算都使用WebGL实现,计算速度快(仅为1~5ms),响应速度快,对大数据的序列的兼容性好(可以支持500~1000张轴位图像);
(3)可以在显示时直接支持窗宽窗位处理等图像操作。
附图说明
图1为本发明所述的基于WebGL的MPR三维重建方法的流程示意图;
图2为本发明所述的基于WebGL的MPR三维重建方法的实施例一的示意图;
图3(a)为本发明所述的基于WebGL的MPR三维重建方法的实施例二中的一个头部轴位图像,图3(b)为本发明所述的基于WebGL的MPR三维重建方法的实施例二中所获得的头部冠状位图像,图3(c)为本发明所述的基于WebGL的MPR三维重建方法的实施例二中所获得的头部矢状位图像;
图4(a)为本发明所述的基于WebGL的MPR三维重建方法的实施例三中的一个胸部轴位图像,图4(b)为本发明所述的基于WebGL的MPR三维重建方法的实施例三中所获得的胸部冠状位图像,图4(c)为本发明所述的基于WebGL的MPR三维重建方法的实施例三中所获得的胸部矢状位图像。
具体实施方式
下面结合附图对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
如图1所示,本发明提供一种基于WebGL的MPR三维重建方法,包括:
步骤一、提供多个轴位图像,通过canvas将各轴位图像绘制到一个指定区域,并缓存成WebGL的texture对象,其中,每个轴位图像上的各点对应有三维坐标和像素值;
步骤二、根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标;再根据所需要显示的切面图像上的各点的三维坐标确定该点所在的轴位图像,进而得到该轴位图像所在的texture对象及该texture对象所在区域,使用GLSL(一种WebGL编程的Shader编程语言)的texture2D函数从该texture对象中取出该点的像素值;
步骤三、根据所需要显示的切面图像上的各点的像素值,将所需要显示的切面图像以二维图像的形式显示出来。
本发明所述方法的原理是根据坐标变化来实现的,即根据轴位图像可得到所有坐标的像素值,当需要显示某一切面图像,根据该切面图像的位置,针对切面图像上的每个点的坐标,可以换算得到该点在由轴面图像构造出来的三维对象中的坐标,从而得到该点的像素值。把切面图像上所有点的像素值都显示出来,就完成一个完整的切面展示了。
此外,步骤二中,由于WebGL对象平面坐标都是小数,当根据所需要显示的切面图像上的各点的三维坐标确定该点所在的轴位图像时,可以通过所需要显示的切面图像上的各点的三维坐标找到一个近似的轴位图像。
如果采用传统Web编程的方法来实现,取图片(即轴位图像)和取图片里某点的像素值都是十分费时的操作。换句话说,采用传统Web编程方法不可能实现MPR的三维重建。但是这种计算方式却正是WebGL所擅长的。
在一个优选的实施例中,所述的基于WebGL的MPR三维重建方法中,所述步骤一中,通过canvas将各轴位图像绘制到一个指定区域,并缓存成WebGL的texture对象,其具体过程为:
构建一个canvas对象,并从该canvas对象取得WebGL context;从WebGL Conext中取出当前设备支持的texture个数的最大值和texture大小的最大值;根据texture大小的最大值生成Canvas对象的最大值,并将需要缓存的轴位图像使用canvas的drawimage方法绘制到指定的区域;当该canvas绘制满后,将该canvas转换成WebGL的texture对象。
即,在该实施例中,通过采用WebGL,极大地提高了计算速度,缩短了响应时间。
在一个优选的实施例中,所述的基于WebGL的MPR三维重建方法中,所述步骤二中,所需要显示的切面图像为冠状位图像或矢状位图像。
在一个优选的实施例中,所述的基于WebGL的MPR三维重建方法中,所述步骤二中,当所需要显示的切面图像为斜位图像时,根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标,其具体过程为:
先根据与所需要显示的切面图像成一定角度的冠状位图像或矢状位图像在由多个轴位图像所构建的三维对象中的位置,确定该冠状位图像或矢状位图像上的各点的三维坐标;再根据该角度,通过三维坐标变换计算出所需要显示的切面图像上的各点的三维坐标。
请查阅图2,所需要显示的切面图像为一个对角线图像(为由右侧三维对象中的虚线表示),以采用冠状位图像(为由左侧三维对象中虚线表示)为例,先根据该冠状位图像在三维对象中的位置,来确定该冠状位图像上各点的三维坐标,再根据冠状位图像与对角线图像所成角度(本实施例中为45°),就可以通过三维坐标变换计算出对角线图像上各点的三维坐标。
在一个优选的实施例中,所述的基于WebGL的MPR三维重建方法中,所述步骤二中,当所需要显示的切面图像为曲面图像时,根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标,其具体过程为:
根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置以及该切面图像的弧度,确定所需要显示的切面图像上的各点的三维坐标。
实施例二
利用本发明所述的基于WebGL的MPR三维重建方法获取切面图像。图3(a)为头部的一个轴位图像;图3(b)为经过运算后获得的头部冠状位图像,图3(c)为经过运算后获得的头部矢状位图像。
实施例三
利用本发明所述的基于WebGL的MPR三维重建方法获取切面图像。图4(a)为胸部的一个轴位图像;图4(b)为经过运算后获得的胸部冠状位图像,图4(c)为经过运算后获得的胸部矢状位图像。
本发明所述的基于WebGL的MPR三维重建方法具有以下优点:
(1)首次实现真正在浏览器中实现了MPR的三维重建;
(2)本发明中,大量的计算都使用WebGL实现,计算速度快(仅为1~5ms),响应速度快,对大数据的序列的兼容性好(可以支持500~1000张轴位图像);
(3)可以在显示时直接支持窗宽窗位处理等图像操作。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。
Claims (5)
1.一种基于WebGL的MPR三维重建方法,其特征在于,包括:
步骤一、提供多个轴位图像,通过canvas将各轴位图像绘制到一个指定区域,并缓存成WebGL的texture对象,其中,每个轴位图像上的各点对应有三维坐标和像素值;
步骤二、根据所需要显示的切面图像在由多个轴位图像所构建出的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标;再根据所需要显示的切面图像上的各点的三维坐标确定该点所在的轴位图像,进而得到该轴位图像所在的texture对象及该texture对象所在区域,使用GLSL的texture2D函数从该texture对象中取出该点的像素值;
步骤三、根据所需要显示的切面图像上的各点的像素值,将所需要显示的切面图像以二维图像的形式显示出来。
2.如权利要求1所述的基于WebGL的MPR三维重建方法,其特征在于,所述步骤一中,通过canvas将各轴位图像绘制到一个指定区域,并缓存成WebGL的texture对象,其具体过程为:
构建一个canvas对象,并从该canvas对象取得WebGL context;从WebGL Conext中取出当前设备支持的texture个数的最大值和texture大小的最大值;根据texture大小的最大值生成Canvas对象的最大值,并将需要缓存的轴位图像使用canvas的drawimage方法绘制到指定的区域;当该canvas绘制满后,将该canvas转换成WebGL的texture对象。
3.如权利要求1或2所述的基于WebGL的MPR三维重建方法,其特征在于,所述步骤二中,所需要显示的切面图像为冠状位图像或矢状位图像。
4.如权利要求1或2所述的基于WebGL的MPR三维重建方法,其特征在于,所述步骤二中,当所需要显示的切面图像为斜位图像时,根据所需要显示的切面图像在由多个轴位图像所构建出的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标,其具体过程为:
先根据与所需要显示的切面图像成一定角度的冠状位图像或矢状位图像在由多个轴位图像所构建的三维对象中的位置,确定该冠状位图像或矢状位图像上的各点的三维坐标;再根据该角度,通过三维坐标变换计算出所需要显示的切面图像上的各点的三维坐标。
5.如权利要求1或2所述的基于WebGL的MPR三维重建方法,其特征在于,所述步骤二中,当所需要显示的切面图像为曲面图像时,根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置,确定所需要显示的切面图像上的各点的三维坐标,其具体过程为:
根据所需要显示的切面图像在由多个轴位图像所构建的三维对象中的位置以及该切面图像的弧度,确定所需要显示的切面图像上的各点的三维坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611121799.5A CN106708981A (zh) | 2016-12-08 | 2016-12-08 | 基于WebGL的MPR三维重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611121799.5A CN106708981A (zh) | 2016-12-08 | 2016-12-08 | 基于WebGL的MPR三维重建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106708981A true CN106708981A (zh) | 2017-05-24 |
Family
ID=58936409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611121799.5A Pending CN106708981A (zh) | 2016-12-08 | 2016-12-08 | 基于WebGL的MPR三维重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106708981A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107424211A (zh) * | 2017-06-15 | 2017-12-01 | 彭志勇 | WebGL 容积重建方法 |
CN109408059A (zh) * | 2018-11-12 | 2019-03-01 | 成都四方伟业软件股份有限公司 | 一种场景生成方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101011266A (zh) * | 2007-02-01 | 2007-08-08 | 上海交通大学 | 三维超声心动图四腔切面图像自动检测的方法 |
CN101034473A (zh) * | 2005-12-06 | 2007-09-12 | 西门子公司 | 计算机辅助识别断层造影中高对比度对象的方法和系统 |
CN101400306A (zh) * | 2006-03-09 | 2009-04-01 | 画像诊断株式会社 | 医用三维图像的显示控制程序和医用三维图像的显示方法 |
CN102411795A (zh) * | 2011-08-16 | 2012-04-11 | 中国人民解放军第三军医大学第一附属医院 | 一种关于冠状动脉的多模图像的联合显示方法 |
-
2016
- 2016-12-08 CN CN201611121799.5A patent/CN106708981A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101034473A (zh) * | 2005-12-06 | 2007-09-12 | 西门子公司 | 计算机辅助识别断层造影中高对比度对象的方法和系统 |
CN101400306A (zh) * | 2006-03-09 | 2009-04-01 | 画像诊断株式会社 | 医用三维图像的显示控制程序和医用三维图像的显示方法 |
CN101011266A (zh) * | 2007-02-01 | 2007-08-08 | 上海交通大学 | 三维超声心动图四腔切面图像自动检测的方法 |
CN102411795A (zh) * | 2011-08-16 | 2012-04-11 | 中国人民解放军第三军医大学第一附属医院 | 一种关于冠状动脉的多模图像的联合显示方法 |
Non-Patent Citations (1)
Title |
---|
陆凌牛: "《HTML5开发精要与实例详解》", 31 January 2012, 机械工业出版社 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107424211A (zh) * | 2017-06-15 | 2017-12-01 | 彭志勇 | WebGL 容积重建方法 |
CN109408059A (zh) * | 2018-11-12 | 2019-03-01 | 成都四方伟业软件股份有限公司 | 一种场景生成方法及装置 |
CN109408059B (zh) * | 2018-11-12 | 2022-02-11 | 成都四方伟业软件股份有限公司 | 一种场景生成方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105223212B (zh) | 安检ct系统及其方法 | |
JP4335817B2 (ja) | 関心領域指定方法、関心領域指定プログラム、関心領域指定装置 | |
CN104331918B (zh) | 基于深度图实时绘制室外地表遮挡剔除以及加速方法 | |
US20100156901A1 (en) | Method and apparatus for reconstructing 3d model | |
CN107369205B (zh) | 一种移动终端城市二三维联动显示方法 | |
US8149237B2 (en) | Information processing apparatus and program | |
CN104574263A (zh) | 一种基于gpu的快速三维超声重建和显示方法 | |
CN102842141A (zh) | 一种旋转x射线造影图像迭代重建方法 | |
CN111009033B (zh) | 一种基于OpenGL的病灶区域的可视化方法和系统 | |
US20040036674A1 (en) | Apparatus and method for associating voxel information with display positions | |
CN1870054A (zh) | 将多分支脉管及其环境同时投影到单个图像上 | |
KR20180023856A (ko) | 그래픽 처리 시스템 및 그래픽 프로세서 | |
JP4122314B2 (ja) | 投影画像処理方法、投影画像処理プログラム、投影画像処理装置 | |
CN106708981A (zh) | 基于WebGL的MPR三维重建方法 | |
CN101853518B (zh) | 基于各向异性体数据的错切变形体绘制方法 | |
CN108510580A (zh) | 一种椎骨ct图像三维可视化方法 | |
JPH0927045A (ja) | 画像処理装置 | |
CN111652876B (zh) | 一种三维盆底超声图像的检测方法 | |
CN104751438A (zh) | 用于交互地对医学三维图像配准的方法和设备 | |
CN105205862B (zh) | 一种图像三维重建方法及系统 | |
CN109829962B (zh) | 一种利用opengl的物体空间消隐线计算加速方法 | |
CN104915992A (zh) | 基于股骨ct图像的实时阴影体绘制方法 | |
CN1667652A (zh) | 一种基于图形处理芯片的视觉凸壳加速绘制方法 | |
CN101937575B (zh) | 一种快速高质量最大密度投影实现方法 | |
CN107330965B (zh) | 一种利用局部保守光栅化方法实现硬阴影反走样的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170524 |