Edge Tracking of subjective contours in Biomedical Imaging
|
|
- Piers Hart
- 6 years ago
- Views:
Transcription
1 Edge Tracking of subjective contours in Biomedical Imaging Giovanni Garibotto(*), Valentina Garibotto (**) * Elsag spa, Genova, Italy ** Un. La Bicocca, Milano, Italy giovanni.garibotto@elsagdatamat.com Abstract The paper describes a method for computer-assisted edge tracking and following by local profile matching. It is used primarily to integrate and support manual selection when complex contours have to be found and identified. Such a requirement is quite common in biomedical applications, including a wide spectrum of image modalities, from histological samples to MRI images. The referred example takes part in the analysis of in vitro receptor auto-radiographic data, for the architectonic characterization of transmitter receptors within a defined cortical region. The practical validation we provide of the proposed tracking model reveals that it is a very promising approach, enabling accurate and fast processing of multiple images. 1. Introduction Edge detection and tracking represents one of the fundamental early processing steps in most imaging and computer vision applications. This is particularly important for tissue segmentation and classification in biomedical imaging. The goal of most imaging applications in the biomedical field is the quantification of different parameters that allow differential characterization of biological tissues and regions. The preliminary step for all such measurements consists in the definition of a region of interest and/or segmentation of the image in its constituting elements. Strategies to address this problem are various tissue classification approaches, whereby single pixel elements are assigned to a tissue class according to their intensities. However, there are no standard tools available yet for 2D images having a high spatial resolution and a relatively poor contrast. These images are commonly found when referring to many different biomedical fields, including histological samples and autoradiographies, electron-microscopic sections, and also in the emerging field of high resolution magnetic resonance imaging (MRI). In these application, the contour definition is currently performed manually, as described in Eickhoff et al., 2007, [1], Sätzler et al, 2002 [2], Walters et al., 2007, [3], and this necessarily increases the analysis time and work burden. In particular, the processing system described in this paper arises from a reference application as described in [1] for the analysis of the concentration of transmitter receptors in the human cerebral cortex. This method is based on the analysis of in vitro receptor auto-radiographic data, aimed to define regional cortical receptor distribution patterns. The cortical contours show a large variability in terms of contrast and are sometime difficult to be detected and identified also by a human operator. Quite often the edge contours must be extrapolated by exploiting visual continuity constraints (subjective contours). On the other hand a complete manual operation is very heavy to be performed, since there are extended areas where edge position can be reasonably predicted and tracked. As such, the best solution seems to be a semi-automatic edge feature tracking approach with human supervision (man in the loop). By the way, this requirement is quite common in many other context of computer vision where edge line features must be detected and fitted with continuous curves. This is the case of most biomedical imaging applications with high resolution images of radiographic records (ribs and other bone features) as well in the identification of relevant areas in the brain images. The technique described in this paper is limited to grey-level images, although it may be easily extended to deal with color pictures too. Main requirements of this edge contour representation are an accurate and regular description of the subjective profile, as well as a mono-pixel representation of the connected edge chain with 4-connectivity constraints, to allow an Proceedings of the 14 th International Conference on Image Analysis & Processing (ICIAP 2007), Sep., Modena, Italy
2 effective implementation of the following steps in the processing chain (region filling with Bresenham line interpolation). 2. Overview of the proposed approach Edge contours are described as a sequence of points that exhibit a similar orthogonal grey level profile. The first step of the process consists in the evaluation of a sequence of interpolation points the density of which is variable and proportional to the degree of curvature of the contour itself. This is a very general model of an edge profile to fit almost all possible situations from high-contrast edges up to highly smoothed transitions. HMI selection of 2 starting points Edge-profile model update HMI correction of last edge position & new edge point selection End of the process Edge-point prediction / tracking by correlation Stop criterion low-correlation Fig.1. flow diagram of the edge tracking process. HMI: Human Machine Interaction The proposed approach is a supervised edge detection and tracking where human intervention is required when the system is unable to decide with a sufficient level of confidence the appropriate position of the next interpolation edge point in the chain. The processing sequence is a combination of local edge chain fitting and tracking, whose length and extension depends on the continuity and regularity of the edge itself. The confidence criterion is based on the matching of the current edge profile with an edge model that is continuously updated along the contour chain. When the accuracy level falls below a suitable threshold, the model function is reset according to the new human selection. The resulting interpolated edge point sequence may be affected by some irregularities as possible intersections, loops, adjacency, that may prevent the following steps of the biomedical application, i.e. region filling and analysis. As such a further step of chain post-processing is required, to achieve a mono-pixel representation of the resulting edge chain sequence. 3. Edge following by local profile matching The starting point of the edge chain must be detected by the human operator in the image. A second manual point is also required, to define the starting local direction v e of the edge contour. Both edge points are used to build a model of the local edge profile. A local smoothing is applied by weighting the edge profile on a local neighborhood across the edge, to reduce noise effects and local irregularities. This regularization is particularly important when dealing with poor contrast images with a low level of signal to noise ratio. Then the automatic process starts, by predicting the next candidate position along the contour direction v e (with a step that is proportional to the local curvature of the edge chain). The new predicted position is corrected, across the edge profile, in the gradient direction v g orthogonal to the edge direction v e to get the maximum of correlation with respect to the current model. The process is repeated by computing the new edge direction v e between the last corrected point and the previous one. This method is summarized in the flow diagram of fig.1. An example of the edge profile matching process is shown in fig.2 where the new edge position s E along the gradient direction s, is shifted of the amount d to get the maximum value of the normalized cross covariance function: C(d) = [P(s-d)-P m ] [M(s)-M m ] ( [P(s-d)-P m ] 2 ) 1/2 ( [M(s)-M m ] 2 ) 1/2 (1) where P(s) is the current edge profile with P m as the average value, and M(s) is the current profile model with M m as the average value. The shift value d corresponds to the best estimation of the new edge position. The process is repeated until the covariance function C(d) has reached a value below a minimum level of acceptance (index of uncertainty). Hence the automatic search process stops and the human intervention is required to confirm the current position (or make a correction).
3 The edge profile model updating is a critical process with a risk of divergence when dealing with ambiguous contours due to the very local analysis of the automated process Fig.2 Matching between the edge profile P(s) ( ) with the current model M(s) ( ). re-projected available contour points and the red line is the results of editing and final edge chain construction as described in the next section. The editing tool provides a series of functions to the human operator, starting from the simple translation and rotation of the full edge chain to correct the main large differences of alignment and registration. Unfortunately a rigid roto-translation is not sufficient since most differences are due to local deformations. The proposed solution is an elastic matching process on the edge points before and after the one that has been corrected by the human operator. He/she may select any part of the contour that is considered too far from ideal position and shift it to the right place. The nearby contour points are modified in such a way to achieve the best fit to the ideal edge profile automatically, by weighting the amount of the edgeshift with the distance along the edge sequence. This is another reason to support the semiautomatic approach and the human supervision to correct any possible error of the system in such variable unconstrained situations. The result of this processing phase is a sequence of edge points with a variable displacement. 4. Contour editing by elastic matching Due to the complexity of the contour shape and variable contrast, it is possible to get a result that needs some additional corrections. In fact some of the detected points may be not properly located at the edge position and the overall shape of the contour sequence may be not fully satisfactory. An editing process is then implemented, to support the human operator in inserting new points, replacing wrong points, and improving the shape of the contour when necessary. Another motivation for the implementation of an interactive editing process is the possibility to transfer some of the already detected contour sequences to other similar images. Such a situation is quite common when dealing with temporal sequences (like heart video recordings) as well with 2D high resolution sectioning of 3D volumes. This is also the case in our reference applications where neighbouring sections are quite similar, considering the individual thickness of 20 µm of each section. Hence, the already available contours, from previous sections, may fit some part of the brain cortical shape, and it is more effective to re-shape the edge chain rather than rebuild it from scratch. An example of such full contour prediction is referred in fig.3. The light blue dots correspond to the Fig.3. Example of re-projection of a sequence of edge points (light-blue) onto neighbouring images and the matching results (red line). The cost function to be minimized is based on the usual cross covariance (1) of the edge profile computed for all neighbouring edge points, with continuity constraints on the first derivative of the interpolated edge contour chain (local curvature constraint). The first order approximation of the edge shift is further optimized for each point (in a local
4 neighbourhood of the corrected edge position) in the direction orthogonal to the edge. 5. Continuous edge chain representation The final step of the process consists in the reconstruction of the full edge chain sequence, starting from the list of the edge tracking result. The continuous chain sequence is obtained by an interpolation function that may be selected (parametric) by the user according to the desired level of curve smoothing. Unfortunately the resulting edge chain sequence may exhibit some local irregularities as multiple touching pixels especially in correspondence to curvature discontinuities or even some local loops. In our application [1] it is required to recover a mono-pixel chain code representation to simplify the following steps of region filling. To achieve such result it was necessary to develop a chain-parsing algorithm based on a local analysis of the 3x3 neighbourhood around the current edge point. The problem can be stated as follows: for any edge point of the chain P 0 =P[ind], only two edge connections are allowed in the local 3x3 neighbourhood as shown below P + P - P 0 P - P 0 P + Case 2. The triplet (P 0,P 1,P 2 ) of the last 3 consecutive points must be considered. With the assumption that the previous connection (P 1,P 2 ) is correct, the only irregular situation may occur when P 0 is found in the same position of P 2 corresponding to an opposite direction along the edge chain. In this case both P 0 and P 1 must be removed from the chain. Case 3: The current edge point P 0 should have only one busy connection in the local neighbourhood (corresponding to the previous edge point in the chain P 1 ) in both 8-connectivity and 4-connectivity models. Otherwise, if there is another edge point P k connected it is necessary to keep P 0, and remove all the intermediate points (from P 1 to P k ) since they would correspond to a local loop. It is worth to remind that these are the only possible situations that may occur and are recursively corrected on a local buffer of edge chain points to satisfy the 1- point connectivity constraint. The length of the buffer determines the length of the loops that can be removed to prevent any edge line crossing effect. 6. Interactive Human Interface A window application has been developed to improve human operation through a user-friendly interface. Fig 5 shows a screen-shot of the HMI (Human Machine Interaction) environment. It provides an overall view of the full image on the left, and the corresponding zoom view on the right, where human corrections and interactions are performed, through mouse interaction. Fig.4 Example of a mono-pixel chain test: 8-connectivity (left) and 4-connectivity (right) The removal of spurious edge points is quite simple, but extremely expensive if performed on the image map. A much more effective solution consists in exploiting the causal sequential representation of the edge chain. The proposed solution is a recursive process that takes into account all possible combination of position of the last edge points of the chain, namely points P 0 =P[ind], P 1 =P[ind-1], P 2 =P[ind-2] P k =P[indk]. This approach can be used both in the case of 4- connectivity and 8-connectivity constraints. A buffer of NP points is fed from the last position P 0 at any new interpolation point. The following situations must be checked, as follows: Case 1: the two last consecutive points (P 0, P 1 ) must satisfy regular continuity. This is achieved by construction and is always valid. Fig.5 Screen-shot of the HMI interactive environment On the right side of the screen the HMI interface contains the list of control keys to perform: Input-output image management
5 Selection of manual / automatic tracking of the edge points The editing menu of previously collected edge chain sequences. Start/stop of the tracking process On-line updating of the process parameters 7. Results The proposed approach has been used for the processing of a large number of images of variable contrast. The reference application has been described in [1] for the analysis of the concentration of transmitter receptors in the human cerebral cortex. cryostat microtome at 20 C. Alternating sections were histological stained, for cell bodies or incubated with sixteen different receptors of all classical transmitter systems. Due to the complexity of the task, the standard method that has been used so far consists in a fully manual segmentation of the contour profiles that are approximated by a series of segments of variable length. As a consequence, the accuracy of the contour representation is manually controlled by the human operator and the time required for each sample image is quite relevant (depending on the complexity of the scene). It is quite clear the great potential advantage of a computer assisted approach, considering the large size of the image samples (approximately in the range of 3500 x 5500 pixels). Such large size images require always a multi-scale control by the human operator, to manage the wide-field view of the overall edge contour shape (lower-resolution) together with a more precise analysis of the edge details (high-resolution requirement). Moreover the database of images to be processed is very huge (more than 700 images per hemisphere are currently available) with a proportional cost of the processing time. Fig.6 Image sample of the cortical shape This method is based on the analysis of in vitro receptor auto-radiographic data, aimed to define regional cortical receptor distribution patterns. Briefly, human hemispheres were obtained at autopsy from subjects with no record of neurological or psychiatric diseases and used for quantitative in vitro receptor auto-radiography. All subjects had given written consent before death and had been included in the body donor program of the Department of Anatomy at the University of Düsseldorf, Germany. Serial coronal sections (thickness 20 µm) were cut using a large-scale Fig.7. Contour tracking result on a small piece of the biomedical image sample The proposed system has been used for the edge reconstruction of the full sequence of brain slices. The obtained results have been found very satisfactory and have been judged comparable to the best quality achieved by full manual operations. In fact the human operator is often constrained to reduce the number of edge points due to time constraints as well as for a poor
6 perception of the local similarity of the edge profiles. The availability of a quantitative matching process is really the most relevant achievement of the proposed edge-tracking process, in terms of accuracy of edge chain reconstruction, more stable area filling within edge tracked contours, and better interpolation of smoothed edge curvatures. The improved accuracy of the automatic process has proved quite relevant for the success of the specific biomedical application. Another great advantage can be measured by the processing time that has been reduced by more than 60% by the use of the semi-automatic process, as compared to the traditional manual data acquisition scheme. 8. Conclusions The paper describes a method for computer-assisted edge tracking and following by local profile matching. It is a very specific solution used primarily to support human operation when complex contours have to be found and identified. It has been motivated by an application of reconstruction and representation of auto-radiographic data on brain sections. Anyway it may be successfully used also in a large variety of other applications including image processing and computer vision and many other biomedical imaging applications, including histological sections, digital radiographies and high resolution MRI imaging. The proposed solution is based on quite standard image processing and pattern recognition technology. Edge tracking is performed with an adaptive control of the local direction and a continuous update of the gradient profile model. Finally a quite significant contribution comes from the chain parsing algorithm that has been developed to achieve an error-free mono-pixel reconstruction of the connected edge chain sequence, an essential requirement for any further processing step of region filling and analysis. Acknowledgment The authors would like to thank the Institute of Medicine of the Jülich Forschungzentrum, and in particular Dr. S.B.Eickhoff, Dr. N. Palomero-Gallagher and Prof. K. Zilles, for the kind support in providing the sample data for the current research. Many thanks are also due to Luisa DeVena (OnAir srl) for the support in the implementation of the Human Interaction application. 9. References [1] Eickhoff S.B., et al. «Analysis of neurotransmitter receptor distribution patterns in the cerebral cortex», Neuroimage. 34(4): (2007). [2] Sätzler K, et al. Three-dimensional reconstruction of a calyx of held and its postsynaptic principal neuron in the medial nucleus of the trapezoid body. J. Neuroscience 22, (2002) [3] Walters N.B., et al Observer-Independent Analysis of High-Resolution MR Images of the Human Cerebral Cortex: In Vivo Delineation of Cortical Areas. Human Brain Mapping 28:1 8(2007) [4] Takashi Fujita et al. Semi-automatic tracking of the diaphragm contour in X-ray image sequences: preliminary results, Proceedings of the Digital Imaging Computing: Techniques and Applications (DICTA 2005). [5] Tae-Yong Kim, Jihun Park, Seong-Whan Lee, Object Boundary Edge Selection using normal direction derivatives of a Contour in a Complex Scene, Proc. Of the 17 th Int. Conf. on Pattern Recognition (ICPR 04), 2004.
Multimodal Elastic Image Matching
Research results based on my diploma thesis supervised by Prof. Witsch 2 and in cooperation with Prof. Mai 3. 1 February 22 nd 2011 1 Karlsruhe Institute of Technology (KIT) 2 Applied Mathematics Department,
More informationCombinatorial optimization and its applications in image Processing. Filip Malmberg
Combinatorial optimization and its applications in image Processing Filip Malmberg Part 1: Optimization in image processing Optimization in image processing Many image processing problems can be formulated
More informationObject Identification in Ultrasound Scans
Object Identification in Ultrasound Scans Wits University Dec 05, 2012 Roadmap Introduction to the problem Motivation Related Work Our approach Expected Results Introduction Nowadays, imaging devices like
More informationMEDICAL IMAGE ANALYSIS
SECOND EDITION MEDICAL IMAGE ANALYSIS ATAM P. DHAWAN g, A B IEEE Engineering in Medicine and Biology Society, Sponsor IEEE Press Series in Biomedical Engineering Metin Akay, Series Editor +IEEE IEEE PRESS
More informationBabu Madhav Institute of Information Technology Years Integrated M.Sc.(IT)(Semester - 7)
5 Years Integrated M.Sc.(IT)(Semester - 7) 060010707 Digital Image Processing UNIT 1 Introduction to Image Processing Q: 1 Answer in short. 1. What is digital image? 1. Define pixel or picture element?
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 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 informationAdaptive active contours (snakes) for the segmentation of complex structures in biological images
Adaptive active contours (snakes) for the segmentation of complex structures in biological images Philippe Andrey a and Thomas Boudier b a Analyse et Modélisation en Imagerie Biologique, Laboratoire Neurobiologie
More informationThe organization of the human cerebral cortex estimated by intrinsic functional connectivity
1 The organization of the human cerebral cortex estimated by intrinsic functional connectivity Journal: Journal of Neurophysiology Author: B. T. Thomas Yeo, et al Link: https://www.ncbi.nlm.nih.gov/pubmed/21653723
More informationEdge Detection (with a sidelight introduction to linear, associative operators). Images
Images (we will, eventually, come back to imaging geometry. But, now that we know how images come from the world, we will examine operations on images). Edge Detection (with a sidelight introduction to
More informationIntroduction to Medical Image Processing
Introduction to Medical Image Processing Δ Essential environments of a medical imaging system Subject Image Analysis Energy Imaging System Images Image Processing Feature Images Image processing may be
More informationDigital Volume Correlation for Materials Characterization
19 th World Conference on Non-Destructive Testing 2016 Digital Volume Correlation for Materials Characterization Enrico QUINTANA, Phillip REU, Edward JIMENEZ, Kyle THOMPSON, Sharlotte KRAMER Sandia National
More informationAnalysis of Functional MRI Timeseries Data Using Signal Processing Techniques
Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques Sea Chen Department of Biomedical Engineering Advisors: Dr. Charles A. Bouman and Dr. Mark J. Lowe S. Chen Final Exam October
More informationImproving Latent Fingerprint Matching Performance by Orientation Field Estimation using Localized Dictionaries
Available Online at www.ijcsmc.com International Journal of Computer Science and Mobile Computing A Monthly Journal of Computer Science and Information Technology IJCSMC, Vol. 3, Issue. 11, November 2014,
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 informationTexture April 17 th, 2018
Texture April 17 th, 2018 Yong Jae Lee UC Davis Announcements PS1 out today Due 5/2 nd, 11:59 pm start early! 2 Review: last time Edge detection: Filter for gradient Threshold gradient magnitude, thin
More informationEECS150 - Digital Design Lecture 14 FIFO 2 and SIFT. Recap and Outline
EECS150 - Digital Design Lecture 14 FIFO 2 and SIFT Oct. 15, 2013 Prof. Ronald Fearing Electrical Engineering and Computer Sciences University of California, Berkeley (slides courtesy of Prof. John Wawrzynek)
More information3D Surface Reconstruction of the Brain based on Level Set Method
3D Surface Reconstruction of the Brain based on Level Set Method Shijun Tang, Bill P. Buckles, and Kamesh Namuduri Department of Computer Science & Engineering Department of Electrical Engineering University
More informationSpatio-Temporal Registration of Biomedical Images by Computational Methods
Spatio-Temporal Registration of Biomedical Images by Computational Methods Francisco P. M. Oliveira, João Manuel R. S. Tavares tavares@fe.up.pt, www.fe.up.pt/~tavares Outline 1. Introduction 2. Spatial
More informationColor Image Segmentation Editor Based on the Integration of Edge-Linking, Region Labeling and Deformable Model
This paper appears in: IEEE International Conference on Systems, Man and Cybernetics, 1999 Color Image Segmentation Editor Based on the Integration of Edge-Linking, Region Labeling and Deformable Model
More informationSURFACE RECONSTRUCTION OF EX-VIVO HUMAN V1 THROUGH IDENTIFICATION OF THE STRIA OF GENNARI USING MRI AT 7T
SURFACE RECONSTRUCTION OF EX-VIVO HUMAN V1 THROUGH IDENTIFICATION OF THE STRIA OF GENNARI USING MRI AT 7T Oliver P. Hinds 1, Jonathan R. Polimeni 2, Megan L. Blackwell 3, Christopher J. Wiggins 3, Graham
More informationComputational Medical Imaging Analysis
Computational Medical Imaging Analysis Chapter 1: Introduction to Imaging Science Jun Zhang Laboratory for Computational Medical Imaging & Data Analysis Department of Computer Science University of Kentucky
More informationCHAPTER-1 INTRODUCTION
CHAPTER-1 INTRODUCTION 1.1 Fuzzy concept, digital image processing and application in medicine With the advancement of digital computers, it has become easy to store large amount of data and carry out
More informationMachine Learning for Medical Image Analysis. A. Criminisi
Machine Learning for Medical Image Analysis A. Criminisi Overview Introduction to machine learning Decision forests Applications in medical image analysis Anatomy localization in CT Scans Spine Detection
More informationLast update: May 4, Vision. CMSC 421: Chapter 24. CMSC 421: Chapter 24 1
Last update: May 4, 200 Vision CMSC 42: Chapter 24 CMSC 42: Chapter 24 Outline Perception generally Image formation Early vision 2D D Object recognition CMSC 42: Chapter 24 2 Perception generally Stimulus
More informationA Survey of Light Source Detection Methods
A Survey of Light Source Detection Methods Nathan Funk University of Alberta Mini-Project for CMPUT 603 November 30, 2003 Abstract This paper provides an overview of the most prominent techniques for light
More informationModule 7 VIDEO CODING AND MOTION ESTIMATION
Module 7 VIDEO CODING AND MOTION ESTIMATION Lesson 20 Basic Building Blocks & Temporal Redundancy Instructional Objectives At the end of this lesson, the students should be able to: 1. Name at least five
More informationMotion Estimation for Video Coding Standards
Motion Estimation for Video Coding Standards Prof. Ja-Ling Wu Department of Computer Science and Information Engineering National Taiwan University Introduction of Motion Estimation The goal of video compression
More informationCAP 5415 Computer Vision Fall 2012
CAP 5415 Computer Vision Fall 01 Dr. Mubarak Shah Univ. of Central Florida Office 47-F HEC Lecture-5 SIFT: David Lowe, UBC SIFT - Key Point Extraction Stands for scale invariant feature transform Patented
More informationCOMP 102: Computers and Computing
COMP 102: Computers and Computing Lecture 23: Computer Vision Instructor: Kaleem Siddiqi (siddiqi@cim.mcgill.ca) Class web page: www.cim.mcgill.ca/~siddiqi/102.html What is computer vision? Broadly speaking,
More informationImage Acquisition Systems
Image Acquisition Systems Goals and Terminology Conventional Radiography Axial Tomography Computer Axial Tomography (CAT) Magnetic Resonance Imaging (MRI) PET, SPECT Ultrasound Microscopy Imaging ITCS
More informationDeformable Registration Using Scale Space Keypoints
Deformable Registration Using Scale Space Keypoints Mehdi Moradi a, Purang Abolmaesoumi a,b and Parvin Mousavi a a School of Computing, Queen s University, Kingston, Ontario, Canada K7L 3N6; b Department
More informationFPGA IMPLEMENTATION FOR REAL TIME SOBEL EDGE DETECTOR BLOCK USING 3-LINE BUFFERS
FPGA IMPLEMENTATION FOR REAL TIME SOBEL EDGE DETECTOR BLOCK USING 3-LINE BUFFERS 1 RONNIE O. SERFA JUAN, 2 CHAN SU PARK, 3 HI SEOK KIM, 4 HYEONG WOO CHA 1,2,3,4 CheongJu University E-maul: 1 engr_serfs@yahoo.com,
More information8/3/2017. Contour Assessment for Quality Assurance and Data Mining. Objective. Outline. Tom Purdie, PhD, MCCPM
Contour Assessment for Quality Assurance and Data Mining Tom Purdie, PhD, MCCPM Objective Understand the state-of-the-art in contour assessment for quality assurance including data mining-based techniques
More informationThe Processing of Form Documents
The Processing of Form Documents David S. Doermann and Azriel Rosenfeld Document Processing Group, Center for Automation Research University of Maryland, College Park 20742 email: doermann@cfar.umd.edu,
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 informationSupplementary Information
Supplementary Information Magnetic resonance imaging reveals functional anatomy and biomechanics of a living dragon tree Linnea Hesse 1,2,*, Tom Masselter 1,2,3, Jochen Leupold 4, Nils Spengler 5, Thomas
More informationGeneration of Hulls Encompassing Neuronal Pathways Based on Tetrahedralization and 3D Alpha Shapes
Generation of Hulls Encompassing Neuronal Pathways Based on Tetrahedralization and 3D Alpha Shapes Dorit Merhof 1,2, Martin Meister 1, Ezgi Bingöl 1, Peter Hastreiter 1,2, Christopher Nimsky 2,3, Günther
More informationAnnouncements. Edge Detection. An Isotropic Gaussian. Filters are templates. Assignment 2 on tracking due this Friday Midterm: Tuesday, May 3.
Announcements Edge Detection Introduction to Computer Vision CSE 152 Lecture 9 Assignment 2 on tracking due this Friday Midterm: Tuesday, May 3. Reading from textbook An Isotropic Gaussian The picture
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 informationNEW CONCEPT FOR JOINT DISPARITY ESTIMATION AND SEGMENTATION FOR REAL-TIME VIDEO PROCESSING
NEW CONCEPT FOR JOINT DISPARITY ESTIMATION AND SEGMENTATION FOR REAL-TIME VIDEO PROCESSING Nicole Atzpadin 1, Serap Askar, Peter Kauff, Oliver Schreer Fraunhofer Institut für Nachrichtentechnik, Heinrich-Hertz-Institut,
More informationTopics to be Covered in the Rest of the Semester. CSci 4968 and 6270 Computational Vision Lecture 15 Overview of Remainder of the Semester
Topics to be Covered in the Rest of the Semester CSci 4968 and 6270 Computational Vision Lecture 15 Overview of Remainder of the Semester Charles Stewart Department of Computer Science Rensselaer Polytechnic
More informationCHAPTER 1 INTRODUCTION
CHAPTER 1 INTRODUCTION 1.1 Introduction Pattern recognition is a set of mathematical, statistical and heuristic techniques used in executing `man-like' tasks on computers. Pattern recognition plays an
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 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 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 informationComputer Vision and Pattern Recognition in Homeland Security Applications 1
Computer Vision and Pattern Recognition in Homeland Security Applications 1 Giovanni B. Garibotto Elsag spa Genova, Italy giovanni.garibotto@elsagdatamat.com Abstract. The tutorial will summarize the status
More informationChapter 9 Conclusions
Chapter 9 Conclusions This dissertation has described a new method for using local medial properties of shape to identify and measure anatomical structures. A bottom up approach based on image properties
More informationBiomedical Image Processing
Biomedical Image Processing Jason Thong Gabriel Grant 1 2 Motivation from the Medical Perspective MRI, CT and other biomedical imaging devices were designed to assist doctors in their diagnosis and treatment
More informationA Generic Lie Group Model for Computer Vision
A Generic Lie Group Model for Computer Vision Within this research track we follow a generic Lie group approach to computer vision based on recent physiological research on how the primary visual cortex
More informationGENERAL AUTOMATED FLAW DETECTION SCHEME FOR NDE X-RAY IMAGES
GENERAL AUTOMATED FLAW DETECTION SCHEME FOR NDE X-RAY IMAGES Karl W. Ulmer and John P. Basart Center for Nondestructive Evaluation Department of Electrical and Computer Engineering Iowa State University
More informationVisualization of cross sectional data for morphogenetic studies
Visualization of cross sectional data for morphogenetic studies G. Brunnett, M. Vančo, Ch. Haller, S. Washausen, H.-J. Kuhn, W. Knabe Abstract: We report on a visualization system that has been implemented
More informationAutomatic segmentation of the cortical grey and white matter in MRI using a Region Growing approach based on anatomical knowledge
Automatic segmentation of the cortical grey and white matter in MRI using a Region Growing approach based on anatomical knowledge Christian Wasserthal 1, Karin Engel 1, Karsten Rink 1 und André Brechmann
More informationComputer Vision I - Filtering and Feature detection
Computer Vision I - Filtering and Feature detection Carsten Rother 30/10/2015 Computer Vision I: Basics of Image Processing Roadmap: Basics of Digital Image Processing Computer Vision I: Basics of Image
More informationQuantifying Three-Dimensional Deformations of Migrating Fibroblasts
45 Chapter 4 Quantifying Three-Dimensional Deformations of Migrating Fibroblasts This chapter presents the full-field displacements and tractions of 3T3 fibroblast cells during migration on polyacrylamide
More informationBasic principles of MR image analysis. Basic principles of MR image analysis. Basic principles of MR image analysis
Basic principles of MR image analysis Basic principles of MR image analysis Julien Milles Leiden University Medical Center Terminology of fmri Brain extraction Registration Linear registration Non-linear
More informationRange Imaging Through Triangulation. Range Imaging Through Triangulation. Range Imaging Through Triangulation. Range Imaging Through Triangulation
Obviously, this is a very slow process and not suitable for dynamic scenes. To speed things up, we can use a laser that projects a vertical line of light onto the scene. This laser rotates around its vertical
More informationA Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation
, pp.162-167 http://dx.doi.org/10.14257/astl.2016.138.33 A Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation Liqiang Hu, Chaofeng He Shijiazhuang Tiedao University,
More informationNew Approaches for EEG Source Localization and Dipole Moment Estimation. Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine
New Approaches for EEG Source Localization and Dipole Moment Estimation Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine Outline Motivation why EEG? Mathematical Model equivalent
More informationJournal of Asian Scientific Research FEATURES COMPOSITION FOR PROFICIENT AND REAL TIME RETRIEVAL IN CBIR SYSTEM. Tohid Sedghi
Journal of Asian Scientific Research, 013, 3(1):68-74 Journal of Asian Scientific Research journal homepage: http://aessweb.com/journal-detail.php?id=5003 FEATURES COMPOSTON FOR PROFCENT AND REAL TME RETREVAL
More informationBiomedical Image Analysis based on Computational Registration Methods. João Manuel R. S. Tavares
Biomedical Image Analysis based on Computational Registration Methods João Manuel R. S. Tavares tavares@fe.up.pt, www.fe.up.pt/~tavares Outline 1. Introduction 2. Methods a) Spatial Registration of (2D
More informationClassification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging
1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant
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 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 informationContent Based Image Retrieval Using Color Quantizes, EDBTC and LBP Features
Content Based Image Retrieval Using Color Quantizes, EDBTC and LBP Features 1 Kum Sharanamma, 2 Krishnapriya Sharma 1,2 SIR MVIT Abstract- To describe the image features the Local binary pattern (LBP)
More informationCS334: Digital Imaging and Multimedia Edges and Contours. Ahmed Elgammal Dept. of Computer Science Rutgers University
CS334: Digital Imaging and Multimedia Edges and Contours Ahmed Elgammal Dept. of Computer Science Rutgers University Outlines What makes an edge? Gradient-based edge detection Edge Operators From Edges
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Review of Motion Modelling and Estimation Introduction to Motion Modelling & Estimation Forward Motion Backward Motion Block Motion Estimation Motion
More informationSemi-Automatic Segmentation of the Patellar Cartilage in MRI
Semi-Automatic Segmentation of the Patellar Cartilage in MRI Lorenz König 1, Martin Groher 1, Andreas Keil 1, Christian Glaser 2, Maximilian Reiser 2, Nassir Navab 1 1 Chair for Computer Aided Medical
More informationPROSTATE CANCER DETECTION USING LABEL IMAGE CONSTRAINED MULTIATLAS SELECTION
PROSTATE CANCER DETECTION USING LABEL IMAGE CONSTRAINED MULTIATLAS SELECTION Ms. Vaibhavi Nandkumar Jagtap 1, Mr. Santosh D. Kale 2 1 PG Scholar, 2 Assistant Professor, Department of Electronics and Telecommunication,
More informationRegistration Techniques
EMBO Practical Course on Light Sheet Microscopy Junior-Prof. Dr. Olaf Ronneberger Computer Science Department and BIOSS Centre for Biological Signalling Studies University of Freiburg Germany O. Ronneberger,
More informationVector Field Visualisation
Vector Field Visualisation Computer Animation and Visualization Lecture 14 Institute for Perception, Action & Behaviour School of Informatics Visualising Vectors Examples of vector data: meteorological
More informationA MORPHOLOGY-BASED FILTER STRUCTURE FOR EDGE-ENHANCING SMOOTHING
Proceedings of the 1994 IEEE International Conference on Image Processing (ICIP-94), pp. 530-534. (Austin, Texas, 13-16 November 1994.) A MORPHOLOGY-BASED FILTER STRUCTURE FOR EDGE-ENHANCING SMOOTHING
More informationImplementation of Advanced Image Guided Radiation Therapy
Image Acquisition Course Outline Principles, characteristics& applications of the available modalities Image Processing in the T x room Image guided treatment delivery What can / can t we do in the room
More information3/27/2012 WHY SPECT / CT? SPECT / CT Basic Principles. Advantages of SPECT. Advantages of CT. Dr John C. Dickson, Principal Physicist UCLH
3/27/212 Advantages of SPECT SPECT / CT Basic Principles Dr John C. Dickson, Principal Physicist UCLH Institute of Nuclear Medicine, University College London Hospitals and University College London john.dickson@uclh.nhs.uk
More informationScene-Based Segmentation of Multiple Muscles from MRI in MITK
Scene-Based Segmentation of Multiple Muscles from MRI in MITK Yan Geng 1, Sebastian Ullrich 2, Oliver Grottke 3, Rolf Rossaint 3, Torsten Kuhlen 2, Thomas M. Deserno 1 1 Department of Medical Informatics,
More informationTexture April 14 th, 2015
Texture April 14 th, 2015 Yong Jae Lee UC Davis Announcements PS1 out today due 4/29 th, 11:59 pm start early! 2 Review: last time Edge detection: Filter for gradient Threshold gradient magnitude, thin
More informationVolumetric Analysis of the Heart from Tagged-MRI. Introduction & Background
Volumetric Analysis of the Heart from Tagged-MRI Dimitris Metaxas Center for Computational Biomedicine, Imaging and Modeling (CBIM) Rutgers University, New Brunswick, NJ Collaboration with Dr. Leon Axel,
More informationBasic fmri Design and Analysis. Preprocessing
Basic fmri Design and Analysis Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial filtering
More informationAn Evaluation of Volumetric Interest Points
An Evaluation of Volumetric Interest Points Tsz-Ho YU Oliver WOODFORD Roberto CIPOLLA Machine Intelligence Lab Department of Engineering, University of Cambridge About this project We conducted the first
More informationCS 229 Final Project Report Learning to Decode Cognitive States of Rat using Functional Magnetic Resonance Imaging Time Series
CS 229 Final Project Report Learning to Decode Cognitive States of Rat using Functional Magnetic Resonance Imaging Time Series Jingyuan Chen //Department of Electrical Engineering, cjy2010@stanford.edu//
More informationStructural MRI analysis
Structural MRI analysis volumetry and voxel-based morphometry cortical thickness measurements structural covariance network mapping Boris Bernhardt, PhD Department of Social Neuroscience, MPI-CBS bernhardt@cbs.mpg.de
More informationCS534: Introduction to Computer Vision Edges and Contours. Ahmed Elgammal Dept. of Computer Science Rutgers University
CS534: Introduction to Computer Vision Edges and Contours Ahmed Elgammal Dept. of Computer Science Rutgers University Outlines What makes an edge? Gradient-based edge detection Edge Operators Laplacian
More informationFunctional MRI in Clinical Research and Practice Preprocessing
Functional MRI in Clinical Research and Practice Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization
More informationStructured Light II. Guido Gerig CS 6320, Spring (thanks: slides Prof. S. Narasimhan, CMU, Marc Pollefeys, UNC)
Structured Light II Guido Gerig CS 6320, Spring 2013 (thanks: slides Prof. S. Narasimhan, CMU, Marc Pollefeys, UNC) http://www.cs.cmu.edu/afs/cs/academic/class/15385- s06/lectures/ppts/lec-17.ppt Variant
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 10 Segmentation 14/02/27 http://www.ee.unlv.edu/~b1morris/ecg782/
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 informationA Model-Independent, Multi-Image Approach to MR Inhomogeneity Correction
Tina Memo No. 2007-003 Published in Proc. MIUA 2007 A Model-Independent, Multi-Image Approach to MR Inhomogeneity Correction P. A. Bromiley and N.A. Thacker Last updated 13 / 4 / 2007 Imaging Science and
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 informationLocating ego-centers in depth for hippocampal place cells
204 5th Joint Symposium on Neural Computation Proceedings UCSD (1998) Locating ego-centers in depth for hippocampal place cells Kechen Zhang,' Terrence J. Sejeowski112 & Bruce L. ~cnau~hton~ 'Howard Hughes
More informationEPI Data Are Acquired Serially. EPI Data Are Acquired Serially 10/23/2011. Functional Connectivity Preprocessing. fmri Preprocessing
Functional Connectivity Preprocessing Geometric distortion Head motion Geometric distortion Head motion EPI Data Are Acquired Serially EPI Data Are Acquired Serially descending 1 EPI Data Are Acquired
More informationRT_Image v0.2β User s Guide
RT_Image v0.2β User s Guide RT_Image is a three-dimensional image display and analysis suite developed in IDL (ITT, Boulder, CO). It offers a range of flexible tools for the visualization and quantitation
More informationExperiments with Edge Detection using One-dimensional Surface Fitting
Experiments with Edge Detection using One-dimensional Surface Fitting Gabor Terei, Jorge Luis Nunes e Silva Brito The Ohio State University, Department of Geodetic Science and Surveying 1958 Neil Avenue,
More informationLow Cost Motion Capture
Low Cost Motion Capture R. Budiman M. Bennamoun D.Q. Huynh School of Computer Science and Software Engineering The University of Western Australia Crawley WA 6009 AUSTRALIA Email: budimr01@tartarus.uwa.edu.au,
More informationSu et al. Shape Descriptors - III
Su et al. Shape Descriptors - III Siddhartha Chaudhuri http://www.cse.iitb.ac.in/~cs749 Funkhouser; Feng, Liu, Gong Recap Global A shape descriptor is a set of numbers that describes a shape in a way that
More informationVolumetry of hypothalamic substructures by multimodal morphological image registration
Volumetry of hypothalamic substructures by multimodal morphological image registration Dominik Löchel 14.09.2011 Institute for Applied and Numerical Mathematics KIT University of the State of Baden-Württemberg
More informationVolume rendering for interactive 3-d segmentation
Volume rendering for interactive 3-d segmentation Klaus D. Toennies a, Claus Derz b a Dept. Neuroradiology, Inst. Diagn. Radiology, Inselspital Bern, CH-3010 Berne, Switzerland b FG Computer Graphics,
More informationSegmentation and Grouping
Segmentation and Grouping How and what do we see? Fundamental Problems ' Focus of attention, or grouping ' What subsets of pixels do we consider as possible objects? ' All connected subsets? ' Representation
More informationA Statistical Consistency Check for the Space Carving Algorithm.
A Statistical Consistency Check for the Space Carving Algorithm. A. Broadhurst and R. Cipolla Dept. of Engineering, Univ. of Cambridge, Cambridge, CB2 1PZ aeb29 cipolla @eng.cam.ac.uk Abstract This paper
More informationDigital Image Processing Fundamentals
Ioannis Pitas Digital Image Processing Fundamentals Chapter 7 Shape Description Answers to the Chapter Questions Thessaloniki 1998 Chapter 7: Shape description 7.1 Introduction 1. Why is invariance to
More informationModern Medical Image Analysis 8DC00 Exam
Parts of answers are inside square brackets [... ]. These parts are optional. Answers can be written in Dutch or in English, as you prefer. You can use drawings and diagrams to support your textual answers.
More information