processamento de imagens mÉdicas -...

37
PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia Rittner [email protected] FEEC/UNICAMP

Upload: truongnhi

Post on 16-Dec-2018

218 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

PROCESSAMENTO DE IMAGENS MÉDICAS

IA369O – Capítulo 2

Leticia Rittner [email protected]

FEEC/UNICAMP

Page 2: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Modalidades existentes •  Radiology •  Sonography •  Mammography •  Genitourology •  DSA •  Digital Fluoroscopy •  MRI •  CT •  Digital Mammography •  Bone Density •  Ultrasounds •  X-Ray •  Breast Thermography •  Angiography

Page 3: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Modalidades existentes

Page 4: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Modalidades ionizantes • Raios-X

•  Convencional •  Fluoroscopia •  Angiografia •  Tomografia Computadorizada

• Raios-γ (Medicina Nuclear) •  Cintilografia •  SPEC •  PET

Page 5: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Modalidades não-ionizantes • Ultrassonografia (US)

•  Convencional •  Doppler

• Ressonância magnética (MRI) •  T1 •  T2 •  Flair •  fMRI •  DTI

Page 6: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

RAIO-X

Page 7: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Raio-X • Detecta fótons de

raios-X que não foram absorvidos durante exposição do paciente ao feixe de raios-X

• 

the scattered x-rays enter the grid at an oblique angle and are absorbed by the lead strips.The grid ratio, h/D, where h is the height of the lead strips and D is their separation,determines that angle and is the primary determinant of the effectiveness of a grid inremoving scattered radiation. Grid ratios vary between about 6:1 and 16:1 in routineradiographic systems, with less scattered radiation reaching the image receptor for largegrid ratios. Since scatter depends on the thickness of the body part, lower grid ratios, or nogrid at all, can be used with thin body parts (e.g. extremities, such as the hand or foot) butlarger grid ratios are used for thick body parts (e.g. abdomen or chest). In mammographicsystems the grid ratio can be as low as 2:1. The simple linear grid has now been largelyreplaced by the focused linear grid, where the lead strips are progressively angled onmoving away from the central axis (Fig. 3.14(ii)). This eliminates the problem of “cut-off” at the periphery of the grid but requires careful placement of the grid so that theangles of the strips match the divergence of the x-rays. In order that the lead strips do notcast a shadow on the image, the grid is usually placed in a special holder (“Bucky”) whichmoves the grid during the exposure to blur this effect.

There is a price to be paid for using grids to reduce scatter in order to preserve contrastin the image. Some primary radiation strikes the lead strips head-on and is absorbed. Inorder to compensate for this loss, the x-ray tube current needs to be increased somewhat,which increases the radiation dose to the patient. This is part of a general underlying

Space between strips

Height oflead strip

LeadAluminum

D

h

! ! ! Patient

Absorption

Grid

Scattering

Absorption

(i) (ii)

! !

Figure 3.14 (i) A parallel grid, shown end on; the lead strips are of height h and separated by a distance D.(ii) A focused grid matches the divergence of the x-ray beam to allow more primary x-rays topass through. (After Wolbarst, 1993, p. 175.)

3.2 Images from x-rays 59

I0 = Ii.e!µ."x

Page 8: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Raio-X Convencional

(i) (ii)

Figure 3.7 Posterior–anterior chest radiographs of (i) a normal patient and (ii) a patient suffering fromtuberculosis.

Focal spot size

Focus–objectdistance

Patient

Object–imagedistance

A B

P

C

Detector plane

Feature

Blur Blur

D

S

Figure 3.8 Blurring at the detector plane results from the x-ray focal spot size and placement of the patient.The closer the patient is to the detectors, the smaller the blurring.

54 Medical images obtained with ionizing radiation

Page 9: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Angiografia

averaged. Frame averaging may work well for static images, at the cost of increasedradiation exposure to the patient, but the increased image lag may be unacceptable forproducing images of a dynamic process. The signal-to-noise ratio is also directlyproportional to the concentration of the contrast medium, but increasing its concentrationinvolves increased risk to the patient also.

The patient should be immobile between the acquisition of the mask and live images,otherwise the images will not be registered properly and motion artifacts, in the form ofwhitish streaks, will be visible in the subtracted image. The solution is to shift the mask

(i) (ii)

(iii)

Figure 3.22 (i) Mask or pre-contrast image; (ii) live or post-contrast images; (iii) the subtracted image(live-mask), clearly showing the blood vessels.

66 Medical images obtained with ionizing radiation

averaged. Frame averaging may work well for static images, at the cost of increasedradiation exposure to the patient, but the increased image lag may be unacceptable forproducing images of a dynamic process. The signal-to-noise ratio is also directlyproportional to the concentration of the contrast medium, but increasing its concentrationinvolves increased risk to the patient also.

The patient should be immobile between the acquisition of the mask and live images,otherwise the images will not be registered properly and motion artifacts, in the form ofwhitish streaks, will be visible in the subtracted image. The solution is to shift the mask

(i) (ii)

(iii)

Figure 3.22 (i) Mask or pre-contrast image; (ii) live or post-contrast images; (iii) the subtracted image(live-mask), clearly showing the blood vessels.

66 Medical images obtained with ionizing radiation

Máscara (sem contraste) Imagem (com contraste) Vasos (subtração)

•  Fluoroscopia com contraste e subtração em tempo-real

Page 10: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Angiografia – Artefato de movimento

image by a few pixels, to obviate the motion, before subtraction (Fig. 3.23). This pixel-shifting tends to be a trial-and-error process, involving a combination of shifts in differentdirections and by differing amounts. Motion artifacts can be a significant problem incardiac studies, resulting from the involuntary motion of the soft tissues.

3.2.5 Computed tomography

Conventional radiographic procedures, those which we have described so far, produceplanar images that are projections of three-dimensional objects onto two-dimensionalplanes. This results in a considerable loss of information. The superpositioning of over-lying organs complicates their identification, unequal magnification effects cause distor-tion, and x-ray scattering can result in poor dynamic range. Tomographic imaging, ofwhich x-ray computed tomography (CT) is an example, is a technique that was developedfor producing transverse images, by scanning a slice of tissue from multiple directionsusing a narrow fan-shaped beam. The data from each direction comprise a one-dimen-sional projection of the object, and a transverse image can be retrospectively recon-structed from multiple projections. The body can be compared to a loaf of sliced bread,and a transverse image can be produced as if it were a selected slice viewed face-on(Fig. 3.24). The slice thickness can be reduced to 1mm or so, so that very little super-positioning occurs. Indeed, if many transverse images are obtained, the data can bepresented as an image in any plane, or even as a three-dimensional composite image.

A computed tomography (CT) scanner looks like a big, square doughnut. The patientis placed within the aperture of the rotating frame or gantry (Fig. 3.25(i)). The geometryof scanning has been changing with time, with each major step forward being called a

(i) (ii)

Figure 3.23 (i) Motion artifacts appearing as white streaks in the image. (ii) Minimization of motion artifactsby pixel-shifting.

3.2 Images from x-rays 67

Page 11: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Tomografia computadorizada

• Dispositivo de radiografia giratório (gantry) + computador

the gantry rotated by a small angle (!1°), and another scan performed. This is repeated togive scans covering 180°. The sampling interval and the angle of rotation determine thepixel size in the reconstructed transverse image, and the collimator width determines theslice thickness. The motion is known as translate-rotate.

A transverse slice of the body is schematically divided into many small volumeelements or voxels (Fig. 3.27). They are displayed on a monitor as an array of pixels,where the third dimension, the slice thickness, has been flattened. Using the first-generation scanning as the model, the measured x-ray intensity depends on the sum of

X-ray tube

Collimator

Detector

1st scan

60thscan

(i) (ii)

Figure 3.26 Operation of a first-generation CTscanner. (i) Translatemotion comprising a single scan. (ii) Differentscans are taken after a rotation of the beam and detector. (After Wolbarst, 1993, p. 322.)

l0

µ1,1 µ1,2

µ2,2µ2,1

µm,1 µm,2 µm,n

µ2,n

µ1,n–! µi, j

ti = 1

nI = l0 [ ]exp( )

–! µi, j t

i = 1

mI = l0 [ ]exp( )

–! µi,1 t

i = 1

mI = l0 [ ]exp( )

Figure 3.27 Matrix of voxels comprising the patient, with x-rays shown entering and leaving atdiffering angles.

3.2 Images from x-rays 69

Page 12: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Tomografia computadorizada • Equaçoes de atenuação

• Métodos de reconstrução •  Método iterativo •  Retroprojeção filtrada

the attenuation coefficients (!i,j) for each of the voxels in that particular path. Substitutingthe sum of the attenuation coefficients along a line of voxels (!i,j) into Equation (3.3) andtaking the natural logarithm gives the ray sum, p, along this particular path:

p ! "ln#I0=I$ ! !! i;jt (3:7)

The reference intensity, I0, is measured for each detector in a calibration step. A completeset of ray sums at a particular gantry angle comprises a projection or profile.Measurements of the x-ray beam intensity after penetrating the patient, I, are recorded;since the output of the x-ray tube, I0, and the size of the voxels, t, are known, the sum ofthe attenuation coefficients (! !i,j) along particular paths can be calculated. These valuesare calculated for many different directions. The task of x-ray computed tomographyimage reconstruction is to solve for the individual attenuation coefficients of each voxel,and to assign a value depending on the attenuation coefficients to each pixel in a two-dimensional array which then describes the transverse image.

Several correction factors need to be implemented prior to the calculation. Thex-ray beam is not mono-energetic and attenuation is known to depend on x-ray energy.The effective energy of the x-ray beam increases as it passes through the patient dueto the greater attenuation of lower-energy x-rays, an effect known as beam hardening. Thisresults in the effective attenuation coefficients of a voxel containing a particular materialdecreasing with its distance into the patient; correction algorithms typically estimate thedistance traveled by the x-rays to each voxel. A further complication is that x-ray beams atan angle to the voxel grid traverse different path lengths through different voxels.

The individual attenuation coefficients for each voxel can be obtained exactly, byusing simultaneous equations, essentially by matrix inversion. However, for a matrix sizeof 512! 512, at least 262144 simultaneous equations would be needed to ensure thatthey were independent, and more would be required to account for noise and patientmotion. Matrix inversion is being used in industrial applications where limited rotationangles are possible. In medical imaging a number of alternative algorithms, such as

µ1 µ2

µ4µ3

P 4

P2

P1

3

21

6

0

1

4

4

5 2P3

Figure 3.28 An object comprising four voxels and its projections.

70 Medical images obtained with ionizing radiation

I1 = Ii.e!(µ1+µ3 ).x

I2 = Ii.e!(µ2+µ4 ).x

I3 = Ii.e!(µ1+µ4 ).x

I6 = Ii.e!(µ2+µ3 ).x

Page 13: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Reconstrução – método iterativo

the attenuation coefficients (!i,j) for each of the voxels in that particular path. Substitutingthe sum of the attenuation coefficients along a line of voxels (!i,j) into Equation (3.3) andtaking the natural logarithm gives the ray sum, p, along this particular path:

p ! "ln#I0=I$ ! !! i;jt (3:7)

The reference intensity, I0, is measured for each detector in a calibration step. A completeset of ray sums at a particular gantry angle comprises a projection or profile.Measurements of the x-ray beam intensity after penetrating the patient, I, are recorded;since the output of the x-ray tube, I0, and the size of the voxels, t, are known, the sum ofthe attenuation coefficients (! !i,j) along particular paths can be calculated. These valuesare calculated for many different directions. The task of x-ray computed tomographyimage reconstruction is to solve for the individual attenuation coefficients of each voxel,and to assign a value depending on the attenuation coefficients to each pixel in a two-dimensional array which then describes the transverse image.

Several correction factors need to be implemented prior to the calculation. Thex-ray beam is not mono-energetic and attenuation is known to depend on x-ray energy.The effective energy of the x-ray beam increases as it passes through the patient dueto the greater attenuation of lower-energy x-rays, an effect known as beam hardening. Thisresults in the effective attenuation coefficients of a voxel containing a particular materialdecreasing with its distance into the patient; correction algorithms typically estimate thedistance traveled by the x-rays to each voxel. A further complication is that x-ray beams atan angle to the voxel grid traverse different path lengths through different voxels.

The individual attenuation coefficients for each voxel can be obtained exactly, byusing simultaneous equations, essentially by matrix inversion. However, for a matrix sizeof 512! 512, at least 262144 simultaneous equations would be needed to ensure thatthey were independent, and more would be required to account for noise and patientmotion. Matrix inversion is being used in industrial applications where limited rotationangles are possible. In medical imaging a number of alternative algorithms, such as

µ1 µ2

µ4µ3

P 4

P2

P1

3

21

6

0

1

4

4

5 2P3

Figure 3.28 An object comprising four voxels and its projections.

70 Medical images obtained with ionizing radiation

Page 14: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Reconstrução – retroprojeção filtrada

Transformada de Fourier unidimensional da projeção

Transformada de Fourier bidimensional da imagem

Page 15: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Tomografia computadorizada

new generation. In the third-generation scanners, an x-ray tube mounted on the gantryrevolves around the patient, a tightly collimated fan-beam of x-rays enters the patient,and an arc of detectors on the opposite side, and rotating synchronously with the x-raytube, records the intensity of emerging radiation (Fig. 3.25(ii)). The x-ray tube is alwayson, and readings from the detectors are taken about a thousand times during the rotation.

The first-generation CT scanner, used for scans of the head, illustrates many of theprinciples of the method, and is easier to visualize (Fig. 3.26). A tightly collimated beamminimizes scatter and radiation dose to adjacent tissue. The beam is swept linearly acrossthe patient’s head, and a single detector, moving synchronously with it, measures thetransmitted radiation at regular intervals; this is known as a scan. The tube is turned off,

Figure 3.24 The body as a loaf of sliced bread; a transverse image is a slice of bread viewed end on.

(i) (ii)

Figure 3.25 (i) Outside view of third-generation CT scanner showing the patient table and gantry aperture.(ii) When the cover is removed, the x-ray tube and the arc of detectors can be seen. The gantryholding the x-ray tube and detector rotate around the patient as the data are gathered.

68 Medical images obtained with ionizing radiation

Page 16: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Tomografia computadorizada

Page 17: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Tomografia computadorizada – Artefatos

as !N. Increasing the number of emitted quanta or the scan time would increase the doseto the patient; more useful would be to increase the quantum efficiency of the detectors. Ifthe voxel size is increased, then signal-to-noise ratio would improve but at the expense ofreduced spatial resolution, either in-plane or axially. Using a reconstruction algorithmthat incorporates smoothing would also improve the signal-to-noise ratio at the expenseof spatial resolution.

Because CT images are digital their contrast can be easily manipulated by the displaylook-up table used after image reconstruction. This makes noise the major limitation indetecting low-contrast details. Typically computed tomography can distinguish differ-ences of " 0.5%, i.e. " 5HU, which is better than film-screen radiography which needsabout a 10% difference.

If the patient moves during the acquisition time, characteristic streaks appear in theimage. Ring artifacts occur if the detectors are poorly calibrated or if several fail(Fig. 3.31(i)). Beam-hardening artifacts can occur if the necessary pre-processing correc-tions were not accurate. Partial volume effects result when voxels are partially filled byhighly attenuating materials; the weighted average attenuation can be close to the fullvalue for that material, resulting in narrow, highly attenuating objects, such as the bonyskull, appearing wider than they should be. Streak artifacts (Fig. 3.31(ii)) are commonaround metal objects, such as dental fillings, and are due to a combination of effectsincluding partial voluming, scatter and incomplete beam-hardening correction.

In conventional CT scanners only a single slice is acquired at a time, with the patienttable stationary. The gantry would spin 360° in one direction and enough readings would

(i) (ii)

Figure 3.31 CT artifacts: (i) ring artifacts; (ii) streak artifacts.

3.2 Images from x-rays 75

Page 18: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Equipamentos de Tomografia

•  1a. Geração = 5 minutos por fatia

•  2a. Geração = 10 segundos por fatia

•  3a. Geração = 2 segundos por fatia

•  4a. Geração = 1 segundo por fatia

•  5a. Geração = 50 milisegundos por fatia

Page 19: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

MEDICINA NUCLEAR

Page 20: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Cintilografia • Mapeia raios-γ

emitidos por radioisótopos (radiofármaco)

• 

radioactivity. The collimator only allows !-photons traveling in a certain range ofdirections to interact with the scintillation crystal, so that the site of origin ofthe radioactive event in the patient can be found. However, (Compton) scattering withinthe patient results in some scattered !-photons entering the field of view and degrading the

!-rays

!-rays

Collimator

Nal(TI)scintillation

crystal

Photomultipliertubes

Analog-to-digital converter

Computer

V1 V2 V3

Figure 3.34 Schematic diagram for obtaining a planar nuclear medicine image, using a gamma camera. (AfterWebb, 2003, p. 58.)

1

5

10

16

23

29

34 35 36 37

30 31 32 33

24 25 26 27 28

17 18 19 20 21 22

11 12 13 14 15

6 7 8 9

2 3 4

coordinate

coordinate

coordinatey –

x – x +

y +

coordinate

Figure 3.35 Position encoding within the gamma camera.

3.3 Images from !-rays 79

Page 21: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Cintilografia • Baixa resolução spacial • Baixo contraste •  Imagem ruidosa

•  Fisiologia • Metabolismo

3.3.2 SPECT imaging

In single-photon emission computed tomography, SPECT, a rotating gamma camera,with one, two or three detector heads, rotates around and as closely as possible to thepatient because spatial resolution decreases with the distance from the collimator.Sensitivity increases and acquisition time decreases with more detector heads.

The different acquired projections are used to reconstruct cross-sectional or three-dimensional images by filtered backprojection. Reconstruction computations are morecomplicated than with x-ray computed tomography because the detected signals dependupon both the spatial distribution of the radioisotopes and the attenuation properties ofthe voxels. As with x-ray computed tomography imaging, the main advantage of SPECT

Figure 3.37 Whole-body bone scan, using 99mTc-MDP. The anterior view is on the left, and the posteriorview is on the right. Increased uptake of the tracer around the right knee indicates the presenceof a tumor.

3.3 Images from !-rays 81

Page 22: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

SPECT • Single photon emission computed tomography • Câmara gama giratória • Reconstrução por

retroprojeção filtrada

Page 23: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

SPECT

Page 24: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

ULTRASSONOGRAFIA

Page 25: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Ultrassonografia (US) • Mede a refração de ondas de

ultrassom quando atingem um determinado tecido

which is smooth compared to the ultrasound wavelength and perpendicular to thedirection of wave propagation, is given by

R ! Z1 " Z2# $2

Z1 % Z2# $2(4:4)

where Z1, Z2 are the acoustic impedances of the materials to either side of the surface.Acoustic impedance is a constant for a specific material (Table 4.1) and is analogous toelectrical impedance: as impedance increases it inhibits velocity (current) for a givenpressure (voltage). Acoustic impedance is given by

Z ! rv !!!!!r!

r(4:5)

The SI unit of acoustic impedance is the rayl (kg m!2 s!1).There is little reflection if the materials are acoustically similar and a lot of reflection

when there is a mis-match of acoustic impedances. For example, at an interface betweensoft tissue and bone a very large reflected signal or echo results, comprising about 40% ofthe incident intensity. This greatly attenuates the transmitted beam and makes theimaging of structures deeper-lying than bone extremely difficult. At a soft tissue/gasinterface, around 99% of the beam intensity is reflected, making it impossible to scandistal structures deeper than the lungs or gas-containing bowel.

Even when the beam is traveling through biological tissue, it loses intensity continu-ously as a result of scattering and absorption (Fig. 4.2). Although the mechanisms arecomplicated, the overall effect is that the beam energy decreases more or less exponen-tially as it penetrates the tissue, similar to the attenuation of x-ray intensity. Thus, theultrasound beam intensity, I, decreases with propagation distance, x, from its startingvalue at x = 0, I0, according to

I#x$ ! I0e#"mx$ (4:6)

Tissue 1 Tissue 2

Attenuation

Attenuation

Some reflection,some transmission

Attenuation

Figure 4.2 The returning echo pulse suffers continuous attenuation along its path, and an abrupt change inintensity on reflection at the interface. (After Wolbarst, 1993, p. 408.)

4.1 Ultrasound imaging 93

Sound waves are longitudinal waves, i.e. the particles of the material move back andforth in the same direction that the wave is traveling. The speed of sound in a material, v,is characteristic of that material and depends on the density of the material, !, and itscompressibility, !. The easier it is to compress a material, the higher is its compressi-bility. Thus:

v ! 1!!!!!!!r

p (4:3)

We can compare bone with soft tissue. Although bone has a larger density, it has a muchsmaller compressibility than soft tissue. The product !! is smaller for bone than softtissue, resulting in a larger velocity for sound waves through bone.

A pulse of ultrasound, which is what is often used in medical ultrasound rather than acontinuous wave, actually comprises a range of frequencies: the shorter the pulse thelarger the range of frequencies comprising it. Luckily, the velocity of sound in a mediumis nearly independent of frequency or wavelength, otherwise the pulse would spread outas it traveled leading to pulse blurring. This behavior is different from that of light: thespeed of light in a medium depends on wavelength, which is why prisms split sunlightinto its constituent colors.

When an ultrasound wave encounters a tissue surface, separating tissues with differentacoustical properties, a fraction of the wave is backscattered and detected by thetransducer on its return. Generally, only those waves that reflect back through about180° can contribute to an ultrasound image. By measuring the delay between pulsetransmission and pulse reception, and knowing the speed of propagation, the depth of thefeature can be calculated. For example, if the time delay is 160 "s and the pulse is passingthrough soft tissue with a speed of 1540m s!1, the round-trip path is 24.6 cm and thetissue depth is 12.3 cm.

The intensity of the echo is used to determine the brightness of the image at thereflecting tissue surface (Fig. 4.2). The intensity reflection coefficient, R, at a boundary

Skin

Transducer

pulse duration

1–10 MHzsound frequency

Sound Velocity1540 m s–1

(in soft tissue)

Pulse repetition time: 1 ms

Pulse repetition rate:1000 Hz

1 µs

Figure 4.1 Schematic diagram of a typical clinical ultrasound beam. (After Wolbarst, 1993, p. 408.)

92 Medical images obtained with non-ionizing radiation

Page 26: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Ultrassonografia (US) • Não-ionizante • Não-invasiva •  Tempo real • Procedimentos

terapêuticos

Page 27: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Doppler • Mede o fluxo

sanguíneo através de ondas sonoras

Page 28: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

RESSONÂNCIA MAGNÉTICA NUCLEAR

Page 29: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Ressonância Magnética Nuclear

of pixels. For the last dimension the signal is encoded in terms of phase. This is not easyto understand: suffice it to say that a number of gradients are needed to create phasechanges from row to row prior to acquisition so that the FT is provided with enoughinformation to encode fully the final image. This is known as phase-encoding.

Three separate gradient coils are required to encode the three spatial dimensionsunambiguously. Multiple projections can be obtained by electronically rotating themagnetic field gradients around the patient. The Fourier transform of the echo signalsreceived produce what are effectively projections of the object along certain directions.These one-dimensional projections can be assembled to give the two-dimensionalFourier transform of the object, which can then be inverse Fourier transformed to givean axial image; this is the direct Fourier reconstruction (DFR) technique (Section 7.8.2),which is almost universally used in MRI reconstruction.

The basic components of a MRI scanner are shown in Figure 4.16. The primarymagnet polarizes the protons in the patient. It can either be an electromagnet usingsuperconducting coils, which require cooling to very low temperatures, and whichresembles the gantry of a CT scanner; or it can be a permanent magnet constructed ofrare-earth alloys, in which case it has a more open structure. The gradient coils producelinear variations on the magnetic field, so that the proton resonant frequencies within thepatient are spatially dependent; they are fixed permanently within the primary magnet.The RF coil produces the oscillating magnetic field necessary for creating phase coher-ence between protons, and receives the FID signal by magnetic induction; it is placedaround the body part of the patient to be imaged, with a geometry optimized for thatspecific part.

B1

RF coil

RFreceiver

RFtransmitter

Computer systemand display

Pulseprogrammer

x y z

Gradient coils

Gradient drivers

Superconductingmagnet

Bx, By, Bz

B0

Figure 4.16 Basic components of a MRI scanner.

106 Medical images obtained with non-ionizing radiation

Figure 4.19 is a spin echo sequence diagram. The bottom line illustrates the evolutionof the MR signal (the FID immediately after the 90° pulse and the echo at time TE). Notethat the repetition time is also labeled. Gradients are illustrated by rectangular blocks, thearea of which represents the amplitude and the sign (i.e. positive or negative) dictated bythe position above or below the “time” axis. In this example the phase encoding is in the ydirection and the phase encoding gradient (Gy) is drawn as multiple lines to illustrate thatthe amplitude of this changes each time the sequence is repeated. In contrast, frequencyencoding (Gx) is performed all at once at the time of signal detection. A de-phasing lobe(negative half of area) compensates for changes in phase, such that at the time of the echoonly a frequency change is exhibited. Lastly, the slice-selection gradient (Gz) has to beapplied at the time of both RF pulses so that only the spins within the slice of interest areexcited and refocused. It also uses a de-phasing lobe.

The acquisition time for the spin echo sequence is given by the product of the TR of thepulse sequence and the number of phase encoding steps (the number of pixels in thephase direction). If multiple acquisitions were done to improve the signal-to-noise ratio(SNR) then the total acquisition time would be multiplied by this factor. By recording theecho more than once the coherent signal is additive but the incoherent noise cancels out.

RF pulses

T2 decay

*T2 decay

90° 180° 180° 180° 180°

Figure 4.18 A spin echo pulse sequence, comprising a 90° pulse followed by carefully timed 180° pulsesproduces a series of spin echoes of decreasing amplitudes. Although each individual echo decays asT !2 , the envelope decays as T2.

RF 90° 180°

Gz

Gy

Gx

Signal

TE

TR

Figure 4.19 Pulse sequence for spin echo imaging.

108 Medical images obtained with non-ionizing radiation

Page 30: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

MRI

Page 31: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

MRI – T1, T2 e FLAIR

Page 32: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

MRI – diferentes contrastes

Page 33: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

MRI – Cirurgia guiada por imagem • Protótipo de Robô guiado por imagens de Ressonância

Magnética (Hitachi)

Page 34: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Imagem de tensores de difusão (DTI)

⎟⎟⎟

⎜⎜⎜

=

zzzyzx

yzyyyx

xzxyxx

DDDDDDDDD

D

Page 35: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

DTI

Page 36: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

DTI - Alzheimer

Page 37: PROCESSAMENTO DE IMAGENS MÉDICAS - Unicampadessowiki.fee.unicamp.br/media/Attachments/courseIA369O1S2011/Cap... · PROCESSAMENTO DE IMAGENS MÉDICAS IA369O – Capítulo 2 Leticia

Resolução espacial

the object. Typical values for the spatial resolution of various medical imaging systemsare given in Table 3.1.

ContrastX-rays interact with a patient’s body in various ways (Fig. 3.12). They can penetrate thebody, moving in a straight line from source to detector; these are the primary x-rays thatproduce a faithful shadowgram. Some x-rays are absorbed within the body by thephotoelectric effect; these x-rays do not contribute to producing an image, and insteadincrease the dose, and therefore the risk, to the patient. Others undergo (Compton)scattering within the patient, producing scattered or secondary x-rays, some of whichreach the detector; they produce random blackening over the film. If the subjectcontrast is defined as log10 (I1/I2), where I1, I2 are the transmitted x-ray intensities

Table 3.1 Typical spatial resolution of medical imaging modalities.

Modality Spatial resolution (mm)

Film-screen radiography !0.05Computed radiography 0.1–0.2Computed tomography (CT) 0.25–0.5Fluoroscopy !0.5Magnetic resonance imaging (MRI) 0.5–1.0Ultrasound 1.0–5.0Nuclear medicine (NM) 3.0–10.0

AbsorptionPatient

Secondary x-rays

Detector

Scattering

Primary x-rays

Figure 3.12 Interaction of x-rays with a patient’s body.

3.2 Images from x-rays 57