Insight into Images Principles and Practice for Segmentation Registration and Image Analysis 1st Edition Terry S. Yoo - Discover the ebook with all chapters in just a few seconds
Insight into Images Principles and Practice for Segmentation Registration and Image Analysis 1st Edition Terry S. Yoo - Discover the ebook with all chapters in just a few seconds
com
https://ebookgate.com/product/insight-into-images-
principles-and-practice-for-segmentation-registration-and-
image-analysis-1st-edition-terry-s-yoo/
OR CLICK BUTTON
DOWLOAD EBOOK
https://ebookgate.com/product/image-registration-for-remote-
sensing-1st-edition-jacqueline-le-moigne/
ebookgate.com
https://ebookgate.com/product/microarray-image-and-data-analysis-
theory-and-practice-1st-edition-luis-rueda/
ebookgate.com
https://ebookgate.com/product/coaching-science-theory-into-practice-
wiley-sport-texts-1st-edition-terry-mcmorris/
ebookgate.com
https://ebookgate.com/product/apa-handbook-of-behavior-analysis-
vol-2-translating-principles-into-practice-gregory-j-madden-editor/
ebookgate.com
Insight into IELTS Updated Edition Vanessa Jakeman
https://ebookgate.com/product/insight-into-ielts-updated-edition-
vanessa-jakeman/
ebookgate.com
https://ebookgate.com/product/deep-learning-for-medical-image-
analysis-1st-edition-s-kevin-zhou/
ebookgate.com
https://ebookgate.com/product/middleton-s-allergy-principles-and-
practice/
ebookgate.com
https://ebookgate.com/product/an-insight-into-chemical-
engineering-1st-edition-m-subbu/
ebookgate.com
https://ebookgate.com/product/image-processing-principles-and-
applications-1st-edition-tinku-acharya/
ebookgate.com
Insight Into Images
Principles and Practice for
Segmentation, Registration,
and Image Analysis
A K Peters
Wellesley, Massachusetts
Editorial, Sales, and Customer Service Office
AK Peters, Ltd.
Weliesey, MA 02482
888 Worcester Street, Suite 230
www.akpeters.com
Ali rights reserved. No part of the material protected by this copyright notice may
be reproduced or utilized in any form, electronic or mechanical, including photo
copying, recording, or by any information storage and retrieval system, without
written permission from the copyright owner.
p. cm.
I ncludes bibliographical references and index.
ISBN 1 -56881 -217-5
1. Diagnostic imaging-Digital techniques. 2. I mage processing.
I. Yoo, Terry S., 1963-
Printed in Canada
08 07 06 05 04 10 9 8 7 6 5 4 3 2 1
To my parents, T.S.Y.
To my wife Geanne, Baby Ada and my parents, D.N.M.
To Konrad, Marcin and Tomasz C.1.
To Stephen P izer, G.D.S.
To the ITK team at the lab, J.C.G.
To my family and parents, JK
To Geoff, L.N.
For Dan and Marlene, R.T.W
To my beloved wife and my parents, Yinpeng Jin
To Daniel and Camilo, L.1.
To the friendships formed while writing ITK, S.R.A.
To my grandmother, PKS.
To Rebecca and her big heart, B.B.A.
To lab friends old and new, T.A.S.
To my daughter, Y.Z.
To my parents and friends in the ITK group, T.C.
The purpose of computing is Insight, not numbers.
-Richard Hamming
Contents
Foreword xi
1 I ntroduction 3
1.1 Medical Image Processing . . . . . . . . . . . 4
1 .2 A Brief Retrospective on 3D Medical Imaging 5
1 .3 Medical Imaging Technology . . . . . . 6
1 .4 Acquisition, Analysis, and Visualization . 16
1 .5 Summary . . . . . . . . . . . . . . . . . 17
vii
viii Contents
3.7 Conclusion . . . . . . . . . . . . . . . . . . 92
3.8 Appendix: Extruded Gaussian Distributions . 93
II Segmentation 119
5 Segmentation Basics 1 21
5. 1 Introduction . . . . · 121
5.2 Statistical Pattern Recognition . · 1 23
5.3 Region Growing . . . . . . . . · 1 24
5.4 Active SurfaceslFront Evolution . 1 26
5.5 Combining Segmentation Techniques 1 27
5.6 Looking Ahead . . . . . . . . . . . . 1 28
389
I ndex
Foreword
xi
xii Foreword
• GE Global Research
• Harvard BWH Surgical Planning Lab
• Kitware, Inc.
• Insightful, Inc.
ITK has met its primary goals for the creation of an archival vehicle for im
age processing algorithms and an established functioning platform to accelerate
new research efforts. New projects have connected ITK to existing software ap
plications for medical and scientific supercomputing visualization such as "An
alyze" from the Mayo Clinic and "SCIRun" from the University of Utah's SCI
Institute. These projects were sponsored to demonstrate that the application pro
grammer's interface (API) for ITK was suitable for supporting large applications
in established medical and scientific research areas. Also, recent projects have
been funded to add algorithms not previously funded to the collection, proving
the versatility of the underlying software architecture to accommodate emerging
ideas in medical imaging research.
There is growing evidence that ITK is beginning to influence the international
research community. In addition to the sponsored growth of the ITK community,
the mailing list of users and developers includes over 300 members in over 30
countries. These facts may be indicators that a nascent research community may
be forming about the core ITK software. The applications for ITK have begun
to diversify including work in liver RF ablation surgical planning, vascular seg
mentation, live-donor hepatic transplant analysis, intra-operative 3D registration
for image-guided intervention, rapid analysis of neurological function in veter
ans, longitudinal studies of Alzheimer's patients using MR!, and even handwrit
ing analysis. The Insight Program of the Visible Human Project is emerging as a
clear success and a significant and permanent contribution to the field of medical
imaging.
In 1 999, the NLM Office of High Performance Computing and Communi
cations, supported by an alliance of NIH ICs and Federal Funding Agencies,
awarded six contracts for the formation of a software development consortium to
create and develop an application programmer interface and first implementation
of a segmentation and registration toolkit, subsequently named the Insight Toolkit
(ITK). The final deliverable product of this group has been a functional collection
of software components, compatible for direct insertion into the public domain
via Internet access through the NLM or its licensed distributors. Ultimately, NLM
hopes to sponsor the creation of a public segmentation and registration software
toolkit as a foundation for future medical image understanding research. The in
tent is to amplify the investment being made through the Visible Human Project
and future programs for medical image analysis by reducing the reinvention of
basic algorithms. We are also hoping to empower young researchers and small
research groups with the kernel of an image analysis system in the public domain.
PaIt to attempt to describe it all. The concept of Finite Element Models (FEM)
is a book-length treatise in itself, as are each of the topics of Level Set Methods
and Statistical Pattern Recognition. Indeed, we reference many valuable texts on
these topics in particular.
We therefore limit this material to methods that have been implemented in the
Insight ToolKit. All methods described herein are reflected in the source code
for ITK. We limit the perspectives in this text to the viewpoints presented in ITK,
and though there are many valid and valuable alternate approaches incorporating
many of these same ideas, we cannot treat them all, and in fairness, we claim
expertise only in those areas that we can deomnstrably express directly through
software development.
In keeping with the principles of the Insight project itself, we concentrate
only on medical image filtering, segmentation, and registration. These were the
primary areas of concentration of the software project, and thus they are the topics
described in this book. As such, we do not directly treat the ares of visualization,
identification, classification, or computer-aided detection and diagnosis. All of
the methods presented in this book are relevant to these topics, and we cannot
show our results without some form of visualization, but we do not include an
in-depth treatise of these topics or explicitly describe the integration of ITK with
visualization environments.
instructors and their staff to create the software infrastructure necessary to run
high-dimensional image processing software. For such a course, this book would
represent theoretical course material, and the ITK Software Guide would be an
invaluable resource for laboratory exercises and course projects.
We are extraordinarily proud of ITK and hope that you find our offering in
structive, valuable, and a useful supporting addition to your software projects. We
humbly and gratefully present this material with hope for your future success in
research and medical software application development.
Acknowledgements
Following the consortium-based nature of this project, this book is the collec
tive effort of many, many people. The contributing authors include in alphabet
ical order: Brian Avants, Stephen Aylward, Celina Imielinska, Jisung Kim, Bill
Lorensen, Dimitris Metaxas, Lydia Ng, Punam K Saha, George Stetten, Tessa
Sundaram, Jay Udupa, Ross Whitaker, and Terry Yoo. Other people have made
essential contributions, including: Luis Ibanez, Will Schroeder, and others. The
Insight ToolKit would not have been possible without the selfless efforts of Bill
Hoffman, Brad King, Josh Cates, and Peter Ratiu.
We are indebted to the coalition of federal funding partners who have made
this project possible: the National Library of Medicine, the National Institute
for Dental and Craniofacial Research, the National Institute for Mental Health,
the National Eye Institute, the National Institute on Neurological Disorders and
Stroke, the National Science Foundation, the DoD Telemedicine and Advanced
Technologies Research Center, the National Institute on Deafness and other Com
munication Disorders, and the National Cancer Institute. We would especially
like to recognize the patronage of the Visible Human Project from the Office of
High Performance Computing and Communications of the Lister Hill National
Center for Biomedical Communications of the National Library of Medicine. In
particular, we'd like to thank Drs. Donald Lindberg and Michael Ackerman for
their vision and leadership in sponsoring this project.
Part One
Introduction and Basics
CHAPTER
ONE
Introduction
Terry S. Yoo
National Library of Medicine, NIH
3
4 Chapter 1 . Introduction
The medical mission differs from the other forms o f image processing arising
from non-medical data. In satelite surveillance analysis, the purpose is largely a
screening and cartographic task, aligning multiple types of data and correspond
ing them to a known map and highlighting possible points of interest. In computer
vision, camera views must be analyzed accounting for the perspective geometry
and photogrammetric distortions associated with the optical systems that are the
basis for robotic sensors. In many of these systems, autonomous navigation, tar
get identification and acquisition, and threat avoidance are the primary tasks. For
the most part, the incoming information arrives as 2D images, data arrays that
can be organized using two cartesian dimensions. In addition, the tasks are to
be performed independently by the machine, relying on the development of ma
chine learning and artificial intelligence algorithms to automatically accomplish
the tasks.
In medicine, the problem as well as the input data stream are usually three
dimensional, and the effort to solve the primary tasks is often a partnership of
human and machine. Medicine is notably a human enterprise, and computers are
merely assistants, not surrogates nor possible replacements for the human expert.
The medical task can often be split into three areas; ( 1 ) data operations of filtering,
noise removal, and contrast and feature enhancement, (2) detection of medical
conditions or events, and (3) quantitative analysis of the lesion or detected event.
Of these subtasks, detection of lesions or other pathologies is often a subjective
and qualitative decision, a type of process ill-suited for execution by a computer.
By contrast, the computer is vastly more capable of both quantitative measurent
of the medical condition (such as tumor volume or the length of a bone fracture)
and the preprocessing tasks of filtering, sharpening, and focusing image detail.
The natural partnership of humans and machines in medicine is to provide the
clinician with powerful tools for image analysis and measurement, while relying
on the magnificent capabilities of the human visual system to detect and screen
for the primary findings.
We divide the problems inherent in medical image processing into three basic
categories;
• Filtering: These are the basic tasks involved in filtering and preprocessing
the data before detection and analysis are performed either by the machine
or the human operator.
• Segmentation: This is the task of partitioning an image (2D array or vol
ume) into contiguous regions with cohesive properties.
• Registration: This is the task of aligning multiple data streams or images,
permitting the fusion of different information creating a more powerful di
agnostic tool than any single image alone.
1 .2. A Brief Retrospective on 3D Medical Imaging 5
These three basic divisions represent the organization of both this book and the
Insight Toolkit for which this book was written. In order to best understand the
approaches and difficulties associated with these tasks, we being with a history of
medical imaging and a brief overview of modem imaging modalities.
In the broad realm of the sciences, medical imaging is a young field. The physics
enabling diagnostic imaging are barely one hundred years old. In 1 895, Wilhelm
Roentgen discovered x-rays while experimenting with a Crookes tube, the pre
cursor to the cathode ray tube common in video applications today. It is worth
noting that he immediately recognized the potential of x-ray radiation in diagnos
tic imaging and that one of his earliest images is of the bones of his wife's hand.
Roentgen's discovery was so noteworthy and revolutionary that he received a No
bel prize within six years of announcing his initial work. Thus began a bountiful
and tightly-coupled one hundred year partnership between physics and medicine.
The meaning and purpose of medical imaging has been to provide clinicians
with the ability to see inside the body, to diagnose the human condition. The
primary focus for much of this development has been to improve the quality of the
images for humans to evaluate. Only recently has computer technology become
sufficiently sophisticated to assist in the process of diagnosis. There is a natural
partnership between computer science and radiology.
Early in the twentieth century, a Czech mathematician named Johann Radon
derived a transform for reconstructing cross-sectional information from a series
of planar projections taken from around an object. While this powerful theory had
been known for over fifty years, the ability to compute the transform on real data
was not possible until digital computers began to mature in the 1 970s.
Imaging in 3D emerged in 1972 when x-ray computed tomography (CT) was
developed independently by Godfrey Hounsfield and Alan Cormack. These inno
vators later shared the 1 979 Nobel Prize in Medicine. Their achievement is note
worthy because it is largely based on engineering, the theoretical mathematics and
the underlying science had been described decades earlier. The contribution is in
the application of mechanical engineering and computer science to complete what
had previously only been conceived on paper. Clinical systems were patented in
1 975 and began service immediately thereafter.
While techniques for using x-rays in medical imaging were being refined, or
ganic chemists had been exploring the uses of nuclear magnetic resonance (NMR)
to analyze chemical samples. Felix Bloch and Edward Purcell were studying
NMR in the mid 1940s. Together, they shared the Nobel Prize in Physics in
1 952. Paul Lauterbur, Peter Mansfield, and Raymond Damadian were the first
to develop imaging applications from NMR phenomena. Lauterbur created tomo-
6 Chapter 1 . Introduction
The majority of medical visualization involves the display of data acquired di
rectly from a patient. The radiologist is trained through long years of education
and practice to read relatively simple presentations of the raw data. The key to
improving a diagnosis is in the careful crafting of the acquisition, applying the
physics of radiology to maximize the contrast among the relevant tissues and sup
pressing noise, fog, and scatter that may obscure the objects of interest. This is no
less true for more complex visualizations that incorporate 3D renderings or vol
ume projections; improving the quality of the acquired data will fundamentally
affect the quality of the resulting visualization.
This section will briefly cover some of the more common sources of three or
higher-dimensional medical imaging data. This treatment is necessarily super
ficial, serving mostly as an annotated glossary for many of the terms taken for
granted in the radiology community. To create effective medical visualization
tools, the computer scientist requires a fundamental understanding of the source
of the image data, the technology involved, and the physical principles from which
the image values are derived. This cursory introduction will be insufficient for in
depth research, but will serve as background for this text; the reader is encouraged
to continue the exploration of this topic. This additional command of the basics
of medical image acquisition will enhance your practice in medical visualization,
1.3. Medical Imaging Technology 7
ease communication between you and the doctor and the technologist, improve
the quality of the source data, and smooth transitions among the many interfaces
from acquisition to display.
An accurately calibrated moving bed to translate the patient through the scan
ner, an x-ray tube mounted in such a way to allow it to revolve about the patient,
and an array of x-ray detectors (gas filled detectors or crystal scintillation detec
tors) comprise the essential system components of a CT machine. The x-ray tube
and detector array are mounted in a gantry that positions the detector assembly
directly across from the x-ray source. The x-ray source is collimated by a pair of
lead jaws so that the x-rays form a flat fan beam with a thickness determined by
the operator. During the acquisition of a "slice" of data, the source-detector ring
is rotated around the patient. The raw output from the detector array is backpro
jected to reconstruct a cross-sectional transaxial image of the patient. By reposi
tioning the patient, a series of slices can be aggregated into a 3D representation of
the patient's anatomy. Figure 1 . 1 is a picture of a CT table and gantry.
Within the last fifteen years there have been significant advances in CT tech
nology, allowing for faster spiral acquisition and reduced dose to the patient as
well as multi slice detector arrays permitting simultaneous acquisition of several
slices at a time. Early CT gantries were constructed with a revolving detector
array positioned directly across the patient from the moving x-ray source. Cable
lengths connecting the moving detector assembly and the x-ray tube permitted
only a single slice to be acquired at one time. The revolving assembly then had
to be "unwound", the patient advanced the distance of one slice, and the process
repeated.
As small affordable detectors have become available, scanners have been de
signed with a fixed array of x-ray detectors. The only remaining revolving part
of the gantry is the x-ray tube. This has simplified the engineering and cabling of
the gantry. The x-ray tube can now be cabled using a slip ring, permitting con
tinuous revolution of the tube about the patient. This type of design is capable of
helical (or spiral) CT acquisition. By simultaneously revolving the x-ray source
about the patient and continuously moving the patient through the bore of the
gantry, the data are acquired via a spiral path. These methods have enabled very
fast image acquisition, improving the patient throughput in a CT facility, reducing
artifacts from patient motion, and reducing absorbed dose by the patient. Com
bined with multiple layers of sensors in the detection ring, these new scanners are
generating datasets of increasingly larger sizes.
A CT scanner is an x-ray modality and is subject to all of the physics as
sociated with the generation, dispersion, absorption, and scatter associated with
all x-ray photons. The process of generating a CT scan is similar to creating a
standard x-ray film; however, while a single x-ray exposure generates a complete
film-based exam, the CT image is not acquired in a complete form. Rather, it must
be reconstructed from multiple views. The advantage of being an x-ray modality
is that laymen and clinicians alike have considerable intuition when dealing with
x-ray-based images. The concepts of dense objects like bone absorbing more
photons relative to less dense tissues like muscle or fat come naturally from our
1 .3. Medical Imaging Technology 9
experience and expectations about x-ray imaging. A typical CT scanner can gen
erally acquire the data for a transaxial slice in a matter of seconds (within 1 to
5 seconds). An exam can include several series of slices, in some cases with
and without pharmaceutical contrast agents injected into the patient to aid in di
agnostic reading. Slices can be spaced such that they are either overlapping or
contiguous, though some protocols call for gaps between the slices. A large study
can include well over 1 00 separate 5 1 2 x 5 1 2 pixel images. The radiation dose
from a CT scan is comparable with that of a series of traditional x-rays.
The concept of "resolution" should be divided in the researcher's mind into
spatial resolution (i.e., how much area in each dimension a voxel covers) and
sampling resolution (how many voxels in each dimension of the slice). Sampling
resolution in current scanners usually creates images that are either 256 x 256
or 512 x 5 1 2 voxels square. Sampling resolution in the longitudinal direction is
limited only by the number of slices acquired. Spatial resolution in the longitu
dinal direction is bound by physical limitations of collimating the photons into
thin planes. The physical lower limit is approximately I mm; narrower collima
tion requires inordinate amounts of x-ray flux to image and also leads to diffrac
tion interference. Sampling resolution in the transaxial dimensions of a voxel is
based on the field of view selected by the operator of the CT scanner and the ma
trix (256 x 256 or 5 1 2 x 5 1 2) yielding pixel dimensions that are generally 0.5 to
2 mm. Attempting to achieve higher spatial resolution will lead to voxels with
too little signal to accurately measure x-ray absorption. With today's diagnos
tic equipment, if pixels smaller than 0.25 mm are attempted, low signal-to-noise
ratios become a problem.
Units of measure for the pixel values of CT imaging are standard across the
industry. Each pixel ideally represents the absorption characteristics of the small
volume within its bounds. By convention, these measurements are normalized
relative to the x-ray absorption characteristics of water. This is a unit of measure
known as Hounsfield units (HU). The Hounsfield unit scale is calibrated upon
the attenuation coefficient for water, with water reading 0 HU. On this scale air is
-1000 HU, fat tissue will be in the range of -300 to -100 HU, muscle tissue 1 0-70
HU and bone above 200 HU.
Every CT image contains artifacts that should be understood and handled
when visualizing the data. The process of reconstruction and sampling leads to
aliasing artifacts, just as any sampling procedure in computer graphics or signal
processing. Another common artifact is partial voluming, the condition where
the contents of a pixel are distributed across multiple tissue types, blending the
absorption characteristics of different materials. Patient motion while scanning
generates a variety of blurring and ring artifacts during reconstruction.
There is another class of artifacts that arises from the fact that CT is an x-ray
modality. Embedded dense objects such as dental fixtures and fillings or bullets
lead to beam shadows and streak artifacts. More common is the partial filtering
10 Chapter 1 . Introduction
of the x-ray beam by dense tissue such as bone which leads to beam hardening,
a condition which causes slight shadows (0 halo the dense features of the image.
When creating visualizations, beam hardening anifacts may cause the researcher
to underestimate the volume of dense objects, making them seem smaller than
they are.
trast and resolution. Figure 1 .2 is a picture of an MRI magnet and table, showing
a head coil, an antenna for imaging the head and neck.
Unlike CT scanners, the bore of an MRI scanner is often up to two meters in
length (6-8 feet). Patients are inserted into the middle of the magnetic field, often
inducing claustrophobia. The environments for MRI scanners must be shielded
for magnetic and radiofrequency interference. Large Faraday cages and substan
tial masses of iron usually surround the magnet and sometimes the entire room
about the scanner.
In MRI scanning, the patient is placed within a high intensity magnetic field.
Field strengths vary from 0.35 Tesla to 1 .5 Tesla for most diagnostic MRI devices
(for reference, 1 Tesla = 1 0,000 Gauss, and the earth's magnetic field, though
variable, is approximate 0.5 Gauss). The induced magnetic field causes the mag
netic moments of the hydrogen atoms within the patient to align along the prin
cipal direction of the superconducting magnet. Low-level radio waves in the mi
crowave frequencies (approximately 1 5 to 60 MHz) are then transmitted through
the patient, causing the magnetic moments of the hydrogen nuclei to resonate and
re-emit microwaves after each pulse. The microwaves emitted by the body are
recorded using a radio frequency antenna, filtered, amplified, and reconstructed
into tomographic slices. While all of the hydrogen nuclei typically resonate at
a frequency fixed by the strength of the induced magnetic field, different tissue
types resonate longer than others, allowing the viewer to discriminate among them
based on the magnitude of the signal from different points in space over time.
Spatial locations can be determined by varying the magnetic field about the
patient in different directions at different times. Linear gradients in the induced
magnetic field are produced by magnets that supplement the main superconduct
ing magnet. These gradient coils are activated in a variety of sequences, altering
the phase and frequency of the microwave pulses that are received and re-emitted
by the hydrogen nuclei. The design and crafting of pulse sequences is a field of
research unto itself. The careful selection of pulse sequences can illuminate a
variety of clinical conditions including the display of swelling and bleeding, even
enhancing the blood vessels deep within the brain. Often, several sequences will
be taken of the patient during an exam to capture the variety of information avail
able in the different pulse sequences. The visualization researcher is encouraged
to learn about these pulse sequences and their properties in order to select the best
ensemble of them to use when crafting visualizations.
The output of an MRI scanner is similar to CT. Slices representing slabs of
the object scanned are produced. However, unlike CT which always produces
transaxial slices, the slices from MRI can be oriented in any plane. The output
values at each image element are not calibrated to any particular scale. Gen
erally they are l O-bit data samples. The values will vary depending upon the
scan parameters, and the patient's size and magnetic characteristics. Addition
ally, the values are not constant over the entire scan space since inhomogeneity
12 Chapter 1. Introduction
in the magnetic field causes pixels that may represent the same tissue, but located
some distance apart to give different signals. This lack of an absolute scale for a
dataset is a cause of much consternation to the researcher attempting to segment
MRI data.
Visualization researchers seldom have good intuition for the meaning of MRI
signals. Unlike x-ray-based imaging, there are no physical analogs to what the
viewer is seeing. An MRI device measures the radio signals emitted by drops of
water over time. Usually skin and fat are much brighter than bone which has vir
tually no signal at all. Segmentation and classification are therefore significantly
harder and are the subject of much research.
As with CT, the concept of resolution must be divided into spatial resolu
tion and sampling resolution. The matrix of the image can typically be selected
to be either 256 x 256 or 5 12 x 5 1 2, depending on the sampling resolution of the
frequency and phase decoders of the receiving equipment. Square pixels in a rect
angular matrix are often selected since reducing the number of rows or columns
can significantly reduce the time required to complete a scan. Spatial resolution,
or field of view, is primarily dependent on the strength of the gradient magnets
and their ability to separate the slices along their gradient directions. The ra
diofrequency receiving equipment must be able to distinguish among frequencies
and phases which have only slight variations. Stronger gradients create greater
separations and improve spatial resolution.
Remember that the MR scanner is attempting to measure the radio signal res
onating from a drop of water that is as small as 1 mm x 1 mm x 2 mm. There
are significant trades to be made when attempting to increase spatial resolution.
Signal to noise will degrade as electronic distortion of the antenna coil and the
amplifiers begin to overcome the small signals involved in diagnostic imaging.
Imaging with larger voxels, thicker slices, or repeating and averaging multiple ac
quisitions are solutions. The relative consequences are increased aliasing, partial
voluming, and possible patient motion artifact.
MRI data is subject to several artifacts. The issues of partial voluming, patient
motion, and aliasing are common with CT. While MRI does not have x-ray-related
artifacts, it has its own host of radio and magnetic artifacts of which the visual
ization expert should be aware. Like any microwave used for cooking, the mass
of the patient (or food) will affect how well the material absorbs the radiofre
quency energies. This leads to "cold" spots in the imaging volume. In addition,
the distribution of the antenna coverage for both transmission and reception and
inhomogeneity in the induced magnetic field lead to inconsistent quantitative val
ues and even inaccurate geometry within an image. There is significant research
being conducted addressing methods to correct for these distortions.
Most MRI scanners can acquire a set of slices (30 to 50) within five to ten
minutes. An entire study of a patient generally represents two to three series of
slices, with a study time of 30 to 45 minutes. Each slice generally represents a
1 .3. Medical Imaging Technology 13
thickness of 2 to 1 0 mm and contains 256 x 256 pixels. As with CT, the pixel
dimensions are set by the image matrix and the field of view parameters.
New work in MR imaging is leading to novel capabilities in this modality.
Studies in perfusion and diffusion of various agents across the body are being
enabled by new capabilities in fast imaging. These studies can illuminate blood
flow, creating high-resolution images of vascular structure. Functional MRl has
been introduced not only to record the patient's anatomy, but also the physiolog
ical functions of the tissues being studied, largely used today to map the cerebral
cortex of the brain. The complex pulse sequences of Diffusion Tensor Imaging (or
DTI) are beginning to reveal the direction and course of nerve bundles deep wi thin
the white matter of the brain, providing new research areas of tractography and
patient specific anatomical analysis. The pace of research in MRl is increasing as
new capabilities are enabled.
Figure 1.3. A two-headed Gamma Camera used in nuclear medicine (circa 1998).
1 .3. Medical Imaging Technology 15
1 .3.4 U ltrasou nd
Unlike the three other modalities described in this chapter, Ultrasonography uses
classical physics to perform its imaging rather than the more esoteric phenomena
of nuclear decay or x-ray quantum mechanics. No photons are used. A piezoelec
tric quartz crystal, similar to those used to clock CPUs, wristwatches, or to start
electronic ignition lighters for barbeque grills, is used to create high frequency
acoustic energy (3 to 10 Megahertz) which is then reflected off surfaces and in
terfaces between organs deep within the body.
The same device that is used to create the acoustic signal, a transducer, is then
used to measure the returning echo information. Partial reflections are created by
interfaces between areas of differing acoustical impedance. The result is imaging
very similar to the SONAR systems used in maritime and undersea naval imag
ing. The sonographer places the probe against the patient and moves it to obtain
images of various parts of the body. The result is usually a 2D representation of
the area under the transducer. Most ultrasound machines consist of a linear ar
ray of transducers, and produce an image representing a pie-shaped slice of the
body. One of the advantages of ultrasound is that it produces images in real time.
Another advantage is the absence of ionizing radiation.
16 Chapter 1 . Introduction
Ultrasound machines are fairly inexpensive compared with the other diagnos
tic 3D modalities in common use in the health care industry. A high-end diag
nostic ultrasound system can be purchased for approximately $250,000, while
functional, more moderately priced units can be acquired for less than $1 00,000
today.
3D images can easily be created by extruding several slices through space.
Most commercial ultrasound equipment allows for 3D imaging; however, accu
rate spatial information is seldom provided. Spatial tracking of the transducer
array is imperative for clinical 3D visualization. Some approaches to correcting
this defect have been to use physical, optical, and electronic devices for locating
the transducer in space. Providing a rigid path for the transducer creates some
compelling results. Other techniques involve rotating the transducer assembly,
sweeping a volume in a cylindrical section similar to aeronautical RADAR imag
ing. Finally, recent advances in transducer design have yielded demonstrable 3D
transducers that can acquire volume images from ultrasound with no mechanically
moving parts.
However, once the position and orientation of the probe are known, the data
are still often sampled in irregular intervals, and defy many image processing
techniques. Ultrasound images typically contain a large amount of noise termed
speckle, that adds to the problem of identifying structures. Object detection in
volume ultrasound is the subject of much advanced medical image processing
research.
ested in medical visualization would be well served to learn where these trades
are being made and learn to cope with their consequences.
An even deeper understanding of the clinical findings of a case is required to
validate a 3D visualization. Knowledge of what the physician is seeking in the
patient will help to direct the acquisition so that the results not only capture the
desired pathology in high detail, but also assure the computer scientist that the
data are in an appropriate form for advanced rendering.
In his 1 993 presentation on 3D Medical Visualization, Derek Ney, an assistant
professor for the Johns Hopkins Department of Radiology, wrote:
Indeed, as with all data processing, "if garbage goes in, only garbage comes
out." Or more precisely, if one beautifully renders garbage data, then one only
has a picture of well-dressed garbage. The result serves no one in particular and
wastes valuable resources that would be better used to serve the patient, to serve
the community, and to serve future generations of researchers and patients alike.
1 .5 S u m mary
Keeping the basic concepts of medical image acquistion in mind, we tum our
attention to the processing and analysis of the generated data. The multiplica
tion of imaging modalities with differing strengths and weaknesses in imaging
anatomy v. physiology and hard tissue v. soft tissue requires careful registration
and alignment of multimodal data, and differences in size and resolution lead to
multiscale methods as well. We refrain from covering the algorithms and methods
for generating visualizations of the data, instead concentrating on the analysis and
refinement of the data stream in preparation for the ultimate creation of rendered
images of the patient anatomy and physiology.
This book is divided into three main parts:
• Basics: These chapters describe techniques for preprocessing the image
data before complex semantic, know lege-based operations are performed
either by the machine or the human operator. Statistical, linear, and more
18 Chapter 1 . Introduction
Terry S. Yoo
National Library of Medicine, NIH
George D. Stetten
University of Pittsburgh
Bill Lorensen
GE Global Research
2.1 I ntroduction
The topics discussed throughout this book fall into the realm of what is usually
called image analysis. Before researchers adopted that term, they often referred
to their field as image processing, an extension of traditional signal processing
to images. An image can be thought of as a signal in space rather than time,
arrayed in two or three dimensions rather than one. This leads to fundamental
differences between conventional signal processing and image processing. While
one generally only moves forward through time, in space, one moves left or right
with equal ease, as well as up or down, forwards or backwards. The extension
to multiple dimensions raises new issues, such as rotation and other interactions
between the axes.
As with conventional signals, images can be either continuous or discrete.
Optical images, consisting of light, call for analysis in the continuous domain, as
they pass through lenses, for example. In contrast, images found on a computer
19
20 Chapter 2. Basic Image Processing and Linear Operators
2.2 I mages
In this book, we often are concerned with 3D, or volume, images. A volume
image, cp, is a mapping from �3 to �n where n = 1 for the typical, scalar-valued
volume. More precisely:
where U is the domain of the volume. The image cp is often written as a function
cp(x,y,z).
Throughout this chapter, many examples will be presented as either ID or 2D,
where they may be more easily examined (at least in print), and then generalized
to higher dimensions. In general, for a specific number of dimensions, we will
use the corresponding function notation (cp(x) , cp(x,y) , or cp(x,y,z) ).
2.3. Point Operators 21
(2.2)
We likewise use notation to distinguish between discrete and continuous differ
ential operators. To find the partial derivative in the x-direction of a discretely
sampled image we can say that
(2.3)
oxFiI,} =
-
. Fi+l,j - Fi-l ,j Fi+ l,j - Fi- l,j (2.4)
Xi+l - Xi-l 2h
to approximate the partial derivative, where h is the grid spacing (nonnally as
sumed to be 1 pixel). The value OxFi,j is an approximation of the instantaneous
partial derivative of <1>(Xi,Yi) at (Xi ,Yi ) in the x-direction. Equation (2.4) is the
method of central differences, one of a number of commonly used methods to
approximate the derivatives of unifonnly, discretely sampled datasets.
2 .3 Po i n t Operators
The simplest operators on an image are those that treat each pixel independently
from its neighbors. Point operators usually ignore infonnation about pixel loca
tion, relying only on pixel intensity. A surprising number of useful things can be
done with such operators. Some examples follow.
2.3.1 Thresholding
Figure 2.1. A simple threshold segments the bone from other tissue in a cr scan of a
fractured tibia.
2.3.3 Windowing
Often it is advantageous to change the range of intensity values for an image to
adjust the brightness and contrast of what will be displayed on a screen. Dis
play is an essential pan of image analysis, because it is the portal 10 the human
visual system, whose capabilities far surpass any computer system likely to be
2.3. Point Operators 23
Figure 2.2. Histogram h( <1» of pixel intensity <I> for tissue types A and B, with optimum
threshold for segmentation (arrow).
built in the near future. The dynamic range of image display systems, between
the brightest and darkest pixel the screen can produce, is quite small compared
to the dynamic range of the human eye. Many medical imaging modalities are
also capable of a greater dynamic range than can be displayed. The display itself
therefore represents a bottleneck in the overall system. To minimize loss of in
formation, the brightness and contrast of the image may be adjusted so that the
brightest and darkest pixels in the image exactly match the extremes of the dis
play, as shown in Figure 2.3. This is a linear function remapping the intensity of
each pixel, with contrast being the slope of the function and brightness being the
offset. In some clinical systems, the process is called windowing.
h ( <I»
max '"
Figure 2.3. Windowing intensity histograms. Left: Histogram h(<1» of pixel intensity
<I> of original image, on a scale from zero to the maximum value of the display; Center:
windowing function 'V 1(<1» to yield a new intensity 'V; Right: new histogram
= h('I') in
which the dynamic range of the image matches that of the display.
24 Chapter 2. Basic Image Processing and Linear Operators
A further step may be taken to optimize the match between the image and the
display system, which ensures that each level of pixel intensity is equally rep
resented in the image. This is called histogram equalization, and although it is
non-linear, we include it here as a useful point operator. As shown in Figure 2.4,
histogram equalization entails finding a monotonic function ", = f ( cj» that remaps
pixel intensity so that the histogram becomes a constant function of intensity ",.
Color video images typically have three separate channels to record the red, green,
and blue levels from the camera. Color displays are capable of communicating
these images to the human visual system, but such displays also can be used to
enhance the display of single-channel gray-scale images through the use of colors,
artificially assigned to particular pixel intensities. Such color maps are often used
to overlay additional information in anatomical images, such as Doppler flow in
ultrasound images. The color map usually follows some perceptually consistent
ordering scheme, such as the rainbow, or especially bright colors may be assigned
to emphasize pixel values above certain levels.
Some point operators combine two images, keeping each pixel separate from its
neighbors, while combining it with the corresponding pixel in another image.
Algebraic point operators may be used to add images together, for example, in
averaging many images acquired from the same sample to reduce noise. Or they
2.4. Linear Filtering 25
may be used to subtract one image from another, such as in Digital Subtraction
Angiography where the difference between sequential fluoroscopic images, be
fore and after contrast, leaves just the contrast without the anatomical structures.
Another use of algebraic point operators is to mask an image by multiplying it
by a second, binary image, containing ones where the first image is to be passed
through, and zeroes where it is to be masked.
All the point operators described so far apply the same operation to every pixel
in the image. In some cases it is useful to make the operation depend upon pixel
location, for example, to correct for spatial inhomogeneity in the image acquisi
tion. Pixels are kept independent of each other, but in a framework where location
within the image makes a difference.
Although many useful things can be done to an image using point operators on
individual pixels, most image processing tasks require the simultaneous consider
ation of multiple pixels. Techniques that combine multiple pixels in a linear and
space-invariant manner are known collectively as linear filtering. Two standard
techniques for linear filtering, convolution and the Fourier transform, comprise
the rest of this chapter.
If you have ever adjusted the focus on a camera you have performed convolution
in the continuous domain. If you have run a blurring function or performed edge
enhancement using a commercial graphics program, it is likely that you were
using discrete convolution. In image processing, convolution filters are used to
make measurements, detect features, smooth noisy signals, and deconvolve the
effects from image acquisition (e.g., deblurring the known optical artifacts from
a telescope). The effect of the particular filter depends on the nature of the input
filter, or kernel.
Convolution, denoted with the operator 181 can be described in I D as a contin
uous operation applying the filter kernel h(x) (itself an image) to an image <j>(x)
using the following integral form:
Note that the expression <j>(x) 181 h(x) itself describes an image mapping. Thus
the following equality holds:
26 Chapler 2. Basic Image Processing and linear Operators
The expressions !j)(x) 0 h(x) and ($ 0 h)(x) are used interchangeably, with the
position index x often omitted to help streamline and clarify the nOlation espe
cially when higher dimensional (e.g., 3D) images are involved. In the practice of
digital image processing, hex) typically does nOI have infinite extent, nor is it in
finitely dense; however, when working in the continuous domain of linear filtering
theory, considering kemels that are infinitely wide (thus avoiding the complica
tions of truncation) is often morc convenient. We will elaborate later on discrete
functions with finite extent in the discussion on sampling.
We can generalize the I D convolution operation to 20 and 3D images. Thus,
in 20, convolution becomes
$(x,y,z) ®h (x,y,z) =
Figure 2.5 shows a noisy 1 D step discontinuity, or edge. The function $(x) is
shown from x = 0 to x = 255 with a step function at x = 128. The signal to noise
ratio is approximately 4: I.
Figures 2.6. 2.7. and 2.8, show a progression of filtered versions of the input
from Figure 2.5 with Gaussian filter kernels of increasing aperture or width a.
The Gaussian function g(x) is defined as
, �
g {x,a) = p;=:: e-2a! (2.9)
v 27ta
0.'"
• ,.
.,
(,) (b) (c)
Figure 2.6. Gaussian filtering of input �(x) where (I = 16. (a) Input 4I(x); (b) a ID
Gaussian kernel, 8(X,(I) where (I = 16; (e) output of �(X)0g(x, 16).
'''1..
:��J��� �
.,
•
�
., ,. ••
Figure 2.7. Gaussian filtering of input 4I(x) where (I = 24. (a) Input 4I(x); (b) a ID
Gaussian kernel, g(x,(I) where (I = 24; (e) output of 4I(x) <8> g(x, 24).
'"
Figure 2.8. Gaussian filtering of input 4I(x) where (I 32. (a) Inpul ¢I(x); (b)
= a lD
Gaussian kernel. g(x,(I) where (I = 32; (e) Output of,(x) <8>g(x, 32).
28 Chapter 2. Basic Image Processing and Linear Operators
Notice how increasing the aperture of the Gaussian decreases the noise but
also blurs the edge. This trade-off of resolution for noise reduction is one of the
many considerations in the design of linear filter systems. An entire approach
to image analysis, called "scale space," depends upon the Gaussian to produce a
broad range of blurred versions of a given image.
The nature of a linear filtering operation depends on the properties of the filter
kernel. For instance, the "shape" of a 2D or 3D kernel will determine whether
the operation remains invariant with respect to rotation. Independent of kernel
shape, however, the convolution operation has many useful properties, including
the following:
Convolution is linear:
Convolution is commutative:
p 0 q = q 0 p. (2. 1 1)
Convolution is associative:
(p 0 q) 0 r = p 0 (q 0 r) . (2. 1 2)
p 0 (q + r) = p 0 q + p 0 r. (2. 1 3)
These combined properties create the justification for using convolution as the
principal operation in linear filtering.
(2. 14)
Equation (2. 14) depicts the differential operator as a kernel by which convolution
accomplishes differentiation. This is hard to illustrate using the ideal differential
2.4. Linear Filtering 29
� r\��� ..•
•
.-
.-
i\�� ..
.. ..,
jk - -
....., .• •
Figure 2.9. Taking derivatives of a noisy input by convolution with the derivative of a
Gaussian kernel. (a) I-D input 41(-,); (b) tg(-" a) where a = 3; (c) 41(-,) 0 tg(-" O").
kemel, because such a kernel has infinitesimal width and infinite height. In the
real world such a kernel is impossible. In any event, because differentiation en
hances high-frequency noise, it is often necessary to regularize (smooth) a noisy
image before computing its derivative, by first convolving the function with a
smoothing kernel h{-,). It follows from the associative and commutative proper
ties of convolution thaI
(2.15)
In other words, the derivative of a function $(x) convolved with a filter kernel
h(x) is equivalent to convolving 1/I(x) with the derivative of h(x). This suggests
that one of the easiest ways to compute the derivative of a function is through con
volution with the derivative of some smoothing kernel. We have already demon
strated the use of the Gaussian g(x) as a smoothing kernel. The Gaussian's in
finitely differentiable propenies make it attractive for taking derivatives as well as
smoothing.
Figure 2.9 depicts a noisy 10 input signal for which derivative infomJation is
desired. When convolved with the derivative of a Gaussian, the resulting output
repons the derivative of a smoothed version of the input. By the commutative
and associative properties of convolution, it can just as well be described as a
smoothed version of the derivative of the input. This technique for differentiating
functions can be extended to higher order derivatives, and to higher numbers of
dimensions.
kernel. In 2D, a discrete convolution of image P with finite kernel Q looks like
domainx[Q] domailly[Q]
Px,), (9 Qx,), = L L Px-i,y-j Qi,j. (2. 16)
i j
The discrete version of convolution shares all the above-mentioned attributes
of its continuous cousin. Convolution of discretely sampled data is linear, commu
tative, associative, and distributive over addition. As in the continuous domain,
convolution in the discrete domain can be used to smooth noisy data or detect
boundaries, depending on the choice of kernel. An example of discrete convo
lution for differentiation in 2D has already been seen in the central differences
[
operation (Equation (2.4)). The kernel in that case would be
0 0
-0.5 0 (2. 17)
o 0
which has the effect of subtracting the pixel to the left from the pixel to the right,
just as Equation (2.4) specifies, to represent the x-component of the gradient.
The following section shows more examples using convolution for smoothing and
taking derivatives in 2D.
One particularly useful kernel for discrete convolution is the binomial kernel, gen
erated by repeated convolutions with a simple box of identical values, as shown
here:
[: :H� �n
[ 11 , (2. 1 8)
With successive iterations, the binomial kernel approaches a Gaussian shape (by
the Central Limit Theorem). The scale of the Gaussian is determined by the
number of iterations. (The kernels shown in Equation (2. 1 8) would be normalized
by 114 with each iteration.)
Boundaries in an image represent areas of high gradient magnitude, i.e., the image
intensity increases or decreases rapidly across the edge. The gradient vector,
(2. 19)
2.4. Linear Filtering 31
is oriented in the direction of the steepest change in image intensity, normal to the
implied boundary. The gradient magnitude,
(2.20)
[ 1
gaussians (DoOO).
-1 -2 0 2 1
-1 o 1 ], -2 -4 0 4 2 . (2.21 )
-1 -2 0 2 1
The two DoOO kernels shown in Equation (2.2 1 ) correspond to iterations 1
and 3 in Equation (2. 18). Each kernel is the difference between two copies of
the corresponding binomial kernels displaced along the x-axis. Convolution with
these DoOO kernels measures the gradient component in the x-direction. Simi
lar kernels can be constructed by displacing copies of the binomial kernel in the
y-direction to measure the y-component of the gradient. DoOO kernels can be
constructed for images in 3D or higher dimensions. The number of iterations of
the underlying binomial kernel determines the scale of the DoOO.
Figure 2.10 shows the results of detecting the individual gradient components
as well as the gradient magnitude of a simple image containing a rectangular
object. Notice the orientations of the edges detected by the individual gradient
components, and the orientation-independence of the gradient magnitude.
We have seen how the gradient vector represents the strength and orientation of
the boundary. But what of higher order differentials? Where the first derivative
(gradient) is represented as a vector, the second derivative is a matrix, known as
the Jacobian
[ � a1
dYdx
a2 e!> a2 e!>
ay
•
The elements of this matrix are useful in a number of ways. For example, the
(2.22)
Laplacian,
(2.23)
32 Chapter 2. Basic Image Processing and linear Operators
Figure 2.10. Gradient components and gradient magnitude: (A) original image; (8) x
component of gradient; (C) y-component of gradient; (D) gradient magnitude.
sums the diagonal tenns of the Jacobian to yield a rotationally invariant repre
[
sentation of the second derivative of image intensity. A common 3 x 3 kemel
representing the Laplacian in 20 is
1
\ -\ \
-\ 8 -\ . (2.24)
-\ -\ -\
Convolution with this matrix yields zero along a boundary, no matter what
orientation the boundary has. For example, convolution with the following 3 x 3
patch of an image containing a diagonal boundary (between regions with intensi
ties of I and 3 respectively),
[: ; n
yields zero for the center pixel (which is directly on the boundary),
(2.25)
af
f
ax
x x
Figure 2.1 1 . ID case of the intensity function across a boundary, its first derivative, and
its second derivative showing a zero crossing at the boundary location.
In two or more dimensions, we use the Laplacian to capture the "total" second
derivative at a pixel. Mathematically, the Laplacian is the divergence of the gra
dient. Divergence is a common concept in fluid dynamics, describing how much
more fluid enters a region than leaves it. In our case, it describes how much more
gradient is "entering" a pixel than "leaving" it from any direction (making it ro
tationally invariant). A boundary exists where there is no net change in gradient
(i.e., the gradient is at a maximum on the boundary), making the Laplacian zero,
just as in the ID case. The difference of Gaussian (DOG) kernel, which we will
discuss in Section 2.4.9, yields results similar to the Laplacian, because the DOG
kernel has a similar shape to the Laplacian.
Let us now consider the terms off the diagonal of the Jacobian matrix (Equa
tion (2.22» . These can be used to measure the curvature of boundaries. The
individual components of the gradient at a curved boundary change as one moves
orthogonally to that particular component. This type of change shows up as partial
second derivatives of intensity off the diagonal of the Jacobian. For illustration,
Exploring the Variety of Random
Documents with Different Content
Clive, I am positive that the man was none other than the original of
this likeness, and——"
He was interrupted by a passionate cry of pain and anger, and
Ralph, snatching the photograph from his hand, stood confronting
him with blazing eyes.
"It is false!" he cried. "You know it is false! I believe that you are
responsible for my father's disappearance!"
CHAPTER III
MR. ST. CLIVE PROVES HIMSELF A TRUE
FRIEND
"He is a fine young fellow, but his past life has been spent amidst
very different scenes, and he is far from having a fitting education.
But he is very intellectual and will acquire knowledge quickly. His
father must have been a gentleman, and he has taught his son to be
one also."
It was Mr. St. Clive who spoke, and his words were addressed to Dr.
Beverly, the principal of Marlthorpe College—the best school in all
the county.
A fine-looking man was the doctor, tall, erect, dignified, with firm
face and piercing eyes—eyes which could look terribly severe when
their owner was angry, but which otherwise were gentle, and even
mirthful.
Dr. Beverly was proud of his school, but prouder still of his work. He
did not labour to make scholars only, but also to build up men—
good, noble men—who should be a credit to the old school, and a
blessing to their country. Work or play, the doctor believed in
everything being done as well as it could be, for his watchword was
"Whatever you do, do it to the glory of God," and nothing can be
done to God's glory that is not done as well as it possibly can be.
Mr. St. Clive had explained how Ralph came to be under his care,
and had told the doctor how much he owed to him; and he finished
by mentioning the cruel statement which Lord Elgert had made, and
the angry way in which Ralph had answered it.
"I tell you this," he said, "that you may know everything. I attach no
weight to Elgert's statement myself—it is too absurd, but you must
exercise your own discretion," and the doctor smiled slightly.
"Lord Elgert is rather prone to make rash statements," he said. "I
shall be quite willing to receive your young friend, and I will do my
best to turn him into a good man."
"That I am sure of," was the hearty reply, "and I am also sure that
you will have good material to work upon. Then I will bring Ralph
over."
"And do you propose that he shall board here entirely, or return to
you every Saturday, as most of the lads do?"
"Oh, come home. That is how I did in my day—you know I want to
watch the boy. Good-day, doctor," and Mr. St. Clive came away.
Marlthorpe College was a splendid old building, with large playing
fields at the back, and a great quadrangle in front, to which entrance
was gained through a pair of great iron gates, against which the
porter's lodge was built.
The school itself was at the other side of the quadrangle, directly
facing the gates—a two storey building, with the hall, in which the
whole school assembled upon special occasions, below, and with the
classrooms above. It had two wings; the one to the right being the
doctor's own residence, and that on the left the undermaster's
quarters.
At the back there were again buildings on the right and left—on the
left junior dormitories, the dining-hall, and matron's rooms; and on
the right senior dormitories and studies.
Mr. St. Clive drove home and told Ralph the result of his visit.
"I am sure that you will like the doctor," he said, "and you will find
your companions a nice lot of fellows. Of course there will be some
unpleasant ones; and Ralph, if things are as they used to be, you
will find that there are two sets of fellows—those who mean to work
honestly, and those who never intend to take pains. I need not ask
which set you will belong to," and Mr. St. Clive smiled. "But now," he
added, "I want you to try and be brave. You have a very terrible
sorrow, I know; and it is hard to put it from my mind——"
"It is never from my mind, sir," interrupted Ralph sadly. "I am always
thinking of it."
"But you must not brood over it. To do that, will unfit you for all else.
Leave it with God, Ralph, and do not let even so great a grief
interfere with life's duties. Will you promise me to try and remember
this?"
"I will indeed, sir," answered Ralph. "If I have lost father, I mean to
try and think that he knows, and just do that which would please
him."
"That is good; but still better is it to remember that we have to try
and do that which shall please our Heavenly Father. Now, Ralph, I
suppose that out where you made your home, blows often were the
only way of settling troubles. I do not say that blows are never
justifiable, for sometimes we are placed in such circumstances as
warrant fighting, but do not be too ready to quarrel, or to avenge
every fancied insult with your fist. But there, I am sure that I can
leave that to you. Now come to lunch, and then we must see about
starting."
"I am so glad that you are coming home every week, Ralph," so said
Irene St. Clive, when she heard of the arrangements which her
father had made. "My own lessons are finished on Friday, and we
can have all Saturday to ourselves. I shall count all the days until
each Saturday comes."
So with kindly words to cheer him on his way, Ralph started off with
Mr. St. Clive, and was introduced to Dr. Beverly; and Ralph felt that
he liked the doctor from the very first moment that he saw him; and
he determined that he would do all that he could to get on and
prove to Mr. St. Clive that he meant to keep his word.
Then when his friend had gone, the doctor questioned Ralph to see
just what he knew; and at the conclusion of the examination he laid
his hand on his shoulder.
"My boy," he said, "it is my desire always to have the fullest
confidence in my scholars, and also to enjoy their confidence. I want
you to remember that I desire to be your friend as well as your
master, and that out of school hours I am always glad to see any of
my boys who want to talk with me. I do not mean who want to
come tale-bearing," he added, and Ralph smiled as he answered—
"Thank you, sir. I think I understand."
"You will have to be in the Fourth Form at first, that is the lowest
Form in the Senior House," the doctor continued. "But if you work
well, you will soon be in the Fifth. Now, if you will come with me I
will introduce you to your master, Mr. Delermain, and I think you will
find him ever ready to help you in any way he can."
Ralph thanked the Head again, and followed him, with more of
curiosity than of nervousness, to make the acquaintance of the boys
with whom he was to study; and twenty pairs of eyes glanced up as
the Head opened the door, and then dropped as quickly when they
saw who had entered.
But the master rose from his seat and came forward to meet the
doctor, who said, patting Ralph on the shoulder—
"I have brought you a new scholar, Mr. Delermain. This is Ralph
Rexworth, and he is the young gentleman of whom you have heard
—the one who saved Mr. St. Clive's daughter." Hereat the eyes were
stealthily raised, and glances of something like respectful awe
followed. Of course every one there had heard of the incident about
the bull, and of the disappearance of Mr. Rexworth.
"Rexworth is rather backward," the Head continued. "His life has
been spent abroad, and he has not had the opportunities for study;
but I believe that he will soon pick up." And with this Dr. Beverly
went, and Mr. Delermain, having spoken a few words of welcome,
beckoned to a boy to come forward.
"Warren, let Rexworth sit beside you this afternoon, and give him a
set of the sums we are doing. If you find them too difficult," he
added to Ralph, "do not hesitate to come to me."
But Ralph did not need to ask for aid, he could do the sums and the
exercises that followed. Indeed, he did better than some who had
been there longer, notably one big lad with a sickly flabby face, who
was seated at the bottom of the class, and who received a
reprimand from his master for his indolence.
"It is shameful, Dobson! Here, a new boy has done better than you
have. Your idleness is disgraceful."
A writing exercise followed; and Ralph was bending over his book,
when flop!—a wad of wet blotting-paper hit him in the cheek. He
looked up, but every one seemed busy with their work, so wiping his
cheek he put the wet mass on one side, and went on with his task.
Flop! A second wad came. Ralph noted the direction, and saw that at
the end of the form Dobson was seated, and Ralph had his
suspicions. Pretending to be absorbed in his work, he kept a covert
watch; and presently he was rewarded by seeing Dobson extract a
third wad from his mouth, where he had been chewing it into a
convenient pellet, and under cover of the boy in front of him prepare
to fire it by a flick of his thumb. Ralph raised his eyes and looked
him full in the face, and, somehow, Dobson seemed confused. He
turned red, and bent over his work hastily; and no more pellets were
fired at Ralph that afternoon.
It seemed rather a wearisome afternoon to the boy, used as he was
to his open-air life, but he worked away with all his might; and
presently the bell rang and work was over; and then Warren, the
boy beside whom he had sat, came to him and held out his hand.
"I am first monitor of our form," he said, "and I hope that we shall
be friends. If you come with me I will take you round the school."
"Rexworth."
Ralph turned as his name was called; his master stood there.
"I want you a few minutes. Warren, you can take him round
afterwards. I want to arrange about his study."
"We have only got one vacant, sir," the monitor said. "Charlton has
that."
"I know," was the quiet answer; and then, when Warren ran off, the
master turned to Ralph.
"Rexworth," he said, "I must explain that in our form every two boys
have one study between them, and as you heard Warren say, we
have only one study that is not fully occupied. A lad named Charlton
has it, and you must chum with him. It is about him I want to speak
to you."
"Yes, sir," said Ralph, wondering why his master spoke so gravely.
"Rexworth, I am sorry to say that Charlton is not quite in favour with
his schoolmates. His father got into some trouble and has
disappeared—it is supposed that he is dead—and the boy managed
to gain a scholarship at another and poorer school, and has come
here. He is a real nice lad, but very weakly and timid, and the others
put upon him, partly on that account, partly because of his father's
disappearance, and partly because he is poor—a sad crime in the
eyes of many. It would have been wiser, I think, if he had not come
here, but Dr. Beverly wished him to do so. I wish, Rexworth, that
you would try to be his friend, for he needs one; some of the lads
are nice enough to him, but he seems so very much alone."
"I would like to help him, sir," was the ready answer. And the master
smiled.
"I thought that I was not mistaken in you," he said. "Look, there the
lad is. Charlton, come here."
The lad came up. He was a pale boy, very delicate in appearance,
and with a sad, wistful face.
"Yes, sir," he said.
"Charlton, there is only one vacancy in our studies, and that is with
you. Rexworth will have to chum with you." The boy cast a startled
glance at Ralph. "Take him and show him where it is, and try to
make him feel at home."
"Yes, sir." The boy beckoned to Ralph. "Please come with me," he
said, in troubled tones, as if he doubted whether Ralph would care
about sharing the study with him.
"Have we got to be chums?" asked Ralph; and the other boy
nodded.
"Yes. That is what we call it. It means sharing studies; but you need
not speak to me if you don't want to, and I will not be in the study
much. I am not as it is, for they are always disturbing me and
spoiling my things."
"They! Who?" demanded Ralph; and the lad answered—
"The other chaps and the Fifths. Dobson, in ours, and Elgert of the
Fifth, are the worst. They go in and spoil my things."
"They have no business to, of course?"
"Go in? No, of course not—only the two who chum have any right in
it. Here we are, and—there, they are in now!"—as a scuffling and
burst of laughter came from the inside of the study before which the
boy had halted. "Oh, what are they doing! Will you stop until they
have gone?"
"Not I," answered Ralph grimly. "That study is mine as well as yours,
and I mean to see that we have it to ourselves, Charlton. Come on,
and we will see what is up." And saying this, Ralph threw open the
door and walked into the little room, followed by his companion.
CHAPTER V
MAKING THINGS STRAIGHT
ebookgate.com