Digital Image Processing in Radiography
Digital Image Processing in Radiography
Historical Development
Discovery of computed tomography (CT) in the 1970s
Introduction of computed tomography scanning created the concept of electronic data
sets for medical images
The explosion of the imaging sciences that took place between the early 1970s and the
1990s was caused by the first wave of digital computerization in CT and diagnostic
medical sonography (DMS)
o Computed tomography captured the imagination of the public
Digital processing has now been successfully applied to magnetic resonance imaging
(MRI), nuclear medicine, cardiovascular imaging, and mammography, and
predominates in diagnostic imaging
Digital radiography imaging systems replace traditional film with a reusable detector
2
One standard that has been central to the development of modern radiological imaging
is the Digital Imaging and Communication in Medicine (DICOM)
The radiologic community has promulgated this digital standard to ensure that all
equipment from all manufacturers who choose to adhere to the standard is
speaking the same computer language that can be understood across
manufacturer and device lines
Each sampled piece of analog data is then computed and assigned a discrete
analog value, a process called quantification
These analog values are then converted to a binary digit, using binary counting
methods
3
Just as English has a 26-letter-symbol alphabet, the binary system operates with a two-
symbol alphabet
Because electrical currents are most easily understood as being either on or off, the
binary system consists of information recorded as either a 0 for off or a 1 for on
Each binary number is called a bit, for binary digit
Nominal numbers function on a base 10 system instead of base 2
An 8-bit word is required to form the 26 letters of the English alphabet
an 8-bit word is called a byte
Computer memory is often rated in terms of the total byte memory (rounded off)
- Example: a 10-megabyte magnetic hard disk will store over 10,000,000
bytes (80,000,000 bits) of information
4
That three-dimensional data point is referred to as a volume element or a voxel
The overall dimension of the image matrix is called the field of view (FOV)
In digital radiography imaging systems, the FOV is determined by the size of
the detector
In other modalities, such as CT and MRI, the operator can select the FOV for a
particular study
If the FOV remains the same, then as matrix size increases, the pixels must get
smaller, creating a sharper image
If the overall dimensions of the FOV are 400 × 400 mm and the matrix size is
200 × 200, then every pixel would be 2 mm in size
If the overall dimensions of the FOV remain
unchanged at 400 × 400 mm and the matrix
size is increased to 400 × 400, then every pixel
would be 1 mm in size
The greater the matrix size for the same FOV,
the better the spatial resolution because a
larger matrix provides smaller pixels
Spatial resolution or image sharpness is a critical image
quality factor that is determined by pixel size
Spatial resolution is measured in line pairs per
millimeter (lp/mm)
In digital imaging it can be no better than the size
of the pixel
Each pixel in the matrix is capable of representing a wide
range of different shades of gray from pure white through
total black
The number of shades of gray is determined by the
grayscale bit depth
Grayscale bit depth ranges from 8 bits to 32 bits
A grayscale bit depth of 8–32 equals a range of 1–4 bytes of storage required per
pixel in the image matrix
A grayscale bit depth of 12 produces 2¹² gray levels
This represents 4,096 different shades of gray that are available
This is also known as density resolution
5
Spatial Domains
Digital images are represented in two domains, based on how the image is acquired
Images in the spatial location domain are based on a matrix that has specific locations
for each pixel
The spatial location is usually described as the pixel location
the x coordinate representing the location of the pixel on a horizontal line from
left to right
followed by the y coordinate representing the location of the pixel on a vertical
line from top to bottom
Images in the spatial frequency domain are based on the number of cycles per unit
length
which is the definition of frequency
The spatial frequency is usually described as the resolution (object size) and the
contrast
Small structures produce high frequency and high contrast
Low frequency is produced by larger structures and low contrast
The Fourier transformation (FT) is the mathematical algorithm that is applied to
change an image from the spatial location domain to the spatial frequency domain
and an inverse FT can be used to return an image back to the spatial location
domain
Physicians and radiographers view images in the spatial location domain
o whereas physicists and engineers can extract information from the
mathematical representations of the image in the spatial frequency domain
There are a number of computer operations that can occur to process the image
These image processing operations are primarily intended to change the input values
of the pixels to improve the diagnostic quality of the output image
They include point processing, local processing, and geometric processing
6
Of these, point processing operations are the most important to understand for digital
radiography imaging systems
Point processing operations are those that are performed between the receipt of the input
image from the image receptor and the output image that is viewed on the monitor
Point processing operations involve
adjusting the value of an input pixel
(point) to the corresponding output
pixel
They are the most common
processing operations in digital
radiography imaging systems
One common point processing
operation is called grayscale
processing, which allows for
adjustments to image brightness and
contrast
This process involves the creation of a
histogram and the application of the
look-up table (LUT) and windowing
7
If the exposure is outside the range from underexposure or overexposure, then the
computer will correct the image by shifting or rescaling the histogram to the correct area
This is referred to as histogram modification or stretching
A wide histogram demonstrates higher contrast
A narrow histogram shows lower contrast
If the histogram values are concentrated in lower comparative values, the image
appears dark or dim
If the histogram values are
concentrated in higher
comparative values, the
image appears bright or
light
Windowing is a point processing operation that changes the contrast and brightness of the
image on the monitor
The brightness
and contrast of the
digital image
depend on the
shades of gray,
which are
controlled by
varying the
numerical values
of each pixel
The human visual
range
encompasses 32 or
fewer shades of
8
gray, whereas the photon beam exiting the patient encompasses over 1,000
variations
A 12-bit depth has 4,096 shades of gray
Because the stored shades exceed the visual range, any digital image is only a
window of the total data
The window width (WW) is the range of shades of gray that will be displayed
A narrow width → few shades, high contrast
A large width → many shades, low contrast
The window level is the center of the window width and controls brightness
A low level → image will be light
A high level → image will be darker
The computer can easily change the level and width by mathematical recalculations
The quantity of information stored for each pixel varies depending on the grayscale (bit
depth) imaging abilities of the system
The radiographer controls which portion of this information will be seen by
manipulating the data prior to sending it to PACS
Radiologists can manipulate the data sent to them
However, the more the technologist changes the original data before it is sent
to PACS, the less information the radiologist has to work with
It is critical that radiographers be familiar with how the computer determines what
portion of the image information will be displayed
Careful image processing is critical to avoid obscuring diagnostically critical
information
Grayscale processing, both image subtraction and temporal averaging are considered
point processing operations
o These functions are used primarily in digital subtraction angiography
Local processing operations are those in which mathematical calculations are applied
to only a small group of pixels, although the operation will often be continued until all
pixels have been subjected to it
o This process is often referred to as a kernel
o A kernel is processing code that is mandatory and common to the computer
system
It is normally applied over and over to the entire set of data being
processed
o Local processing is also called area or group processing
9
Spatial frequency filtering is considered a type of local processing operation
o It is used to sharpen, smooth, blur, reduce noise, or pull elements of interest
from an image
o It can occur in the spatial location domain or in the spatial frequency domain
o When spatial frequency filtering is done in the spatial frequency domain, the
Fourier transform is used
o When it is done in the spatial location domain, the pixel values themselves are
used
10
subtracts a low-pass filtered (blurred) image from the original image, thus
producing a new subtracted and sharper image
Spatial location filtering is another type of local processing operation, and is also
known as convolution
A kernel is applied repeatedly to each pixel in a matrix in order to weigh the
values or apply a coefficient across the matrix
Application of convolution kernels usually works better when array processors
can be used to rapidly apply them
Geometric processing operations are used to change the position or orientation of the
pixels in the image
This allows rotation, magnification, and other operations similar to those used in
digital photographs
Spatial Resolution
Spatial resolution is controlled by the matrix size and how many pixels can be
displayed by the monitor
There is a direct relationship between matrix and spatial resolution
When the matrix size is increased and pixel size goes down, spatial resolution
increases
Matrix size is not typically variable within a system
11
It is controlled by the design of image detectors and the accompanying
electronics built into the system
12
Noise
13
Detective Quantum Efficiency (DQE)
A measure of the sensitivity and accuracy by which the image receptor converts the
incoming data to the output viewing device is known as the detective quantum
efficiency (DQE)
A perfect device would perform this task with 100 percent efficiency (no loss of
information) or with a DQE of 1
Digital radiography imaging systems have DQEs of 0.3–0.7 or 30–70 percent
DQE is affected by the:
SNR
quantum noise
system noise
DQE can be described as the ratio of the squared
output signal-to-noise ratio (SNRo)² to the squared
input signal-to-noise ratio (SNRI)² of the imaging
detector
It can be calculated using the formula
DQE = (SNRo)² / (SNRI)²
In addition to considering SNR, DQE is dependent on
modulation transfer function (MTF)
incident x-ray energy
detector material
spatial frequency
Modulation transfer function (MTF) is used to measure the capacity of the detector to
pass its spatial resolution characteristics to the final image
A perfect detector would have an MTF of 1
DQE and MTF of a digital detector are directly proportional
DQE and SNR have an inverse relationship
High DQE values indicate that less radiation is necessary to achieve identical image
quality that a lower DQE detector would produce
o a digital receptor with a higher DQE does not always achieve improved image
quality at a lower radiation dose
o Many other factors affect image quality:
14
incident beam energy
signal sampling and processing
quality of image display
Receptor Exposure
Exposure Indicators
Exposure indicators reflect the amount of x-ray exposure to an image receptor, which
looks at the quantity of photons that strike the detector
Exposure indicator numbers can be used as a
quality control tool for image quality and
radiation protection
Departments have established target exposure
indicators (EIt), which reflect optimum
exposure and acceptable noise levels
The correct exposure factors will deliver an
exposure indicator within an acceptable range and
result in the best image quality
It is important for the radiographer to understand
the exposure indicator numbers for each image
and how to use this data to assess image quality
The exposure sensitivity of imaging detectors ranges from 0.1 mR to 100 mR, a range of
approximately 1,000:1
If soft tissue and bone densities are imaged, one exposure will give the necessary
information due to the ability to change window levels
To evaluate exposure, most manufacturers have an exposure indicator that provides
information on the average amount of radiation used to produce an ideal diagnostic
image
15
Different manufacturers use different systems and acceptable noise levels
16
The most important concept is the deviation index (DI)
Defined as an indicator of whether the detector response agrees with KTGT
KIND = indicated equivalent air Kerma (radiation dose at detector surface)
KTGT = target equivalent air Kerma (radiation dose for proper exposure)
DI compares actual exposure vs proper exposure
Summary
Digital radiography imaging systems replace traditional film with a reusable detector
These systems are divided into computed radiography (CR) and digital
radiography (DR)
A digital image is one that has been converted into numerical values for transmission
or processing
A matrix is a series of boxes laid out in rows and columns that gives form to the
image
The individual matrix boxes are known as picture elements or pixels
Each pixel can represent a wide range of different shades of gray
The number of shades of gray is determined by the grayscale bit depth
17
There are a number of computer operations that process the image
These operations change the input values of the pixels to improve diagnostic
quality
They include point processing, local processing, and geometric processing
One common processing operation is grayscale processing
Allows adjustments of image brightness and contrast
Involves creation of a histogram, use of a look-up table (LUT), and windowing
A histogram is generated during initial processing
It allows the digital system to find the useful signal by locating the minimum and
maximum signal within the anatomical region of interest
An LUT is applied to the data that has the standard contrast for that exam
Gives the desired image contrast for display
Windowing changes the contrast and brightness of the image on the monitor
Digital image quality factors include
Spatial resolution
Noise
DQE
Artifacts
When matrix size is increased and pixel size goes down, spatial resolution increases
Low-contrast resolution is a type of contrast resolution
Deals with the ability to visualize subtle energy differences in soft tissues
Noise can be classified as
Electronic system noise
Ambient noise
Quantum mottle noise
Detective quantum efficiency (DQE) is a measure of the sensitivity and accuracy by
which the image receptor converts incoming data to the output viewing device
Exposure indicators reflect the amount of x-ray exposure to an image receptor
Look at the quantity of photons that strike the detector
Exposure indicator numbers can be used as a quality control tool for image
quality and radiation protection
18