Elasticity Imaging with Ultrasound JEE 4980 Final Report. George Michaels and Mary Watts
|
|
- Maude Ross
- 5 years ago
- Views:
Transcription
1 Elasticity Imaging with Ultrasound JEE 4980 Final Report George Michaels and Mary Watts University of Missouri, St. Louis Washington University Joint Engineering Undergraduate Program St. Louis, Missouri April 29, 2008
2 Abstract Gaining elasticity information of tissues from ultra-sound images would be useful in diagnosing cancers in patients without invasive biopsies. The elasticity would be determined by generating a strain image that tells you the tissues response to compression. Two tofu phantoms were used to simulate tissue with differing elasticity. Ultrasound images were then taken of the phantoms being compressed. Two images at different stages of compression were selected to be analyzed. These images were compared using cross-correlation to get the displacement of the second image with respect to the first. Then a displacement image was created. This image was then averaged to filter out noise and errors. The differential was then taken of the displacement image to create the strain image. Using this method, strain images were produced of the tofu phantoms that accurately portrayed the elasticity of the two types of tofu in the phantoms. Introduction Some cancers are hard to detect with current technology. Breast cancer mammography struggles with small cancers and with people with dense breast tissue. Ultrasound can detect different bodies in tissue but cannot determine whether those bodies are cancerous. Palpation is the current method used to evaluate the stiffness of these tumors; the relative stiffness can indicate the type of tumor or lump present. However, palpation is not possible with some of the lumps or tumors that are found in people, and even when practical, it is a highly subjective measure. Biopsy is then needed to achieve a diagnosis. A tool that could reliably determine the stiffness or elasticity of the tissue is needed to create a faster and more accurate diagnosis of cancer. We have developed a method of creating elasticity images. Elasticity imaging is a byproduct of the standard ultrasound. [A1] The images are constructed from the displacement and strain of ultrasound images of tissue in various stages of compression. Elasticity imaging has emerged as a sensitive and accurate method that more clearly demonstrates relative tissue stiffness or hardness, infiltration, and cystic region. [A1] This technology, if implemented correctly, could be crucial in the early detection of cancer and become a non-invasive replacement for biopsies. We have created an algorithm in MatLab that takes ultrasound images as the input and calculates and displays the strain image of the ultrasound signal.
3 Methodology A large part of our methodology was to take the project and break it down into smaller parts or stages. Then attack the smaller problems one at a time, building upon each completed part. The first step in the process was to create a simple algorithm that takes two signals as the input, cross-correlates them, and outputs the displacement. To accomplish this task we used simulation data that was provided to us. We used the simulation data because we knew that the displacement curve should look like a sinc function, so we would know if the algorithm was working properly based on the output. We developed xcorrfun.m, a function that takes a segment from each of the two signals and cross-correlates them. The function breaks the two signals into segments (windows) of different sizes, one larger than the other. The smaller of the two windows scans along the larger and cross-correlation matches up the data. The function then computes the amplitude difference of the two signals and outputs that number as the result. We found that up-sampling was very useful in the cross-correlation process. Next, we developed For1.m, a script that calls our xcorr function and plots the results. This script is where we obtained our displacement curve. The For1.m script allowed us to have variable parameters, for the window sizes and the upsampling rate. Having these adjustable parameters lead to a process of trial and error that assisted in reducing the error in our displacement image. After we were satisfied with our results for the displacement curve, we decided to apply our cross-correlation function to more complicated data. We used the Forearm Reduced 4980.mat file to test out function. Before we began to run the data on our function, we used group four s movie script to find frames that seemed to have the least amount of lateral movement. We wanted as little lateral movement as possible because we felt it would help our function have fewer errors as it runs only in one dimension, axially. We decided that frames 24 to 28 were the best frames for what we needed. Then we selected frames 24 and 25 to use in our script. We then selected a small square from the two frames that seemed to have some parts that compressed and some that did not. Groups from last semester also chose this general area and that aided in our decision as well.
4 Below are frames 24 and 25. These squares are 375 by 60 in size.
5 Magnified image of areas of interest in frame 24 and 25 We then wrote a script called ForearmFor.m i that used two for loops to implement our cross correlation function on the two squares. Our first attempts had many errors. To combat these errors we increased the upsampling rate from 6 to 9. We increased our small window size from 36 to 60 and we increased our maxlags to three times the sampling rate. These changes decreased our errors and produced a reasonable displacement image. We feel that it was a reasonable displacement image, despite the fact there were still errors, because it had linear increase in displacement as the compression occurred. After we decided that we were satisfied with our displacement image script, we moved on to work on the strain image. Obtaining the strain image was our main goal. All the work up to this point was the basis for the main part of the project. Before we attempted to get the strain image from the data files, we created test data to make sure we had our code working correctly. We made a simple program that loads and plots a sine wave. We used a sine wave because we knew that differentiating the sine function would produce a cosine function. Knowing this, we could judge how well the code worked based on the output. After we had the sine wave, we used the diff function to obtain the cosine wave. We then went one step farther and added random noise to the test data. We added the noise to simulate a more realistic situation because mostly like any practical application would not have perfectly smooth data. The noise presented a problem because when we applied the diff function to the noisy sine wave, it did not produce a clean cosine wave. We tried two different approaches to solve this problem.
6 The first was using the polyfit function. This function fits a smooth line to the data. Applying this function to our noisy sine wave produced a good result, but in the end, we opted not to use this method. Although this method gave good results, we concluded that it might be problematic on complicated data. The next approach we tried was to average the noisy sine wave and then use the diff function. We used an averaging function that we found in the MatLab help section. The function takes a single point, creates a window around it, and then averages all the points within the window. Then the window slides to a new point and does the same process. This continually happens until the end of the data. The averaging function essentially smoothes out the curve, eliminating any large spikes. This is the method we used to help generate a smoother displacement image, essentially eliminating noise. We then created a function, Strainpic.m, to generate the actual strain images from the data files. The function takes the displacement image as the input and applies the averaging function on it. Then runs the diff function on the smooth, averaged displacement and outputs the strain image. This is the function we used to generate our final results.
7 Results Preliminary Results: Displacement Curve for simulation data from our For1.m ii script The ideal displacement curve should look like a smooth sinc function. Our curve has the general outline of the sinc function, but has undesirable blips. We determined that the blips were acceptable for us to continue with our method.
8 Displacement image for ForearmFor.m ii script using the forearm data. We are satisfied with this displacement image. Although it may look a little crude, we feel it was a good result. A good displacement image will have a linearly increasing displacement as it moves down the data due to the nature of compression from the top. As you can see our picture does have this displacement characteristic so we felt our method was working.
9 Final Results: Below are results from Tofu1 phantom created by Dr. Trobaugh. This phantom has softer tofu inside of the firmer tofu. Tofu1: Displacement image (created with TofuFor.m Frames 40 and 50 Size 100:1100 and 1:256) This is a good displacement image. Again we have a linearly increasing displacement from top to bottom. The very top we expect some errors due to the ultrasound wand. The specs are errors from the cross-correlation function and are acceptable for our purposes.
10 Tofu1: Strain image (created with Strain.m) This is a good result. There is a clear area where the tofu changes density. We would expect the softer tissue inside the circle to compress more and the color inside the tofu is about.4 whereas the outside, which we would expect to have less compression, is.1 to.2. The red and orange on the top left of the cylinder might be an air bubble in the tofu 1 data. The area at the top is a result of our averaging function, the area is the same length as our averaging window, up until that point the averaging function is using a lot of zeros to average each point. Below are results from Tofu1 phantom created by Dr. Trobaugh. This phantom has a firmer tofu cylinder inside of the softer tofu.
11 Tofu2: Displacement image (created with TofuFor.m using Frames 40 and 50 Size 100:1200 and 1:256) Once again this displacement image is giving good results. There is a linearly increasing displacement which is disrupted slightly by the firmer tofu. There are still errors in this image but they are manageable. We have a larger amount of errors, which show up as red dots under the cylinder of tofu in this image. There are more errors in this region than any other. We think our cross-correlation function struggled with the displacement in this area due to the firmer tofu causing the displacement to change more drastically in this area.
12 Tofu2: Strain image (created with Strain.m) This is the best result we obtained. We would expect from the phantom composition that the cylinder would have less compression and the surrounding tofu to have more compression. Using the colorbar, the cylinder of firmer tofu has a strain of very close to 0 and the softer tofu has a strain of.2. The errors from the displacement picture transferred over to the strain image but it does not seem to affect the results very much.
13 Conclusion We believe this project was a success. Based on the results we feel that we have successfully created a software program that is capable of taking ultrasound images and outputting an accurate strain image. We feel this way because our Tofu2 results closely resemble the results from the Righetti paper, which used an almost identical phantom. It is clear from our strain images that the elasticity is different in the two sections of tofu, just as it is in the published paper. Our images may not be as clear as those published, that is due to the errors produced in our code, which are a result of our limited coding ability. Overall, we feel our results are good enough to demonstrate that we have a clear understanding of the problem and have engineered an adequate solution. Recommendations Based on our project so far, we believe we have a solid foundation for future work. There are things that could be done to improve the function of our algorithm and the quality of the results. The first step for improvement would be to refine the code. The number of errors would need to be reduced. This could be addressed by adding spike prevention which other groups have successfully done. This would eliminate the visible blips from our images and provide clearer results. Two-dimensional tracking should also eventually be added. This addition could fully compensate for any lateral movement which our current code has great difficulty handling. Another basic area for improvement would be to make the algorithm run faster so it could work in a real-life situation. Our code takes a long time to run on large pieces of data, this would prevent real-time strain image production. Our algorithm is comprised of very basic coding techniques and more advanced code would most likely help speed up the process.
14 MatLab Scripts and Functions Xcorrfun.m %% George Michaels & Mary Watts -- Senior Design %% takes two signals and xcorr's them, plots the xcorr and also plots the signals against each %% other. function [shift] = XcorrFun(y1,y2,samplerate) %% Cross-correlation of y1 and y2 maxlags = 6*samplerate; [c,lags] = xcorr(y1,y2,maxlags); %% Getting Max location of cross-correlation [Y, I] = max(c); L = lags(i); %% gave us displacement of the lags max %figure(1); %clf %hold on %%Temp comment %%plot(lags, c) %% plotting cross-correlation vs. lags %%Temp Comment shift = L; figure(1) clf hold on %%plot signal one plot(y1) %%plot signal two plot([1:length(y2)] + L, y2, 'r') title ('Cross-Correlation Match-up')
15 TofuFor.m %% Mary Watts & George Michaels %% For loop code takes two signals and using xcorr function to get a %% displacement image. %% Data Input load tfu2.mat y1 = b_data040; y2 = b_data050; %% load appropriate data into workspace %% Enter first frame file name %% Enter second frame file name %% Selecting portion of data to cross-correlate rowindex = 100:1200; %% This is the vertical range can be changed to zoom in on areas of interest columnindex = 1:256; %% This is the horizontal range can be changed to zoom in on areas of interest shift = [1:1100,1:1100]; %% This creates a shift matrix to store the final portion in, the size is at least the size of the data. %% Loop to breakdown image matrix to a vector for upsampling and cross-correlating for x = columnindex; y1vector = y1(rowindex, x); y2vector = y2(rowindex, x); %% Selects out a vertical vector of data %% upsampling samplerate = 9; %% Sample rate can be adjusted to get best results, we found 5-10 a good range. y1up = resample(y1vector,samplerate,1); y2up = resample(y2vector,samplerate,1); [row, column] = size(y1up); %% getting size of data %% Window size parameters size1 = 30; %% These can be adjusted as needed. size2 = 50; diffofsize = size2 - size1; for i = 1+size2*samplerate:2*samplerate:row-size2*samplerate central point to move around %%starting i as a Start1 = i-size1*samplerate; Stop1 = i+size1*samplerate; %%small window start %%small window end
16 Start2 = i-size2*samplerate; %%big window start Stop2 = i+size2*samplerate; %%big window end y1shift = [zeros(diffofsize*samplerate,1); y1up(start1:stop1,1)]; %% Padding zeros to small window to make them start at %% beginning point of big window shift(i, x) = XcorrFun(y1shift, y2up(start2:stop2), samplerate); end figure(2) plot(shift(1+size2*samplerate:2*samplerate:row-size2*samplerate, x)); %% This function allows us to watch the function and see how it is working, %% will help gauge if maxlags is set to correct amount. drawnow end [row2, column2] = size(shift); shiftcut = shift(1+size2*samplerate:2*samplerate:row-size2*samplerate, columnindex); is the final displacement matrix. %%This figure(3) imagesc([0 4], [0 4],shiftcut); figure(gcf) %% Display displacement image colorbar caxis([15 40]) %% Colorbar settings might need adjustment for different data title ('Displacement Picture', 'FontSize', 12) xlabel('cm','fontsize',12) ylabel('cm','fontsize',12) axis equal
17 StrainPic.m %% Mary Watts and George Michaels %% JEE 4980 Spring Semester 2008 %% Strain Image %% This script takes the output of the TofuFor.m file, shiftcut and creates a strain image. windowsize = 25; %% Set size of window to average around average = filter(ones(1,windowsize)/windowsize,1,shiftcut); %% Averaging output from TofuFor.m DiffFit = diff(average); %% Take difference of averaged data clf figure(1) imagesc([0 4], [0 4],average); %% Display averaged displacement image title ('Tofu2 Phantom Averaged Displacement Image') xlabel('cm','fontsize',12) ylabel('cm','fontsize',12) axis equal clf figure(2) imagesc([0 4], [0 4],DiffFit); figure(gcf) colorbar caxis([0 1]) title ('Tofu2 Phantom Strain Image') xlabel('cm','fontsize',12) ylabel('cm','fontsize',12) axis equal %% Display strain image
18 References [A1] Revolutionary Elasticity Imaging Siemens. 27 Feb < AAPM/RSNA Physics Tutorial for Residents: Topics in US. Timothy J. Hall, PhD. RadioGraphics The feasibility of estimating and imaging the mechanical behavior of poroelastic materials using axial strain elastgraphy. Righetti, Rafaella, Mariapaola Righetti, Jonathan Ophir and Thomas A Krouskop, Physics in Medicine and Biology 15 May 2007 Group 1 Final Project Groups 2 Final Project Group 3 Final Project Group 4 Final Project
19
JEE 4980 Senior Design
WASHINGTON UNIVERSITY JEE 498 Senior Design Ultrasound Elasticity Imaging Final Report Steven Goodwin Mark Green 12/16/28 ABSTRACT Within this report, we have developed a useful way of retrieving a strain
More informationSenior Design Project
Senior Design Project Strain Imaging Josh Tischler and Ken Nguyen JEE 4980, Fall 2008 UMSL/Washington University Joint Undergraduate Engineering Program December 16, 2008 1 Introduction The purpose of
More informationDATA COMPRESSION USING THE FFT
EEE 407/591 PROJECT DUE: NOVEMBER 21, 2001 DATA COMPRESSION USING THE FFT INSTRUCTOR: DR. ANDREAS SPANIAS TEAM MEMBERS: IMTIAZ NIZAMI - 993 21 6600 HASSAN MANSOOR - 993 69 3137 Contents TECHNICAL BACKGROUND...
More informationChapter 7. Scanner Controls
Chapter 7 Scanner Controls Gain Compensation Echoes created by similar acoustic mismatches at interfaces deeper in the body return to the transducer with weaker amplitude than those closer because of the
More informationDoubletalk Detection
ELEN-E4810 Digital Signal Processing Fall 2004 Doubletalk Detection Adam Dolin David Klaver Abstract: When processing a particular voice signal it is often assumed that the signal contains only one speaker,
More informationSpectrum Analyser Basics
Hands-On Learning Spectrum Analyser Basics Peter D. Hiscocks Syscomp Electronic Design Limited Email: phiscock@ee.ryerson.ca June 28, 2014 Introduction Figure 1: GUI Startup Screen In a previous exercise,
More informationAutomatic LP Digitalization Spring Group 6: Michael Sibley, Alexander Su, Daphne Tsatsoulis {msibley, ahs1,
Automatic LP Digitalization 18-551 Spring 2011 Group 6: Michael Sibley, Alexander Su, Daphne Tsatsoulis {msibley, ahs1, ptsatsou}@andrew.cmu.edu Introduction This project was originated from our interest
More informationGetting Started. Connect green audio output of SpikerBox/SpikerShield using green cable to your headphones input on iphone/ipad.
Getting Started First thing you should do is to connect your iphone or ipad to SpikerBox with a green smartphone cable. Green cable comes with designators on each end of the cable ( Smartphone and SpikerBox
More informationThe Cocktail Party Effect. Binaural Masking. The Precedence Effect. Music 175: Time and Space
The Cocktail Party Effect Music 175: Time and Space Tamara Smyth, trsmyth@ucsd.edu Department of Music, University of California, San Diego (UCSD) April 20, 2017 Cocktail Party Effect: ability to follow
More informationLabView Exercises: Part II
Physics 3100 Electronics, Fall 2008, Digital Circuits 1 LabView Exercises: Part II The working VIs should be handed in to the TA at the end of the lab. Using LabView for Calculations and Simulations LabView
More informationPredicting the immediate future with Recurrent Neural Networks: Pre-training and Applications
Predicting the immediate future with Recurrent Neural Networks: Pre-training and Applications Introduction Brandon Richardson December 16, 2011 Research preformed from the last 5 years has shown that the
More informationNENS 230 Assignment #2 Data Import, Manipulation, and Basic Plotting
NENS 230 Assignment #2 Data Import, Manipulation, and Basic Plotting Compound Action Potential Due: Tuesday, October 6th, 2015 Goals Become comfortable reading data into Matlab from several common formats
More informationDELTA MODULATION AND DPCM CODING OF COLOR SIGNALS
DELTA MODULATION AND DPCM CODING OF COLOR SIGNALS Item Type text; Proceedings Authors Habibi, A. Publisher International Foundation for Telemetering Journal International Telemetering Conference Proceedings
More informationSupplemental Material for Gamma-band Synchronization in the Macaque Hippocampus and Memory Formation
Supplemental Material for Gamma-band Synchronization in the Macaque Hippocampus and Memory Formation Michael J. Jutras, Pascal Fries, Elizabeth A. Buffalo * *To whom correspondence should be addressed.
More informationDistortion Analysis Of Tamil Language Characters Recognition
www.ijcsi.org 390 Distortion Analysis Of Tamil Language Characters Recognition Gowri.N 1, R. Bhaskaran 2, 1. T.B.A.K. College for Women, Kilakarai, 2. School Of Mathematics, Madurai Kamaraj University,
More informationStaMPS Persistent Scatterer Practical
StaMPS Persistent Scatterer Practical ESA Land Training Course, Leicester, 10-14 th September, 2018 Andy Hooper, University of Leeds a.hooper@leeds.ac.uk This practical exercise consists of working through
More informationA Novel Approach towards Video Compression for Mobile Internet using Transform Domain Technique
A Novel Approach towards Video Compression for Mobile Internet using Transform Domain Technique Dhaval R. Bhojani Research Scholar, Shri JJT University, Jhunjunu, Rajasthan, India Ved Vyas Dwivedi, PhD.
More informationStaMPS Persistent Scatterer Exercise
StaMPS Persistent Scatterer Exercise ESA Land Training Course, Bucharest, 14-18 th September, 2015 Andy Hooper, University of Leeds a.hooper@leeds.ac.uk This exercise consists of working through an example
More informationHow to Obtain a Good Stereo Sound Stage in Cars
Page 1 How to Obtain a Good Stereo Sound Stage in Cars Author: Lars-Johan Brännmark, Chief Scientist, Dirac Research First Published: November 2017 Latest Update: November 2017 Designing a sound system
More informationPlease feel free to download the Demo application software from analogarts.com to help you follow this seminar.
Hello, welcome to Analog Arts spectrum analyzer tutorial. Please feel free to download the Demo application software from analogarts.com to help you follow this seminar. For this presentation, we use a
More informationhit), and assume that longer incidental sounds (forest noise, water, wind noise) resemble a Gaussian noise distribution.
CS 229 FINAL PROJECT A SOUNDHOUND FOR THE SOUNDS OF HOUNDS WEAKLY SUPERVISED MODELING OF ANIMAL SOUNDS ROBERT COLCORD, ETHAN GELLER, MATTHEW HORTON Abstract: We propose a hybrid approach to generating
More informationSimilarity Measurement of Biological Signals Using Dynamic Time Warping Algorithm
Similarity Measurement of Biological Signals Using Dynamic Time Warping Algorithm Ivan Luzianin 1, Bernd Krause 2 1,2 Anhalt University of Applied Sciences Computer Science and Languages Department Lohmannstr.
More informationThe Effect of Time-Domain Interpolation on Response Spectral Calculations. David M. Boore
The Effect of Time-Domain Interpolation on Response Spectral Calculations David M. Boore This note confirms Norm Abrahamson s finding that the straight line interpolation between sampled points used in
More informationLecture 2 Video Formation and Representation
2013 Spring Term 1 Lecture 2 Video Formation and Representation Wen-Hsiao Peng ( 彭文孝 ) Multimedia Architecture and Processing Lab (MAPL) Department of Computer Science National Chiao Tung University 1
More informationBroken Wires Diagnosis Method Numerical Simulation Based on Smart Cable Structure
PHOTONIC SENSORS / Vol. 4, No. 4, 2014: 366 372 Broken Wires Diagnosis Method Numerical Simulation Based on Smart Cable Structure Sheng LI 1*, Min ZHOU 2, and Yan YANG 3 1 National Engineering Laboratory
More informationApplication Note AN-708 Vibration Measurements with the Vibration Synchronization Module
Application Note AN-708 Vibration Measurements with the Vibration Synchronization Module Introduction The vibration module allows complete analysis of cyclical events using low-speed cameras. This is accomplished
More informationPole Zero Correction using OBSPY and PSN Data
Pole Zero Correction using OBSPY and PSN Data Obspy provides the possibility of instrument response correction. WinSDR and WinQuake already have capability to embed the required information into the event
More informationGraphing Your Motion
Name Date Graphing Your Motion Experiment 33 Graphs made using a Motion Detector can be used to study motion. In this experiment, you will use a Motion Detector to make graphs of your own motion. OBJECTIVES
More informationECE438 - Laboratory 1: Discrete and Continuous-Time Signals
Purdue University: ECE438 - Digital Signal Processing with Applications 1 ECE438 - Laboratory 1: Discrete and Continuous-Time Signals By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015 1 Introduction
More informationCHARACTERIZATION OF END-TO-END DELAYS IN HEAD-MOUNTED DISPLAY SYSTEMS
CHARACTERIZATION OF END-TO-END S IN HEAD-MOUNTED DISPLAY SYSTEMS Mark R. Mine University of North Carolina at Chapel Hill 3/23/93 1. 0 INTRODUCTION This technical report presents the results of measurements
More informationCURIE Day 3: Frequency Domain Images
CURIE Day 3: Frequency Domain Images Curie Academy, July 15, 2015 NAME: NAME: TA SIGN-OFFS Exercise 7 Exercise 13 Exercise 17 Making 8x8 pictures Compressing a grayscale image Satellite image debanding
More informationEE373B Project Report Can we predict general public s response by studying published sales data? A Statistical and adaptive approach
EE373B Project Report Can we predict general public s response by studying published sales data? A Statistical and adaptive approach Song Hui Chon Stanford University Everyone has different musical taste,
More informationImplementation of an MPEG Codec on the Tilera TM 64 Processor
1 Implementation of an MPEG Codec on the Tilera TM 64 Processor Whitney Flohr Supervisor: Mark Franklin, Ed Richter Department of Electrical and Systems Engineering Washington University in St. Louis Fall
More informationCS229 Project Report Polyphonic Piano Transcription
CS229 Project Report Polyphonic Piano Transcription Mohammad Sadegh Ebrahimi Stanford University Jean-Baptiste Boin Stanford University sadegh@stanford.edu jbboin@stanford.edu 1. Introduction In this project
More informationLab 5 Linear Predictive Coding
Lab 5 Linear Predictive Coding 1 of 1 Idea When plain speech audio is recorded and needs to be transmitted over a channel with limited bandwidth it is often necessary to either compress or encode the audio
More informationReconstruction of Ca 2+ dynamics from low frame rate Ca 2+ imaging data CS229 final project. Submitted by: Limor Bursztyn
Reconstruction of Ca 2+ dynamics from low frame rate Ca 2+ imaging data CS229 final project. Submitted by: Limor Bursztyn Introduction Active neurons communicate by action potential firing (spikes), accompanied
More informationAudio Compression Technology for Voice Transmission
Audio Compression Technology for Voice Transmission 1 SUBRATA SAHA, 2 VIKRAM REDDY 1 Department of Electrical and Computer Engineering 2 Department of Computer Science University of Manitoba Winnipeg,
More informationCSE 166: Image Processing. Overview. Representing an image. What is an image? History. What is image processing? Today. Image Processing CSE 166
CSE 166: Image Processing Overview Image Processing CSE 166 Today Course overview Logistics Some mathematics MATLAB Lectures will be boardwork and slides Take written notes or take pictures of the board
More informationKoester Performance Research Koester Performance Research Heidi Koester, Ph.D. Rich Simpson, Ph.D., ATP
Scanning Wizard software for optimizing configuration of switch scanning systems Heidi Koester, Ph.D. hhk@kpronline.com, Ann Arbor, MI www.kpronline.com Rich Simpson, Ph.D., ATP rsimps04@nyit.edu New York
More informationA Matlab toolbox for. Characterisation Of Recorded Underwater Sound (CHORUS) USER S GUIDE
Centre for Marine Science and Technology A Matlab toolbox for Characterisation Of Recorded Underwater Sound (CHORUS) USER S GUIDE Version 5.0b Prepared for: Centre for Marine Science and Technology Prepared
More informationECG Denoising Using Singular Value Decomposition
Australian Journal of Basic and Applied Sciences, 4(7): 2109-2113, 2010 ISSN 1991-8178 ECG Denoising Using Singular Value Decomposition 1 Mojtaba Bandarabadi, 2 MohammadReza Karami-Mollaei, 3 Amard Afzalian,
More informationA Parametric Autoregressive Model for the Extraction of Electric Network Frequency Fluctuations in Audio Forensic Authentication
Journal of Energy and Power Engineering 10 (2016) 504-512 doi: 10.17265/1934-8975/2016.08.007 D DAVID PUBLISHING A Parametric Autoregressive Model for the Extraction of Electric Network Frequency Fluctuations
More informationCharacterization and improvement of unpatterned wafer defect review on SEMs
Characterization and improvement of unpatterned wafer defect review on SEMs Alan S. Parkes *, Zane Marek ** JEOL USA, Inc. 11 Dearborn Road, Peabody, MA 01960 ABSTRACT Defect Scatter Analysis (DSA) provides
More informationPS User Guide Series Seismic-Data Display
PS User Guide Series 2015 Seismic-Data Display Prepared By Choon B. Park, Ph.D. January 2015 Table of Contents Page 1. File 2 2. Data 2 2.1 Resample 3 3. Edit 4 3.1 Export Data 4 3.2 Cut/Append Records
More informationDigital Image and Fourier Transform
Lab 5 Numerical Methods TNCG17 Digital Image and Fourier Transform Sasan Gooran (Autumn 2009) Before starting this lab you are supposed to do the preparation assignments of this lab. All functions and
More informationLaboratory Assignment 3. Digital Music Synthesis: Beethoven s Fifth Symphony Using MATLAB
Laboratory Assignment 3 Digital Music Synthesis: Beethoven s Fifth Symphony Using MATLAB PURPOSE In this laboratory assignment, you will use MATLAB to synthesize the audio tones that make up a well-known
More informationAudacity Tips and Tricks for Podcasters
Audacity Tips and Tricks for Podcasters Common Challenges in Podcast Recording Pops and Clicks Sometimes audio recordings contain pops or clicks caused by a too hard p, t, or k sound, by just a little
More informationMotion Video Compression
7 Motion Video Compression 7.1 Motion video Motion video contains massive amounts of redundant information. This is because each image has redundant information and also because there are very few changes
More informationSingle-switch Scanning Example. Learning Objectives. Enhancing Efficiency for People who Use Switch Scanning. Overview. Part 1. Single-switch Scanning
Enhancing Efficiency for People who Use Switch Scanning Heidi Koester, Ph.D. hhk@kpronline.com, Ann Arbor, MI www.kpronline.com Rich Simpson, Ph.D., ATP richard.c.simpson@gmail.com Duquesne University
More informationA Parametric Autoregressive Model for the Extraction of Electric Network Frequency Fluctuations in Audio Forensic Authentication
Proceedings of the 3 rd International Conference on Control, Dynamic Systems, and Robotics (CDSR 16) Ottawa, Canada May 9 10, 2016 Paper No. 110 DOI: 10.11159/cdsr16.110 A Parametric Autoregressive Model
More informationMIE 402: WORKSHOP ON DATA ACQUISITION AND SIGNAL PROCESSING Spring 2003
MIE 402: WORKSHOP ON DATA ACQUISITION AND SIGNAL PROCESSING Spring 2003 OBJECTIVE To become familiar with state-of-the-art digital data acquisition hardware and software. To explore common data acquisition
More informationMarc I. Johnson, Texture Technologies Corp. 6 Patton Drive, Hamilton, MA Tel
Abstract Novel Automated Method for Analyzing Peel Adhesion Ben Senning, Territory Manager, Texture Technologies Corp, Hamilton, MA Marc Johnson, President, Texture Technologies Corp, Hamilton, MA Most
More informationAnalysis of WFS Measurements from first half of 2004
Analysis of WFS Measurements from first half of 24 (Report4) Graham Cox August 19, 24 1 Abstract Described in this report is the results of wavefront sensor measurements taken during the first seven months
More informationMITOCW ocw f08-lec19_300k
MITOCW ocw-18-085-f08-lec19_300k The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free.
More informationUsing Multiple DMs for Increased Spatial Frequency Response
AN: Multiple DMs for Increased Spatial Frequency Response Using Multiple DMs for Increased Spatial Frequency Response AN Author: Justin Mansell Revision: // Abstract Some researchers have come to us suggesting
More informationAn Overview of Video Coding Algorithms
An Overview of Video Coding Algorithms Prof. Ja-Ling Wu Department of Computer Science and Information Engineering National Taiwan University Video coding can be viewed as image compression with a temporal
More informationCommunication Theory and Engineering
Communication Theory and Engineering Master's Degree in Electronic Engineering Sapienza University of Rome A.A. 2018-2019 Practice work 14 Image signals Example 1 Calculate the aspect ratio for an image
More informationSTAT 113: Statistics and Society Ellen Gundlach, Purdue University. (Chapters refer to Moore and Notz, Statistics: Concepts and Controversies, 8e)
STAT 113: Statistics and Society Ellen Gundlach, Purdue University (Chapters refer to Moore and Notz, Statistics: Concepts and Controversies, 8e) Learning Objectives for Exam 1: Unit 1, Part 1: Population
More informationONE SENSOR MICROPHONE ARRAY APPLICATION IN SOURCE LOCALIZATION. Hsin-Chu, Taiwan
ICSV14 Cairns Australia 9-12 July, 2007 ONE SENSOR MICROPHONE ARRAY APPLICATION IN SOURCE LOCALIZATION Percy F. Wang 1 and Mingsian R. Bai 2 1 Southern Research Institute/University of Alabama at Birmingham
More informationType-2 Fuzzy Logic Sensor Fusion for Fire Detection Robots
Proceedings of the 2 nd International Conference of Control, Dynamic Systems, and Robotics Ottawa, Ontario, Canada, May 7 8, 2015 Paper No. 187 Type-2 Fuzzy Logic Sensor Fusion for Fire Detection Robots
More informationWhat to consider when choosing a mammography display
What to consider when choosing a mammography display Screen size and resolution In digital breast imaging, the quality of the medical display has a direct impact on the decisions you make. Next to display
More informationDetection and demodulation of non-cooperative burst signal Feng Yue 1, Wu Guangzhi 1, Tao Min 1
International Conference on Applied Science and Engineering Innovation (ASEI 2015) Detection and demodulation of non-cooperative burst signal Feng Yue 1, Wu Guangzhi 1, Tao Min 1 1 China Satellite Maritime
More informationMulti-Shaped E-Beam Technology for Mask Writing
Multi-Shaped E-Beam Technology for Mask Writing Juergen Gramss a, Arnd Stoeckel a, Ulf Weidenmueller a, Hans-Joachim Doering a, Martin Bloecker b, Martin Sczyrba b, Michael Finken b, Timo Wandel b, Detlef
More informationSupplemental Information. Dynamic Theta Networks in the Human Medial. Temporal Lobe Support Episodic Memory
Current Biology, Volume 29 Supplemental Information Dynamic Theta Networks in the Human Medial Temporal Lobe Support Episodic Memory Ethan A. Solomon, Joel M. Stein, Sandhitsu Das, Richard Gorniak, Michael
More informationEDDY CURRENT IMAGE PROCESSING FOR CRACK SIZE CHARACTERIZATION
EDDY CURRENT MAGE PROCESSNG FOR CRACK SZE CHARACTERZATON R.O. McCary General Electric Co., Corporate Research and Development P. 0. Box 8 Schenectady, N. Y. 12309 NTRODUCTON Estimation of crack length
More informationWhy Engineers Ignore Cable Loss
Why Engineers Ignore Cable Loss By Brig Asay, Agilent Technologies Companies spend large amounts of money on test and measurement equipment. One of the largest purchases for high speed designers is a real
More informationSystem Identification
System Identification Arun K. Tangirala Department of Chemical Engineering IIT Madras July 26, 2013 Module 9 Lecture 2 Arun K. Tangirala System Identification July 26, 2013 16 Contents of Lecture 2 In
More informationSpeech and Speaker Recognition for the Command of an Industrial Robot
Speech and Speaker Recognition for the Command of an Industrial Robot CLAUDIA MOISA*, HELGA SILAGHI*, ANDREI SILAGHI** *Dept. of Electric Drives and Automation University of Oradea University Street, nr.
More informationModule 8 VIDEO CODING STANDARDS. Version 2 ECE IIT, Kharagpur
Module 8 VIDEO CODING STANDARDS Lesson 27 H.264 standard Lesson Objectives At the end of this lesson, the students should be able to: 1. State the broad objectives of the H.264 standard. 2. List the improved
More informationFilm Sequence Detection and Removal in DTV Format and Standards Conversion
TeraNex Technical Presentation Film Sequence Detection and Removal in DTV Format and Standards Conversion 142nd SMPTE Technical Conference & Exhibition October 20, 2000 Scott Ackerman DTV Product Manager
More informationSmoothing Techniques For More Accurate Signals
INDICATORS Smoothing Techniques For More Accurate Signals More sophisticated smoothing techniques can be used to determine market trend. Better trend recognition can lead to more accurate trading signals.
More informationON THE INTERPOLATION OF ULTRASONIC GUIDED WAVE SIGNALS
ON THE INTERPOLATION OF ULTRASONIC GUIDED WAVE SIGNALS Jennifer E. Michaels 1, Ren-Jean Liou 2, Jason P. Zutty 1, and Thomas E. Michaels 1 1 School of Electrical & Computer Engineering, Georgia Institute
More informationADJUSTABLE RANGE OF PARAMETERS [End of the Book]
Page Introduction FEATURE OF AG-HPX250 S BUILT-IN SENSOR 2 Chapter 1. SCENE FILE 4 Chapter 2. DETAIL / V DETAIL 6 Chapter 3. KNEE 8 Chapter 4. GAMMA 10 Chapter 5. DRS 12 Chapter 6. DETAIL CORING /SKIN
More informationTime Domain Simulations
Accuracy of the Computational Experiments Called Mike Steinberger Lead Architect Serial Channel Products SiSoft Time Domain Simulations Evaluation vs. Experimentation We re used to thinking of results
More informationWEAVE: Web-based Educational Framework for Analysis, Visualization, and Experimentation. Steven M. Lattanzio II 1
WEAVE: Web-based Educational Framework for Analysis, Visualization, and Experimentation Steven M. Lattanzio II 1 Abstract WEAVE (Web-Based Educational Framework for Analysis, Visualization, and Experimentation)
More informationy POWER USER MUSIC PRODUCTION and PERFORMANCE With the MOTIF ES Mastering the Sample SLICE function
y POWER USER MUSIC PRODUCTION and PERFORMANCE With the MOTIF ES Mastering the Sample SLICE function Phil Clendeninn Senior Product Specialist Technology Products Yamaha Corporation of America Working with
More informationAppendix D. UW DigiScope User s Manual. Willis J. Tompkins and Annie Foong
Appendix D UW DigiScope User s Manual Willis J. Tompkins and Annie Foong UW DigiScope is a program that gives the user a range of basic functions typical of a digital oscilloscope. Included are such features
More informationGetting Started with the LabVIEW Sound and Vibration Toolkit
1 Getting Started with the LabVIEW Sound and Vibration Toolkit This tutorial is designed to introduce you to some of the sound and vibration analysis capabilities in the industry-leading software tool
More informationRemoving the Pattern Noise from all STIS Side-2 CCD data
The 2010 STScI Calibration Workshop Space Telescope Science Institute, 2010 Susana Deustua and Cristina Oliveira, eds. Removing the Pattern Noise from all STIS Side-2 CCD data Rolf A. Jansen, Rogier Windhorst,
More informationHugo Technology. An introduction into Rob Watts' technology
Hugo Technology An introduction into Rob Watts' technology Copyright Rob Watts 2014 About Rob Watts Audio chip designer both analogue and digital Consultant to silicon chip manufacturers Designer of Chord
More informationME EN 363 ELEMENTARY INSTRUMENTATION Lab: Basic Lab Instruments and Data Acquisition
ME EN 363 ELEMENTARY INSTRUMENTATION Lab: Basic Lab Instruments and Data Acquisition INTRODUCTION Many sensors produce continuous voltage signals. In this lab, you will learn about some common methods
More informationNormalization Methods for Two-Color Microarray Data
Normalization Methods for Two-Color Microarray Data 1/13/2009 Copyright 2009 Dan Nettleton What is Normalization? Normalization describes the process of removing (or minimizing) non-biological variation
More informationCHAPTER 2 SUBCHANNEL POWER CONTROL THROUGH WEIGHTING COEFFICIENT METHOD
CHAPTER 2 SUBCHANNEL POWER CONTROL THROUGH WEIGHTING COEFFICIENT METHOD 2.1 INTRODUCTION MC-CDMA systems transmit data over several orthogonal subcarriers. The capacity of MC-CDMA cellular system is mainly
More informationReduced complexity MPEG2 video post-processing for HD display
Downloaded from orbit.dtu.dk on: Dec 17, 2017 Reduced complexity MPEG2 video post-processing for HD display Virk, Kamran; Li, Huiying; Forchhammer, Søren Published in: IEEE International Conference on
More informationAltiumLive 2017: Effective Methods for Advanced Routing
AltiumLive 2017: Effective Methods for Advanced Routing Charles Pfeil Senior Product Manager Dave Cousineau Sr. Field Applications Engineer Charles Pfeil Senior Product Manager Over 50 years of experience
More information2. AN INTROSPECTION OF THE MORPHING PROCESS
1. INTRODUCTION Voice morphing means the transition of one speech signal into another. Like image morphing, speech morphing aims to preserve the shared characteristics of the starting and final signals,
More informationMODE FIELD DIAMETER AND EFFECTIVE AREA MEASUREMENT OF DISPERSION COMPENSATION OPTICAL DEVICES
MODE FIELD DIAMETER AND EFFECTIVE AREA MEASUREMENT OF DISPERSION COMPENSATION OPTICAL DEVICES Hale R. Farley, Jeffrey L. Guttman, Razvan Chirita and Carmen D. Pâlsan Photon inc. 6860 Santa Teresa Blvd
More informationMusical Hit Detection
Musical Hit Detection CS 229 Project Milestone Report Eleanor Crane Sarah Houts Kiran Murthy December 12, 2008 1 Problem Statement Musical visualizers are programs that process audio input in order to
More informationPeak Dynamic Power Estimation of FPGA-mapped Digital Designs
Peak Dynamic Power Estimation of FPGA-mapped Digital Designs Abstract The Peak Dynamic Power Estimation (P DP E) problem involves finding input vector pairs that cause maximum power dissipation (maximum
More informationSCANNER TUNING TUTORIAL Author: Adam Burns
SCANNER TUNING TUTORIAL Author: Adam Burns Let me say first of all that nearly all the techniques mentioned in this tutorial were gleaned from watching (and listening) to Bill Benner (president of Pangolin
More informationUNIVERSITY OF BAHRAIN COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING
UNIVERSITY OF BAHRAIN COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EENG 373: DIGITAL COMMUNICATIONS EXPERIMENT NO. 3 BASEBAND DIGITAL TRANSMISSION Objective This experiment
More informationTHE DIGITAL DELAY ADVANTAGE A guide to using Digital Delays. Synchronize loudspeakers Eliminate comb filter distortion Align acoustic image.
THE DIGITAL DELAY ADVANTAGE A guide to using Digital Delays Synchronize loudspeakers Eliminate comb filter distortion Align acoustic image Contents THE DIGITAL DELAY ADVANTAGE...1 - Why Digital Delays?...
More informationUnderstanding Compression Technologies for HD and Megapixel Surveillance
When the security industry began the transition from using VHS tapes to hard disks for video surveillance storage, the question of how to compress and store video became a top consideration for video surveillance
More informationFor the SIA. Applications of Propagation Delay & Skew tool. Introduction. Theory of Operation. Propagation Delay & Skew Tool
For the SIA Applications of Propagation Delay & Skew tool Determine signal propagation delay time Detect skewing between channels on rising or falling edges Create histograms of different edge relationships
More informationUC San Diego UC San Diego Previously Published Works
UC San Diego UC San Diego Previously Published Works Title Classification of MPEG-2 Transport Stream Packet Loss Visibility Permalink https://escholarship.org/uc/item/9wk791h Authors Shin, J Cosman, P
More informationExample: compressing black and white images 2 Say we are trying to compress an image of black and white pixels: CSC310 Information Theory.
CSC310 Information Theory Lecture 1: Basics of Information Theory September 11, 2006 Sam Roweis Example: compressing black and white images 2 Say we are trying to compress an image of black and white pixels:
More informationjamaseis Guide for Displaying Seismic Data
jamaseis Guide for Displaying Seismic Data jamaseis is free cross-platform software that allows users to collect and view near real-time seismic data from a local sensor or from around the world. STEP
More informationMore Info at Open Access Database Process Control for Computed Tomography using Digital Detector Arrays
Digital Industrial Radiology and Computed Tomography (DIR 2015) 22-25 June 2015, Belgium, Ghent - www.ndt.net/app.dir2015 More Info at Open Access Database www.ndt.net/?id=18082 Process Control for Computed
More informationCOSC3213W04 Exercise Set 2 - Solutions
COSC313W04 Exercise Set - Solutions Encoding 1. Encode the bit-pattern 1010000101 using the following digital encoding schemes. Be sure to write down any assumptions you need to make: a. NRZ-I Need to
More informationAssessing and Measuring VCR Playback Image Quality, Part 1. Leo Backman/DigiOmmel & Co.
Assessing and Measuring VCR Playback Image Quality, Part 1. Leo Backman/DigiOmmel & Co. Assessing analog VCR image quality and stability requires dedicated measuring instruments. Still, standard metrics
More information