1 / 51

Medical Imaging

Simultaneous measurements on a spatial grid. Many modalities: mainly EM radiation and sound. Medical Imaging “To invent you need a good imagination and a pile of junk.” Thomas Edison 1879 Electron rapidly decelerates at heavy metal target, giving off X-Rays. Bremsstrahlung 1896

adamdaniel
Download Presentation

Medical Imaging

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Simultaneous measurements on a spatial grid. Many modalities: mainly EM radiation and sound. Medical Imaging

  2. “To invent you need a good imagination and a pile of junk.” Thomas Edison 1879

  3. Electron rapidly decelerates at heavy metal target, giving off X-Rays. Bremsstrahlung

  4. 1896

  5. Projection of X-Ray silhouette onto a piece of film or detector array, with intervening fluorescent screen. X-Ray and Fluoroscopic Images

  6. From a series of projections, a tomographic image is reconstructed using Filtered Back Projection. Computerized Tomography

  7. Radioactive isotope separated by difference in inertia while bending in magnetic field. Mass Spectrometer

  8. Gamma camera for creating image of radioactive target. Camera is rotated around patient in SPECT (Single Photon Emission Computed Tomography). Nuclear Medicine

  9. Ultrasound beam formed and steered by controlling the delay between the elements of the transducer array. Phased Array Ultrasound

  10. Real Time 3D Ultrasound

  11. Positron-emitting organic compounds create pairs of high energy photons that are detected synchronously. Positron Emission Tomography

  12. MRI (Magnetic Resonance Imaging) OCT (Optical Coherence Tomography) Other Imaging Modalities

  13. 3D Higher speed Greater resolution Measure function as well as structure Combining modalities (including direct vision) Current Trends in Imaging

  14. Dissection: Medical School, Day 1: Meet the Cadaver. From Vesalius to the Visible Human The Gold Standard

  15. Local Operators and Global Transforms

  16. Some things work in n dimensions, some don’t. It is often easier to present a concept in 2D. I will use the word “pixel” for n dimensions. Images are n dimensional signals.

  17. Geometric (rigid body) n translations and rotations. Similarity Add 1 scale (isometric). Affine Add n scales (combined with rotation => skew). Parallel lines remain parallel. Projection Global Transforms in n dimensions

  18. Capable of geometric, similarity, or affine. Homogeneous coordinates. Multiply in reverse order to combine SGI “graphics engine” 1982, now standard. Orthographic Transform Matrix

  19. Translation by (tx , ty) Scale x by sx and y by sy

  20. Rotation in 2D • 2 x 2 rotation portion is orthogonal (orthonormal vectors). • Therefore only 1 degree of freedom, .

  21. Rotation in 3D • 3 x 3 rotation portion is orthogonal (orthonormal vectors). • 3 degree of freedom (dotted circled), , as expected.

  22. Non-Orthographic Projection in 3D • For X-ray or direct vision, projects onto the (x,y) plane. • Rescales x and y for “perspective” by changing the “1” in the homogeneous coordinates, as a function of z.

  23. f is usually monotonic, and shift invariant. Inverse may not exist due to discrete values of intensity. Brightness/contrast, “windowing”. Thresholding. Color Maps. f may vary with pixel location, eg., correcting for inhomogeneity of RF field strength in MRI. Point Operators

  24. A pixel-wise intensity mapping is found that produces a uniform density of pixel intensity across the dynamic range. Histogram Equalization

  25. Assumes bimodal distribution. Trough represents boundary points between homogenous areas. Adaptive Thresholding from Histogram

  26. Assumes registration. Averaging multiple acquisitions for noise reduction. Subtracting sequential images for motion detection, or other changes (eg. Digital Subtractive Angiography). Masking. Algebraic Operators

  27. Can result in denser or sparser pixels. Two general approaches: Forward Mapping (Splatting) Backward Mapping (Interpolation) Nearest Neighbor Bilinear Cubic 2D and 3D texture mapping hardware acceleration. Re-Sampling on a New Lattice

  28. Template matching uses correlation, the primordial form of image analysis. Kernels are mostly used for “convolution” although with symmetrical kernels equivalent to correlation. Convolution flips the kernel and does not normalize. Correlation subtracts the mean and generally does normalize. Convolution and Correlation

  29. Discrete images always requires a specific scale. “Inner scale” is the original pixel grid. Size of the kernel determines scale. Concept of Scale Space, Course-to-Fine. Neighborhood PDE Operators

  30. Vector Direction of maximum change of scalar intensity I. Normal to the boundary. Nicely n-dimensional. Intensity Gradient

  31. Scalar Maximum at the boundary Orientation-invariant. Intensity Gradient Magnitude

  32. Classic Edge Detection Kernel (Sobel)

  33. 100% opaque watertight surface Fast, 28 = 256 combinations, pre-computed Isosurface, Marching Cubes (Lorensen)

  34. Marching cubes works well with raw CT data. • Hounsfield units (attenuation). • Threshold calcium density.

  35. Ixy = Iyx= curvature Orientation-invariant. What about in 3D? Jacobian of the Intensity Gradient

  36. Divergence of the Gradient. Zero at the inflection point of the intensity curve. Laplacian of the Intensity I Ix Ixx

  37. Repeated averaging of neighbors => Gaussian by Central Limit Theorem. Binomial Kernel

  38. Not the conventional concentric DOG Subtracting pixels displaced along the x axis after repeated blurring with binomial kernel yields Ix Binomial Difference of Offset Gaussian (DooG)

  39. Two regions with the same intensity but differentiated by texture are easily discriminated by the human visual system. Texture Boundaries

  40. 2D Fourier Transform analysis or synthesis

  41. Most of the usual properties, such as linearity, etc. Shift-invariant, rather than Time-invariant Parsevals relation becoms Rayleigh’s Theorem Also, Separability, Rotational Invariance, and Projection (see below) Properties

  42. Separability

  43. Rotation Invariance

  44. Projection Combine with rotation, have arbitrary projection.

  45. Gaussian seperable Since the Fourier Transform is also separable, the spectra of the 1D Gaussians are, themselves, separable.

  46. Hankel Transform For radially symmetrical functions

  47. Elliptical Fourier Series for 2D Shape Parametric function, usually with constant velocity. Truncate harmonics to smooth.

  48. Fourier surface of 3D shapes (parameterized on surface). Spherical Harmonics (parameterized in spherical coordinates). Both require coordinate system relative to the object. How to choose? Moments? Problem of poles: sigularities cannot be avoided Fourier shape in 3D

  49. Quaternions – 3D phasors Product is defined such that rotation by arbitrary angles from arbitrary starting points become simple multiplication.

More Related