Enhancement of Optical Coherence Tomography Images of the Retina by Normalization and Fusion
|
|
- Hilary Gaines
- 6 years ago
- Views:
Transcription
1 86 F. ROSSANT, F. AMIEL, T. EA, C. MARAVER MARTINEZ, M. PAQUES, ENHANCEMENT OF OPTICAL COHERENCE... Enhancement o Optical Coherence Tomography Images o the Retina by Normalization and Fusion Florence ROSSANT, Frederic AMIEL, Thomas EA, C. MARAVER MARTINEZ, Michel PÂQUES 2 Institut Supérieur d Electronique de Paris, ISEP, 2 rue d Assas Paris, France 2 Centre d'investigation Clinique - Centre Hospitalier National d'ophtalmologie des Quinze-Vingts, 28 rue de Charenton 7502 Paris, France lorence.rossant rederic.amiel michel.paques@gmail.com Abstract. This paper describes an image processing method applied to Optical Coherence Tomography (OCT) images o the retina. The aim is to achieve improved OCT images rom the usion o sequential OCT scans obtained at identical retinal locations. The method is based on the normalization o the acquired images and their usion. As a result, a noise reduction and an image enhancement are reached. Thans to the resulting improvement in retinal imaging, clinical specialists are able to evaluate more eiciently eyes pathologies and anomalies. This paper presents the proposed method and gives some evaluation results. Keywords Optical coherence tomography, denoising, iltering, correlation, usion.. Introduction Optical coherence tomography (OCT) is a cross-sectional imaging technique allowing micrometric-scale resolution o retinal structures. It enables reliable demonstration o changes in overall retinal thicness, detection o luid in and behind the neurosensory retina, and identiication o the retinal nerve iber and photoreceptor layers. It is routinely used or diagnosis o retinal diseases. The OCT principle is to obtain a histology image o the retina by measuring the relection intensity o a low coherence inra-red light beam on the retina. The ovea o the eye, region where the analysis is done on, has a length between 0.5 and mm. The cut taes place on a dimension o 3 to 0 mm and on 52 points. The best resolution is, thereore, obtained or a length o 3 mm. The depth o 20 μm is measured by the tomograph on 024 points []. The obtained relection coeicients are displayed as an image (52x024 pixels) in alse colors. Although these images give an excellent idea about the retina, the excessive presence o noise maes the intervening layers o the neurosensory retina be only vaguely discernible. A primary goal o our study was to obtain a better characterization o the outer retina, that is, the interace between the retinal pigment epithelium and the photoreceptors. With ultra high resolution instruments, based on titanium sapphire lasers [3], urther improvements would be possible. However, this instrument is not commercialized yet. This is the reason why the research o a method that increases the signal to noise ratio and enhance the image has become an issue to study [4]. Image usion generally results in a better quality, since it leads to an increased signal-to-noise ratio. In OCT imaging o the retina, a series o scans (I I N ) are acquired rom identical undus locations. A irst pre-processing step is applied to each one, in order to normalize the images. This step is called alignment" in what ollows. It allows a better readability o each image taen separately, and is necessary to achieve the usion. The usion step itsel is based on an average, but other unctions could be considered. The bloc diagram o the global system is shown in Fig.. I Max intensity I 2 relexion Normalisation A 2... layer process :... I N detection alignment A N Fig.. Global system. A Fusion I used We have explored this idea using data rom the conventional StratusOCT instrument [], [2]. We have wored directly with the relection coeicients instead o the processed alse color image. Indeed, the relection coeicients provide much more inormation (about 500 dierent values), whereas the given output images are quantized on about 245 dierent values. So, the inormation obtained by the StratusOCT tomography, to be processed by our method, is a matrix o relection coeicients (024 rows per 52 columns). The developed algorithm treats this matrix as an image, where each value corresponds to a pixel. The paper is organized as ollows. Sections 2 and 3 present the normalization method, which is based on the detection o the maximal intensity relection layer (Section
2 RADIOENGINEERING, VOL. 7, NO. 4, DECEMBER ), ollowed by an image transormation, that results in the alignment o this layer (Section 3). Section 4 describes the usion method. Finally, some evaluation results are analyzed in Section 5. The last section presents a conclusion and introduces uture wor. 2. Maximum Intensity Relection Layer Detection 2. Region o Interest Fig. 2 shows two scans o retinal OCT images, I and I 2, captured consecutively rom the same patient. The images are displayed in alse colors. They have been also cropped and resized or a better readability o the igure. Our coordinate system is deined by the origin at the upper let corner o the image, the vertical x-axis and the horizontal y-axis. It is obviously not possible to use directly these images: some distortions can be observed between both, since they do not exactly represent the same retina cut (because o ocular movements). So, an essential preprocessing step consists in detecting a common region o interest, the maximal intensity relection layer (Fig. 2), and aligning all images, so that the pixels belonging to this layer are located at the same vertical position (Fig. 5). Thus, the retinal structure will be at the same position in all images, and the usion can be achieved. Another important goal o this transormation is that it allows a better clinical interpretation o the retina OCT images. (c) Fig. 2. Two scans taen rom the same patient (I, I 2 ), (c) a zoom o igure. The maximal intensity relection layer (in red) is the common region o interest used as reerence to achieve the normalization. The major diiculty or this detection is the high noise presence (Fig. 2(c)). The dominating noise source in OCT images is usually specle noise arising rom intererence between coherent waves bacscattered rom nearby scatters in the measuring volume o the retina. Some o the classic methods used or denoising are averaging each pixel with its neighbors (spatial low-pass iltering), applying a median iltering or a low pass iltering in the requency domain [5]. These three techniques are used in dierent steps all along the method that is described below. 2.2 Detection Method The aim is to ind an internal line on the retinal layer. In the OCT image, this layer corresponds to a roughly horizontal set o pixels taing the highest values (displayed in red). Because o the low signal-to-noise ratio, such simple methods lie looing or the maximal pixel value column by column, do not wor. In addition, nown algorithms about boundary detection, as or example active contours [6], do not wor properly either, as the noisy pixels bloc the evolution o the contour. Consequently, more complex treatments are proposed. Firstly, a pixel belonging to the region o interest, called "internal point", is ound with certainty. This pixel is then used to initialize an iterative algorithm that allows to deduce the median line o the maximal intensity layer. This algorithm uses correlation inormation between adjacent pixel columns, in order to reach a more robust and representative result. It is also applied on images that have been beorehand low-pass iltered, in order to reduce the noise. The bloc diagram o the detection method, applied to each source image I n ( n N ), is shown in Fig. 3. I n Internal Point Detection LP FIR Filter J n (x 0,y 0 ) Median Line Detection ( ( x,, 0 y < j ) opt 52 Fig. 3. Maximal relection intensity layer detection method. Searching or an internal point: The image I n is irst smoothed, using a large convolution mas (0*40 coeicients equal to /400). The pixel reaching the highest value, at the (x 0,y 0 ) coordinates, is certain to belong to the maximal intensity relection layer, and consequently, it is taen as the internal point. In our experimentations, it is proved that a rectangular mas leads to better results, since it is more representative o the globally horizontal structure we want to detect. Low-pass iltering: The initial image I n is iltered in the requency domain, using a FIR ilter o order 6, with a normalized cut-o requency equal to As a result, the signal to noise ratio is increased while the boundaries o the region o interest are not signiicantly delocalized. The output image, denoted by J n, is passed to the median line detection algorithm. Median line detection: Starting rom the detected internal point (x 0,y 0 ), the algorithm deduces iteratively, column by column, the median line o the maximal intensity relection layer. Let consider the next column at the horizontal coordinate y=y 0 +. The pixels o this column are low-pass iltered, using a convolution mas o K = 2 coeicients, all equal to /K: S ( = K K / 2 J n x +, = K / 2 (. ()
3 88 F. ROSSANT, F. AMIEL, T. EA, C. MARAVER MARTINEZ, M. PAQUES, ENHANCEMENT OF OPTICAL COHERENCE... The highest outputs o this ilter are located inside the maximal intensity relection layer. In order to avoid alse detections and smooth the searched median line, a recursive low-pass ilter is also applied to the S( coeicients, along the y coordinate. Let denote by C( the output o this second ilter (Eq. 2). The decision at column y is taen by retaining the x ( opt vertical coordinate corresponding to the maximum output value (Eq. 3): C( = ( α ) S( + αc( y ), (2) max x ( { C( } C( x, =. (3) opt In this method, the S( coeicients are continuously integrated to provide the decision at the column y. The parameter α o this recursive ilter expresses the relative importance between the current local results, and the previous ilter outputs. In our experiments, we use α=0.8. This value results rom a compromise: with a greater value, the median line is not accurately ollowed; with a smaller value, the algorithm is too sensitive to bright noise pixels. The same algorithm is applied in the other direction, or decreasing y coordinates, starting also rom the irst detected internal point (x 0,y 0 ). The result obtained ater applying the method to the image I (Fig. 2) is presented below, superimposed on the iltered image J. The detected median line (in dar blue) ollows accurately the maximal intensity relection layer. estimated by averaging the pixels on a small area around the median line. Fig. 5 shows the superior boundary ound with this method. Fig. 5. Detection o the superior boundary o the maximum relection intensity layer, and alignment against this line. Finally, the alignment is carried out. It consists on a column by column simple shit (Fig. 5). The resulting images are denoted by A n in what ollows. Fig. 4. Maximal intensity relection layer detection. The dar line represents the median position o this layer. 3. Alignment The aim o the alignment is to normalize the dierent source images, so that they can be used, and also to improve their readability or the clinical diagnostic. The alignment is not carried out against the median internal line ound by the method previously described, but against the superior boundary o the maximum intensity relection layer. This boundary can be deduced through a region growing algorithm, whose seed is the internal median line. In order to suppress the noise, a median ilter is applied beorehand column per column. Its size is equal to the typical width o the maximum intensity relection layer. The region growing stops when a pixel value is below a threshold T that is dynamically set to the hal o the mean pixel value M o the region o interest (T = 0.5M). M is 4. Fusion The aim o this section is to obtain an enhanced image in order to satisy the medical interest. The idea is to generate rom several aligned images a new one, with a higher signal-to-noise ratio. Beore using, the correct superposition between the images has to be ound. Translation, rotation and homothety transormation could be considered. But only translation is studied, because a rotation will give a dierent cut and thus a dierent image. Homothety transormation is not signiicant because the patient does not move his head during snapshots. The images are correlated two by two, in order to ind the optimal superposition between all couples. The higher is the correlation score, the higher is the similarity between both images. So, the maximum correlation value corresponds to the optimal translation shit. The maximum correlation scores obtained over all the aligned images allow also to choose a reerence image, rom which all the translations will be achieved, and to reject images that are not enough similar to it. The reerence image A R is chosen as ollows: it is the one that leads, in average, to the two highest correlation scores with the other images. Images that get a correlation score below S=0.5 (threshold experimen-
4 RADIOENGINEERING, VOL. 7, NO. 4, DECEMBER tally ound) are considered to be erroneous ones, and will not be included in the usion. A A 2... A N Inter correlation Reerence Image A R Other accepted images A Rejected C( AR, A ) > S images A C ( A,A ) S R A R A v u Fig. 6. Fusion process. Fusion Dierent usion unctions can be envisaged. In our irst experiments, we use a simple average unction that applies spatially on each retained image A and also between all these images: A ( x = = i= j= I used K, A ( x + u + i, y + v + j) (4) 9K where K (K N), the number o accepted images, indexed by, including the reerence image A R, (u,v ) the horizontal and vertical shits between the image A and the reerence image A R. The results can be still improved, by excluding rom the average the pixels that are liely to be noise pixels. A ( represents the mean pixel value, calculated over the K retained images, on a small 3x3 neighborhood around (. In the same way we compute an estimation o the standard deviation σ (. The pixels o the neighborhood that dier too much rom the mean value A ( are excluded rom the averaging (5). The resulting image I used is consequently improved, compared to the mean image A. A ( A ( x + u + i, y + v + j) > 0.5σ ( x, A ( x + u + i, y + v j) excluded rom the usion. (5) + and the variance calculated on the used image (6) provides an estimation o the signal to noise ratio improvement. This study has been conducted only on images o healthy retinas, otherwise the three regions cannot be all correctly deined. Tab. indicates the number o images used or the usion and the gains (in db) obtained or each region. G ( i) σ 2 K = ( i) ( i 0log, σ m = σ ( i) 2 K = σ m ( i) Number o images (K) Tab.. The signal to noise ratio gain measured on three layers, or 7 cases o healthy eye. The signal to noise ratio is increased in all cases, with a gain between 3 and 2 db. The greater is the number o used images, the greater is the SNR gain. The Fovea and the ONL regions are clearly denoised (gain over 6 db) while the layer has been signiicantly smoothed. ONL G () (db) G (2) (db) ) G (3) (db) ONL (6) An example o usion is shown in Fig. 8. Four images A were used. This process results in a higher signal to noise ratio and an enhanced contrast (see next section). 5. Results Our study has been done over 0 patients, healthy or presenting dierent pathologies. Three to iteen images were captured per patient. The proposed method has been evaluated and quantiied based on the ollowing procedure. We have considered three homogenous areas where the grey levels should be almost constant without the presence o noise: the ovea (), the ganglion cell layer and the inner plexiorm layer (), and the outer nuclear layer (ONL) (Fig. 7). We have estimated the noise power by calculating the image variance on these regions (manually delimitated). Let us denote by σ (i) the standard deviation o the region i in the image, and by σ (i) the corresponding measure in the used image. Then, the ratio between the mean variance calculated on the original images ONL Fig. 7. Fusion o 0 images: one o the original images, the result image. Images are represented in grey levels and normalized between 0 and. The results were also submitted to clinical specialists rom the Quinze-Vingts hospital (Paris). Doctors assessed in all cases that the normalization and usion method leads to an enhanced image that maes the clinical interpretation easier and more accurate. It is clearly apparent rom Fig. 7 and 8 that the contrast between the retinal layers is improved. In particular, the relective layer attributed to the external
5 90 F. ROSSANT, F. AMIEL, T. EA, C. MARAVER MARTINEZ, M. PAQUES, ENHANCEMENT OF OPTICAL COHERENCE... limiting membrane, that is, the rontier between the inner and outer segment o the photoreceptors, becomes clear, while it was slightly apparent in non-processed images. (c) Fig. 8. Fusion o 4 aligned images (c)(d) to provide the enhanced output image (e). In the example shown in Fig. 9, the interruption o the photoreceptor line is evident. Alterations o the intraretinal structure are also observed, consisting o collection o oedema within the retina. (c) (d) (e) (d) Fig. 9. Fusion o 6 aligned images (three o them are represented in igures (c)) to provide the enhanced output image (d). 6. Conclusion This paper describes a method o optical coherence tomography imaging enhancement o the retinal area o the eye. The methods and algorithms proposed in this article, i.e. relection layer detection, image alignment and usion, have proved to increase considerably the quality and the signal to noise ratio. In act, usion o a collection o image OCT scans rom the same retinal area enhances enough the quality o imaging to reveal new details o the retina. To date, only total retinal thicness was considered in most studies o OCT; alterations o retinal layering, which is an indicator o retinal diseases, could not be reliably deined. Averaging OCT scans thus provides additional inormation about the intraretinal structures. This is in agreement with the results presented by [4]. The higher deinition o the retinal layering oers the opportunity to better deine the alterations o retinal structures. Such improved precision contributes to a better diagnosis o retinal alterations and hence to the cause o visual impairment. Future wor will include investigation o new usion techniques, in order to improve the processing speed and the quality o the result. The automatic detection o the dierent physiological layers o the retina, in order to calculate some measures useul or the clinical diagnostic, is also a subject o urther study. Reerences [] StratusOCT tomography, Zeiss documentation, Frame/2BF7095D5578B4D F9AD. [2] Carl Zeiss Meditec AG, Goeschwitzer Str 5-52, Jena, Germany, [3] KENT, C. Taing OCT technology to the next level. Review o Ophthalmology, 2006, vol. 3:06, Issue 6/7/2006. [4] SANDER, B., LARSEN, M., THRANE, L., HOUGAARD, J. L., JØRGENSEN, T. M. Enhanced optical coherence tomography imaging by multiple scan averaging. Br. J. Ophthalmol. 2005; vol. 89; p doi:0.36/bjo , December [5] GONZALES, R. C., WOODS, R. E. Digital Image Processing. New Jersey: Prentice Hall, [6] SETHIAN, J. A. Level Set Method and Fast Matching Methods. Cambridge University Press, 999.
Blood Vessel Tree Reconstruction in Retinal OCT Data
Blood Vessel Tree Reconstruction in Retinal OCT Data Gazárek J, Kolář R, Jan J, Odstrčilík J, Taševský P Department of Biomedical Engineering, FEEC, Brno University of Technology xgazar03@stud.feec.vutbr.cz
More informationNoise Removal from ECG Signal and Performance Analysis Using Different Filter
International Journal o Innovative Research in Electronics and Communication (IJIREC) Volume. 1, Issue 2, May 214, PP.32-39 ISSN 2349-442 (Print) & ISSN 2349-45 (Online) www.arcjournal.org Noise Removal
More informationDARK CURRENT ELIMINATION IN CHARGED COUPLE DEVICES
DARK CURRENT ELIMINATION IN CHARGED COUPLE DEVICES L. Kňazovická, J. Švihlík Department o Computing and Control Engineering, ICT Prague Abstract Charged Couple Devices can be ound all around us. They are
More informationPublic Watermarking Surviving General Scaling and Cropping: An Application for Print-and-Scan Process
Public Watermarking urviving General caling and ropping: An Application or Print-and-can Process hing-yung Lin Department o Electrical Engineering olumbia University 500 W0th t. #3 New York, NY 007, UA
More informationMeasuring the Speed of Light
Physics Teaching Laboratory Measuring the peed o Light Introduction: The goal o this experiment is to measure the speed o light, c. The experiment relies on the technique o heterodyning, a very useul tool
More informationLength-Sensing OpLevs for KAGRA
Length-Sensing OpLevs or KAGRA Simon Zeidler Basics Length-Sensing Optical Levers are needed in order to measure the shit o mirrors along the optical path o the incident main-laser beam with time. The
More informationImage Enhancement II: Neighborhood Operations
Image Enhancement II: Neighborhood Operations Image Enhancement:Spatial Filtering Operation Idea: Use a mask to alter piel values according to local operation Aim: De)-Emphasize some spatial requencies
More informationImage Enhancement. Image Enhancement
SPATIAL FILTERING g h * h g FREQUENCY DOMAIN FILTERING G H. F F H G Copright RMR / RDL - 999. PEE53 - Processamento Digital de Imagens LOW PASS FILTERING attenuate or eliminate high-requenc components
More informationDetection and direction-finding of spread spectrum signals using correlation and narrowband interference rejection
Detection and direction-inding o spread spectrum signals using correlation and narrowband intererence rejection Ulrika Ahnström,2,JohanFalk,3, Peter Händel,3, Maria Wikström Department o Electronic Warare
More informationPLANNING AND DESIGN OF FRONT-END FILTERS
PLANNING AND DESIGN OF FRONT-END FILTERS AND DIPLEXERS FOR RADIO LINK APPLICATIONS Kjetil Folgerø and Jan Kocba Nera Networks AS, N-52 Bergen, NORWAY. Email: ko@nera.no, jko@nera.no Abstract High capacity
More informationCHAPTER 4 LOCATING THE CENTER OF THE OPTIC DISC AND MACULA
90 CHAPTER 4 LOCATING THE CENTER OF THE OPTIC DISC AND MACULA The objective in this chapter is to locate the centre and boundary of OD and macula in retinal images. In Diabetic Retinopathy, location of
More informationArtefact Characterisation for JPEG and JPEG 2000 Image Codecs: Edge Blur and Ringing
I'.NCINEER- Vol. XXXX, No. 3, pp. 25-3, 27
More informationOptimized Color Transforms for Image Demosaicing
International Journal O Computational Engineering Research (ceronlinecom) Vol Issue Optimized Color ransorms or Image Demosaicing Evgen Gershiov Department o Electrical Engineering, Ort Braude Academic
More informationStudy on Visual Style of Twist Silk Knitted Fabrics
2nd International Conerence on Inormation, Electronics and Computer (ICIEAC 2014) Study on Visual Style o Twist Sil Knitted Fabrics Yuan Fang, Ya-qun Zhang Department o Textile Engineering, College o Textile
More informationA MATLAB Model of Hybrid Active Filter Based on SVPWM Technique
International Journal o Electrical Engineering. ISSN 0974-2158 olume 5, Number 5 (2012), pp. 557-569 International Research Publication House http://www.irphouse.com A MATLAB Model o Hybrid Active Filter
More informationNew metallic mesh designing with high electromagnetic shielding
MATEC Web o Conerences 189, 01003 (018) MEAMT 018 https://doi.org/10.1051/mateccon/01818901003 New metallic mesh designing with high electromagnetic shielding Longjia Qiu 1,,*, Li Li 1,, Zhieng Pan 1,,
More informationTime distributed update of the NLMS algorithm coefficients for Acoustic Echo Cancellers
Time distributed update o the NLMS algorithm coeicients or Acoustic Echo Cancellers Fotis E. Andritsopoulos, Yannis M. Mitsos, Christos N. Charopoulos, Gregory A. Doumenis, Constantin N. Papaodysseus Abstract
More informationFatigue Life Assessment Using Signal Processing Techniques
Fatigue Lie Assessment Using Signal Processing Techniques S. ABDULLAH 1, M. Z. NUAWI, C. K. E. NIZWAN, A. ZAHARIM, Z. M. NOPIAH Engineering Faculty, Universiti Kebangsaan Malaysia 43600 UKM Bangi, Selangor,
More informationAN APPROXIMATION-WEIGHTED DETAIL CONTRAST ENHANCEMENT FILTER FOR LESION DETECTION ON MAMMOGRAMS
AN APPROXIMATION-WEIGHTED DETAIL CONTRAST ENHANCEMENT FILTER FOR LESION DETECTION ON MAMMOGRAMS Zhuangzhi Yan, Xuan He, Shupeng Liu, and Donghui Lu Department of Biomedical Engineering, Shanghai University,
More informationTraditional Analog Modulation Techniques
Chapter 5 Traditional Analog Modulation Techniques Mikael Olosson 2002 2007 Modulation techniques are mainly used to transmit inormation in a given requency band. The reason or that may be that the channel
More informationECE5984 Orthogonal Frequency Division Multiplexing and Related Technologies Fall Mohamed Essam Khedr. Channel Estimation
ECE5984 Orthogonal Frequency Division Multiplexing and Related Technologies Fall 2007 Mohamed Essam Khedr Channel Estimation Matlab Assignment # Thursday 4 October 2007 Develop an OFDM system with the
More informationGoing beyond the surface of your retina OCT-HS100 OPTICAL COHERENCE TOMOGRAPHY
Going beyond the surface of your retina OCT-HS100 OPTICAL COHERENCE TOMOGRAPHY Automatic functions make examinations short and simple. Perform the examination with only two simple mouse clicks! 1. START
More informationGoing beyond the surface of your retina OCT-HS100 OPTICAL COHERENCE TOMOGRAPHY
Going beyond the surface of your retina OCT-HS100 OPTICAL COHERENCE TOMOGRAPHY Full Auto OCT High specifications in a very compact design Automatic functions make examinations short and simple. Perform
More informationBlood Vessel Tracking Technique for Optic Nerve Localisation for Field 1-3 Color Fundus Images
Blood Tracing Technique for Optic Nerve Localisation for Field 1-3 Color Fundus Images Hwee Keong Lam, Opas Chutatape School of Electrical and Electronic Engineering Nanyang Technological University, Nanyang
More informationSECTION I - CHAPTER 2 DIGITAL IMAGING PROCESSING CONCEPTS
RADT 3463 - COMPUTERIZED IMAGING Section I: Chapter 2 RADT 3463 Computerized Imaging 1 SECTION I - CHAPTER 2 DIGITAL IMAGING PROCESSING CONCEPTS RADT 3463 COMPUTERIZED IMAGING Section I: Chapter 2 RADT
More informationVISULAS Trion. Treatment flexibility to the power of three. Multicolor Photocoagulation Laser
VISULAS Trion Treatment flexibility to the power of three Multicolor Photocoagulation Laser Carl Zeiss: A pioneer in retinal therapy For many years, Carl Zeiss has fostered a culture of highest precision,
More informationPerformance of LTE Linear MIMO Detectors: Achievable Data Rates and Complexity
Perormance o LTE Linear MIMO Detectors: Achievable Data Rates and Complexity Dragan Samardzija, Milos Pilipovic, Dusica Marijan, Jaroslav Farkas, Miodrag Temerinac University o Novi Sad Novi Sad, Serbia
More informationFurther developments on gear transmission monitoring
Further developments on gear transmission monitoring Niola V., Quaremba G., Avagliano V. Department o Mechanical Engineering or Energetics University o Naples Federico II Via Claudio 21, 80125, Napoli,
More informationFourier Domain (Spectral) OCT OCT: HISTORY. Could OCT be a Game Maker OCT in Optometric Practice: A THE TECHNOLOGY BEHIND OCT
Could OCT be a Game Maker OCT in Optometric Practice: A Hands On Guide Murray Fingeret, OD Nick Rumney, MSCOptom Fourier Domain (Spectral) OCT New imaging method greatly improves resolution and speed of
More informationThe TRC-NW8F Plus: As a multi-function retinal camera, the TRC- NW8F Plus captures color, red free, fluorescein
The TRC-NW8F Plus: By Dr. Beth Carlock, OD Medical Writer Color Retinal Imaging, Fundus Auto-Fluorescence with exclusive Spaide* Filters and Optional Fluorescein Angiography in One Single Instrument W
More informationAN EFFICIENT SET OF FEATURES FOR PULSE REPETITION INTERVAL MODULATION RECOGNITION
AN EFFICIENT SET OF FEATURES FOR PULSE REPETITION INTERVAL MODULATION RECOGNITION J-P. Kauppi, K.S. Martikainen Patria Aviation Oy, Naulakatu 3, 33100 Tampere, Finland, ax +358204692696 jukka-pekka.kauppi@patria.i,
More informationTIME-FREQUENCY ANALYSIS OF NON-STATIONARY THREE PHASE SIGNALS. Z. Leonowicz T. Lobos
Copyright IFAC 15th Triennial World Congress, Barcelona, Spain TIME-FREQUENCY ANALYSIS OF NON-STATIONARY THREE PHASE SIGNALS Z. Leonowicz T. Lobos Wroclaw University o Technology Pl. Grunwaldzki 13, 537
More informationBode Plot based Auto-Tuning Enhanced Solution for High Performance Servo Drives
Bode lot based Auto-Tuning Enhanced Solution or High erormance Servo Drives. O. Krah Danaher otion GmbH Wachholder Str. 4-4 4489 Düsseldor Germany Email: j.krah@danaher-motion.de Tel. +49 3 9979 133 Fax.
More informationIntroduction to OFDM. Characteristics of OFDM (Orthogonal Frequency Division Multiplexing)
Introduction to OFDM Characteristics o OFDM (Orthogonal Frequency Division Multiplexing Parallel data transmission with very long symbol duration - Robust under multi-path channels Transormation o a requency-selective
More informationSinusoidal signal. Arbitrary signal. Periodic rectangular pulse. Sampling function. Sampled sinusoidal signal. Sampled arbitrary signal
Techniques o Physics Worksheet 4 Digital Signal Processing 1 Introduction to Digital Signal Processing The ield o digital signal processing (DSP) is concerned with the processing o signals that have been
More informationFovea and Optic Disc Detection in Retinal Images with Visible Lesions
Fovea and Optic Disc Detection in Retinal Images with Visible Lesions José Pinão 1, Carlos Manta Oliveira 2 1 University of Coimbra, Palácio dos Grilos, Rua da Ilha, 3000-214 Coimbra, Portugal 2 Critical
More informationSpeckle Noise Reduction for the Enhancement of Retinal Layers in Optical Coherence Tomography Images
Iranian Journal of Medical Physics Vol. 12, No. 3, Summer 2015, 167-177 Received: February 25, 2015; Accepted: July 8, 2015 Original Article Speckle Noise Reduction for the Enhancement of Retinal Layers
More informationSoftware Defined Radio Forum Contribution
Committee: Technical Sotware Deined Radio Forum Contribution Title: VITA-49 Drat Speciication Appendices Source Lee Pucker SDR Forum 604-828-9846 Lee.Pucker@sdrorum.org Date: 7 March 2007 Distribution:
More informationOptimizing Reception Performance of new UWB Pulse shape over Multipath Channel using MMSE Adaptive Algorithm
IOSR Journal o Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vol. 05, Issue 01 (January. 2015), V1 PP 44-57 www.iosrjen.org Optimizing Reception Perormance o new UWB Pulse shape over Multipath
More informationISSUE: April Fig. 1. Simplified block diagram of power supply voltage loop.
ISSUE: April 200 Why Struggle with Loop ompensation? by Michael O Loughlin, Texas Instruments, Dallas, TX In the power supply design industry, engineers sometimes have trouble compensating the control
More informationCOMPENSATION OF CURRENT TRANSFORMERS BY MEANS OF FIELD PROGRAMMABLE GATE ARRAY
METROLOGY AD MEASUREMET SYSTEMS Index 330930, ISS 0860-89 www.metrology.pg.gda.pl COMPESATIO OF CURRET TRASFORMERS BY MEAS OF FIELD PROGRAMMABLE GATE ARRAY Daniele Gallo, Carmine Landi, Mario Luiso Seconda
More informationWith the proposed technique, those two problems will be overcome. reduction is to eliminate the specific harmonics, which are the lowest orders.
CHAPTER 3 OPTIMIZED HARMONIC TEPPED-WAVEFORM TECHNIQUE (OHW The obective o the proposed optimized harmonic stepped-waveorm technique is to reduce, as much as possible, the harmonic distortion in the load
More informationSEG/San Antonio 2007 Annual Meeting. Summary. Morlet wavelet transform
Xiaogui Miao*, CGGVeritas, Calgary, Canada, Xiao-gui_miao@cggveritas.com Dragana Todorovic-Marinic and Tyler Klatt, Encana, Calgary Canada Summary Most geologic changes have a seismic response but sometimes
More informationARRAY PROCESSING FOR INTERSECTING CIRCLE RETRIEVAL
16th European Signal Processing Conference (EUSIPCO 28), Lausanne, Switzerland, August 25-29, 28, copyright by EURASIP ARRAY PROCESSING FOR INTERSECTING CIRCLE RETRIEVAL Julien Marot and Salah Bourennane
More informationOCT mini-symposium. Presenters. Donald Miller, Indiana Univ. Joseph Izatt, Duke Univ. Thomas Milner, Univ. of Texas at Austin Jay Wei, Zeiss Meditec
OCT mini-symposium Presenters Donald Miller, Indiana Univ. Joseph Izatt, Duke Univ. Thomas Milner, Univ. of Texas at Austin Jay Wei, Zeiss Meditec Starlight, eyebright Canberra Times, Australia Combining
More informationLousy Processing Increases Energy Efficiency in Massive MIMO Systems
1 Lousy Processing Increases Energy Eiciency in Massive MIMO Systems Sara Gunnarsson, Micaela Bortas, Yanxiang Huang, Cheng-Ming Chen, Liesbet Van der Perre and Ove Edors Department o EIT, Lund University,
More informationLock-In Amplifiers SR510 and SR530 Analog lock-in amplifiers
Lock-In Ampliiers SR510 and SR530 Analog lock-in ampliiers SR510/SR530 Lock-In Ampliiers 0.5 Hz to 100 khz requency range Current and voltage inputs Up to 80 db dynamic reserve Tracking band-pass and line
More informationOn the Impact of Fading and Inter-piconet Interference on Bluetooth Performance
On the Impact o Fading and Inter-piconet Intererence on Bluetooth Perormance Andrea Zanella Dept. o Inormation Engineering University o Padova, Padova, Italy zanella@dei.unipd.it Andrea Tonello Bell Labs,
More informationECC419 IMAGE PROCESSING
ECC419 IMAGE PROCESSING INTRODUCTION Image Processing Image processing is a subclass of signal processing concerned specifically with pictures. Digital Image Processing, process digital images by means
More informationDetermination of Pitch Range Based on Onset and Offset Analysis in Modulation Frequency Domain
Determination o Pitch Range Based on Onset and Oset Analysis in Modulation Frequency Domain A. Mahmoodzadeh Speech Proc. Research Lab ECE Dept. Yazd University Yazd, Iran H. R. Abutalebi Speech Proc. Research
More informationWhite-light interferometry, Hilbert transform, and noise
White-light interferometry, Hilbert transform, and noise Pavel Pavlíček *a, Václav Michálek a a Institute of Physics of Academy of Science of the Czech Republic, Joint Laboratory of Optics, 17. listopadu
More informationIntroduction Approach Work Performed and Results
Algorithm for Morphological Cancer Detection Carmalyn Lubawy Melissa Skala ECE 533 Fall 2004 Project Introduction Over half of all human cancers occur in stratified squamous epithelia. Approximately one
More information(b) By measuring the image height for various image distances (adjusted by sliding the tubes together or apart) a relationship can be determined.
(c) The image is smaller, upright, virtual, ann the same side o the lens. Applying Inquiry Skills 7. (a) (b) By measuring the image height or various image distances (adjusted by sliding the tubes together
More informationAutomated Detection of Early Lung Cancer and Tuberculosis Based on X- Ray Image Analysis
Proceedings of the 6th WSEAS International Conference on Signal, Speech and Image Processing, Lisbon, Portugal, September 22-24, 2006 110 Automated Detection of Early Lung Cancer and Tuberculosis Based
More information287. The Transient behavior of rails used in electromagnetic railguns: numerical investigations at constant loading velocities
287. The Transient behavior o rails used in electromagnetic railguns: numerical investigations at constant loading velocities L. Tumonis 1, a, R. Kačianauskas 1,b, A. Kačeniauskas 2,c, M. Schneider 3,d
More informationRadionuclide Imaging MII Single Photon Emission Computed Tomography (SPECT)
Radionuclide Imaging MII 3073 Single Photon Emission Computed Tomography (SPECT) Single Photon Emission Computed Tomography (SPECT) The successful application of computer algorithms to x-ray imaging in
More informationEstimation and Compensation of IQ-Imbalances in Direct Down Converters
Estimation and Compensation o IQ-Imbalances in irect own Converters NRES PSCHT, THOMS BITZER and THOMS BOHN lcatel SEL G, Holderaeckerstrasse 35, 7499 Stuttgart GERMNY bstract: - In this paper, a new method
More informationSPEECH PARAMETERIZATION FOR AUTOMATIC SPEECH RECOGNITION IN NOISY CONDITIONS
SPEECH PARAMETERIZATION FOR AUTOMATIC SPEECH RECOGNITION IN NOISY CONDITIONS Bojana Gajić Department o Telecommunications, Norwegian University o Science and Technology 7491 Trondheim, Norway gajic@tele.ntnu.no
More informationRecommendation ITU-R P (06/2017)
Recommendation ITU-R P.1238-9 (06/2017) Propagation data and prediction methods or the planning o indoor radiocommunication systems and radio local area networks in the requency range 300 MHz to 100 GHz
More informationCorrection of Clipped Pixels in Color Images
Correction of Clipped Pixels in Color Images IEEE Transaction on Visualization and Computer Graphics, Vol. 17, No. 3, 2011 Di Xu, Colin Doutre, and Panos Nasiopoulos Presented by In-Yong Song School of
More informationSignals and Systems II
1 To appear in IEEE Potentials Signals and Systems II Part III: Analytic signals and QAM data transmission Jerey O. Coleman Naval Research Laboratory, Radar Division This six-part series is a mini-course,
More informationNoise. Interference Noise
Noise David Johns and Ken Martin University o Toronto (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) University o Toronto 1 o 55 Intererence Noise Unwanted interaction between circuit and outside world
More informationCooperative Wideband Spectrum Sensing Based on Joint Sparsity
Virginia Commonwealth University VCU Scholars Compass Theses and Dissertations Graduate School 2017 Cooperative Wideband Spectrum Sensing Based on Joint Sparsity ghazaleh jowkar Follow this and additional
More informationPower Optimization in Stratix IV FPGAs
Power Optimization in Stratix IV FPGAs May 2008, ver.1.0 Application Note 514 Introduction The Stratix IV amily o devices rom Altera is based on 0.9 V, 40 nm Process technology. Stratix IV FPGAs deliver
More informationComputing for Engineers in Python
Computing for Engineers in Python Lecture 10: Signal (Image) Processing Autumn 2011-12 Some slides incorporated from Benny Chor s course 1 Lecture 9: Highlights Sorting, searching and time complexity Preprocessing
More informationPreprocessing and Segregating Offline Gujarati Handwritten Datasheet for Character Recognition
Preprocessing and Segregating Offline Gujarati Handwritten Datasheet for Character Recognition Hetal R. Thaker Atmiya Institute of Technology & science, Kalawad Road, Rajkot Gujarat, India C. K. Kumbharana,
More informationStudy on 3D CFBG Vibration Sensor and Its Application
PHOTONIC SENSORS / Vol. 6, No. 1, 2016: 90 96 Study on 3D CFBG Vibration Sensor and Its Application Qiuming NAN 1,2* and Sheng LI 1,2 1 National Engineering Laboratory or Fiber Optic Sensing Technology,
More informationImage Processing Of Oct Glaucoma Images And Information Theory Analysis
University of Denver Digital Commons @ DU Electronic Theses and Dissertations Graduate Studies 1-1-2009 Image Processing Of Oct Glaucoma Images And Information Theory Analysis Shuting Wang University of
More information3.6 Intersymbol interference. 1 Your site here
3.6 Intersymbol intererence 1 3.6 Intersymbol intererence what is intersymbol intererence and what cause ISI 1. The absolute bandwidth o rectangular multilevel pulses is ininite. The channels bandwidth
More informationPLL AND NUMBER OF SAMPLE SYNCHRONISATION TECHNIQUES FOR ELECTRICAL POWER QUALITY MEASURMENTS
XX IMEKO World Congress Metrology or Green Growth September 9 14, 2012, Busan, Republic o Korea PLL AND NUMBER OF SAMPLE SYNCHRONISATION TECHNIQUES FOR ELECTRICAL POWER QUALITY MEASURMENTS Richárd Bátori
More informationThe Research of Electric Energy Measurement Algorithm Based on S-Transform
International Conerence on Energy, Power and Electrical Engineering (EPEE 16 The Research o Electric Energy Measurement Algorithm Based on S-Transorm Xiyang Ou1,*, Bei He, Xiang Du1, Jin Zhang1, Ling Feng1,
More informationMicropulse Duty Cycle. # of eyes (20 ms) Total spots (200 ms)
Micropulse Duty Cycle Total spots (2 ms) # of eyes (2 ms) Total spots (2 ms) % 269 44 3 47% 9 4 4 25% 3 5 4 4 5% 2 4 3 5 2% 5 2 NA NA 9% 2 4 6% NA NA 57 2 5% 4 5 6 3 3% 39 5 35 5 # of eyes (2 ms) Supplemental
More informationValidation of a crystal detector model for the calibration of the Large Signal Network Analyzer.
Instrumentation and Measurement Technology Conerence IMTC 2007 Warsaw, Poland, May 1-3, 2007 Validation o a crystal detector model or the calibration o the Large Signal Network Analyzer. Liesbeth Gommé,
More informationPiecewise Mapping in HEVC Lossless Intraprediction
Piecewise Mapping in HEVC Lossless Intraprediction Coding Victor Sanchez Member IEEE Francesc Aulí-Llinàs Senior Member IEEE and Joan Serra-Sagristà Senior Member IEEE Abstract The lossless intra-prediction
More informationAmplifiers. Department of Computer Science and Engineering
Department o Computer Science and Engineering 2--8 Power ampliiers and the use o pulse modulation Switching ampliiers, somewhat incorrectly named digital ampliiers, have been growing in popularity when
More informationA new method for segmentation of retinal blood vessels using morphological image processing technique
A new method for segmentation of retinal blood vessels using morphological image processing technique Roya Aramesh Faculty of Computer and Information Technology Engineering,Qazvin Branch,Islamic Azad
More informationTopic 3 - Image Enhancement. (Part 2) Spatial Filtering
Topic 3 - Image Enhancement Part Spatial Filtering Spatial iltering - iltering operations that are perormed directl on the piels o an image Use o spatial mask Spatial iltering Operation moing the mask
More informationImage Database and Preprocessing
Chapter 3 Image Database and Preprocessing 3.1 Introduction The digital colour retinal images required for the development of automatic system for maculopathy detection are provided by the Department of
More informationEXPLOITING RMS TIME-FREQUENCY STRUCTURE FOR DATA COMPRESSION IN EMITTER LOCATION SYSTEMS
NAECON : National Aerospace & Electronics Conerence, October -,, Dayton, Ohio 7 EXPLOITING RMS TIME-FREQUENCY STRUCTURE FOR DATA COMPRESSION IN EMITTER LOCATION SYSTEMS MARK L. FOWLER Department o Electrical
More information24 Bit Image Noise Reduction with Median Filtering Algorithm
24 Bit Image Noise Reduction with Median iltering Algorithm Harmayani 1, Robbi Rahim 2* 1 Departement o Computer Engineering, AMIK Intel Com Global Indo Kisaran, Jl. Graha Abdi Satya No 4-48, Kisaran,
More informationNumerical Derivatives See also T&V, Appendix A.2 Gradient = vector of partial derivatives of image I(x,y) = [di(x,y)/dx, di(x,y)/dy]
I have put some Matlab image tutorials on Angel. Please take a look i you are unamiliar with Matlab or the image toolbox. Lecture : Linear Operators Administrivia I have posted Homework on Angel. It is
More informationSimulation Results for Permutation Trellis Codes using M-ary FSK
Simulation Results or Permutation Trellis Codes using M-ary FSK T.G. Swart, I. de Beer, H.C. Ferreira Department o Electrical and Electronic Engineering University o Johannesburg Auckland Park, South Arica
More informationDigital Image Processing Chapter 3: Image Enhancement in the Spatial Domain
Digital Image Processing Chapter 3: Image Enhancement in the Spatial Domain Principle Objective o Enhancement Process an image so that the result will be more suitable than the original image or a speciic
More informationPhase contrast microscopy with full numerical aperture illumination
Phase contrast microscopy with ull numerical aperture illumination Christian Maurer, Alexander Jesacher, Stean Bernet and Monika Ritsch-Marte Division or Biomedical Physics, Innsbruck Medical University,
More informationCLARUS 500 from ZEISS HD ultra-widefield fundus imaging
CLARUS 500 from ZEISS HD ultra-widefield fundus imaging Imaging ultra-wide without compromise. ZEISS CLARUS 500 // INNOVATION MADE BY ZEISS Compromising image quality may leave some pathology unseen. Signs
More informationIEEE C802.16h-05/022r1. IEEE Broadband Wireless Access Working Group <
Project IEEE 802.16 Broadband Wireless Access Working Group Title Cognitive radio concepts or 802.16h Date Submitted 2005-07-11 Source(s) Mariana Goldhamer Alvarion Tel Aviv, 21
More informationA Wavelet Approach to Wideband Spectrum Sensing for Cognitive Radios
A Wavelet Approach to Wideband Spectrum Sensing or Cognitive Radios Zhi Tian Department o Electrical & Computer Engineering Michigan Technological University Houghton, MI 4993 USA ztian@mtu.edu Georgios
More informationStudy Guide for the First Exam
Study Guide or the First Exam Chemistry 838 Fall 27 T V Atkinson Department o Chemistry Michigan State Uniersity East Lansing, MI 48824 Table o Contents Table o Contents...1 Table o Tables...1 Table o
More informationCLARUS 500 from ZEISS HD ultra-widefield fundus imaging
CLARUS 500 from ZEISS HD ultra-widefield fundus imaging Imaging ultra-wide without compromise. ZEISS CLARUS 500 // INNOVATION MADE BY ZEISS Compromising image quality may leave some pathology unseen. Signs
More informationConsumers are looking to wireless
Phase Noise Eects on OFDM Wireless LAN Perormance This article quantiies the eects o phase noise on bit-error rate and oers guidelines or noise reduction By John R. Pelliccio, Heinz Bachmann and Bruce
More informationSPEECH ENHANCEMENT BASED ON ITERATIVE WIENER FILTER USING COMPLEX SPEECH ANALYSIS
SPEECH ENHANCEMENT BASED ON TERATVE WENER FLTER USNG COMPLEX SPEECH ANALYSS Keiichi Funaki Computing & Networking Center, Univ. o the Ryukyus Senbaru, Nishihara, Okinawa, 93-3, Japan phone: +(8)98-895-8946,
More informationRequest Request Request Request Request Request Request
TITLE: DATE: March, 0 AFFECTED DOCUMENT: OCuLink.0 SPONSOR: Part I:. Summary o the Functional Changes PCI-SIG ENGINEERING CHANGE REQUEST OCuLink Cable Spec ECR Rev. Alex Haser (Molex), Jay Neer (Molex)
More informationClinical evaluation and management of glaucoma is largely
Macular Segmentation with Optical Coherence Tomography Hiroshi Ishikawa, 1,2 Daniel M. Stein, 1 Gadi Wollstein, 1,2 Siobahn Beaton, 1,2 James G. Fujimoto, 3 and Joel S. Schuman 1,2 PURPOSE. To develop
More informationOptimization of pupil design for point-scanning and line-scanning confocal microscopy
Optimization o pupil design or point-scanning and line-scanning conocal microscopy Yogesh G. Patel, 1,* Milind Rajadhyaksha, 2 and Charles A. DiMarzio 1,3 1 Electrical & Computer Engineering Department,
More informationCOMPARITIVE STUDY OF IMAGE DENOISING ALGORITHMS IN MEDICAL AND SATELLITE IMAGES
COMPARITIVE STUDY OF IMAGE DENOISING ALGORITHMS IN MEDICAL AND SATELLITE IMAGES Jyotsana Rastogi, Diksha Mittal, Deepanshu Singh ---------------------------------------------------------------------------------------------------------------------------------
More informationThe fourier spectrum analysis of optical feedback self-mixing signal under weak and moderate feedback
University o Wollongong Research Online Faculty o Inormatics - Papers (Archive) Faculty o Engineering and Inormation Sciences 8 The ourier spectrum analysis o optical eedback sel-mixing signal under weak
More informationA new zoom algorithm and its use in frequency estimation
Waves Wavelets Fractals Adv. Anal. 5; :7 Research Article Open Access Manuel D. Ortigueira, António S. Serralheiro, and J. A. Tenreiro Machado A new zoom algorithm and its use in requency estimation DOI.55/wwaa-5-
More informationSUPPLEMENTARY INFORMATION
Computational high-resolution optical imaging of the living human retina Nathan D. Shemonski 1,2, Fredrick A. South 1,2, Yuan-Zhi Liu 1,2, Steven G. Adie 3, P. Scott Carney 1,2, Stephen A. Boppart 1,2,4,5,*
More informationPhysics 142 Lenses and Mirrors Page 1. Lenses and Mirrors. Now for the sequence of events, in no particular order. Dan Rather
Physics 142 Lenses and Mirrors Page 1 Lenses and Mirrors Now or the sequence o events, in no particular order. Dan Rather Overview: making use o the laws o relection and reraction We will now study ormation
More informationMultiband Joint Detection with Correlated Spectral Occupancy in Wideband Cognitive Radios
Multiband Joint Detection with Correlated Spectral Occupancy in Wideband Cognitive Radios Khalid Hossain, Ayman Assra, and Benoît Champagne, Senior Member, IEEE Department o Electrical and Computer Engineering,
More information