X-Ray Contrast Review

1. X-Ray Contrast Review Photoelectric Effect and Compton Scatter are main causes of attenuation. µ = pNg {f() + Cp (z3.8/3.2)} where Ng = NA/2 for most elements; NA for hydrogen and f() ≈.6 x 10-24 exp [.0028 (-30)] Cp = 0.8*10-24

2. r  d Source Intensity Review Io = KN/(4πd2) Ii = Io cos3  Source Detector Plane cos() - Obliquity cos2  - Inverse Square Law

3. Depth Dependent Magnification d dr xd,yd µ (x,y,z) z An incremental path of the x-ray, dr, can be described by its x, y, and z components. X-ray path is longer off axis

4. To study, let’s make intensity expression parametric in z. Each point in the body (x,y) can be defined in terms of the detector coordinates it will be imaged at. We can write the response in terms of detector plane coordinates Id (xd,yd) = Ii exp [ - ∫ µo ((xd/d)z, (yd/d)z, z) dz] Putting it all together gives, Id = Io /(((1 + rd2)/d2)3/2) exp [ - ∫ µo ((xd/d)z, (yd/d)z, z) dz] source obliquity object obliquity

5. Two examples Id (xd,yd) y Source z x L z0 Object µ (x,y,z) = µo rect ((z - zo)/L) Object is not a function of x or y, just z. Id (xd,yd) = Ii exp [ - µo L]

6. Id (xd,yd) = Ii exp [ - µo L] If we assume detector is entirely in the near axis, rd2 << d2 (paraxial approximation) Then, simplication results, Id = Io e -µ0L

7. W Id (xd,yd) µ0 0 L x out of plane Infinite in x z0 d Example 2 For µo (x,y,z) = µoP (y/L) P((z - z0)/w) Find the intensity on the detector plane Three cases: 1. Blue Line: X-Ray goes through entire object 2. Red Line: X-ray misses object completely 3. Orange Line: X-ray partially goes through object Id (xd,yd) = Ii exp [ - ∫ µo ((xd/d)z, (yd/d)z, z) dz]

8. W Id (x,y) µ0 0 L x out of plane Infinite in x z0 d Example 2 Id (xd,yd) = Ii exp [ - ∫ µo ((xd/d)z, (yd/d)z, z) dz] For the blue line, we don’t have to worry that the path length through the object will increase as rd increases. That is taken care of by the obliquity term For the red path, Id (xd,yd) = Ii For the orange path, the obliquity term will still help describe the lengthened path. But we need to know the limits in z to integrate

9. µo (x,y,z) = µoP (y/L) P((z - z0)/w) w Id (xd,yd) 0 L x out of plane Infinite in x z0 d Example 2 If we think of thin planes along z, each plane will form a rect in yd of width dL/z. Instead of seeing this as a P in y, let’s mathematically consider it as a P in z that varies in width according to the detector coordinate yd. Then we have integration only in the variable z. The P define limits of integration

10. P in y dimension P in z dimension z z0 z0 + w/2 z0 - w/2 -dL/2|yd| dL/2|yd| X-ray misses object completely As yd grows,first P contracts and no overlap exists between the P , functions. No overlap case when dL/2yd <z0 – w/2 |yd | > dL/(2 zo -w ) Id = Ii

11. P in y dimension P in z dimension z z0 z0 + w/2 z0 - w/2 -dL/2|yd| dL/2|yd| X-ray goes completely through object As yd 0, x-ray goes completely through object This is true for dL/2yd >z0 + w/2 |yd |< dL/ (2 zo + w) Id = Ii exp {-µo w}

12. P in y dimension P in z dimension z z0 z0 + w/2 z0 - w/2 -dL/2|yd| dL/2|yd| • Partial overlap case. Picture above gives the limits of integration. min { zo + w/2 dL/2|yd| Id = Ii exp {-µo ∫ - dz} • min {zo-w/2 • dL/2|yd|

13. 3) Id = Ii exp {-µo ((dL/ 2|yd | ) + w/2 - zo)} X-rays miss object X-rays partially travel through object X-rays go entirely through object The above diagram ignores effects of source obliquity and the factor in the exponential How would curve look differently if we accounted for both of these?

14. Thin section Analysis See object as an array of planes µ (x,y,z) =  (x,y)  (z - z0) Analysis simplifies since only one z plane Id = Ii exp {-  (xd/M , yd/M) } where M = d/z0 If we ignore obliquity, Id = I0 exp {-  (xd/M , yd/M) } Or in terms of the notation for transmission, t, Id = I0 t (xd/M , yd/M) Note: No resolution loss yet. Point remains a point.