Three-dimensional synthetic blood vessel generation using stochastic L-systems
|
|
- Darcy Griffith
- 5 years ago
- Views:
Transcription
1 Three-dimensional synthetic blood vessel generation using stochastic L-systems Miguel A. Galarreta-Valverde a, Maysa M. G. Macedo a, Choukri Mekkaoui b and Marcel P. Jackowski* a a Department of Computer Science, Institute of Mathematics and Statistics, University of São Paulo, São Paulo, Brazil; b Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA. ABSTRACT Segmentation of blood vessels from magnetic resonance angiography (MRA) or computed tomography angiography (CTA) images is a complex process that usually takes a lot of computational resources. Also, most vascular segmentation and detection algorithms do not work properly due to the wide architectural variability of the blood vessels. Thus, the construction of convincing synthetic vascular trees makes it possible to validate new segmentation methodologies. In this work, an extension to the traditional Lindenmayer system (L-system) that generates synthetic 3D blood vessels by adding stochastic rules and parameters to the grammar is proposed. Towards this aim, we implement a parser and a generator of L-systems whose grammars simulate natural features of real vessels such as the bifurcation angle, average length and diameter, as well as vascular anomalies, such as aneurysms and stenoses. The resulting expressions are then used to create synthetic vessel images that mimic MRA and CTA images. In addition, this methodology allows for vessel growth to be limited by arbitrary 3D surfaces, and the vessel intensity profile can be tailored to match real angiographic intensities. Keywords: L -system, stochastic L-system, synthetic vessel, 3D blood vessel. 1. DESCRIPTION OF PURPOSE This paper aims to improve previous work 1, 2 concerning the generation of synthetic blood vessels by L-systems to produce more realistic synthetic blood vessels in 3D by adding stochastic and parametric rules to the L-system grammar. We show that convincing angiographic volumetric images can be obtained using our methodology. 2. INTRODUCTION Magnetic resonance angiography (MRA) or computed tomography angiography (CTA) images allows for a detailed analysis of blood vessels. However, due to its three-dimensional nature, powerful tools are needed to support radiologists in making the right diagnoses. For this purpose, in some cases, it is necessary to segment vessels from angiography images. Although much research in recent years has focused on solving this task, it still remains an open subject. In order to validate these techniques, a comparison between the segmentation results and ground-truth data must be performed. This information is easier to obtain from a synthetic vessel image than from a physical phantom, which may not be straightforward to build, possibly requiring specific acquisition parameters. The synthesis of artificial vessels using a mathematical model offers advantages such as acquisition time, cost and the flexibility to create vessels with different growth configurations and the possibility to incorporate anomalies. Synthetic blood vessels could also help in the simulation of surgeries, decreasing the risks involved with invasive real procedures. 1 While Lindenmayer systems 3 (L-systems) were originally developed to model cellular interactions, they have been used to generate synthetic blood vessels. 1, 2 However, these synthetic vessels were confined to 2D and their properties did not resemble those measured from real MR or CT angiographic images. Extensions to 3D have been proposed, but have been specific to modeling plant growth. 4 In this work, we propose a methodology that extends the traditional L-system grammar with stochastic and parametric rules which can be used to synthesize angiographic images that simulate real CT or MR angiography datasets. * mjack@ime.usp.br, Telephone: +55 (11) Medical Imaging 2013: Image Processing, edited by Sebastien Ourselin, David R. Haynor, Proc. of SPIE Vol. 8669, 86691I 2013 SPIE CCC code: /13/$18 doi: / Proc. of SPIE Vol I-1
2 3. METHODS Our methodology is divided into three steps: (i) string generation, (ii) synthetic vessel generation and (iii) discretization. The first step requires the grammar as input, and generates a string that represents the execution of the grammar in a given iteration. This string represents a sequence of instructions, where each character has an associated action. These actions are executed by the second step in order to generate a sequence of points with their respective properties, such as vessel diameter. This set of points and their properties compose our synthetic blood vessel meta-data. The discretization step finally creates the synthetic image by connecting the previously generated points while adding volumetric information to mimic a real angiographic image. 3.1 Grammar Definition Following formal L-system definitions, 4, 5 we define a rule as Name successor. An L-system will replace symbols according to its rules. For example, the rule F af b, will generate afb in the first iteration and aafbb in the second iteration. A stochastic L-system is one that contains at least one rule with an associated probability. For example, a grammar with the next two rules, F : 0.5 af b and F : 0.5 a, will generate afb or a in the first iteration, each with a probability of 0.5. Subsequently, one could get aab, aafbb or a in the second iteration. Our rules are composed by symbols whose respective actions are: f(length, diam): goes forward by length units in the direction vector and record vessel diameter diam at that point; +(angle) and -(angle): rotate the direction vector clockwise or counterclockwise by angle degrees on the plane defined by the current direction vector and a perpendicular vector to it, /(angle) and *(angle): rotate the current plane about the direction vector by angle degrees, [ : stores the current state formed by the diameter and position on the stack, ] : restores the last saved state from the stack and allows the generation of a new branch. Finally, { and } define the limits of a vessel segment, which is used in the discretization step. 3.2 String Generation In order to produce the final string whose actions will generate a synthetic vessel tree, intermediate functions are generated for each grammatical rule using a lexical and syntactic analyzer, this sequence is shown in the Figure 1. (String generation Grammar syntactic analyzer String of instructions Python functions intermediate functions generator Figure 1: Vessel string generation sequence Lexical and syntactic analyzer The task of a lexical analyzer is to produce tokens from a string given as an input. 6 These tokens will be used by the syntactic analyzer to generate an object code for each rule. An example for a very simple rule using the grammar is F(d0) f[+(th1)f(d1)]-(th2)f(d2). The corresponding object code will be generated as a function in the Python language, as follows: 1 def F(n, d0 ) : 2 i f n>0: 3 params=c a l c u l a t e B i f u r c a t i o n ( d0 ) 4 return f + [ + +( +s t r ( params [ th1 ] )+ ) +F(n 1,params [ d1 ] )+ ] + ( + 5 s t r ( params [ th2 ] )+ ) +F(n 1,params [ d2 ] ) 6 else : return F Proc. of SPIE Vol I-2
3 In this function, the parameter n represents the number of prescribed iterations and d0 is the diameter in the main vessel segment. The function calculatebifurcation uses d0 and vessel diameter definitions from previous studies 2 to define values for parameters d1 and d2, which are the diameters of the segments after a bifurcation, and th1 and th2, which represent the corresponding rotation angles, respectively. In the case of stochastic grammars, rules contain probabilities associated with their occurrence. Hence, an additional function is created for each group of rules. As an example, for a grammar formed by the two rules F(d0):0.3 f and F(d0):0.7 ff, the additional function created is: 1 def F(n, d0 ) : 2 r=random. random ( ) 3 i f r >=0.0 and r <0.3: return F1(n, d0 ) 4 i f r >=0.3 and r <1.0: return F2(n, d0 ) Once these functions are generated, the first rule known as the axiom, is executed. It expands each of the non-terminal character by the rule that defines it. The resulting string obtained after execution of a sample rule given a fixed number of iterations is shown the table below. Iteration Resulting string for the rule F(d0) f(d0)[+(th1)f(d1)]-(th2)f(d2) 1 f[+(43.47)f][-(31.58)f] 2 f[+(53.57)f[+(39.98)f][-(34.98)f]][-(22.24)f[+(65.38)f][-(12.39)f]] 3 f[+(37.46)f[+(38.94)f[+(18.5)f][-(57.87)f]][-(36.0)f[+(37.47)f][-(37.47)f]]][- (37.47)f[+(37.47)f[+(40.32)F][-(34.64)F]][-(37.47)f[+(55.6)F][-(20.46)F]]] 3.3 Synthetic vessel generation To obtain the synthetic vessels, we use the string generated from the previous step by splitting it into tokens and applying the respective actions. The instructions and their respective parameters are: -(angle), +(angle), *(angle), /(angle) and f(length, diameter). If one of the parameters is omitted, a default value is obtained by using definitions from previous studies. 2 For example, the string, +(25) is a token that represents a clockwise rotation by 25 degrees. After that, each token is executed iteratively according to the function that defines its behavior. However, the execution of this string may exceed the physical space where the vessel is allowed to grow. Therefore, before creating each point, a verification is performed to confirm whether or not the point is contained within the allotted space. If a point happens to fall outside this space, which is represented by a bounding surface, a heuristic approach is used to determine the closest point that falls inside the surface. 3.4 Discretization To generate a synthetic angiographic image, the final image size and spacing must be known. In addition, synthetic vessel points created by the previous step need to be scaled according to the desired image size. If the generated vessel trajectories are enclosed by a box defined by points A 0 and A 1, the size of the final image is ( X, Y, Z), and the desired spacing is (E x, E y, E z ) for each point P 1 (x, y, z) the correspoding scaled point P 2 is given by: ( P1x A 0x P 2 = X P1y A 0y, Y P1z A 0z, Z ). A 1x A 0x E x A 1y A 0y E y A 1z A 0z E z Once that all the points are normalized, a vessel skeleton is initially created by using the points combined with a modification of Bresenham s algorithm for 3D. When a vessel segment is formed by sub-segments, a B-Spline interpolation function is used to obtain intermediate points, producing smoother curves while avoiding abrupt angles. Proc. of SPIE Vol I-3
4 3.4.1 Intensity distributions A spherical kernel with a sigmoid or Gaussian intensity distribution is used to convolve each one of the vessel skeletons. Note that this intensity distribution is user selectable and mimic the vessel intensity profiles of real CT and MR angiographies. The kernel size is proportional to the vessel diamater and are parametrized as follows: I p = e d2 2σ 2, with σ = radius/2, radius is the value of the radius in the vessel segment. distribution is 1 P (d) = 1 + e, γ(d β) the values for γ and β with better experimental results are γ = 12/radius and β = 0.7 radius. The equation for the sigmoid (a) (b) (c) (d) (e) (f) (g) (h) (i) (j) (k) (l) Figure 2: (a)-(l) MIP projections of generated synthetic images. (g) simulation of the basic structure of kidneys, (j) simulation of aneurysms and (k) and (l) simulation of the estenosis. 4. RESULTS We have generated a synthetic angiographic image database of more than 40 images along with their respective grammars. A few of the images from this database are depicted by the maximum intensity projection (MIP) and are shown in the Figure 2. Figure 3 shows two perspectives of two datasets generated with spatial constraints, a liver, and a thigh segment, respectively. To create the constraint surfaces, a liver, 7, 8 a thigh and a femoral bone were segmented from real CT images. The time required to generate the images showed in Figure 2 using 10 iterations and without constraints was less than two minutes for each image, on a machine with a dual-core Proc. of SPIE Vol I-4
5 processor of 1596 Mhz and 8 GB of RAM. When constraints are on, the computational cost depends on the number of points of the bounding the surfaces. Thus, the vessel tree created within the liver surface containing a total of took 55 minutes for 11 iterations; and the vascular tree within the thigh surface took 10 minutes for 9 iterations. Resulting image sizes were 512x512x212, and 270x230x230 pixels, respectively. (a) (b) (c) (d) (e) (f) (g) (h) Figure 3: Different perspectives of synthetic blood vessels generated using the liver and thigh as physical constraints. (a), (b), (c) and (d) show the trajectories of the vessel segments (red) and the limiting surfaces (orange). (d), (e), (f) and (g) show two MIP projections of the resulting 3D angiographic images. 5. NEW OR BREAKTHROUGH WORK TO BE PRESENTED To the best of our knowledge, this is the first methodology that provides three-dimensional image results that simulate intensity profiles of real angiography images determined from a stochastic L-system grammar. Previous work on synthetic blood vessel generation 1 using L-systems were limited to 2D and did not produce volumetric angiographic images as a result. In addition, the grammar can be parameterized to generate organ-specific blood vessels and also allows for physical constraints to limit their growth. 6. DISCUSSION The main motivation for this work was to generate volumetric images of vessels with known properties (i.e. diameters) to help validate blood vessel segmentation algorithms. The synthetic images presented here show the flexibility of this methodology, and the diversity of the artificial vessel trees that can be generated by simply altering grammar rules and associated parameters. It is important to point out that, in order to achieve the realistic-looking vessels in this work, distinct grammars and parameters had to be used. This often required tedious and intricate handwork and parameter tuning to match the properties of real vascular trees. Hence, future work will examine the application of machine learning techniques to infer rules and parameters automatically from real angiographic images. Proc. of SPIE Vol I-5
6 7. CONCLUSIONS We have developed a methodology that generates realistic-looking synthetic blood vessels in three dimensions, taking into account arbitrary surfaces as physical constraints. The resulting volumetric images are visually convincing and the methodology is flexible enough to synthesize a wide diversity of vascular trees, and associated anomalies. We believe that this work may be of help in validating vascular segmentation algorithms using CT and MR images and may be useful as a tool for virtual surgery. REFERENCES [1] Liu, X., Liu, H., Hao, A., and Zhao, Q., Simulation of blood vessels for surgery simulators, in [2010 International Conference on Machine Vision and Human-machine Interface], , IEEE (2010). [2] Zamir, M., Arterial branching within the confines of fractal L-system formalism, The Journal of general physiology 118(3), 267 (2001). [3] Lindenmayer, A., Mathematical models for cellular interactions in development I. filaments with one-sided inputs, Journal of Theoretical Biology 18(3), (1968). [4] Qi, H., Qiu, R., and Jia, J., L-system based interactive and lightweight web3d tree modeling, in [Proceedings of the 10th International Conference on Virtual Reality Continuum and Its Applications in Industry], , ACM (2011). [5] Samal, A., Peterson, B., and Holliday, D., Recognition of plants using a stochastic L-system model, Journal of Electronic Imaging 11, 50 (2002). [6] Xiao, X. and Xu, Y., The design and implementation of c-like language interpreter, in [Intelligence Information Processing and Trusted Computing (IPTC), nd International Symposium on], , IEEE (2011). [7] Heimann, T. and et al., I. W., Comparison and Evaluation of Methods for Liver Segmentation from CT datasets, IEEE Transactions on Medical Imaging 28(8), (2009). [8] Dário, O., Raul, F., and Mauro, C., Segmentation of liver, its vessels and lesions from CT images for surgical planning, BioMedical Engineering OnLine 10(1), 1 23 (2011). Proc. of SPIE Vol I-6
Vessel Centerline Tracking in CTA and MRA Images Using Hough Transform
Vessel Centerline Tracking in CTA and MRA Images Using Hough Transform Maysa M.G. Macedo 1, Choukri Mekkaoui 2, and Marcel P. Jackowski 1 1 University of São Paulo, Department of Computer Science, Rua
More informationProbabilistic Tracking and Model-based Segmentation of 3D Tubular Structures
Probabilistic Tracking and Model-based Segmentation of 3D Tubular Structures Stefan Wörz, William J. Godinez, Karl Rohr University of Heidelberg, BIOQUANT, IPMB, and DKFZ Heidelberg, Dept. Bioinformatics
More informationComparison of Vessel Segmentations Using STAPLE
Comparison of Vessel Segmentations Using STAPLE Julien Jomier, Vincent LeDigarcher, and Stephen R. Aylward Computer-Aided Diagnosis and Display Lab, The University of North Carolina at Chapel Hill, Department
More informationComparison of Vessel Segmentations using STAPLE
Comparison of Vessel Segmentations using STAPLE Julien Jomier, Vincent LeDigarcher, and Stephen R. Aylward Computer-Aided Diagnosis and Display Lab The University of North Carolina at Chapel Hill, Department
More informationCS 223B Computer Vision Problem Set 3
CS 223B Computer Vision Problem Set 3 Due: Feb. 22 nd, 2011 1 Probabilistic Recursion for Tracking In this problem you will derive a method for tracking a point of interest through a sequence of images.
More informationTopics. Recursive tree models. Procedural approach L-systems. Image-based approach. Billboarding
Plant Modeling Topics Recursive tree models Billboarding Procedural approach L-systems Image-based approach Tree Model The structure of a tree Trunk (linkage) Branches (linkage, child of trunk node) Leaves/Buds/flowers/fruit
More informationAutomatic Cerebral Aneurysm Detection in Multimodal Angiographic Images
Automatic Cerebral Aneurysm Detection in Multimodal Angiographic Images Clemens M. Hentschke, Oliver Beuing, Rosa Nickl and Klaus D. Tönnies Abstract We propose a system to automatically detect cerebral
More informationCS 231A Computer Vision (Fall 2012) Problem Set 3
CS 231A Computer Vision (Fall 2012) Problem Set 3 Due: Nov. 13 th, 2012 (2:15pm) 1 Probabilistic Recursion for Tracking (20 points) In this problem you will derive a method for tracking a point of interest
More informationABSTRACT 1. INTRODUCTION 2. METHODS
Finding Seeds for Segmentation Using Statistical Fusion Fangxu Xing *a, Andrew J. Asman b, Jerry L. Prince a,c, Bennett A. Landman b,c,d a Department of Electrical and Computer Engineering, Johns Hopkins
More information3D Guide Wire Navigation from Single Plane Fluoroscopic Images in Abdominal Catheterizations
3D Guide Wire Navigation from Single Plane Fluoroscopic Images in Abdominal Catheterizations Martin Groher 2, Frederik Bender 1, Ali Khamene 3, Wolfgang Wein 3, Tim Hauke Heibel 2, Nassir Navab 2 1 Siemens
More informationRigid and Deformable Vasculature-to-Image Registration : a Hierarchical Approach
Rigid and Deformable Vasculature-to-Image Registration : a Hierarchical Approach Julien Jomier and Stephen R. Aylward Computer-Aided Diagnosis and Display Lab The University of North Carolina at Chapel
More informationNovel Approach of Modeling Self Similar Objects using Parallel String Rewriting Methods through Combined L-System Techniques
International Journal of Information & Computation Technology. ISSN 0974-2255 Volume 2, Number 1 (2012), pp. 1-12 International Research Publications House http://www. ripublication.com Novel Approach
More informationDevelopmental Systems
Developmental Systems 1 Biological systems Early development of the Drosophila fly http://flybase.bio.indiana.edu dorsal view lateral view 2 Biological systems Early development of Drosophila [Slack 2006]
More informationAutomatic Ascending Aorta Detection in CTA Datasets
Automatic Ascending Aorta Detection in CTA Datasets Stefan C. Saur 1, Caroline Kühnel 2, Tobias Boskamp 2, Gábor Székely 1, Philippe Cattin 1,3 1 Computer Vision Laboratory, ETH Zurich, 8092 Zurich, Switzerland
More informationModeling and preoperative planning for kidney surgery
Modeling and preoperative planning for kidney surgery Refael Vivanti Computer Aided Surgery and Medical Image Processing Lab Hebrew University of Jerusalem, Israel Advisor: Prof. Leo Joskowicz Clinical
More information82 REGISTRATION OF RETINOGRAPHIES
82 REGISTRATION OF RETINOGRAPHIES 3.3 Our method Our method resembles the human approach to image matching in the sense that we also employ as guidelines features common to both images. It seems natural
More informationA Study of Medical Image Analysis System
Indian Journal of Science and Technology, Vol 8(25), DOI: 10.17485/ijst/2015/v8i25/80492, October 2015 ISSN (Print) : 0974-6846 ISSN (Online) : 0974-5645 A Study of Medical Image Analysis System Kim Tae-Eun
More informationVolume visualization. Volume visualization. Volume visualization methods. Sources of volume visualization. Sources of volume visualization
Volume visualization Volume visualization Volumes are special cases of scalar data: regular 3D grids of scalars, typically interpreted as density values. Each data value is assumed to describe a cubic
More informationVolume 2, Issue 9, September 2014 ISSN
Fingerprint Verification of the Digital Images by Using the Discrete Cosine Transformation, Run length Encoding, Fourier transformation and Correlation. Palvee Sharma 1, Dr. Rajeev Mahajan 2 1M.Tech Student
More informationVessel Explorer: a tool for quantitative measurements in CT and MR angiography
Clinical applications Vessel Explorer: a tool for quantitative measurements in CT and MR angiography J. Oliván Bescós J. Sonnemans R. Habets J. Peters H. van den Bosch T. Leiner Healthcare Informatics/Patient
More informationGenetic L-System Programming: Breeding and Evolving Artificial Flowers with Mathematica
Genetic L-System Programming: Breeding and Evolving Artificial Flowers with Mathematica C. Jacob, jacob@informatik.uni-erlangen.de, Chair of Programming Languages, Department of Computer Science, University
More informationImproving Positron Emission Tomography Imaging with Machine Learning David Fan-Chung Hsu CS 229 Fall
Improving Positron Emission Tomography Imaging with Machine Learning David Fan-Chung Hsu (fcdh@stanford.edu), CS 229 Fall 2014-15 1. Introduction and Motivation High- resolution Positron Emission Tomography
More informationApplication of level set based method for segmentation of blood vessels in angiography images
Lodz University of Technology Faculty of Electrical, Electronic, Computer and Control Engineering Institute of Electronics PhD Thesis Application of level set based method for segmentation of blood vessels
More informationChapter 87 Real-Time Rendering of Forest Scenes Based on LOD
Chapter 87 Real-Time Rendering of Forest Scenes Based on LOD Hao Li, Fuyan Liu and Shibo Yu Abstract Using the stochastic L-system for modeling the trees. Modeling the trees includes two sides, the trunk
More informationdoi: /
Yiting Xie ; Anthony P. Reeves; Single 3D cell segmentation from optical CT microscope images. Proc. SPIE 934, Medical Imaging 214: Image Processing, 9343B (March 21, 214); doi:1.1117/12.243852. (214)
More informationComputational Radiology Lab, Children s Hospital, Harvard Medical School, Boston, MA.
Shape prior integration in discrete optimization segmentation algorithms M. Freiman Computational Radiology Lab, Children s Hospital, Harvard Medical School, Boston, MA. Email: moti.freiman@childrens.harvard.edu
More informationCREATION OF THE TREE MODEL ADAPTING TO ENVIRONMENT
CREATION OF THE TREE MODEL ADAPTING TO ENVIRONMENT Ryota Ueno Yoshio Ohno {ryota ohno}@on.ics.keio.ac.jp Graduate School of Science and Technology, Keio University 3-14-1 Hiyoshi, Kohoku-ku Yokohama 223-8522
More informationMedical Image Processing: Image Reconstruction and 3D Renderings
Medical Image Processing: Image Reconstruction and 3D Renderings 김보형 서울대학교컴퓨터공학부 Computer Graphics and Image Processing Lab. 2011. 3. 23 1 Computer Graphics & Image Processing Computer Graphics : Create,
More informationBioimage Informatics
Bioimage Informatics Lecture 13, Spring 2012 Bioimage Data Analysis (IV) Image Segmentation (part 2) Lecture 13 February 29, 2012 1 Outline Review: Steger s line/curve detection algorithm Intensity thresholding
More informationcoding of various parts showing different features, the possibility of rotation or of hiding covering parts of the object's surface to gain an insight
Three-Dimensional Object Reconstruction from Layered Spatial Data Michael Dangl and Robert Sablatnig Vienna University of Technology, Institute of Computer Aided Automation, Pattern Recognition and Image
More informationJournal of Applied Mathematics and Computation (JAMC), 2018, 2(1), 13-20
Journal of Applied Mathematics and Computation (JAMC), 2018, 2(1), 13-20 http://www.hillpublisher.org/journal/jamc ISSN Online:2576-0645 ISSN Print:2576-0653 Generation of Fractal Vessel Structure Functions
More informationPearling: Medical Image Segmentation with Pearl Strings
Pearling: Medical Image Segmentation with Pearl Strings Jarek Rossignac 1, Brian Whited 1, Greg Slabaugh 2, Tong Fang 2, Gozde Unal 2 1 Georgia Institute of Technology Graphics, Visualization, and Usability
More informationSemi-Automatic Detection of Cervical Vertebrae in X-ray Images Using Generalized Hough Transform
Semi-Automatic Detection of Cervical Vertebrae in X-ray Images Using Generalized Hough Transform Mohamed Amine LARHMAM, Saïd MAHMOUDI and Mohammed BENJELLOUN Faculty of Engineering, University of Mons,
More informationEdge Detection in Angiogram Images Using Modified Classical Image Processing Technique
Edge Detection in Angiogram Images Using Modified Classical Image Processing Technique S. Deepak Raj 1 Harisha D S 2 1,2 Asst. Prof, Dept Of ISE, Sai Vidya Institute of Technology, Bangalore, India Deepak
More informationThin Plate Spline Feature Point Matching for Organ Surfaces in Minimally Invasive Surgery Imaging
Thin Plate Spline Feature Point Matching for Organ Surfaces in Minimally Invasive Surgery Imaging Bingxiong Lin, Yu Sun and Xiaoning Qian University of South Florida, Tampa, FL., U.S.A. ABSTRACT Robust
More informationMore modelling. Ruth Aylett
More modelling Ruth Aylett Overview The great outdoors: fractals L-systems City Models Procedural generation of models How to make models.. Interactive modeling Human with a software modeling tool Scanning
More informationAutomatic Vascular Tree Formation Using the Mahalanobis Distance
Automatic Vascular Tree Formation Using the Mahalanobis Distance Julien Jomier, Vincent LeDigarcher, and Stephen R. Aylward Computer-Aided Diagnosis and Display Lab, Department of Radiology The University
More informationCalculating the Distance Map for Binary Sampled Data
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Calculating the Distance Map for Binary Sampled Data Sarah F. Frisken Gibson TR99-6 December 999 Abstract High quality rendering and physics-based
More informationAn Acquisition Geometry-Independent Calibration Tool for Industrial Computed Tomography
4th International Symposium on NDT in Aerospace 2012 - Tu.3.A.3 An Acquisition Geometry-Independent Calibration Tool for Industrial Computed Tomography Jonathan HESS *, Patrick KUEHNLEIN *, Steven OECKL
More information9 Tumor Blood Vessels: 3-D Tubular Network Analysis
Chap. 9 2015/8/20 14:02 page 1 1 9 Tumor Blood Vessels: 3-D Tubular Network Analysis Chap. 9 2015/8/20 14:02 page 2 2 Chap. 9 2015/8/20 14:02 page 3 3 Contents 9 Tumor Blood Vessels: 3-D Tubular Network
More informationIMPROVEMENT AND RESTORATION OF BIOMEDICAL IMAGES BLURRED BY BODY MOVEMENT
Proceedings 3rd Annual Conference IEEE/EBS Oct.5-8, 00, Istanbul, TURKEY IPROVEET AD RESTORATIO OF BIOEDICA IAGES BURRED BY BODY OVEET J. Hori, Y. Saitoh, T. Kiryu, K. Okamoto 3, K. Sakai 3 Faculty of
More informationMedical Image Segmentation
Medical Image Segmentation Xin Yang, HUST *Collaborated with UCLA Medical School and UCSB Segmentation to Contouring ROI Aorta & Kidney 3D Brain MR Image 3D Abdominal CT Image Liver & Spleen Caudate Nucleus
More informationVASCULAR TREE CHARACTERISTIC TABLE BUILDING FROM 3D MR BRAIN ANGIOGRAPHY IMAGES
VASCULAR TREE CHARACTERISTIC TABLE BUILDING FROM 3D MR BRAIN ANGIOGRAPHY IMAGES D.V. Sanko 1), A.V. Tuzikov 2), P.V. Vasiliev 2) 1) Department of Discrete Mathematics and Algorithmics, Belarusian State
More informationNIH Public Access Author Manuscript Proc Soc Photo Opt Instrum Eng. Author manuscript; available in PMC 2014 October 07.
NIH Public Access Author Manuscript Published in final edited form as: Proc Soc Photo Opt Instrum Eng. 2014 March 21; 9034: 903442. doi:10.1117/12.2042915. MRI Brain Tumor Segmentation and Necrosis Detection
More informationAnimating Plant Growth in L-System By Parametric Functional Symbols
Animating Plant Growth in L-System By Parametric Functional Symbols Somporn Chuai-aree, Suchada Siripant, and Chidchanok Lursinsap Advanced Virtual and Intelligent Computing Center (AVIC) Department of
More informationThe SIFT (Scale Invariant Feature
The SIFT (Scale Invariant Feature Transform) Detector and Descriptor developed by David Lowe University of British Columbia Initial paper ICCV 1999 Newer journal paper IJCV 2004 Review: Matt Brown s Canonical
More informationCS 4300 Computer Graphics. Prof. Harriet Fell Fall 2012 Lecture 28 November 8, 2012
CS 4300 Computer Graphics Prof. Harriet Fell Fall 2012 Lecture 28 November 8, 2012 1 Today s Topics Fractals Mandelbrot Set Julia Sets L-Systems 2 Fractals The term fractal was coined in 1975 by Benoît
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB NO. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationExtraction and recognition of the thoracic organs based on 3D CT images and its application
1 Extraction and recognition of the thoracic organs based on 3D CT images and its application Xiangrong Zhou, PhD a, Takeshi Hara, PhD b, Hiroshi Fujita, PhD b, Yoshihiro Ida, RT c, Kazuhiro Katada, MD
More information1238 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 29, NO. 6, JUNE /$ IEEE
1238 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 29, NO. 6, JUNE 2010 Spatio-Temporal Data Fusion for 3D+T Image Reconstruction in Cerebral Angiography Andrew D. Copeland, Rami S. Mangoubi*, Mukund N. Desai,
More informationAutomated Lesion Detection Methods for 2D and 3D Chest X-Ray Images
Automated Lesion Detection Methods for 2D and 3D Chest X-Ray Images Takeshi Hara, Hiroshi Fujita,Yongbum Lee, Hitoshi Yoshimura* and Shoji Kido** Department of Information Science, Gifu University Yanagido
More informationHuman Heart Coronary Arteries Segmentation
Human Heart Coronary Arteries Segmentation Qian Huang Wright State University, Computer Science Department Abstract The volume information extracted from computed tomography angiogram (CTA) datasets makes
More informationEngineering Drawings Recognition Using a Case-based Approach
Engineering Drawings Recognition Using a Case-based Approach Luo Yan Department of Computer Science City University of Hong Kong luoyan@cs.cityu.edu.hk Liu Wenyin Department of Computer Science City University
More informationUniversiteit Leiden Computer Science
Universiteit Leiden Computer Science Dynamically evolving L-system generated plant visualizations Name: Sander Ruijter S1415212 Date: 25/08/2016 1st supervisor: Dr. M.T.M. Emmerich 2nd supervisor: Dr.
More informationSCIENCE & TECHNOLOGY
Pertanika J. Sci. & Technol. 26 (1): 309-316 (2018) SCIENCE & TECHNOLOGY Journal homepage: http://www.pertanika.upm.edu.my/ Application of Active Contours Driven by Local Gaussian Distribution Fitting
More informationAutomatic Detection of Multiple Organs Using Convolutional Neural Networks
Automatic Detection of Multiple Organs Using Convolutional Neural Networks Elizabeth Cole University of Massachusetts Amherst Amherst, MA ekcole@umass.edu Sarfaraz Hussein University of Central Florida
More informationCSC Computer Graphics
// CSC. Computer Graphics Lecture Kasun@dscs.sjp.ac.lk Department of Computer Science University of Sri Jayewardanepura Polygon Filling Scan-Line Polygon Fill Algorithm Span Flood-Fill Algorithm Inside-outside
More informationA client-server architecture for semi-automatic segmentation of peripheral vessels in CTA data
A client-server architecture for semi-automatic segmentation of peripheral vessels in CTA data Poster No.: C-2174 Congress: ECR 2013 Type: Authors: Keywords: DOI: Scientific Exhibit A. Grünauer, E. Vuçini,
More informationBlood Vessel Visualization on CT Data
WDS'12 Proceedings of Contributed Papers, Part I, 88 93, 2012. ISBN 978-80-7378-224-5 MATFYZPRESS Blood Vessel Visualization on CT Data J. Dupej Charles University Prague, Faculty of Mathematics and Physics,
More informationComputer Science 173 Fall, Project 3: Lindenmayer Systems Due in class on Monday, September 28
Computer Science 173 all, 2015 Project 3: Lindenmayer Systems Due in class on Monday, September 28 Work on this project with a partner. Create a new directory for this project, copy the Makefile from the
More informationConstrained Diffusion Limited Aggregation in 3 Dimensions
Constrained Diffusion Limited Aggregation in 3 Dimensions Paul Bourke Swinburne University of Technology P. O. Box 218, Hawthorn Melbourne, Vic 3122, Australia. Email: pdb@swin.edu.au Abstract Diffusion
More informationThe Research of the Improved 3D L- System and Its Application in Plant Modeling
The Research of the Improved D L- System and Its Application in Plant Modeling Ke Guan Key Laboratory of Biomimetic Sensing and Advanced Robot Technology of Anhui Province Institute of Intelligent machines,
More informationFrom medical imaging to numerical simulations
From medical imaging to numerical simulations Christophe Prud Homme, Vincent Chabannes, Marcela Szopos, Alexandre Ancel, Julien Jomier To cite this version: Christophe Prud Homme, Vincent Chabannes, Marcela
More informationNorbert Schuff VA Medical Center and UCSF
Norbert Schuff Medical Center and UCSF Norbert.schuff@ucsf.edu Medical Imaging Informatics N.Schuff Course # 170.03 Slide 1/67 Objective Learn the principle segmentation techniques Understand the role
More informationAutomated segmentation methods for liver analysis in oncology applications
University of Szeged Department of Image Processing and Computer Graphics Automated segmentation methods for liver analysis in oncology applications Ph. D. Thesis László Ruskó Thesis Advisor Dr. Antal
More informationPattern Recognition Using Graph Theory
ISSN: 2278 0211 (Online) Pattern Recognition Using Graph Theory Aditya Doshi Department of Computer Science and Engineering, Vellore Institute of Technology, Vellore, India Manmohan Jangid Department of
More informationAPPLIED MECHANICS 2018
2018 APPLIED MECHANICS 2018 April 9-11, 2018, Myslovice, Czech republic IDENTIFICATION OF MULTI-COMPARTMENT DARCY FLOW MODEL MATERIAL PARAMETERS J. Brašnová 1, V. Lukeš 2, E. Rohan 3 Abstract: This papper
More informationAn Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy
An Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy Chenyang Xu 1, Siemens Corporate Research, Inc., Princeton, NJ, USA Xiaolei Huang,
More informationCT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION
CT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION Frank Dong, PhD, DABR Diagnostic Physicist, Imaging Institute Cleveland Clinic Foundation and Associate Professor of Radiology
More informationSTATISTICS AND ANALYSIS OF SHAPE
Control and Cybernetics vol. 36 (2007) No. 2 Book review: STATISTICS AND ANALYSIS OF SHAPE by H. Krim, A. Yezzi, Jr., eds. There are numerous definitions of a notion of shape of an object. These definitions
More informationL-Systems and Affine Transformations
L-Systems and Affine Transformations Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna http://www.moreno.marzolla.name/ Copyright 2014, Moreno Marzolla, Università di
More informationPlant synthesis Due: Monday, November 17
CMSC 23700 Introduction to Computer Graphics Project 2 Fall 2003 October 30 Plant synthesis Due: Monday, November 17 1 Introduction One of the xenobiology groups at NASA has finally found alien DNA hanging
More information4D Motion Modeling of the Coronary Arteries from CT Images for Robotic Assisted Minimally Invasive Surgery
4D Motion Modeling of the Coronary Arteries from CT Images for Robotic Assisted Minimally Invasive Surgery Dong Ping Zhang 1, Eddie Edwards 1,2, Lin Mei 1,2, Daniel Rueckert 1 1 Department of Computing,
More informationCreating an Automated Blood Vessel. Diameter Tracking Tool
Medical Biophysics 3970Z 6 Week Project: Creating an Automated Blood Vessel Diameter Tracking Tool Peter McLachlan - 250068036 April 2, 2013 Introduction In order to meet the demands of tissues the body
More informationHybrid Approach for MRI Human Head Scans Classification using HTT based SFTA Texture Feature Extraction Technique
Volume 118 No. 17 2018, 691-701 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu Hybrid Approach for MRI Human Head Scans Classification using HTT
More informationImage Segmentation and Registration
Image Segmentation and Registration Dr. Christine Tanner (tanner@vision.ee.ethz.ch) Computer Vision Laboratory, ETH Zürich Dr. Verena Kaynig, Machine Learning Laboratory, ETH Zürich Outline Segmentation
More informationEfficient Rendering of Glossy Reflection Using Graphics Hardware
Efficient Rendering of Glossy Reflection Using Graphics Hardware Yoshinori Dobashi Yuki Yamada Tsuyoshi Yamamoto Hokkaido University Kita-ku Kita 14, Nishi 9, Sapporo 060-0814, Japan Phone: +81.11.706.6530,
More informationImage Enhancement Techniques for Fingerprint Identification
March 2013 1 Image Enhancement Techniques for Fingerprint Identification Pankaj Deshmukh, Siraj Pathan, Riyaz Pathan Abstract The aim of this paper is to propose a new method in fingerprint enhancement
More information2 Michael E. Leventon and Sarah F. F. Gibson a b c d Fig. 1. (a, b) Two MR scans of a person's knee. Both images have high resolution in-plane, but ha
Model Generation from Multiple Volumes using Constrained Elastic SurfaceNets Michael E. Leventon and Sarah F. F. Gibson 1 MIT Artificial Intelligence Laboratory, Cambridge, MA 02139, USA leventon@ai.mit.edu
More informationLandmark-based 3D Elastic Registration of Pre- and Postoperative Liver CT Data
Landmark-based 3D Elastic Registration of Pre- and Postoperative Liver CT Data An Experimental Comparison Thomas Lange 1, Stefan Wörz 2, Karl Rohr 2, Peter M. Schlag 3 1 Experimental and Clinical Research
More informationVolume Illumination and Segmentation
Volume Illumination and Segmentation Computer Animation and Visualisation Lecture 13 Institute for Perception, Action & Behaviour School of Informatics Overview Volume illumination Segmentation Volume
More informationConformal flattening maps for the visualization of vessels
Conformal flattening maps for the visualization of vessels Lei Zhua, Steven Hakerb, and Allen Tannenbauma adept of Biomedical Engineering, Georgia Tech, Atlanta bdept of Radiology, Brigham and Women's
More informationMedical Images Analysis and Processing
Medical Images Analysis and Processing - 25642 Emad Course Introduction Course Information: Type: Graduated Credits: 3 Prerequisites: Digital Image Processing Course Introduction Reference(s): Insight
More informationEdge Detection for Dental X-ray Image Segmentation using Neural Network approach
Volume 1, No. 7, September 2012 ISSN 2278-1080 The International Journal of Computer Science & Applications (TIJCSA) RESEARCH PAPER Available Online at http://www.journalofcomputerscience.com/ Edge Detection
More informationAC : EDGE DETECTORS IN IMAGE PROCESSING
AC 11-79: EDGE DETECTORS IN IMAGE PROCESSING John Schmeelk, Virginia Commonwealth University/Qatar Dr. John Schmeelk is a Professor of mathematics at Virginia Commonwealth University teaching mathematics
More informationAdaptive Zoom Distance Measuring System of Camera Based on the Ranging of Binocular Vision
Adaptive Zoom Distance Measuring System of Camera Based on the Ranging of Binocular Vision Zhiyan Zhang 1, Wei Qian 1, Lei Pan 1 & Yanjun Li 1 1 University of Shanghai for Science and Technology, China
More informationCopyright 2009 Society of Photo Optical Instrumentation Engineers. This paper was published in Proceedings of SPIE, vol. 7260, Medical Imaging 2009:
Copyright 2009 Society of Photo Optical Instrumentation Engineers. This paper was published in Proceedings of SPIE, vol. 7260, Medical Imaging 2009: Computer Aided Diagnosis and is made available as an
More informationWhole Body MRI Intensity Standardization
Whole Body MRI Intensity Standardization Florian Jäger 1, László Nyúl 1, Bernd Frericks 2, Frank Wacker 2 and Joachim Hornegger 1 1 Institute of Pattern Recognition, University of Erlangen, {jaeger,nyul,hornegger}@informatik.uni-erlangen.de
More informationarxiv: v1 [cs.cv] 6 Jun 2017
Volume Calculation of CT lung Lesions based on Halton Low-discrepancy Sequences Liansheng Wang a, Shusheng Li a, and Shuo Li b a Department of Computer Science, Xiamen University, Xiamen, China b Dept.
More information2D Vessel Segmentation Using Local Adaptive Contrast Enhancement
2D Vessel Segmentation Using Local Adaptive Contrast Enhancement Dominik Schuldhaus 1,2, Martin Spiegel 1,2,3,4, Thomas Redel 3, Maria Polyanskaya 1,3, Tobias Struffert 2, Joachim Hornegger 1,4, Arnd Doerfler
More informationCS 231A Computer Vision (Winter 2014) Problem Set 3
CS 231A Computer Vision (Winter 2014) Problem Set 3 Due: Feb. 18 th, 2015 (11:59pm) 1 Single Object Recognition Via SIFT (45 points) In his 2004 SIFT paper, David Lowe demonstrates impressive object recognition
More informationContours & Implicit Modelling 4
Brief Recap Contouring & Implicit Modelling Contouring Implicit Functions Visualisation Lecture 8 lecture 6 Marching Cubes lecture 3 visualisation of a Quadric toby.breckon@ed.ac.uk Computer Vision Lab.
More information12.2 Plants. CS Dept, UK
1 12.2 Plants - modeling and animation of plants represents an interesting and challenging area - exhibit arbitrary complexity while possessing a constrained branching structure - grow from a single source
More informationDesign by Subdivision
Bridges 2010: Mathematics, Music, Art, Architecture, Culture Design by Subdivision Michael Hansmeyer Department for CAAD - Institute for Technology in Architecture Swiss Federal Institute of Technology
More informationCSE 167: Lecture #17: Procedural Modeling. Jürgen P. Schulze, Ph.D. University of California, San Diego Fall Quarter 2011
CSE 167: Introduction to Computer Graphics Lecture #17: Procedural Modeling Jürgen P. Schulze, Ph.D. University of California, San Diego Fall Quarter 2011 Announcements Important dates: Final project outline
More informationDecision Algorithm for Pool Using Fuzzy System
Decision Algorithm for Pool Using Fuzzy System S.C.Chua 1 W.C.Tan 2 E.K.Wong 3 V.C.Koo 4 1 Faculty of Engineering & Technology Tel: +60-06-252-3007, Fax: +60-06-231-6552, E-mail: scchua@mmu.edu.my 2 Faculty
More informationComputer-Aided Diagnosis in Abdominal and Cardiac Radiology Using Neural Networks
Computer-Aided Diagnosis in Abdominal and Cardiac Radiology Using Neural Networks Du-Yih Tsai, Masaru Sekiya and Yongbum Lee Department of Radiological Technology, School of Health Sciences, Faculty of
More informationRADIOMICS: potential role in the clinics and challenges
27 giugno 2018 Dipartimento di Fisica Università degli Studi di Milano RADIOMICS: potential role in the clinics and challenges Dr. Francesca Botta Medical Physicist Istituto Europeo di Oncologia (Milano)
More informationNovel evaluation method of low contrast resolution performance of dimensional X-ray CT
More Info at Open Access Database www.ndt.net/?id=18754 Novel evaluation method of low contrast resolution performance of dimensional X-ray CT Makoto Abe 1, Hiroyuki Fujimoto 1, Osamu Sato 1, Kazuya Matsuzaki
More informationTexture Segmentation by Windowed Projection
Texture Segmentation by Windowed Projection 1, 2 Fan-Chen Tseng, 2 Ching-Chi Hsu, 2 Chiou-Shann Fuh 1 Department of Electronic Engineering National I-Lan Institute of Technology e-mail : fctseng@ccmail.ilantech.edu.tw
More information