COMPUTATIONAL IMAGING. Berthold K.P. Horn
|
|
- Marvin Heath
- 5 years ago
- Views:
Transcription
1 COMPUTATIONAL IMAGING Berthold K.P. Horn
2 What is Computational Imaging? Computation inherent in image formation
3 What is Computational Imaging? Computation inherent in image formation (1) Computing is getting faster and cheaper precision physical apparatus is not
4 What is Computational Imaging? Computation inherent in image formation (1) Computing is getting faster and cheaper precision physical apparatus is not (2) Can t refract or reflect some radiation
5 What is Computational Imaging? Computation inherent in image formation (1) Computing is getting faster and cheaper precision physical apparatus is not (2) Can t refract or reflect some radiation (3) Detection is at times inherently coded
6 Computational Imaging System
7 Examples of Computational Imaging: (1) Synthetic Aperture Imaging (2) Coded Aperture Imaging (3) Exact Cone Beam Reconstruction (4) Diaphanography DiffuseTomography
8 (1) SYNTHETIC APERTURE IMAGING Traditional approach: Coupling of resolution, DOF, FOV to NA Precision imaging flat illumination with: Michael Mermelstein, Jekwan Ryu, Stanley Hong, and Dennis Freeman
9 Objective Lens Parameter Coupling
10 Synthetic Aperture Imaging Traditional approach: Coupling of resolution, DOF, FOV to NA Precision imaging flat illumination New approach: Precision illumination Simple imaging Multiple images Textured illumination
11
12
13
14
15 Synthetic Aperture Imaging Precision illumination Simple imaging Multiple images Textured illumination Image detail in response to textures Non-uniform samples in FT space
16 SAM M6
17 Creating Interference Pattern
18 Reflective Optics M6
19 Creating Interference Pattern
20 Fourier Transform of Texture Pattern
21 Interference Pattern Texture
22 Synthetic Aperture Microscopy Interference of many Coherent Beams Amplitude and Phase Control of Beams
23 Amplitude and Phase Control
24 Synthetic Aperture Microscopy Interference of many Coherent Beams Amplitude and Phase Control of Beams On the fly calibration Non-uniform inverse FT Least Squares
25 Wavenumber Calibration using FT
26 Hough Transform Calibration
27 Least Squares Match in FT
28 Polystyrene Micro Beads (1µm)
29
30 (2) CODED APERTURE IMAGING Can t refract or reflect gamma rays Pinhole tradeoff resolution and SNR with: Richard Lanza, Roberto Accorsi, Klaus Ziock, and Lorenzo Fabris.
31 Coded Aperture Imaging Can t refract or reflect gamma rays Pinhole tradeoff resolution and SNR Multiple pinholes Complex masks can cast shadows
32 Coded Aperture Principle
33 Decoding Method Rationale
34 Coded Aperture Imaging Can t refract or reflect gamma rays Pinhole tradeoff resolution and SNR Complex masks can cast shadows Decoding by Correlation Special Masks with Flat Power Spectrum
35 Mask Design Inverse Systems
36 Masks XRT Coarse
37 Mask Design 1D Definition: q is a quadratic residue (mod p) if n s.t. n 2 q(mod p) Legendre symbol ( a ) p = 1 { 1 if a is quadratic residue otherwise Correlation with zero shift (p 1)/2 Correlation with non-zero shift (p 1)/4
38 Mask Design Auto Correlation a(i) = (p 1) 4 (1 + δ(i)) Power Spectrum A(j) = (p 1) (δ(j) + 1) 4
39 Coded Aperture Extensions Imaging Nearby Objects Mask / Countermask Combination * Dynamic Reconstruction
40 Coded Aperture Application Detection of Fissile Material Large Area Detector Myth Signal and Background Amplified
41 Large Area Alone Doesn t Help
42 Imaging and Large Area Do!
43 Coded Aperture Example Imaging 1/R instead of 1/R 2
44 Coded Aperture Detector Array
45 Computational Imaging System
46
47
48 Dynamic Reconstruction Three weak, distant radioactive sources Reconstruction Animation
49 Coded Aperture Applications Detection of Fissile Material Imaging 1/R instead of 1/R 2 Increasing Gamma Camera Resolution Replacing Rats with Mice.
50 (4) EXACT CONE BEAM ALGORITHM Faster Scanning Fewer Motion Artifacts Lower Exposure Uniform Resolution with: Xiaochun Yang
51 Exact Cone Beam Reconstruction Faster Scanning Fewer Motion Artifacts Lower Exposure Uniform Resolution Parallel Beam Fan Beam Planar Fan Cone Beam
52 Parallel Beam to Fan Beam Coordinate Transform in 2D Radon Space
53 Cone Beam Geometry 3D
54 Radon s Formula In 2D: ~ derivatives of line integrals In 3D: derivatives of plane integrals Can t get plane integrals from projections ( ) f (r, θ)dr dθ 1 f (x, y) dx dy r
55 Radon s Formula in 3D f(x) = 1 8π 2 2 R f (l, β) S 2 l 2 l=x β dβ where R f (l, β) = f(x) δ(x β l)dv
56 Grangeat s Trick z f (x, y, z) dx dy = f (r, φ, θ) dr dφ θ
57 Exact Cone Beam Reconstruction Data Sufficiency Condition Good Orbit for Radiation Source
58 Radon Space 2D
59 Circular Orbit is Inadequate (3D)
60 Data Insufficiency
61 Good Source Orbit
62 Exact Cone Beam Reconstruction Data Sufficiency Condition Good Orbit for Radiation Source Practical Issue: Spiral CT Scanners Practical Issue: Long Body Problem.
63 (3) DIAPHANOGRAPHY (Diffuse Optical Imaging) Highly Scattering Low Absorption Many Sources Many Detectors with: Xiaochun Yang, Richard Lanza, David Boas and Anna Custo.
64
65
66
67 Diaphanography Randomization of Direction Scalar Flux Density
68 30 Head coronal slice scalp mm skull CSF gray mm white
69 Total fluence [CW] Log 10 Fluence source detector separation [mm]
70 scalp skull Normalized PPF brain source detector separation [mm]
71 mm 20 Total fluence [TD] Log 10 TPSF mm mm time steps [ns]
72 Diaphanography Approximation: Diffusion Equation v(x,y) + ρ(x,y)c(x,y) = 0 v(x,y) flux density ρ(x, y) scattering coefficient c(x,y) absorption coefficient Forward: given c(x,y) find v(x,y)
73
74 3 3 node grid example Here we have three input nodes (1, 2, 3), three output nodes (7, 8, 9), and three interior nodes (4, 5, 6). In three experiments we apply currents to each of the input nodes in turn, each time reading out all of the output nodes, yielding a total of nine measurements. We try and recover the nine leakage conductances to ground from each of the nine nodes. It is natural to partition the conductance matrix as follows given that I 4, I 5, I 6, I 7, I 8, and I 9 are always zero, and that we do not measure V 1, V 2, V 3, V 4, V 5, and V 6.. G... V 11 G12.. G13 3. G... V 21 G G23 V 5 = V. G G32.. G33 V 7 V 8 V 9 V 1 V 2 I 1 I 2 I 3 I 4 I 5 I 6 I 7 I 8 I 9
75 . g 1 g 12 g.. 13 g g 12 g 2 g g g 13 g 23 g g g g 4 g 45 g.. 46 g g g 45 g 5 g g g.. 36 g46 g 56 g g g g 7 g 78 g g g 78 g 8 g g 69. g 79 g 89 g 9 V 1 V 2 V 3 V 4 V 5 = V 6 V 7 V 8 V 9 I 1 I 2 I 3 I 4 I 5 I 6 I 7 I 8 I 9
76 We can write the inverse as follows:. C C12.. C13. C C22.. C23. C C32.. C33 I 1 I 2 I 3 I 4 I 5 I 6 I 7 I 8 I 9 V 1 V 2 V 3 V 4 = V 5 V 6 V 7 V 8 V 9
77
78 c(k)(r) FORWARD SOLUTION I(r 1,r 2 ) I (r 1,r 2 ) UPDATE RULE - "ASSIGNING BLAME"
79 Diaphanography Invert Diffusion Equation Regions of Influence.
80 COMPUTATIONAL IMAGING (1) Synthetic Aperture Imaging (2) Coded Aperture Imaging (3) Exact Cone Beam Reconstruction (4) Diaphanography DiffuseTomography
81 COMPUTATIONAL IMAGING
82 Synthetic Aperture Lithography Create pattern controlled interference Example: Two Dots Example: Straight Line Destructive interference safe zone Example: Bessel Ring.
83 SAM M4
84 Amplitude and Phase Control
85
86 Interference Pattern Texture
87 Fourier Transform of Texture Pattern
88 Uneven Fourier Sampling
89 Resolution Enhancement Reflective Optics Illumination Vaccum UV Short Wavelength
90 Resolution Enhancement Reflective Optics Illumination Vaccum UV Short Wavelength Fluorescence Mode Resolution Determined by Illumination
91 Masks Fresnel Camera
92 Masks Legri URA
93 Masks XRT Fine
94 Masks Hexagonal
95 Maximizing SNR n n min w 2 i subject to w i = 1 i=1 i=1 yields w i = 1 n
96 Spatially Varying Background
97 Coded Aperture Extensions Artifacts due to Finite Distance Mask / Countermask Combination Multiple Detector Array Positions Synthetic Aperture Radiography
98 Dynamic Reconstruction
99 Diaphanography Approximation: Diffusion Equation Leaky Resistive Sheet Analog (2D)
100
101
102 mm mm 40mm 0.8 Normalized PPF brain scalp skull time steps [ns] (a) Deviation of Sensitivity to Scalp Skull (PPF PPF o )/PPF o 0 20 mm mm mm time steps [ns] (b) (PPF PPF o )/PPF o mm mm Deviation of Sensitivity to Brain mm time steps [ns] (c)
103 Sample Resistive Grid Inversion Berthold K.P. Horn 1996 February 14th Here we consider about the simplest possible case of the two-d resistive sheet to get some insight into the more general problem. (draw your own picture here :) There are four nodes, arranged in N = 2 rows and M = 2 columns. On the left are the two nodes used for input ((1) and (2)) and on the right are the two output nodes ((3) and (4)). Four horizontal resistors with conductance g 12, g 13, g 24 and g 34 connect these four nodes. These resistors represent scattering, and are assumed to be of known value. There is also a vertical leakage path from each of the four nodes to ground with conductances g 1, g 2, g 3 and g 4. These represent absorption, and are the unknowns. We are to recover the values of the unknown leakage resistors. We perform two experiments: First we inject current at node (1) and measure the potentials on nodes (3) and (4). Call the trans-impedance (ratio of output potential to injected current) observed this way R 3,1 and R 4,1. Then we inject instead current at node (2) and again measure the potential on nodes (3) and (4). Call the trans-impedance observed this way R 3,2 and R 4,2. If the grid was N M instead of 2 2 then we would have performed N experiments, each time injecting current on one of the N input nodes and reading out the potential on each of the N output nodes. We then try to recover the N M unknown leakage conductances to ground. Clearly there is not enough constraint if M>Nsince there are then more unknowns than measurements. Conversely if M<N, we have redundant information and may want to use a least squares approach to obtain the best possible answer. Here we deal with the simple case where M = N = 2. The node equations in this case are: I 1 = g 1 V 1 + g 13 (V 1 V 3 )+g 12 (V 1 V 2 ) I 2 = g 2 V 2 + g 12 (V 2 V 1 )+g 24 (V 2 V 4 ) I 3 = g 3 V 3 + g 13 (V 3 V 1 )+g 34 (V 3 V 4 ) I 4 = g 4 V 4 + g 24 (V 4 V 2 )+g 34 (V 4 V 3 ) or. (g 1 + g 13 + g 12 ) g.. 12 g13 0 g 12 (g 2 + g 12 + g 24 ). 0 g 24 g (g 3 + g 13 + g 34 ) g 34 0 g g34 (g 4 + g 24 + g 34 ) V 1 V 2 = V 3 V 4 Note that all off diagonal elements are negative, and that the unknown leakage conductances all appear on the diagonal. Also, the matrix becomes singular if all leakages conductances are set to zero, since then each row adds up to zero. Making use of the partitioning indicated above, we can write V 1 G 11. G 12 V 2 = V G 21. G 3 22 V 4 This partitioning is convenient, since in the experiments we always have I 3 = 0 and I 4 = 0, and since V 1 and V 2 are not known. If we invert this set of equations we get: C 11. C 12 C 21. C 22 1 I 1 I 2 I 3 I 4 I 1 I 2 I 3 I 4 V 1 V = 2 V 3 V 4 I 1 I 2 I 3 I 4
104 where Upon multiplying out we obtain C 11 =(G 11 G 12 G 1 22 G 21) 1 C 21 = G 1 22 G 21C 11 C 22 =(G 22 G 21 G 1 11 G 12) 1 C 12 = G 1 11 G 12C 22 ( ) I1 C 11 = I 2 where the term involving C 12 drops out because I 3 = 0 and I 4 = 0. This equation is of no interest since we don t know V 1 and V 2. But we also obtain ( ) I1 C 21 = I 2 where the term involving C 22 drops out because I 3 = 0 and I 4 = 0. The four unkowns (g 1,g 2,g 3,g 4 ) occur in the matrix C 21, but we obviously can t solve for them using a single set of measurements. We can however, combine two sets of measurements and obtain: ( ) ( ) I1,1 I C 1,2 V3,1 V 21 = 3,2 I 2,1 I 2,2 V 4,1 V 4,2 where typically we would choose I 1,1 =1,I 2,1 = 0 for the first experiment and I 1,2 =0,I 2,2 = 1 for the second. We then obtain ( ) R3,1 R C 21 = 3,2 R 4,1 R 4,2 where, provided I 2,1 = 0 and I 1,2 =0,R 3,1 = V 3,1 /I 1,1, R 4,1 = V 4,1 /I 1,1, and R 3,2 = V 3,2 /I 1,2, R 4,2 = V 4,2 /I 1,2. So from image measurements we can recover the matrix C 21, and we know that or then, since we obtain ( V1 V 2 ( V3 V 4 ) ) C 21 = G 1 22 G 21C 11 C 21 C 1 11 = G 1 22 G 21 C 11 =(G 11 G 12 G 1 22 G 21) 1 C 21 (G 11 G 12 G 1 22 G 21) = G 1 22 G 21 We need to manipulate this some more to try and isolate the two matrices G 11 and G 22, which contain the unknowns g 1, g 2, g 3, and g 4. We see that or or finally C 21 G 11 = C 21 G 12 G 1 22 G 21 G 1 22 G 21 C 21 G 11 G 1 21 G 22 = C 21 G 12 I G 11 G 1 21 G 22 = G 12 C21 1 Here C 1 21 is obtained from experimental measurements, while G 11 and G 22 contain the unknown leakage conductances. In our simple example, G 12 and G 21 are diagonal. 2
105 Numerical Example Suppose that the horizontal resistors have conductance as follows: g 13 =1,g 12 =2,g 24 = 1, and g 34 =2. Next assum that the leakage or vertical resistors have conductance g 1 =1,g 2 =2,g 3 = 2, and g 4 =1. Then the conductance matrix is and so and so while so and So in the inverse we have and Finally G 1 11 = 1 16 C 11 = C 21 = ( ) and G 12 G 1 22 G 21 = G 1 22 = 1 16 G 21 G 1 11 G 12 = G 1 11 = 1 16 G 11 G 12 G 1 22 G 21 = 1 16 G 22 G 21 G 1 11 G 12 = 1 16 ( 75 ) ( 23 ) G 1 = C = G 1 22 = 1 16 ( ) ( ) ( ) ( 60 ) ( 75 ) and C 22 = and C 12 = ( 60 ) ( 23 ) In the first experiment we have I 1 = 1 and the other node currents are zero so V 1 75 V 2 = V V 4 20 In the second experiment we have I 2 = 1 and the other node currents are zero so V 1 34 V 2 = V V
106 Note that we can only measure V 3 and V 4 in each case. This is the end of the forward problem (finding trans-impedance given leakage conductances). The inverse task is to recover the unknown leakage conductances. Extracting the relevant parts from the above experimental data we see that C 21 = ( ) so so ( ) C = ( ) G 12 C = and G 11 = ( g ) 2 g 2 +3 and G 22 = ( ) g g 4 +3 So ( )( ) G 11 G 1 21 G g g =. 2 g g 4 +3 So that we get the following equations in the unknown leakage conductances: (g 1 + 3)(g 3 +3)+4=24 2(g 1 +3)+2(g 4 +3)=16 2(g 3 +3)+2(g 2 +3)=20 (g 2 + 3)(g 4 +3)+4=24 or ḡ 1 ḡ 3 =20 ḡ 1 +ḡ 4 =8 ḡ 2 +ḡ 3 =10 ḡ 2 ḡ 4 =20 where ḡ 1 = g 1 +3, ḡ 2 = g 2 +3, ḡ 3 = g 3 + 3, and ḡ 4 = g These equations have only one solution: ḡ 1 =4, ḡ 2 =5,ḡ 3 = 5, and ḡ 4 = 4, that is g 1 =1,g 2 =2,g 3 = 2, and g 4 =1. Summary While this shows a solution method for a 2 2 grid, some of the points noted here also apply to the more general case, although an explicit solution cannot be expected then. In general, the matrix would have to be partitioned into a 3 3 arrangement corresponding to the fact that in addition to input nodes and output nodes there are then also interior nodes. 4
107 We can write the inverse as follows:. C C12.. C13. C C22.. C23. C C32.. C33 I 1 I 2 I 3 I 4 I 5 I 6 I 7 I 8 I 9 V 1 V 2 V 3 V 4 = V 5 V 6 V 7 V 8 V 9 We can use the formula for the inverse of matrix partitioned into four parts twice on this matrix partitioned into nine parts. But it may be a bit much to expect to easily obtain explicit formulae the way we did for the 2 2 case... Note that we are only really interested in the bottom left corner (C 31 ) of the inverse, given that I 4, I 5, I 6, I 7, I 8, and I 9 are always zero, and that we do not measure V 1, V 2, V 3, V 4, V 5, and V 6. Each experiment yields three measurements and thus three equations of the form C 31 I 1 I 2 = V 7 V 8. I 3 V 9 By performing three experiments we can find all nine elements of the matrix C 31. Each of these is a polynomial in the unknown leakage conductances g 1, g 2, g 3, g 4, g 5, g 6, g 7, g 8, and g 9 (or rather, we can cross-multiply to obtain nine such polynomials). The part of the inverse of this conductance matrix that we need is the lower left corner, C 31. Using the decomposition rule for partitioned 2 2 matrices twice, we get C 31 = G 1 33 G 32(G 22 G 23 G 1 33 G 21) 1 G 21 ( G11 G 12 (G 22 G 23 G 1 33 G 21) 1 G 21 ) 1 Note that the term (G 22 G 23 G 1 33 G 21) 1 G 21 appears twice. This can be exploited to save on computation. 6
University of Cyprus. Reflectance and Diffuse Spectroscopy
University of Cyprus Biomedical Imaging and Applied Optics Reflectance and Diffuse Spectroscopy Spectroscopy What is it? from the Greek: spectro = color + scope = look at or observe = measuring/recording
More informationHigh-Resolution. Transmission. Electron Microscopy
Part 4 High-Resolution Transmission Electron Microscopy 186 Significance high-resolution transmission electron microscopy (HRTEM): resolve object details smaller than 1nm (10 9 m) image the interior of
More informationInterference, Diffraction and Fourier Theory. ATI 2014 Lecture 02! Keller and Kenworthy
Interference, Diffraction and Fourier Theory ATI 2014 Lecture 02! Keller and Kenworthy The three major branches of optics Geometrical Optics Light travels as straight rays Physical Optics Light can be
More informationWhat is Image Deblurring?
What is Image Deblurring? When we use a camera, we want the recorded image to be a faithful representation of the scene that we see but every image is more or less blurry, depending on the circumstances.
More informationLecture 9: Introduction to Diffraction of Light
Lecture 9: Introduction to Diffraction of Light Lecture aims to explain: 1. Diffraction of waves in everyday life and applications 2. Interference of two one dimensional electromagnetic waves 3. Typical
More informationMaximal Entropy for Reconstruction of Back Projection Images
Maximal Entropy for Reconstruction of Back Projection Images Tryphon Georgiou Department of Electrical and Computer Engineering University of Minnesota Minneapolis, MN 5545 Peter J Olver Department of
More informationA Brief Introduction to Medical Imaging. Outline
A Brief Introduction to Medical Imaging Outline General Goals Linear Imaging Systems An Example, The Pin Hole Camera Radiations and Their Interactions with Matter Coherent vs. Incoherent Imaging Length
More informationGetting Started with Communications Engineering
1 Linear algebra is the algebra of linear equations: the term linear being used in the same sense as in linear functions, such as: which is the equation of a straight line. y ax c (0.1) Of course, if we
More informationTHE SIMULATION OF THE CODED NEUTRON SOURCE PHASE CONTRAST IMAGING
The 12 th International Conference of the Slovenian Society for Non-Destructive Testing»Application of Contemporary Non-Destructive Testing in Engineering«September 4-6, 2013, Portorož, Slovenia More info
More informationPhase Contrast. Zernike Phase Contrast (PC) Nomarski Differential Interference Contrast (DIC) J Mertz Boston University
Phase Contrast Zernike Phase Contrast (PC) Nomarski Differential Interference Contrast (DIC) J Mertz Boston University Absorption Scattering / n wave-number k 2 n = index of refraction / n wave-number
More informationLecture 11: Introduction to diffraction of light
Lecture 11: Introduction to diffraction of light Diffraction of waves in everyday life and applications Diffraction in everyday life Diffraction in applications Spectroscopy: physics, chemistry, medicine,
More informationTomography and Reconstruction
Tomography and Reconstruction Lecture Overview Applications Background/history of tomography Radon Transform Fourier Slice Theorem Filtered Back Projection Algebraic techniques Measurement of Projection
More informationON A CLASS OF GENERALIZED RADON TRANSFORMS AND ITS APPLICATION IN IMAGING SCIENCE
ON A CLASS OF GENERALIZED RADON TRANSFORMS AND ITS APPLICATION IN IMAGING SCIENCE T.T. TRUONG 1 AND M.K. NGUYEN 2 1 University of Cergy-Pontoise, LPTM CNRS UMR 889, F-9532, France e-mail: truong@u-cergy.fr
More informationProblem Solving. radians. 180 radians Stars & Elementary Astrophysics: Introduction Press F1 for Help 41. f s. picture. equation.
Problem Solving picture θ f = 10 m s =1 cm equation rearrange numbers with units θ factors to change units s θ = = f sinθ fθ = s / cm 10 m f 1 m 100 cm check dimensions 1 3 π 180 radians = 10 60 arcmin
More informationPOLARIZATION OF LIGHT
POLARIZATION OF LIGHT OVERALL GOALS The Polarization of Light lab strongly emphasizes connecting mathematical formalism with measurable results. It is not your job to understand every aspect of the theory,
More informationNature of Light Part 2
Nature of Light Part 2 Fresnel Coefficients From Helmholts equation see imaging conditions for Single lens 4F system Diffraction ranges Rayleigh Range Diffraction limited resolution Interference Newton
More informationProbing the orbital angular momentum of light with a multipoint interferometer
CHAPTER 2 Probing the orbital angular momentum of light with a multipoint interferometer We present an efficient method for probing the orbital angular momentum of optical vortices of arbitrary sizes.
More informationTransmission Electron Microscopy
L. Reimer H. Kohl Transmission Electron Microscopy Physics of Image Formation Fifth Edition el Springer Contents 1 Introduction... 1 1.1 Transmission Electron Microscopy... 1 1.1.1 Conventional Transmission
More informationMath (P)Review Part I:
Lecture 1: Math (P)Review Part I: Linear Algebra Computer Graphics CMU 15-462/15-662, Fall 2017 Homework 0.0 (Due Monday!) Exercises will be a bit harder / more rigorous than what you will do for the rest
More informationSOFT X-RAYS AND EXTREME ULTRAVIOLET RADIATION
SOFT X-RAYS AND EXTREME ULTRAVIOLET RADIATION Principles and Applications DAVID ATTWOOD UNIVERSITY OF CALIFORNIA, BERKELEY AND LAWRENCE BERKELEY NATIONAL LABORATORY CAMBRIDGE UNIVERSITY PRESS Contents
More informationIntroduction to FT-IR Spectroscopy
Introduction to FT-IR Spectroscopy An FT-IR Spectrometer is an instrument which acquires broadband NIR to FIR spectra. Unlike a dispersive instrument, i.e. grating monochromator or spectrograph, an FT-IR
More informationInitial Results on the Feasibility of Hybrid X-Ray Microscopy
CHINESE JOURNAL OF PHYSICS VOL. 43, NO. 5 OCTOBER 2005 Initial Results on the Feasibility of Hybrid X-Ray Microscopy P. K. Tseng, 1 W. F. Pong, 1 C. L. Chang, 1 C. P. Hsu, 1 F. Y. Lin, 2 C. S. Hwang, 2
More informationMidterm Review. Yao Wang Polytechnic University, Brooklyn, NY 11201
Midterm Review Yao Wang Polytechnic University, Brooklyn, NY 11201 Based on J. L. Prince and J. M. Links, Medical maging Signals and Systems, and lecture notes by Prince. Figures are from the textbook.
More information30.5. Iterative Methods for Systems of Equations. Introduction. Prerequisites. Learning Outcomes
Iterative Methods for Systems of Equations 0.5 Introduction There are occasions when direct methods (like Gaussian elimination or the use of an LU decomposition) are not the best way to solve a system
More informationMERIS US Workshop. Instrument Characterization Overview. Steven Delwart
MERIS US Workshop Instrument Characterization Overview Steven Delwart Presentation Overview On-Ground Characterisation 1. Diffuser characterisation 2. Polarization sensitivity 3. Optical Transmission 4.
More informationDesigning Information Devices and Systems I Spring 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way
EECS 16A Designing Information Devices and Systems I Spring 018 Lecture Notes Note 1 1.1 Introduction to Linear Algebra the EECS Way In this note, we will teach the basics of linear algebra and relate
More informationRich Tomography. Bill Lionheart, School of Mathematics, University of Manchester and DTU Compute. July 2014
Rich Tomography Bill Lionheart, School of Mathematics, University of Manchester and DTU Compute July 2014 What do we mean by Rich Tomography? Conventional tomography reconstructs one scalar image from
More informationFitting functions to data
1 Fitting functions to data 1.1 Exact fitting 1.1.1 Introduction Suppose we have a set of real-number data pairs x i, y i, i = 1, 2,, N. These can be considered to be a set of points in the xy-plane. They
More informationFoundations of Image Science
Foundations of Image Science Harrison H. Barrett Kyle J. Myers 2004 by John Wiley & Sons,, Hoboken, 0-471-15300-1 1 VECTORS AND OPERATORS 1 1.1 LINEAR VECTOR SPACES 2 1.1.1 Vector addition and scalar multiplication
More informationNetwork Theory and the Array Overlap Integral Formulation
Chapter 7 Network Theory and the Array Overlap Integral Formulation Classical array antenna theory focuses on the problem of pattern synthesis. There is a vast body of work in the literature on methods
More informationQ2 (Michelson) - solution here
The TA is still preparing the solutions for PS#4 and they should be ready on Sunday or early Monday. Meanwhile here are some materials and comments from me. -RSM Q (Michelson) - solution here some notes/comments
More informationSurface Sensitivity & Surface Specificity
Surface Sensitivity & Surface Specificity The problems of sensitivity and detection limits are common to all forms of spectroscopy. In its simplest form, the question of sensitivity boils down to whether
More informationDesigning Information Devices and Systems I Fall 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way
EECS 16A Designing Information Devices and Systems I Fall 018 Lecture Notes Note 1 1.1 Introduction to Linear Algebra the EECS Way In this note, we will teach the basics of linear algebra and relate it
More informationTopic 4: Waves 4.3 Wave characteristics
Guidance: Students will be expected to calculate the resultant of two waves or pulses both graphically and algebraically Methods of polarization will be restricted to the use of polarizing filters and
More informationUNIVERSITY OF ALABAMA Department of Physics and Astronomy. PH / LeClair Fall Circuits Exercises
UNIVERSITY OF ALABAMA Department of Physics and Astronomy PH 106-4 / LeClair Fall 008 Circuits Exercises 1. Are the two headlights of a car wired in series or in parallel? How can you tell? Have you ever
More informationLecture 20 Optical Characterization 2
Lecture 20 Optical Characterization 2 Schroder: Chapters 2, 7, 10 1/68 Announcements Homework 5/6: Is online now. Due Wednesday May 30th at 10:00am. I will return it the following Wednesday (6 th June).
More informationX-RAY MICRO-TOMOGRAPHY OF PORE-SCALE FLOW AND TRANSPORT. University of California Davis. Dorthe Wildenschild & Annette Mortensen
X-RAY MICRO-TOMOGRAPHY OF PORE-SCALE FLOW AND TRANSPORT Jan W. Hopmans University of California Davis Volker Clausnitzer Dorthe Wildenschild & Annette Mortensen ISSUES: Measurements and modeling of water
More informationMid-IR Methods. Matti Hotokka
Mid-IR Methods Matti Hotokka Traditional methods KBr disk Liquid cell Gas cell Films, fibres, latexes Oil suspensions GC-IR PAS KBr Disk 1 mg sample, 200 mg KBr May be used for quantitative analysis Mix
More information1 cm b. 4.4 mm c. 2.2 cm d. 4.4 cm v
PHY 112: General Physics M. F. Thorpe T, Th 7:40-8:55am Fall 2006 Department of Physics Arizona State University Tempe AZ Final, Friday 8 December from 7:40am -> 9.30am All questions carry equal weight.
More informationImage Reconstruction from Projection
Image Reconstruction from Projection Reconstruct an image from a series of projections X-ray computed tomography (CT) Computed tomography is a medical imaging method employing tomography where digital
More informationVisualization of Xe and Sn Atoms Generated from Laser-Produced Plasma for EUV Light Source
3rd International EUVL Symposium NOVEMBER 1-4, 2004 Miyazaki, Japan Visualization of Xe and Sn Atoms Generated from Laser-Produced Plasma for EUV Light Source H. Tanaka, A. Matsumoto, K. Akinaga, A. Takahashi
More informationDepartment of Physics, Colorado State University PH 425 Advanced Physics Laboratory The Zeeman Effect. 1 Introduction. 2 Origin of the Zeeman Effect
Department of Physics, Colorado State University PH 425 Advanced Physics Laboratory The Zeeman Effect (a) CAUTION: Do not look directly at the mercury light source. It is contained in a quartz tube. The
More informationHow do you make an image of an object?
How do you make an image of an object? Use a camera to take a picture! But what if the object is hidden?...or invisible to the human eye?...or too far away to see enough detail? Build instruments that
More informationA super-high angular resolution principle for coded-mask X-ray imaging beyond the diffraction limit of a single pinhole
Research in Astron. Astrophys. 29 Vol. 9 No. 3, 333 34 http://www.raa-journal.org http://www.iop.org/journals/raa Research in Astronomy and Astrophysics A super-high angular resolution principle for coded-mask
More informationComplex number review
Midterm Review Problems Physics 8B Fall 009 Complex number review AC circuits are usually handled with one of two techniques: phasors and complex numbers. We ll be using the complex number approach, so
More informationA beam of coherent monochromatic light from a distant galaxy is used in an optics experiment on Earth.
Waves_P2 [152 marks] A beam of coherent monochromatic light from a distant galaxy is used in an optics experiment on Earth. The beam is incident normally on a double slit. The distance between the slits
More informationThe Plasma Phase. Chapter 1. An experiment - measure and understand transport processes in a plasma. Chapter 2. An introduction to plasma physics
The Plasma Phase Chapter 1. An experiment - measure and understand transport processes in a plasma Three important vugraphs What we have just talked about The diagnostics Chapter 2. An introduction to
More information홀로그램저장재료. National Creative Research Center for Active Plasmonics Applications Systems
홀로그램저장재료 Holographic materials Material Reusable Processing Type of Exposure Spectral Resol. Max. diff. hologram (J/m2) sensitivity (lim./mm) efficiency Photographic emulsion Dichromated gelatin Photoresists
More informationCarbon Dating The decay of radioactive nuclei can be used to measure the age of artifacts, fossils, and rocks. The half-life of C 14 is 5730 years.
Carbon Dating The decay of radioactive nuclei can be used to measure the age of artifacts, fossils, and rocks. The half-life of C 14 is 5730 years. a) If a sample shows only one-fourth of its estimated
More information5. LIGHT MICROSCOPY Abbe s theory of imaging
5. LIGHT MICROSCOPY. We use Fourier optics to describe coherent image formation, imaging obtained by illuminating the specimen with spatially coherent light. We define resolution, contrast, and phase-sensitive
More informationLecture notes 5: Diffraction
Lecture notes 5: Diffraction Let us now consider how light reacts to being confined to a given aperture. The resolution of an aperture is restricted due to the wave nature of light: as light passes through
More informationActivity: Derive a matrix from input-output pairs
Activity: Derive a matrix from input-output pairs [ ] a b The 2 2 matrix A = satisfies the following equations: c d Calculate the entries of the matrix. [ ] [ ] a b 5 c d 10 [ ] [ ] a b 2 c d 1 = = [ ]
More informationOptics for Engineers Chapter 9
Optics for Engineers Chapter 9 Charles A. DiMarzio Northeastern University Nov. 202 Gaussian Beams Applications Many Laser Beams Minimum Uncertainty Simple Equations Good Approximation Extensible (e.g.
More informationThe graph shows how the distance from their home changes with time. Which part of the graph, A, B, C or D, shows them walking the fastest? ...
YR11 Half Term PHYSICS Paper 2 HIGHER Name:... Attempt ALL Questions and bring the completed test to your first Physics lesson after the holiday. Your Exam Equation Sheet is on the last page. Mark /75
More informationobject objective lens eyepiece lens
Advancing Physics G495 June 2015 SET #1 ANSWERS Field and Particle Pictures Seeing with electrons The compound optical microscope Q1. Before attempting this question it may be helpful to review ray diagram
More informationMEDICAL EQUIPMENT: NUCLEAR MEDICINE. Prof. Yasser Mostafa Kadah
MEDICAL EQUIPMENT: NUCLEAR MEDICINE Prof. Yasser Mostafa Kadah www.k-space.org Recommended Textbook Introduction to Medical Imaging: Physics, Engineering and Clinical Applications, by Nadine Barrie Smith
More informationOptical Tweezers. The Useful Micro-Manipulation Tool in Research
Optical Tweezers The Useful Micro-Manipulation Tool in Research Student: Nikki Barron Class: Modern Physics/ Biophysics laboratory Advisor: Grant Allen Instructor: David Kleinfeld Date: June 15, 2012 Introduction
More informationCopyright 2000 by the Society of Photo-Optical Instrumentation Engineers.
Copyright 2 by the Society of Photo-Optical Instrumentation Engineers. This paper was published in the proceedings of Photomask and X-Ray Mask Technology VII SPIE Vol. 466, pp. 172-179. It is made available
More informationResponse Surface Methods
Response Surface Methods 3.12.2014 Goals of Today s Lecture See how a sequence of experiments can be performed to optimize a response variable. Understand the difference between first-order and second-order
More informationNEON Archive School 2006
NEON Archive School 2006 Introduction to Spectroscopic Techniques (low dispersion) M. Dennefeld (IAP-Paris) Outline Basic optics of gratings and spectrographs (with emphasis on long-slit spectroscopy)
More informationToday. MIT 2.71/2.710 Optics 11/10/04 wk10-b-1
Today Review of spatial filtering with coherent illumination Derivation of the lens law using wave optics Point-spread function of a system with incoherent illumination The Modulation Transfer Function
More informationChapter 4 Nadir looking UV measurement. Part-I: Theory and algorithm
Chapter 4 Nadir looking UV measurement. Part-I: Theory and algorithm -Aerosol and tropospheric ozone retrieval method using continuous UV spectra- Atmospheric composition measurements from satellites are
More informationAstronomy 114. Lecture 26: Telescopes. Martin D. Weinberg. UMass/Astronomy Department
Astronomy 114 Lecture 26: Telescopes Martin D. Weinberg weinberg@astro.umass.edu UMass/Astronomy Department A114: Lecture 26 17 Apr 2007 Read: Ch. 6,26 Astronomy 114 1/17 Announcements Quiz #2: we re aiming
More informationHadamard speckle contrast reduction
Published in OPTICS LETTERS 29, -3 (2004) Hadamard speckle contrast reduction Jahja I. Trisnadi Silicon Light Machines, 385 Moffett Park Drive, Sunnyvale, California 94089 The conditon for a diffuser to
More informationMODERN OPTICS. P47 Optics: Unit 9
MODERN OPTICS P47 Optics: Unit 9 Course Outline Unit 1: Electromagnetic Waves Unit 2: Interaction with Matter Unit 3: Geometric Optics Unit 4: Superposition of Waves Unit 5: Polarization Unit 6: Interference
More informationFXA UNIT G485 Module X-Rays. Candidates should be able to : I = I 0 e -μx
1 Candidates should be able to : HISTORY Describe the nature of X-rays. Describe in simple terms how X-rays are produced. X-rays were discovered by Wilhelm Röntgen in 1865, when he found that a fluorescent
More informationLet us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a).
7.1. Low-Coherence Interferometry (LCI) Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a). The light is split by the beam splitter (BS) and
More informationContact resistance and TLM measurements
Contact resistance and TLM measurements In measuring resistance with the four-point-probe or van der Pauw methods, we used 4 contacts (2 for current, 2 for voltage) to determine the sheet resistance of
More informationMath Methods for Polymer Physics Lecture 1: Series Representations of Functions
Math Methods for Polymer Physics ecture 1: Series Representations of Functions Series analysis is an essential tool in polymer physics and physical sciences, in general. Though other broadly speaking,
More informationFar-field radiation pattern in Coherent Anti-stokes Raman Scattering (CARS) Microscopy.
Far-field radiation pattern in Coherent Anti-stokes Raman Scattering (CARS) Microscopy. David Gachet, Nicolas Sandeau, Hervé Rigneault * Institut Fresnel, Mosaic team, Domaine Univ. St Jérôme, 13397 Marseille
More informationDesigning Information Devices and Systems I Fall 2018 Homework 3
Last Updated: 28-9-5 :8 EECS 6A Designing Information Devices and Systems I Fall 28 Homework 3 This homework is due September 4, 28, at 23:59. Self-grades are due September 8, 28, at 23:59. Submission
More informationSometimes the domains X and Z will be the same, so this might be written:
II. MULTIVARIATE CALCULUS The first lecture covered functions where a single input goes in, and a single output comes out. Most economic applications aren t so simple. In most cases, a number of variables
More informationOptics for Engineers Chapter 9
Optics for Engineers Chapter 9 Charles A. DiMarzio Northeastern University Mar. 204 Gaussian Beams Applications Many Laser Beams Minimum Uncertainty Simple Equations Good Approximation Extensible (e.g.
More informationDiffraction of Electrons
Diffraction of Electrons Object: Apparatus: Verify that electrons are waves; i.e., that they diffract just like light waves. This lab is then used to measure their wavelength or, alternatively, measure
More informationTechnical University of Denmark
Technical University of Denmark Page 1 of 11 pages Written test, 9 December 2010 Course name: Introduction to medical imaging Course no. 31540 Aids allowed: none. "Weighting": All problems weight equally.
More informationEL-GY 6813/BE-GY 6203 Medical Imaging, Fall 2016 Final Exam
EL-GY 6813/BE-GY 6203 Medical Imaging, Fall 2016 Final Exam (closed book, 1 sheets of notes double sided allowed, no calculator or other electronic devices allowed) 1. Ultrasound Physics (15 pt) A) (9
More informationNeutron and Gamma Ray Imaging for Nuclear Materials Identification
Neutron and Gamma Ray Imaging for Nuclear Materials Identification James A. Mullens John Mihalczo Philip Bingham Oak Ridge National Laboratory Oak Ridge, Tennessee 37831-6010 865-574-5564 Abstract This
More informationTechnical University of Denmark
Technical University of Denmark Page 1 of 10 pages Written test, 12 December 2012 Course name: Introduction to medical imaging Course no. 31540 Aids allowed: None. Pocket calculator not allowed "Weighting":
More informationFabry-Perot Interferometer for atmospheric monitoring useful for EAS detection E.Fokitis 1, K. Patrinos 1, Z. Nikitaki 1
Fabry-Perot Interferometer for atmospheric monitoring useful for EAS detection E.Fokitis 1, K. Patrinos 1, Z. Nikitaki 1 ABSTRACT A piezotunable Fabry-Perot interferometer is studied as a candidate Doppler
More information5. 3P PIV Measurements
Micro PIV Last Class: 1. Data Validation 2. Vector Field Operator (Differentials & Integrals) 3. Standard Differential Scheme 4. Implementation of Differential & Integral quantities with PIV data 5. 3P
More informationB 2 P 2, which implies that g B should be
Enhanced Summary of G.P. Agrawal Nonlinear Fiber Optics (3rd ed) Chapter 9 on SBS Stimulated Brillouin scattering is a nonlinear three-wave interaction between a forward-going laser pump beam P, a forward-going
More informationNeutron collimator design of neutron radiography based on the BNCT facility arxiv: v1 [physics.ins-det] 3 May 2013
CPC(HEP & NP), 2009, 33(X): 1 5 Chinese Physics C Vol. 33, No. X, Xxx, 2009 Neutron collimator design of neutron radiography based on the BNCT facility arxiv:1305.0672v1 [physics.ins-det] 3 May 2013 YANG
More information1.1 GRAPHS AND LINEAR FUNCTIONS
MATHEMATICS EXTENSION 4 UNIT MATHEMATICS TOPIC 1: GRAPHS 1.1 GRAPHS AND LINEAR FUNCTIONS FUNCTIONS The concept of a function is already familiar to you. Since this concept is fundamental to mathematics,
More informationAPPLICATION INFORMATION
A-1995A APPLICATION INFORMATION Particle Characterization THE SIZING OF NON-SPHERICAL,SUB-MICRON PARTICLES USING POLARIZATION INTENSITY DIFFERENTIAL SCATTERING (PIDS ) Introduction Nearly all laser-based
More informationName Final Exam May 1, 2017
Name Final Exam May 1, 217 This test consists of five parts. Please note that in parts II through V, you can skip one question of those offered. Some possibly useful formulas appear below. Constants, etc.
More informationCS123 INTRODUCTION TO COMPUTER GRAPHICS. Linear Algebra 1/33
Linear Algebra 1/33 Vectors A vector is a magnitude and a direction Magnitude = v Direction Also known as norm, length Represented by unit vectors (vectors with a length of 1 that point along distinct
More information1. Which of the following statements is true about Bremsstrahlung and Characteristic Radiation?
BioE 1330 - Review Chapters 4, 5, and 6 (X-ray and CT) 9/27/2018 Instructions: On the Answer Sheet, enter your 2-digit ID number (with a leading 0 if needed) in the boxes of the ID section. Fill in the
More informationTHE ZEEMAN EFFECT PHYSICS 359E
THE ZEEMAN EFFECT PHYSICS 359E INTRODUCTION The Zeeman effect is a demonstration of spatial quantization of angular momentum in atomic physics. Since an electron circling a nucleus is analogous to a current
More informationThe Imaging Chain for X-Ray Astronomy
The Imaging Chain for X-Ray Astronomy Pop quiz (1): Which is the X-ray Image? B. A. Answer: B!!! (But You Knew That) Pop quiz (2): Which of These is the X-Ray Image? A. B. C. B. A. The dying star ( planetary
More informationLarbert High School. Quanta and Waves. Homework Exercises ADVANCED HIGHER PHYSICS
Larbert High School ADVANCED HIGHER PHYSICS Quanta and Waves Homework Exercises 3.1 3.6 3.1 Intro to Quantum Theory HW 1. (a) Explain what is meant by term black body. (1) (b) State two observations that
More informationDiffraction I. Physics 2415 Lecture 37. Michael Fowler, UVa
Diffraction I Physics 2415 Lecture 37 Michael Fowler, UVa Today s Topics Michelson s interferometer The Michelson Morley experiment Single-slit diffraction Eye of a fly Angular resolution Michelson Interferometer
More informationPhysics 202 Final Exam May 14, 2012
ID CODE: D Physics 202 Final Exam May 14, 2012 Name:... Student ID:... Section:... TA (please circle): Deepak Agarwal Nicholas Brewer Raghvendra Chaubey Todd Garon Yutao Gong Andrew Loveridge Abhishek
More informationDesign and Correction of optical Systems
Design and Correction of optical Systems Part 10: Performance criteria 1 Summer term 01 Herbert Gross Overview 1. Basics 01-04-18. Materials 01-04-5 3. Components 01-05-0 4. Paraxial optics 01-05-09 5.
More informationLC circuit: Energy stored. This lecture reviews some but not all of the material that will be on the final exam that covers in Chapters
Disclaimer: Chapter 29 Alternating-Current Circuits (1) This lecture reviews some but not all of the material that will be on the final exam that covers in Chapters 29-33. LC circuit: Energy stored LC
More informationSection 22. Radiative Transfer
OPTI-01/0 Geometrical and Instrumental Optics Copyright 018 John E. Greivenkamp -1 Section Radiative Transfer Radiometry Radiometry characterizes the propagation of radiant energy through an optical system.
More informationModern physics ideas are strange! L 36 Modern Physics [2] The Photon Concept. How are x-rays produced? The uncertainty principle
L 36 Modern Physics [2] X-rays & gamma rays How lasers work Medical applications of lasers Applications of high power lasers Medical imaging techniques CAT scans MRI s Modern physics ideas are strange!
More informationElementary Linear Algebra
Matrices J MUSCAT Elementary Linear Algebra Matrices Definition Dr J Muscat 2002 A matrix is a rectangular array of numbers, arranged in rows and columns a a 2 a 3 a n a 2 a 22 a 23 a 2n A = a m a mn We
More informationRadionuclide Imaging MII Positron Emission Tomography (PET)
Radionuclide Imaging MII 3073 Positron Emission Tomography (PET) Positron (β + ) emission Positron is an electron with positive charge. Positron-emitting radionuclides are most commonly produced in cyclotron
More informationScattering Lecture. February 24, 2014
Scattering Lecture February 24, 2014 Structure Determination by Scattering Waves of radiation scattered by different objects interfere to give rise to an observable pattern! The wavelength needs to close
More informationLecture 6. Numerical methods. Approximation of functions
Lecture 6 Numerical methods Approximation of functions Lecture 6 OUTLINE 1. Approximation and interpolation 2. Least-square method basis functions design matrix residual weighted least squares normal equation
More information