Fast Retinal Vessel Analysis

Michael Krause, Ralph Maria Alles, Bernhard Burgeth, Joachim Weickert
Institute of Microelectronics, Building A5.1, Saarland University, 66041 Saarbrucken, Germany
Technical Report No. 320, Department of Mathematics, Saarland University, Saarbrucken, Germany, 2012

   title={Fast Retinal Vessel Analysis},

   author={Krause, Michael and Alles, Ralph Maria and Burgeth, Bernhard and Weickert, Joachim},



Download Download (PDF)   View View   Source Source   



Due to the increasing availability of so called "Non-Mydriatic" cameras, digital imaging has become a very important part of the ophthalmologist’s work. This has created large databases of retinal images. It would be desirable to have a fast image processing tool that allows to analyse such databases in a short time, and to process the image in situ while the patient is examined. In this paper we contribute such a system for fast retinal image analysis. While it achieves a comparable quality as state-of-the-art methods, it differs from most of them by the fact that it is extremely fast. Retinal blood vessels are assumed to be line-like structures and can therefore be enhanced via convolution with suitable, elongated kernels. For this task we use the second derivative of the local Radon kernel. It is rotated at different angles and adapts via a maximisation procedure to the directions of the vessels. We combine smoothing along vessel directions with contrast enhancement across them. Afterwards our algorithm detects vessels as connected structures with very few interruptions. A subsequent skeletonisation allows a higher-level description of the vessel tree. In order to end up with a very fast system, we combine efficient algorithms for numerical integration, numerical differentiation and interpolation, and we propose an automatic parameter selection strategy. Our convolution kernels are precomputed and stored into cached constant memory. All essential components in our algorithms are intrinsically parallel, and the resulting system is implemented on GPUs using the CUDA programming language. Our qualitative evaluations with the publicly available DRIVE database and our own clinical database shows that the system achieves competitive performance. We demonstrate that it is possible to process images of size 4288 x 2848 pixels in 1.2 seconds on a NVIDIA Geforce GTX680, including the time for reading from and writing to disk.
VN:F [1.9.22_1171]
Rating: 0.0/5 (0 votes cast)

* * *

* * *

Follow us on Twitter

HGPU group

1584 peoples are following HGPU @twitter

Like us on Facebook

HGPU group

302 people like HGPU on Facebook

* * *

Free GPU computing nodes at hgpu.org

Registered users can now run their OpenCL application at hgpu.org. We provide 1 minute of computer time per each run on two nodes with two AMD and one nVidia graphics processing units, correspondingly. There are no restrictions on the number of starts.

The platforms are

Node 1
  • GPU device 0: nVidia GeForce GTX 560 Ti 2GB, 822MHz
  • GPU device 1: AMD/ATI Radeon HD 6970 2GB, 880MHz
  • CPU: AMD Phenom II X6 @ 2.8GHz 1055T
  • RAM: 12GB
  • OS: OpenSUSE 13.1
  • SDK: nVidia CUDA Toolkit 6.5.14, AMD APP SDK 3.0
Node 2
  • GPU device 0: AMD/ATI Radeon HD 7970 3GB, 1000MHz
  • GPU device 1: AMD/ATI Radeon HD 5870 2GB, 850MHz
  • CPU: Intel Core i7-2600 @ 3.4GHz
  • RAM: 16GB
  • OS: OpenSUSE 12.3
  • SDK: AMD APP SDK 3.0

Completed OpenCL project should be uploaded via User dashboard (see instructions and example there), compilation and execution terminal output logs will be provided to the user.

The information send to hgpu.org will be treated according to our Privacy Policy

HGPU group © 2010-2015 hgpu.org

All rights belong to the respective authors

Contact us: