CN113409436B - 一种菱形像素排列的体绘制方法 - Google Patents

一种菱形像素排列的体绘制方法 Download PDF

Info

Publication number
CN113409436B
CN113409436B CN202110615894.5A CN202110615894A CN113409436B CN 113409436 B CN113409436 B CN 113409436B CN 202110615894 A CN202110615894 A CN 202110615894A CN 113409436 B CN113409436 B CN 113409436B
Authority
CN
China
Prior art keywords
diamond
pixel
shaped
projection
volume rendering
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
Application number
CN202110615894.5A
Other languages
English (en)
Other versions
CN113409436A (zh
Inventor
周庆华
覃文军
于鲲
张子明
杨金柱
赵大哲
Original Assignee
东北大学
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 东北大学 filed Critical 东北大学
Priority to CN202110615894.5A priority Critical patent/CN113409436B/zh
Publication of CN113409436A publication Critical patent/CN113409436A/zh
Application granted granted Critical
Publication of CN113409436B publication Critical patent/CN113409436B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/08Volume rendering
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Generation (AREA)

Abstract

本发明涉及一种菱形像素排列的体绘制方法,包括如下步骤:S1、采用光线投射法对三维数据场进行投射获得投影成像平面;S2、对投影成像平面绕其法向量方向旋转45度,得到一个菱形的投影平面空间;S3、将投影平面空间中菱形排列的像素通过转换为方形排列的像素,得到最终的绘制图像。本发明提供的菱形像素排列的体绘制方法通过像素排列方式提升显示效果的体绘制加速绘制方案,相同分别率下效果会有很大提升,进而在体绘制应用中,在保证绘制精度的前提下,通过降低绘制分辨率提高绘制效率。

Description

一种菱形像素排列的体绘制方法
技术领域
本发明属于体绘制技术领域,尤其涉及一种菱形像素排列的体绘制方法。
背景技术
体绘制是实现三维数据场可视化的核心算法之一,它可以直接而有效地展示数据场内部的物理现象,具有非常强的表现力。三维体绘制的作用就是将离散分布的三维数据场投影显示到二维屏幕上,然而,体绘制中采样计算与图像合成都是非常费时的,且两者的计算开销都与采样点总数成正比。
对体绘制加速的算法目前多采用针对体数据空间跳跃的方式,在光线传播过程中,通过对数据场进行编码,跳过背景及针对特定传输函数不可见的区域。对体数据场的编码有基于八叉树,k-d树等层次编码方式。
基于体数据编码的空间跳跃方式可以有效的提高体绘制过程的计算效率,减少绘制时间,但同时需要保存编码信息,增加额外的系统开销,而且编码依赖于传输函数信息,当传输函数发生变更时,需要从新调整编码进行,保证有效的绘制结果。
发明内容
(一)要解决的技术问题
针对现有存在的技术问题,本发明提供一种菱形像素排列的体绘制方法,通过像素排列方式提升显示效果的体绘制加速绘制方案,相同分别率下效果会有很大提升,进而在体绘制应用中,在保证绘制精度的前提下,通过降低绘制分辨率提高绘制效率。
(二)技术方案
为了达到上述目的,本发明采用的主要技术方案如下:
一种菱形像素排列的体绘制方法,包括如下步骤:
S1、采用光线投射法对三维数据场进行投射获得投影成像平面;
S2、对投影成像平面绕其法向量方向旋转45度,得到一个菱形的投影平面空间;
S3、将投影平面空间中菱形排列的像素通过转换为方形排列的像素,得到最终的绘制图像。
优选地,所述步骤S1还包括:
从图像空间的每一个像素点f(x,y)出发,按视线方向投射多条射线,并在每条射线上均选取采样点;
根据投射的射线上的采样点数据计算得到每一个像素点f(x,y)处的颜色值,将得到的每一个像素点f(x,y)的颜色值合成得到投影成像平面。
优选地,所述方法还包括:
每条射线均按照设定的步长穿过三维数据场;
沿着其中一条射线选择K个等距的采样点;
由距离采样点最近的8个数据点的颜色值和不透明值做三次线性插值,求出该采样点的颜色值和不透明度;
求出射线上所有采样点的颜色值和不透明度;
再将射线上各采样点的颜色值和不透明度值由前向后或由后向前加以合成,得到发出该射线的像素点f(x,y)处的颜色值。
优选地,所述射线的数量等于所述图像空间像素点f(x,y)数量。
优选地,所述步骤S3中将投影平面空间中菱形排列的像素转换为方形排列的像素的方法还包括:
在投影平面空间中沿着水平方向和垂直方向进行二次采样,获得二次采样点;
最后根据所述二次采样点的像素值绘制得到最终的图像。
优选地,所述二次采样点包括:菱形像素点和菱形区域的中间像素点;
所述菱形像素点的像素值为原像素值;
所述中间像素点的像素值为相邻四个顶点的像素值的均值。
优选地,当所述步骤S3中的菱形投影平面的分辨率为n×n时,将投影平面空间中菱形排列的像素通过转换为方形排列的像素,得到最终的绘制图像的分辨率为2n×2n。
(三)有益效果
本发明提供的一种菱形像素排列的体绘制方法,具有以下有益效果:在相同计算量下,提升体绘制的显示效果。
具体地,在相同分别率下效果会有很大提升,进而在体绘制应用中,在保证绘制精度的前提下,通过降低绘制分辨率提高绘制效率。
附图说明
图1为本发明一种菱形像素排列的体绘制方法的流程框图;
图2为本发明一种菱形像素排列的体绘制方法中的光线投射示意图;
图3为本发明一种菱形像素排列的体绘制方法中的分辨率为7*7的投影平面示意图;
图4为本发明一种菱形像素排列的体绘制方法分辨率为7*7的投影平面的二次采样示意图;
图5为本发明一种菱形像素排列的体绘制方法的实施例中的相同分辨率下的正方形排列a和菱形排列b的对比示意图。
具体实施方式
为了更好的解释本发明,以便于理解,下面结合附图,通过具体实施方式,对本发明作详细描述。
如图1和图5所示所示:本实施例中提供了一种菱形像素排列的体绘制方法,具体包括如下步骤:
S1、采用光线投射法对三维数据场进行投射获得投影成像平面;
S2、对投影成像平面绕其法向量方向旋转45度,得到一个菱形的投影平面空间;
S3、将投影平面空间中菱形排列的像素通过转换为方形排列的像素,得到最终的绘制图像。
如图2所示:本实施例中所述的步骤S1还包括:
从图像空间的每一个像素点f(x,y)出发,按视线方向投射多条射线,并在每条射线上均选取采样点;
根据投射的射线上的采样点数据计算得到每一个像素点f(x,y)处的颜色值,将得到的每一个像素点f(x,y)的颜色值合成得到投影成像平面。
本实施例中所述方法还包括:
每条射线均按照设定的步长穿过三维数据场;
沿着其中一条射线选择K个等距的采样点;
由距离采样点最近的8个数据点的颜色值和不透明值做三次线性插值,求出该采样点的颜色值和不透明度;
求出射线上所有采样点的颜色值和不透明度;
再将射线上各采样点的颜色值和不透明度值由前向后或由后向前加以合成,得到发出该射线的像素点f(x,y)处的颜色值。
应说明的是本实施例中所述射线的数量等于所述图像空间像素点f(x,y)数量。
本实施例中所述的步骤S3中将投影平面空间中菱形排列的像素转换为方形排列的像素的方法还包括:
在投影平面空间中沿着水平方向和垂直方向进行二次采样,获得二次采样点;
最后根据所述二次采样点的像素值绘制得到最终的图像。
本实施例中所述二次采样点包括:菱形像素点和菱形区域的中间像素点;
所述菱形像素点的像素值为原像素值;
所述中间像素点的像素值为相邻四个顶点的像素值的均值。
优选地,当所述步骤S3中的菱形投影平面的分辨率为n×n时,将投影平面空间中菱形排列的像素通过转换为方形排列的像素,得到最终的绘制图像的分辨率为2n×2n。
如图3所示,以一个7*7的投影平面为例,体绘制得到的投影结果按照菱形的像素排列方式,由于显示需要采用方形像素,因此,在进行绘制之前,需要将菱形排列的像素转换为方形排列的像素。
如图4所示,沿着水平方向和垂直方向进行采样,采样点有两类,第一类为菱形像素点,直接取值,第二类为菱形区域的中间像素点,该点的像素值为其相邻的四个顶点的像素值的均值。最后得到有效区域分辨率为14*14的绘制结果。
综上,如图5所示,本实施例中在相同分辨率下得到的两个图像的对比结果,采用菱形采样的结果要比正方形采样的结果具有更加丰富的细节信息。
最后应说明的是,以上所述的各实施例仅限于说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解;其依然可以对前述实施例记载的技术方案进行修改,或者对其中部分或全部技术特征进行等同替换;而这些修改或替换,并不会使相应技术方案的本质脱离本发明个实施例技术方案的范围。

Claims (5)

1.一种菱形像素排列的体绘制方法,其特征在于,包括如下步骤:
S1、采用光线投射法对三维数据场进行投射获得投影成像平面;
S2、对投影成像平面绕其法向量方向旋转45度,得到一个菱形的投影平面空间;
S3、将投影平面空间中菱形排列的像素通过转换为方形排列的像素,得到最终的绘制图像;
所述步骤S3中将投影平面空间中菱形排列的像素转换为方形排列的像素的方法还包括:
在投影平面空间中沿着水平方向和垂直方向进行二次采样,获得二次采样点;
最后根据所述二次采样点的像素值绘制得到最终的图像;
所述二次采样点包括:菱形像素点和菱形区域的中间像素点;
所述菱形像素点的像素值为原像素值;
所述中间像素点的像素值为相邻四个顶点的像素值的均值。
2.根据权利要求1所述的菱形像素排列的体绘制方法,其特征在于,所述步骤S1还包括:
从图像空间的每一个像素点f(x,y)出发,按视线方向投射多条射线,并在每条射线上均选取采样点;
根据投射的射线上的采样点数据计算得到每一个像素点f(x,y)处的颜色值,将得到的每一个像素点f(x,y)的颜色值合成得到投影成像平面。
3.根据权利要求2所述的菱形像素排列的体绘制方法,其特征在于,所述方法还包括:
每条射线均按照设定的步长穿过三维数据场;
沿着其中一条射线选择K个等距的采样点;
由距离采样点最近的8个数据点的颜色值和不透明值做三次线性插值,求出该采样点的颜色值和不透明度;
求出射线上所有采样点的颜色值和不透明度;
再将射线上各采样点的颜色值和不透明度值由前向后或由后向前加以合成,得到发出该射线的像素点f(x,y)处的颜色值。
4.根据权利要求2所述的菱形像素排列的体绘制方法,其特征在于,所述射线的数量等于所述图像空间像素点f(x,y)数量。
5.根据权利要求1所述的菱形像素排列的体绘制方法,其特征在于,
当所述步骤S3中的菱形投影平面的分辨率为n×n时,将投影平面空间中菱形排列的像素通过转换为方形排列的像素,得到最终的绘制图像的分辨率为2n×2n。
CN202110615894.5A 2021-06-02 2021-06-02 一种菱形像素排列的体绘制方法 Active CN113409436B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110615894.5A CN113409436B (zh) 2021-06-02 2021-06-02 一种菱形像素排列的体绘制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110615894.5A CN113409436B (zh) 2021-06-02 2021-06-02 一种菱形像素排列的体绘制方法

Publications (2)

Publication Number Publication Date
CN113409436A CN113409436A (zh) 2021-09-17
CN113409436B true CN113409436B (zh) 2023-10-13

Family

ID=77676005

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110615894.5A Active CN113409436B (zh) 2021-06-02 2021-06-02 一种菱形像素排列的体绘制方法

Country Status (1)

Country Link
CN (1) CN113409436B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5594842A (en) * 1994-09-06 1997-01-14 The Research Foundation Of State University Of New York Apparatus and method for real-time volume visualization
WO1998047105A1 (en) * 1997-04-15 1998-10-22 The Research Foundation Of State University Of New York Apparatus and method for parallel and perspective real-time volume visualization
EP1093085A2 (en) * 1999-10-01 2001-04-18 Mitsubishi Denki Kabushiki Kaisha Controller for a parallel pipelined volume rendering system
US6583787B1 (en) * 2000-02-28 2003-06-24 Mitsubishi Electric Research Laboratories, Inc. Rendering pipeline for surface elements
CN101739661A (zh) * 2009-12-09 2010-06-16 中国科学院自动化研究所 一种快速高保真地构造血管内超声长轴影像的方法
CN106104302A (zh) * 2014-01-07 2016-11-09 三星电子株式会社 辐射检测器及其断层扫描成像设备和辐射检测设备

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5594842A (en) * 1994-09-06 1997-01-14 The Research Foundation Of State University Of New York Apparatus and method for real-time volume visualization
WO1998047105A1 (en) * 1997-04-15 1998-10-22 The Research Foundation Of State University Of New York Apparatus and method for parallel and perspective real-time volume visualization
EP1093085A2 (en) * 1999-10-01 2001-04-18 Mitsubishi Denki Kabushiki Kaisha Controller for a parallel pipelined volume rendering system
US6583787B1 (en) * 2000-02-28 2003-06-24 Mitsubishi Electric Research Laboratories, Inc. Rendering pipeline for surface elements
CN101739661A (zh) * 2009-12-09 2010-06-16 中国科学院自动化研究所 一种快速高保真地构造血管内超声长轴影像的方法
CN106104302A (zh) * 2014-01-07 2016-11-09 三星电子株式会社 辐射检测器及其断层扫描成像设备和辐射检测设备

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
三维数据场体绘制研究进展;何丽君;王晓强;云健;包书哲;;大连民族学院学报;第14卷(第05期);第486-491页 *
交互式动态体绘制及其加速算法;刘保权, 周明全, 耿国华;中国图象图形学报;第7卷(第03期);第223-228页 *
基于CT影像的肺组织分割及其功能定量分析;耿欢;覃文军;杨金柱;边子健;赵大哲;;小型微型计算机系统;第37卷(第03期);第581-587页 *
基于光线投影的医学影像体绘制加速技术;陶玲;钱志余;陈春晓;;光电工程;第35卷(第01期);第89-93页 *
面向移动设备的各向异性纹理映射方法;杨柏林;潘志庚;;计算机辅助设计与图形学学报;第19卷(第05期);第569-574页 *

Also Published As

Publication number Publication date
CN113409436A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
US7443392B2 (en) Image processing program for 3D display, image processing apparatus, and 3D display system
CN106815807B (zh) 一种基于gpu‑cpu协同的无人机影像快速镶嵌方法
JP3229162B2 (ja) 圧縮ビデオ・ビットストリーム合成のシステムおよび方法
CN103077549B (zh) 一种基于kd树的实时大规模地形可视化实现方法
CN101571954B (zh) 亚像素配准
JP2018534705A (ja) 3次元画像のための三角形メッシュの生成
KR20030040920A (ko) 스테레오 카메라의 두영상과 양안차도를 이용한 다시점영상 합성 장치
CN102447925B (zh) 一种虚拟视点图像合成方法及装置
CN101312540A (zh) 基于深度和遮挡信息的虚拟视点合成方法
CN102834849A (zh) 进行立体视图像的描绘的图像描绘装置、图像描绘方法、图像描绘程序
CN101908216B (zh) 一种实现矢量字体的方法及装置
CN104718755A (zh) 减小含有相互类似的信息的多个图像的数据量的方法、程序及装置,以及表现含有相互类似的信息的多个图像的数据结构
CN109643462B (zh) 基于渲染引擎的实时图像处理方法以及显示设备
KR20080101998A (ko) 임의 시점 영상 합성 시 영상 보정 방법 및 장치
CN110913201A (zh) 一种光场显示结构及合成图像编码方法
EP2803041B1 (en) Method for multi-view mesh texturing and corresponding device
CN108696692A (zh) 光场成像的数据预处理方法
CN104737539A (zh) 减小含有相互类似的信息的多个图像的数据量的方法、程序及装置
CN113409436B (zh) 一种菱形像素排列的体绘制方法
CN102111562A (zh) 三维模型的投影变换方法及装置
Battiato et al. SVG rendering of real images using data dependent triangulation
CN102110300A (zh) 模拟镜头畸变的三维模型投影方法及装置
Lin et al. Fast multi-view image rendering method based on reverse search for matching
CN103916652B (zh) 视差矢量生成方法及装置
CN102110298A (zh) 虚拟演播室系统中的三维模型投影方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant