NIST Technical Note 1877

Quantitative Imaging of Clinker and 

Cement Microstructure 

Paul E. Stutzman

Jeffrey W. Bullard

Pan Feng

This publication is available free of charge from:

http://dx.doi.org/10.6028/NIST.TN.1877
NIST Technical Note 1877

Quantitative Imaging of Clinker and 

Cement Microstructure

This publication is available free of charge from:

http://dx.doi.org/10.6028/NIST.TN.1877

April 2015 

U.S. Department of Commerce 

Penny Pritzker, Secretary 

National Institute of Standards and Technology 

Willie May, Acting Under Secretary of Commerce for Standards and Technology and Acting Director

Paul E. Stutzman1

Jeffrey W. Bullard1

Pan Feng1,2

1

Materials and Structural Systems Division

Engineering Laboratory

2Southeast University, Nanjing, Jiangsu 210096, China 
Certain commercial entities, equipment, or materials may be identified in this

document in order to describe an experimental procedure or concept adequately.

Such identification is not intended to imply recommendation or endorsement by the

National Institute of Standards and Technology, nor is it intended to imply that the

entities, materials, or equipment are necessarily the best available for the purpose. 

National Institute of Standards and Technology Technical Note 1877

Natl. Inst. Stand. Technol. Tech. Note 1877, 46 pages (April 2015) 

CODEN: NTNOEF

This publication is available free of charge from: 

http://dx.doi.org/10.6028/NIST.TN.1877
Abstract

e combination of scanning electron microscopy (SEM) with X-ray microanalysis and image

processing provides a powerful ability to image and quantify microstructural features of construction materials. is document provides guidance for collecting backscaered electron

images and X-ray element maps of polished sections of Portland cement clinker. It furnishes

step-by-step procedures for processing image data to produce segmentation of an image eld

into its constituent mineral phases, and describes ways use the segmented image to measure

the abundance, surface area, and spatial distribution of phases in the image eld. Inherent

heterogeneity of the microstructure implies that dierent image elds of the same material

will contain dierent phase abundances; the document suggests ways to estimate and compensate for the heterogeneity length scale through statistical analysis of multiple image elds.

and analysis using a scanning electron microscope will be illustrated using an example from

the development of Standard Reference Materials for portland cement clinker.

i
Contents

List of Tables iii

List of Figures iv

1 Introduction 1

2 Specimen Preparation 2

2.1 Carbon Coating for Imaging . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

2.2 Sample Storage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

2.3 SEM Imaging . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

2.4 X-Ray Microanalysis and Imaging . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

3 Application to SRM 2686a 6

3.1 Image Processing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

3.2 Visualizing Phase Occurrence and Distribution by Color Merge . . . . . . . . . . 14

3.3 Image Processing with MultiSpec . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

3.4 Cleaning up the segmented image . . . . . . . . . . . . . . . . . . . . . . . . . . . 21

3.5 antitative Image Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23

3.5.1 Installing MicroChar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24

3.5.2 Using MicroChar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24

3.5.3 Two-point correlation functions . . . . . . . . . . . . . . . . . . . . . . . . 27

3.6 Larger Image Field . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

3.7 Multiple Image Statistics for SRM 2686a . . . . . . . . . . . . . . . . . . . . . . . . 32

Appendix A Step-by Step Summary for Clinker Analysis 34

References 38

ii
List of Tables

1 Clinker phases, compositions, densities, average atomic number, and backscattered electron coecient ranked according to relative SEM BE image brightness [17]. 5

2 Textural and chemical characteristics useful for distinguishing clinker phases. . . 6

3 Comparison of clinker phase mass fractions, expressed as a percentage of total

solids, estimated for SRM 2686a from the larger eld of view in Fig. 30 to the certied values [15]. Measurements are based on one image eld, so no uncertainty

estimate is given. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31

4 SRM 2686a phase fractions expressed as a mass percent by eld and a running

average. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33

iii
List of Figures

1 An early thin section micrograph (1917) of belite crystals bounded by dark gray

tricalcium aluminate domains (90×, crossed Nichols [2]). . . . . . . . . . . . . . . 2

2 Clinker microscopists relate textural characteristics to raw mix proportioning,

homogenization, and clinkering processes, as from Homannner [9]. . . . . . . . 3  ̈

3 Energy-dispersive X-ray microanalysis shows spectra from the primary clinker

phases exhibit characteristic spectra on a continuous background and standardbased quantitative analysis of chemical composition expressed as oxides. . . . . . 7

4 SEM backscaered electron image of SRM 2686a, showing a medium-grained,

heterogeneous silicate distribution and a coarse-grained matrix with phase identications. Field width is 256 μm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

5 X-ray microanalysis images of the spatial distribution of elements for the eld of

view shown in Figure 4. e image in the upper le is the BE image. Field width

is 256 μm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

6 Reducing noise associated with the continuous spectrum background improves

the delineation of aphthitalite from the raw Na XR image (le) and histogram, as

shown in the processed subtracted image (right). Field width is 256 μm. . . . . . . 10

7 Mg XR map associated with low concentrations in the silicates and matrix phases

may confound segmentations. Subtracting low-level signal associated with background and low concentrations of Mg presents a clearer rendering of periclase

location. Field width is 256 μm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11

8 resholded blur lter applied to the BE image removes minor scratch marks and

makes the gray levels for each phase more uniform as seen in the images and gray

level histograms. Field width is 256 μm. . . . . . . . . . . . . . . . . . . . . . . . . 12

9 resholded Blur algorithm applied to the Ca X-ray image improves the distinction between phases by comparing gray level histograms from the initial (A) and

ltered (B) images, which now shows distinct peaks for alite, belite and matrix.

Field width is 256 μm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

10 Composite image with the BE as red, Mg as green, Al as blue and S as yellow

illustrates how the addition of the X-ray data improves phase identication. Alite

is a brighter red than belite, the matrix phases are bluish (aluminate) and pink

(ferrite), periclase is green, and alkali sulfates are yellow. Field width is 256 μm. . 13

11 Loading the SEM BE image in MultiSpec interprets the image le details with

an enhancement option to trim a percentage of the high- and low-end signal, if

necessary. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

12 Linking the BE image to a set of XR images from the same eld of view is accomplished in MultiSpec by selecting each XR image, specifying “Link selected le(S)

to active image window” in the drop-down menu, and completed by selecting

Cancel. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

iv
13 Saving the linked BE and X-ray images using the Processor→Reformat→Rectify

provides a chance to subsample the set and repair any image registration issues

(rectication). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

14 Selecting the linked, rectied multispectral image set, indicated by le extension

.LAN will allow display of the image set and establishing training elds. . . . . . 18

15 ree-color (RGB) display by channel will provide a rendering of phase distribution for the next step of establishing the training elds. . . . . . . . . . . . . . . . 18

16 MultiSpec Image→Display allows three channels to be assigned to Red-GreenBlue color planes to facilitate identication of phase type and distribution. e

upper image is BE (red), Mg (green), and Al (blue) while the lower is BE (red), K

(green), and Na (blue). Field width is 256 μm. . . . . . . . . . . . . . . . . . . . . . 19

17 Training elds are dened by boxing regions typical for each phase (class) and

assigning a class name. Additional training elds are created by selecting a phase,

boxing a new region, and adding it to the class. . . . . . . . . . . . . . . . . . . . . 19

18 Classication is initiated by selection of a procedure, the use or exclusion of some

image channels, and writing the classication results to a disk le. . . . . . . . . . 20

19 MultiSpec classication table, or accuracy matrix, shows the classication set assignment accuracy, incorrectly assigned pixel counts, an overall accuracy assessment by the Kappa statistic, and the area fraction for each class [12]. Data shown

here are based on one image, so no uncertainty statement is possible. e “Reliability Accuracy” value for each phase refers to the percentage of training set

pixels that are correctly reclassied by the segmentation as the phase to which

they were originally assigned. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

20 Comparing the segmented image to the original BE image (Fig. 4) indicates a successful segmentation. Some incorrectly assigned pixels occur at grain boundaries

and within belite, which will be repaired in the nal processing step. Field width

is 256 μm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

21 Original SEM BE image (le) and the processed, segmented indexed image (right)

for comparison. Segmentation used the Fisher Linear Likelihood algorithm. Field

width is 256 μm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23

22 Installation of MicroChar on Mac OS X computers. . . . . . . . . . . . . . . . . . 24

23 MicroChar main window. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

24 Phase identication is accomplished in MicroChar by accepting the default settings or by selecting the phases from the series of pull down menus. . . . . . . . . 26

25 Files produced by MicroChar store the results of its calculations. . . . . . . . . . . 26

26 Contents of the pfc le for the segmented image of SRM 2686a in Fig. 20. Here

the column heading “Volume” corresponds to the phase area fractions and “Area”

to the perimeter fraction. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

27 (a) False color PNG image of SRM 2686 microstructure, and (b) the phase boundaries. Field width is 256 μm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

v
28 (a) Contents of the alite two-point correlation function data le for SRM 2686a

image in Fig. 20 with distance in the rst column and probability density in the

second column; (b) plot of the two-point correlation function S2 for alite and for

belite in the SRM 2686a image (right). . . . . . . . . . . . . . . . . . . . . . . . . . 28

29 Two-point correlation functions for periclase and arcanite in SRM 2686a. Note

the two dierent scales on the y-axes. . . . . . . . . . . . . . . . . . . . . . . . . . 29

30 Larger image eld of SRM 2686a. Magnication 300×, 1024×768 pixels. (a) False

color PNG image, and (b) the phase mass, volume, and area fractions. Field width

is 1.536 mm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

31 Two-point correlation functions calculated for alite and belite using the large image eld of SRM 2686a shown in Fig. 30. . . . . . . . . . . . . . . . . . . . . . . . 31

32 Running average plots for SEM/IA analysis of SRM2686a shows values beginning

to stabilize aer about ve image elds. . . . . . . . . . . . . . . . . . . . . . . . . 34

33 Box plots by phase for SEM analysis of SRM 2686a. Stars indicate certied values. 35

vi
1 INTRODUCTION

1 Introduction

Microstructural analysis of cementitious materials by microscopy dates back over 100 years when

Henri Le Chatelier [13] investigated portland cement clinker using what was, at the time, a new ˆ

polarized light microscopy method. He documented the optical properties of the main phases

and identied tricalcium silicate as the principal constituent of portland cement clinker. Subsequent work by A.E. Tornebohm, Bates and Klein [2] (see Fig. 1), and others provided a more

detailed description of the phases and textures of clinker and also developed the techniques for

light microscopy as it applies to portland cement clinker [11]. More recent resources for clinker

microscopy include the work of Campbell [6] and of Hofmanner [9]. antitative microscopy  ̈

methods began with Delesse [7], who recognized the relationship between the area fraction of a

mineral in a 2D section and its 3D volume fraction in the material. e rst quantitative procedures involved sketching on paper the observed microstructure, tracing the grain outlines onto

a metal foil, and nally cuing and weighing the mineral constituents from the foil traces. ese

quantitative methods evolved by subsampling the microstructure through linear traverse and

point sampling to achieve quantitative estimates of mineral abundance. ese laer methodologies, although more than a century old and quite time consuming, are still used today in many

disciplines for quantitative phase analysis.

e combination of the scanning electron microscope (SEM) with image analysis has in some

ways brought quantitative microscopy back to its roots. Traditional point count and linear traverse methods may still be applied to images generated by the SEM, but the SEM can be further

exploited to produce a set of high-resolution image elds from a single material. ese image sets

are suciently uniform that feature extraction and quantitative measurement by digital image

processing are possible. In this case, the mineral constituents are digitally traced and extracted

from the bulk microstructure by a set of operations termed image processing and image analysis,

which supplant the older weighing operations with semi-automated pixel counting to determine

phase abundance. e potential for automation in the data collection and processing promises

to provide for a much greater sampling of any microstructure and quantitative estimates of microstructural features that were not easily accomplished in the past.

is report describes imaging and analysis procedures for portland cement clinker using SEM

with backscaered electron (BE), X-ray microanalysis (XR), image processing and image analysis. e procedures are used to uniquely identify and separate a microstructure image into its

constituent phases, a process known as segmentation. Segmented images can be used to make

quantitative measurements such as phase abundance and surface area fractions.

It should be stated clearly at the outset that, for several reasons, microstructure imaging and

image analysis have always had—and likely will continue to have—a partially subjective aspect.

First, microstructure is not a well-dened material property. It is oen characterized by signi-

cant chemical and spatial heterogeneity that can be measured only in a statistical sense, and the

principles of uncertainty estimates for image analysis are largely undeveloped. Second, both an

electron beam’s interaction with the specimen and the digitized nature of a microstructure image

causes a level of coarse graining in which volume is homogenized into small elements that cannot

be resolved more nely than a single pixel. But neither chemical nor phase composition of the ac1
2 SPECIMEN PREPARATION

tual microstructure is constrained to the pixel scale; solid solutions and nely divided phases can

cause variations in contrast that, when averaged over the pixel volume, can make phase assignments uncertain for some pixels, especially at interphase boundaries and within the interstitial

phase regions. e fraction of such “borderline” pixels is usually small within any image, but no

single objective method is currently available for infallibly assigning their phase identity. Instead,

for these cases the user must appeal to several types of possibly inconsistent data about a pixel

or pixel grouping, such as its backscaered contrast, X-ray element signal, and location within

the microstructure, and arrive at and informed judgment based on experience. erefore, phrases

used in this report, such as “accurate phase assignment,” or “incorrectly assigned pixels” should

be qualied in terms of the best judgment of an experienced user.

Figure 1: An early thin section micrograph (1917) of belite crystals bounded by dark gray tricalcium aluminate domains (90×, crossed Nichols [2]).

2 Specimen Preparation

e sample taken for analysis should be as representative of the bulk material as possible. However, aside from ASTM C183 [1], no guide or standardized procedure is available for clinker or

cements. Sample spliers or the cone-and-quarter method may be used to reduce large-volume

samples to representative subsamples, and down to a volume suitable for making an epoxyembedded specimen mount. A medium-viscosity epoxy12 is suitable for the mount because it

1Epotek 301, Epoxy Technology, Inc., 14 Fortune Drive, Billerica, MA 01821

2Certain commercial entities, equipment, or materials are identied in this document in order to describe an

experimental procedure or concept adequately. Such identication is not intended to imply recommendation or

endorsement by the National Institute of Standards and Technology, nor is it intended to imply that the entities,

2
2 SPECIMEN PREPARATION

Figure 2: Clinker microscopists relate textural characteristics to raw mix proportioning, homogenization, and clinkering processes, as from Homannner [9].  ̈

can inltrate any connected pores within clinker. e clinker is embedded in the epoxy to preserve the microstructure, to support the pore and fracture network during cuing and polishing,

and to ll the voids to assist in pore identication

Aer curing the epoxy, the specimen is cut using a low-speed diamond saw with ethanol as

a cuing lubricant to expose a cross section of the material. Surface damage caused by sawing

is removed using successively ner grades of silicon carbide paper (400, 600, 1200 grit) leaving

a planar surface devoid of cuing and grinding fractures. is surface is then polished to remove surface scratches using diamond pastes of successively ner particle size (6 μm, 3 μm, 1 μm,

0.25 μm). Aer the nal polish, the specimen is cleaned and carbon coated in preparation for the

SEM. A more detailed procedure for specimen preparation may be found in Ref. [16, 18].

2.1 Carbon Coating for Imaging

A thin coating of carbon is necessary to provide an SEM specimen with a conductive coating.

ere are circumstances in which this is unnecessary (e.g., if using a variable pressure microscope), but coating will provide the conductive path to ground that will reduce the potential for

charging-related problems like beam deection and associated image dri, and reduction in incident electron beam accelerating voltage. Carbon evaporators made for electron microscope

specimen preparation are available for this purpose. Metal coatings are not recommended bematerials, or equipment are necessarily the best available for the purpose.

3
2 SPECIMEN PREPARATION

cause they reduce imaging signal and contrast, and they also interfere with X-ray imaging.

2.2 Sample Storage

Clinker specimens are moisture-sensitive so storage in a desiccator when not in use is necessary.

Some hydration of the more reactive constituents (free lime) will inevitably occur, necessitating

sample examination as soon as practical aer preparation is complete. Hydration of free lime

may also reect sample storage-related prehydration and may be a useful observation relevant in

understanding anomalous cement behavior. e prehydration features may be distinguished by

their surface being level with the specimen whereas hydration aer specimen preparation will

stand in high relief due to the substantial volume increase.

2.3 SEM Imaging

High-quality images for quantitative analysis can be obtained with the following SEM operating

conditions:

1. 12 kV accelerating voltage,

2. 2 nA to 4 nA probe current, adjusted to keep a relatively high X-ray count rate with minimal

dead time,

3. a scan rate of 5 min/frame for an image size of 1024 pixels × 768 pixels to minimize BE

noise,

4. an accumulation of about 15 frames to improve the denition of the X-ray images

5. A magnication of 300×, which provides a spatial resolution of approximately 1.0 μm/pixel,

which is typically on the order of the width of the beam’s interaction with the specimen.

ese seings are a useful starting point to obtain useful images, and are judged to represent a

balance between fast collection times, a match in resolution of BE images with X-ray maps, and

instrument stability.

e local brightness in a BE compositional image is proportional to the average atomic number, Z, of the phase in that location, and the XR image reects element spatial distribution over

the same eld of view. e backscaer coecient, η, is a measure of the backscaered electron

fraction and, following Goldstein et al. [8], is estimated using the mass fractions and η values for

each constituent. Table 1 lists phases found in clinker and cement in descending order of their

backscaered coecient and gray intensity. Caution is necessary here because the brightness

ranking may change depending upon chemistry variability due to impurities; the exact values

are not as important as the general ranking of the phases in the table. e identication of the

individual phases also includes supplemental information such as domain shape (e.g., angular or

rounded), position within the microstructure (e.g., framework grain, matrix, dispersed phase),

and bulk chemistry as described in Table 2.

4
2 SPECIMEN PREPARATION

Table 1: Clinker phases, compositions, densities, average atomic number, and backscaered electron coecient ranked according to relative SEM BE image brightness [17].

Phase Composition Notation Density (Mg m−3) Z η

Ferrite Ca2(AlxFe1-x)

2O5 C4AF 3.77 16.65 0.1860

Free lime CaO C 3.32 16.58 0.1882

Alite Ca3SiO5 C3S 3.13 to 3.22 15.06 0.1716

Belite Ca2SiO4 C2S 3.28 to 3.31 14.56 0.1662

Arcanite K2SO4 KS 2.67 14.41 0.1652

Aluminate-ca Ca3Al2O6 C3A 3.04 14.34 0.1639

Aluminate-ob Ca3Al2O6 C3A 2.56 13.87 0.1588

Aphthitalite NaCaAl3O9 NCA 2.7 13.69 0.1577

enardite Na2SO4 NS 2.66 10.77 0.1249

Periclase MgO M 3.58 10.41 0.1213

a

cubic polymorph

b orthorhombic polymorph

e BE contrast between alite and belite is relatively strong, but that between belite and

cubic tricalcium aluminate is generally too low to rely on contrast to distinguish them. For such

phase pairs that have similar BE coecients, or for phases that have BE coecients too low to be

resolved in the BE image, XR imaging is used to distinguish them based on their composition. For

example, periclase and some alkali sulfates oen appear black, like the epoxy-lled voids, because

of the brightness and contrast seings necessary to make the other phases distinct. Combining

the XR and BE images facilitates their discrimination much beer than relying on either image

alone. Information redundancies also occur that should be reduced or eliminated to improve the

phase distinction. For example, the Ca and Si XR image, or a ratio of the two, may be considered

for distinguishing the calcium silicate phases. e BE image, however, provides a much beer

distinction between the calcium silicates and is preferred because of it is less “grainy” than the

XR image. Further details and examples of how SEM images may be combined to render these

phases distinct will be shown in a later section.

2.4 X-Ray Microanalysis and Imaging

X-Rays are produced as a result of the interaction of the electron beam with the specimen, and

they are typically displayed as a set of peaks with a continuous background (Fig. 3). e energydispersive X-ray analyzer (EDX) ranks and counts X-rays based upon their energy. e spectrum

from alite shown in Fig. 3 has the characteristic peaks for calcium, silicon, aluminum, magnesium,

sodium, oxygen, and carbon (the conductive coating), and a background. e concentrations of

the elements can be quantied if suitable standards are available, and are typically expressed as

mass percent of the oxide.

5
3 APPLICATION TO SRM 2686A

Table 2: Textural and chemical characteristics useful for distinguishing clinker phases.

Phase Identifying Characteristics

Free Lime Bright BE, rounded crystals, sometimes clustered,

strong Ca and no Fe

Ferrite Bright BE, dendritic to prismatic matrix phase, high Ca,

medium Fe and Al

Alite Medium-high BE, euhedral to anhedral crystals, framework grain, strong Ca and Si, minor Mg and Al.

Belite Medium BE, rounded crystals, framework grain, strong

Ca and Si, minor Mg, Al, and weak S. Sometimes clustered.

Aluminate Medium BE, ne to lath-like crystals, matrix phase between ferrite crystals, Ca, Al, Mg, K, Fe and sometimes

Na.

Arcanite Medium BE, crystals dispersed, typically along boundaries of matrix and framework grains, high K, S

Aphthitalite Low BE, crystals dispersed, typically along boundaries

of matrix and framework grains, high K, Na, S

enardite Low BE, crystals dispersed, typically along boundaries

of matrix and framework grains, high Na, S

Periclase Low BE, equant to dendritic grains, clustered to dispersed

3 Application to SRM 2686a

is section demonstrates the use of the foregoing procedures for analyzing the NIST Standard

Reference Material (SRM) clinker 2686a [15]. A representative BE image of this clinker is shown

in Fig. 4, and all subsequent images in this section will use the same image eld for illustration

purposes. e phase domains in this clinker have an intermediate size compared to the other

SRM clinkers. e spatial distribution of phases is heterogeneous, and lile free lime is detected.

Alite occurs as subhedral to anhedral crystals approximately 30 μm in size. Belite occurs in large

clusters with an approximate crystal size of 20 μm. e matrix is dierentiated with a mediumto ne-grained lath-like ferrite, and ne-grained aluminate lling the inter-lath voids. Equant

periclase crystals up to 15 μm are common throughout the microstructure and the alkali sulfate

phases aphthitalite and arcanite are distributed throughout the microstructure. According to the

Hofmanner diagnostic chart in Fig. 2, the belite clustering represents heterogeneity or over-sized  ̈

silica grains in the raw mix. e ne-grained matrix indicates rapid cooling from 1300 ◦C, and

the slightly elevated belite content may reect a slightly low CaO content.

An annotated BE image of this clinker is shown in Fig. 4, illustrating the distinct contrasts be6
3 APPLICATION TO SRM 2686A

Figure 3: Energy-dispersive X-ray microanalysis shows spectra from the primary clinker phases

exhibit characteristic spectra on a continuous background and standard-based quantitative analysis of chemical composition expressed as oxides.

tween most of the phases. Reference to Tables 1 and 2 should allow one to identify the constituent

phases. While the contrast between belite and aluminate is not distinguishable, as mentioned in

the previous section, they are still texturally distinct and can be distinguished on this basis even

within the BE image, the belite being a large-grained framework grain and the aluminate being

a ne-grained matrix. Periclase and alkali sulfates appear dark, like the voids, although they

would be easier to detect with a brighter image seing. e ancillary XR images for the major

elements will be used to identify the locations of these phases as well as isolate them in the image

processing stage.

is image set was collected at 1200× magnication using a 1024 pixel × 768 pixel image,

producing a spatial resolution of 2 μm/pixel. is is in excess of that used for routine imaging

and was done here to beer illustrate the phases. e set of XR images associated with this

backscaered electron image is shown in Fig. 5. ese images were collected simultaneously

7
3 APPLICATION TO SRM 2686A

Figure 4: SEM backscaered electron image of SRM 2686a, showing a medium-grained, heterogeneous silicate distribution and a coarse-grained matrix with phase identications. Field width

is 256 μm.

so they are registered, with the same pixels on each image representing the same region. e

brightness is proportional to X-ray counts, and roughly proportional to element concentration.

While the gure is displayed as a gray-scale image, pseudo-coloring can sometimes improve the

visualization of the element distribution. e images in Fig. 5 have been modied to remove

contributions from the continuous spectrum background using a simple subtraction routine that

will be described more fully in the next section.

e energy-dispersive X-Ray signal consists of the characteristic spectrum for each element

and a continuous spectrum, oen called the background (Fig. 3). is background will contribute

image noise that may be signicant, especially for 8-bit image formats that can resolve only 256

gray levels. A 16-bit image format, with 65 536 gray levels, can oer a beer signal-to-noise ratio.

In other cases, low concentrations of an element are expressed as a very weak and variable signal,

given the relatively short X-ray imaging dwell time. While these low-concentration elements are

helpful for quantitative microanalysis of clinker mineral compositions, they are not as useful for

phase segmentation, and are another example of redundancy.

e dierences between the silicates, alite and belite, are shown in the BE image (Fig. 4). e

dierences in concentration of Ca and Si in the associated XR images (Fig. 5) reects the dierent

Ca and Si concentrations of alite and belite.

e aluminum and iron images in Fig. 5 highlight the interstitial (matrix) phases, aluminate

and ferrite3. ese elements also exist in low concentrations in some of the other phases, and

will therefore have low counts throughout other portions of their XR images. e iron image is

another example of an image redundancy because it highlights mainly the ferrite and, to a lesser

3For brevity, tricalcium aluminate is routinely designated simply as “aluminate,” distinguished as cubic or orthorhombic where appropriate, and tetracalcium aluminoferrite is designated as “ferrite”.

8
3 APPLICATION TO SRM 2686A

Figure 5: X-ray microanalysis images of the spatial distribution of elements for the eld of view

shown in Figure 4. e image in the upper le is the BE image. Field width is 256 μm.

degree, the aluminate phase. ese phases are beer distinguished in the BE image (Fig. 4) and,

for aluminate, the BE and the Al image, so the iron image is not used for the image processing.

e aluminum XR image, however, is important because it can be used together with the BE

image to distinguish belite from aluminate based on the high concentration of Al the laer phase.

Further features that stand out upon a cursory scan of Fig. 5 are periclase crystals dispersed

throughout the image (regions of high Mg), arcanite (areas with strong signals from both S and

K), and aphthitalite (regions with S, K, and Na).

For a more systematic exploration of the phase constituents and their distribution in the

clinker, the rst step is to examine the individual and then combinations of the image set. is is

accomplished using ImageJ [10], or a similar image viewer, and merging dierent image combinations. Two combinations can generally serve to highlight all the phases. ese combinations

may also be used to see if a particular image improves or confounds a phase distinction, allowing

the user to decide on adding that image to the set that will be used for the image processing used

9
3 APPLICATION TO SRM 2686A

to uniquely identify the clinker phases. ImageJ can merge up to seven images while most other

codes display three, as a red, a green, and a blue color plane (RGB). Displaying the BE image, Mg,

and Al as a RGB image allows distinction of alite, belite, ferrite, aluminate, periclase, and voids.

Displaying the BE image with K and Na provides distinction of the alkali sulfates. is process

will be discussed in more detail in a later section.

Figure 6: Reducing noise associated with the continuous spectrum background improves the delineation of aphthitalite from the raw Na XR image (le) and histogram, as shown in the processed

subtracted image (right). Field width is 256 μm.

3.1 Image Processing

e next four sections provide detailed guidance for processing the experimentally obtained images and transforming them to a single phase segemented image suitable for quantitative stereological analysis. e Appendix provides a condensed set of instructions distilled from this discussion.

As an example applied to SRM 2686a, Fig. 6 shows the Na XR map associated with one eld

of view. In this image, the continuous background and inherently low counts associated with

aphthitalite are problematic for segmentation decisions, because they reduce the signal-to-noise

ratio relative to that in an image with a greater element concentration. e background count

level can be determined from the histogram (see the insets of Fig. 6), or by reading the counts

while scrolling the pointer across the background; the Na image in Fig. 6 registers about 10

counts. Once the background count level is determined, the background can be subtracted. To

do so within ImageJ, choose the Process→Math→Subtract operation. is operation opens a preview option window where the subtraction threshold value can be entered (10, in this case). e

resulting image, shown as the right image in Fig. 6, more clearly distinguishes the aphthitalite

domains.

e Mg XR map for the same eld of SRM 2686a shows the spatial distribution of Mg in

the silicates, the matrix, and, in particular, periclase (Fig. 7). e silicates and matrix phases

10
3 APPLICATION TO SRM 2686A

are more clearly distinguished using a combination of the BE image and Al XR image, so the Mg

distribution information is not needed for their segmentation. e low concentration distribution

of Mg becomes apparent in the Mg image by expanding the brightness and contrast (Fig. 7, le).

Just as for the Na distribution, removal of the low-concentration and background counts from

the Mg map produces a clearer delineation of the Mg distribution for periclase. In this case,

scrolling the pointer over the low-Mg phases, such as alite, in the upper image of Fig. 7 indicates

the maximum background level to be about 20 counts, so 20 can be used as the threshold value

in ImageJ’s Process→Math→Subtract operation to produce the right image in Fig. 7.

ImageJ has other operations in itsImage→Noise menu for removing additional image artifacts

(e.g., speckling) that can complicate phase segmentation decisions. Examples include Despeckle

and Remove Outliers.

Figure 7: Mg XR map associated with low concentrations in the silicates and matrix phases

may confound segmentations. Subtracting low-level signal associated with background and low

concentrations of Mg presents a clearer rendering of periclase location. Field width is 256 μm.

Smoothing lters, several of which are available in ImageJ can improve the distinction between phases, but selecting one that preserves edge detail is important to preserve the texture of

the clinker microstructure. A useful ImageJ lter is called resholded Blur [20], which features

an interactive, selective mean lter and displays an image preview to assess the eect of the lter’s seings prior to commiing it. Used carefully, this lter will enhance distinction between

phases by reducing noise and retaining edge denition. It will later be important for cleaning up

incorrectly assigned pixels aer the segmentation is complete. Figs. 8 and 9 show examples of

the application of the resholded Blur algorithm to the BE image and the Ca XR image, respectively. e BE image exhibits the lowest noise of the set, although there can be some variability

in gray intensity and minor scratches that may appear slightly darker than the surroundings.

ese minor imperfections can be removed by careful selection of smoothing seings. In Fig. 8,

image set (A) is the original BE image, including ne scratches; a magnied portion is shown at

the boom of the gure. e smoothed image (B) in Fig. 8 shows a much more uniform gray for

11
3 APPLICATION TO SRM 2686A

Figure 8: resholded blur lter applied to the BE image removes minor scratch marks and makes

the gray levels for each phase more uniform as seen in the images and gray level histograms. Field

width is 256 μm.

12
3 APPLICATION TO SRM 2686A

Figure 9: resholded Blur algorithm applied to the Ca X-ray image improves the distinction

between phases by comparing gray level histograms from the initial (A) and ltered (B) images,

which now shows distinct peaks for alite, belite and matrix. Field width is 256 μm.

Figure 10: Composite image with the BE as red, Mg as green, Al as blue and S as yellow illustrates

how the addition of the X-ray data improves phase identication. Alite is a brighter red than

belite, the matrix phases are bluish (aluminate) and pink (ferrite), periclase is green, and alkali

sulfates are yellow. Field width is 256 μm.

13
3 APPLICATION TO SRM 2686A

each phase, with the ner scratches removed while retaining the grain boundaries, and providing an improved distinction between the silicate histogram peaks. e smoothing seings that

produced Fig. 8 (B) included a kernel radius of 2 pixels, an intensity threshold of 10, a soness of

0, and strength of 2. ese values were selected based on the description of their functions [10]

and the eects on smoothing the image without loss of phase denition.

3.2 Visualizing Phase Occurrence and Distribution by Color Merge

e Color Merge operation in ImageJ is a useful tool to visualize and identify the phases and their

spatial distribution. is operation combines, or stacks, BE and XR images while assigning unique

colors to each image. e colors blend and, depending upon the image combination chosen,

highlight the dierent phases. is is best done with a subset of the XR images or with XR

images and the BE image. e important images to combine include the BE, Mg, Al, K, Na, and

S. Other images either do not further aid in phase discrimination (e.g., Ca, Si) or have such low

concentrations that their noise confounds the phase segmentation (e.g., Ti, Mn).

e color merge is accomplished in ImageJ using the Image→Color→Merge Channels operation. Depending on the version of ImageJ, up to seven images can be merged at once. Fig. 10

shows the merging of the BE (red), Mg (green), Al (blue), and S (yellow). Arcanite and aphthitalite could be distinguished from each other by adding the K and Na images to the merge. e

aluminum XR image makes the distinction between belite and aluminate clear, even though their

contrast in the BE image is nearly the same. e Mg image reveals the locations of periclase,

while the S image highlights the alkali sulfates. Alkali sulfate phases in Fig. 10 are evident at

grain boundaries and within fractures of the belite crystals, which are typical locations for these

phases. Simpler image combinations than those in Fig. 10, such as BE-Mg-Al and BE-Na-K, can

be useful for identifying phases, especially if the image display code can show only three color

planes.

With the visualization of the phases and their distribution as in Fig. 10, the next step in quantitative image analysis is to train the segmentation routine to recognize the individual phases

and then generate an indexed image in which each phase is uniquely identied and coded. is

process involves classication of each pixel as the phase to which it most likely belongs, based

upon a user-created training set and multispectral image processing.

Based on the exploratory analysis using ImageJ and the color merging operations illustrated in

Fig. 10, a minimum set of images selected to discriminate each phase are included in the segmentation operations. is set typically includes BE, Al, Mg, K, Na, and S. As mentioned previously,

other element images should not be included because they do not add any new information to

the image set, and may even confound the segmentation. More details on the data collection, the

multi-spectral image processing, and analysis are given elsewhere [12, 19].

3.3 Image Processing with MultiSpec

MultiSpec, a soware package developed at Purdue University, is a processing system for interactively analyzing multispectral image data [14]. It is used here to generate the nal segmented

14
3 APPLICATION TO SRM 2686A

image from a small subset of images described in the previous section. A condensed set of instructions for using ImageJ and MultiSpec can be found in the Appendix. is section provides

a summary of the steps required to compile a data set, establish the training elds, classify, and

generate the segmented, indexed image for nal analysis.

MultiSpec has tools to help establish an optimal data set for analysis, but the relative simplicity

of the phase chemistry and image sizes make this a relatively straightforward task that can be

eciently completed with ImageJ’s Color Merge operation described earlier. e less useful Xray images can optionally be included, but they should be placed as the last of the images so that

they can be easily excluded from the classication if necessary.

MultiSpec identies and reports the image size, bit depth, and number of planes (RGB color

or single gray). In most cases, the default parameters suggested by the soware can be accepted.

In some cases when processing is performed, an option is available to write a text le with the

details and a processed image le; such cases will be noted in the following description.

Processing clinker or cement microstructure images with MultiSpec involves seven basic

steps:

1. Load the BE image

2. Link the BE image to the XR image set

3. Rectify and save the linked les

4. Identify phases (class assignment)

5. Set up training elds

6. Perform phase segmentation and check its quality

7. Save the MultiSpec project

Step 1: Load the BE image. Fig. 11 shows how to load the BE image as a multispectral type

using the File→Open Image menu item and selecting the multispectral option. Under the “Multispectral Display Specications” section of the dialog, the defaults can all be accepted with the

exception of the “Enhancement” section; here one will usually opt to leave the image unchanged

since background noise was removed using the ImageJ Process→Math→Subtract function, but

the image can also be clipped on the low and high end, or user selected clipping based upon a

percentage of the gray level.

Step 2: Link BE and XR image sets. XR images that previously have been deemed useful in

making a segmentation must be linked to the loaded BE image using the File→Open Image menu

item, as shown in Fig. 12. Select the (already-modied and optimized) XR images by clicking on

them individually, and check the box or drop-down menu for each new le to indicate “Link to

active image window”. Deciding on a standardized image sequence is helpful to later recall the

15
3 APPLICATION TO SRM 2686A

Figure 11: Loading the SEM BE image in MultiSpec interprets the image le details with an

enhancement option to trim a percentage of the high- and low-end signal, if necessary.

exact image source in the virtual linked les (BE, Al, Mg, K, Na, S). ese linked images can be

named individually later if necessary. Any image in the sequence can also be excluded if it is later

determined to hinder the classication of phases. As already stated, image exclusion is easier if

the potentially problematic XR images are placed at the end of the linked image sequence. Click

on the Cancel buon to complete the linking.

Step 3: Rectify and save linked les. Use the Processor→Reformat→Rectify Image menu

item, shown in Fig. 13, to rectify the linked les. Options are provided to subsample, translate,

and rotate the images, although this is rarely needed because the SEM image set should already

be suitably registered. However, if beam charging caused the BE image to be misaligned relative

to the XR image set, these options can help resolve the discrepancy. e default values of the

other parameters can be accepted. e Processor→Reformat→Add Channel Descriptions option

allows one to add a descriptive name for each channel, which can be useful for keeping track of

image sets that use nonstandard image sequences. Provide a le name with a .LAN extension

and, before saving, check that the path is correct because it may not default to the original folder.

A good practice is to use the original le name and keep everything together together in a folder

of the same name. Note: Using a .LAN extension on the le name is important; without the .LAN

extension MultiSpec will not automatically recognize the le as valid for any subsequent loading

in the next steps.

Step 4: Assign classes. Reload the .LAN image set using the File→Open Image menu item

(Fig. 14). If it is not highlighted as a readable le, designate the le type as “Multispectral Image

Files” section of the dialog window.

Once the image set is opened, a window labeled “Multispectral Display Specications” will

appear (Fig. 15) where channel colors can be assigned (R G B), image enhancements such as

clipping or transformation can be made, and image sets can be subsampled by lines and rows.

is display tool functions much like the Color Merge operation in ImageJ, but it allows only

16
3 APPLICATION TO SRM 2686A

MultiSpecLink

Figure 12: Linking the BE image to a set of XR images from the same eld of view is accomplished

in MultiSpec by selecting each XR image, specifying “Link selected le(S) to active image window”

in the drop-down menu, and completed by selecting Cancel.

Figure 13: Saving the linked BE and X-ray images using the Processor→Reformat→Rectify provides a chance to subsample the set and repair any image registration issues (rectication).

17
3 APPLICATION TO SRM 2686A

Figure 14: Selecting the linked, rectied multispectral image set, indicated by le extension .LAN

will allow display of the image set and establishing training elds.

three color planes. e three-color limitation may seem restrictive, but it reduces the likelihood

of confusion in interpreting the displayed images.

Figure 15: ree-color (RGB) display by channel will provide a rendering of phase distribution

for the next step of establishing the training elds.

When assigning channel colors in the “Multispectral Display Specications” window, choose

channel/color combinations that will provide a clear rendering of the phase distribution. e most

useful for clinker is the combination of BE (Red), Mg (Green), Al (Blue), which can clearly reveal

alite, belite, aluminate, ferrite, periclase, lime, and void, as shown in the le image of Fig. 16. To

complement such a color scheme, alkali sulfates can be readily identied with the combination

of BE (Red), K or Na (Green) and S (Blue), as shown in the right image of Fig. 16. Oen the Na

image is redundant with K, but in this example there are dierences representing two distinct

alkali sulfate phases. is second color combination is useful because alkali sulfates are dicult

to perceive from the BE image alone. e BE coecient for arcanite is close to that of belite, and

aphthitalite’s BE coecient is low enough that it may be mistaken for void space without the

additional X-ray information from the Na image.

18
3 APPLICATION TO SRM 2686A

Figure 16: MultiSpec Image→Display allows three channels to be assigned to Red-Green-Blue

color planes to facilitate identication of phase type and distribution. e upper image is BE

(red), Mg (green), and Al (blue) while the lower is BE (red), K (green), and Na (blue). Field width

is 256 μm.

Figure 17: Training elds are dened by boxing regions typical for each phase (class) and assigning a class name. Additional training elds are created by selecting a phase, boxing a new region,

and adding it to the class.

19
3 APPLICATION TO SRM 2686A

Figure 18: Classication is initiated by selection of a procedure, the use or exclusion of some

image channels, and writing the classication results to a disk le.

Step 5: Establish training elds. When all the phases are identied as in Fig. 16, training

elds must be assigned to condition MultiSpec for the image segmentation step.4 A powerful

feature of MultiSpec is its ability to establish training classes by interactively highlighting regions of the image that are associated with specic phases. To set up the training classes, use

Processor→Statistics, and accept the defaults in the “Set Project Options” dialog window. With

the recommended image combinations to highlight the phases from the previous step, regions

for each phase (or class) can be designated as shown by the white outlined boxes in Fig. 17.

e “New Project” window shows a Class Designation at the top, with ‘New’ as the default.

A list of classes can be initiated here by boxing in a region in the image (mouse click and drag)

across each phase individually. For each phase, the class is then added by clicking the “Add to List”

buon and providing a name (e.g., alite, belite). e denition of class characteristics improved

as more pixels are added to the training sets. erefore, the largest possible box ordinarily should

be selected for each phase. If the box accidently covers an unwanted area, recreate the box before

selecting “Add to List”. An option is available for selecting polygonal regions instead of boxes,

by clicking on line segments until the polygon is closed, but the box approach works quickly

and is usually sucient for dening classes. If a eld appears to be incorrectly assigned aer

being added, it can be deleted before continuing by selecting class and then eld in the project

window to get a list of individual elds. Select the eld you wish to delete and from Edit in the

top command bar, select “Cut Field” to delete.

Classes will have an associated index value that is assigned based on the order in which

the elds were added to the list in the Project window. A given class also optionally can be

rened further into subclasses. For example, a generic alkali sulfate class could contain subclasses

for arcanite and apthitalite, although a generic class is rarely helpful in phase segmentation.

e important point is to use a consistent sequence to facilitate phase identication and further

quantitative analysis of the indexed image. e following convention is highly recommended:

4MultiSpec has a built-in automated segmentation utility, which ordinarily would eliminate the need for training classes described here. However, we have found that MultiSpec’s automated segementation utility yields less

thorough, less accurate segmentations than can be obtained by establishing training classes.

20
3 APPLICATION TO SRM 2686A

1 = alite 2 = belite 3 = aluminate 4 = ferrite

5 = periclase 6 = arcanite 7 = aphthitalite 8 = void

e index assignments should be recorded for later use. If phases are absent, the index values

shi down accordingly. For example, if a clinker has no identiable periclase, then the indices

for arcanite, apthitalite, and void would all decrease by one.

Step 6: Segment phases. Use the Processor→Classify command to begin the image segmentation algorithm. Seven classication procedures are available in MultiSpec, and it is relatively

easy to test them all to determine which one yields the best segmentation. e classication

window is shown in Fig. 18. In general, both the “Minimum Euclidean Distance” and the “Fisher

Linear Likelihood” algorithm work well for segmenting clinker or cement images. e results are

saved by checking the “Disk File” radio buon under the “Write Classication Results” section

of the window. Prior to saving, the le path should be checked and the .GIS extension should be

retained for the image.

e relative merit of a classication should be based on how well the classier is able to reassign the training pixels to their designated class (Fig. 19) and, most importantly, on the agreement between the classied image (.GIS) and the operator’s interpretation of the BE image. e

.GIS image is read into MultiSpec as a “Multispectral” type image, and the class assignments are

shown with distinct colors, as in Fig. 20. e color of any phase can be changed by double clicking

it in the class list and seing its color with the color picker. If the segmentation is acceptable, this

image can be saved as a TIFF le (File→Save Image to TIFF). Note: e le must be saved in TIFF

format because that is the only format that retains the index values that are needed for the nal

analysis steps.

Step 7: Save the project. Select the File→Save Project As menu item to which will retain all

the eorts and results of the segmentation process. Even if a segmentation aempt is deemed

unsuccessful, the project can be reloaded later and modied to add training elds, delete elds,

run alternate classication algorithms, or even exclude specic images to see if the classication

can be improved.

3.4 Cleaning up the segmented image

Even aer the image segmentation in MultiSpec, some incorrectly assigned pixels are commonly

found. Especially along grain boundaries or particle edges, contrast artifacts in the original SEM

BE image can cause pixels to be incorrectly classied. ese eects are ideally minimized beforehand by careful adjustment of brightness and contrast seings during acquisition of the SEM

BE image. In this example (Figs. 20 and 21), the edge eects created a problem for the Maximum

Likelihood classier based upon inspection of the segmented image. e Fisher Linear Likelihood

classier generated a more accurate grouping of pixels but still had some incorrectly assigned pixels along the edges and within the belite grains. If not too severe, incorrectly assigned pixels may

be removed aer segmentation by using some of the smoothing lters in ImageJ. reshold Blur

21
3 APPLICATION TO SRM 2686A

Figure 19: MultiSpec classication table, or accuracy matrix, shows the classication set assignment accuracy, incorrectly assigned pixel counts, an overall accuracy assessment by the Kappa

statistic, and the area fraction for each class [12]. Data shown here are based on one image, so no

uncertainty statement is possible. e “Reliability Accuracy” value for each phase refers to the

percentage of training set pixels that are correctly reclassied by the segmentation as the phase

to which they were originally assigned.

Figure 20: Comparing the segmented image to the original BE image (Fig. 4) indicates a successful segmentation. Some incorrectly assigned pixels occur at grain boundaries and within belite,

which will be repaired in the nal processing step. Field width is 256 μm.

22
3 APPLICATION TO SRM 2686A

Figure 21: Original SEM BE image (le) and the processed, segmented indexed image (right) for

comparison. Segmentation used the Fisher Linear Likelihood algorithm. Field width is 256 μm.

works well for minor edge eects by ignoring pixels that deviate from the mean and by allowing

one to specify the kernel radius (start with 2 or 3), threshold (1, 2, or larger because this is an indexed image and the values dier by increments of 1), soness (aects the edge blur) and strength

(number of times to apply the lter). Aer the edge eects have been reduced or eliminated in

the .GIS le, save this and a TIFF and also write the ASCII Indexed le: File→Save As→Text Image

that will be used in the image analysis. If needed, the .GIS le in Multispec or the *MS.TIF le in

ImageJ can be recalled to see which index value corresponds to each phase.

e completed segmentation with edge eects removed and original BE image are presented

in Fig. 21 to illustrate the nal segmentation results. e image is color-coded and indexed with

alite (blue, 1), belite (brown, 2), aluminate (white, 3), ferrite (red, 4), periclase (pink, 5), arcanite

(blue-green, 6), aphthitalite (yellow, 7), and porosity (black, 8).

3.5 antitative Image Analysis

e nal result of phase segmentation is a 2D color image that identies each pixel with an

integer index associated with its assigned phase. antitative analysis can be performed on such

an image to calculate, for example

1. the area fraction of each phase (on a total solids basis)

2. the perimeter fraction of each phase (on a total solid boundary length basis)

3. the spatial distribution of each phase, quantied by correlation functions

4. the possible percolation, or connectivity, of each phase

Many, if not all, of these operations can be performed by built-in functions, such as generating

a histogram and list from the indexed image or augmenting ImageJ with Java plugins such as the

correlation function. However, most of these plugins have been developed and optimized for

23
3 APPLICATION TO SRM 2686A

image processing in other elds, especially biology, and are not well-documented or validated

for use with clinker or cement microstructure images. As an alternative, NIST has developed a

simple, fast computer program specically for performing quantitative image analysis on clinker

or cement microstructures. e program, called MicroChar, is available for Windows or Mac

OS X computers, performs all these calculations with minimal intervention by the user. Use of

MicroChar is described in the following section.

3.5.1 Installing MicroChar

MicroChar can be downloaded as disk image (on Mac OS X) or installer (on Windows).

Installation on Mac OS X. Double click the disk image to mount it. is will open a window

as shown in Fig. 22. e application can be installed by dragging its icon to the Applications

folder.

Figure 22: Installation of MicroChar on Mac OS X computers.

Installation on Windows. Double click the installer. If the user accepts the disclaimer, then

choose a folder in which to install the application folder. Aer installation, a shortcut icon will

be placed on the Desktop and in the Start Menu.

3.5.2 Using MicroChar

Launch MicroChar by double-clicking its icon. A window will open as shown in Fig. 23. Mousing

over many of the items in the window will display a tooltip. For more context-specic information, click on the blue question mark icon in the toolbar and then click on the window item about

which you want more information.

To begin a session, click on the blue folder icon in the upper le toolbar (or type Ctrl-O, to

open the indexed image resulting from phase segmentation. e application works equally well

on indexed images either of clinker or of powdered cement embedded in epoxy, but it anlayzes

the void space and perimeter identication dierently for the two kinds of images. erefore,

aer le selection, check the “Clinker image” checkbox only if the image is a clinker nodule;

otherwise leave it unchecked. Specify the image resolution in terms of the width of a single pixel

in micrometer (μm) units. e pixel dimension should be known at the time that the SEM image

24
3 APPLICATION TO SRM 2686A

Figure 23: MicroChar main window.

is acquired, and is necessary to associate a length scale with the correlation functions that will

be calculated if the corresponding checkbox is activated.

e remaining items in the window, which are grayed out by default, are related to data

needed to create a cement data package for the Virtual Cement and Concrete Testing Laboratory.

ese items, which can be activated only for cement images by clicking the “VCCTL Package”

checkbox, are beyond the scope of this document. Readers who are interested in using VCCTL

soware should consult the VCCTL soware documentation [21].

Calculations begin by clicking the Execute buon at the boom of the window. e program

reads the image le and records all the phase index numbers. When these have been collected,

the user is prompted to associate each index with a known phase, as shown in Fig. 24. A default

collection is oered, but each phase can be manually assigned by selecting the desired phase

from the pulldown menu next to each phase id found in the image. Click the OK buon when

the desired selections have been made.

MicroChar proceeds automatically to perform the following sequence of operations:

1. e area of each phase is calculated and normalized by the total solid area to yield its area

fraction on a total solids basis.

2. e mass fraction of each phase, on a total solids basis, is calculated from the area fractions

and knowledge of the density of each phase.

3. e perimeter of each phase is calculated. For clinker images, the perimeter is dened

as the boundary of the given phase with any other phase, including porosity; for cement

images, the perimeter is dened as the boundary of the given phase only with interparticle

porosity, because such boundaries determine the reactive surface area with water.

4. False color PNG images are created of (a) the microstructure and (b) the phase perimeters.

25
3 APPLICATION TO SRM 2686A

Figure 24: Phase identication is accomplished in MicroChar by accepting the default seings

or by selecting the phases from the series of pull down menus.

5. Two-point correlation functions, also known as autocorrelation functions in the eld of

signal processing, are calculated for each phase (only if the checkbox was activated).

6. e program halts and prompts you to save your work if you wish, which you may do by

clicking on the disk icon in the tool bar or by typing Ctrl-S.

As it proceeds through this sequence, MicroChar writes a collection of les to catalog the

results of its calculations. ese les, wrien to the same folder as the original indexed image, all

have the same root name as the indexed image but dierent le extensions, as shown in Fig. 25.

Figure 25: Files produced by MicroChar store the results of its calculations.

ree types of les are produced for clinker images:

1. a le ending with ’pfc’5, which contain the mass, area, and perimeter fractions of each

phase in the microstructure on a total solids basis, as shown in Fig. 25;

2. les ending with ’png’ or ’surf.png’, which are the color images of the microstructure and

perimeter, respectively, examples of which are shown in Fig. 26; and

5e extension ’pfc’ is an acronym for phase fractions of cement (or clinker).

26
3 APPLICATION TO SRM 2686A

3. les ending with a phase name, which store the two-point correlation function results for

that phase (see example in Fig. 27).

Figure 26: Contents of the pfc le for the segmented image of SRM 2686a in Fig. 20. Here the

column heading “Volume” corresponds to the phase area fractions and “Area” to the perimeter

fraction.

(a) (b)

Figure 27: (a) False color PNG image of SRM 2686 microstructure, and (b) the phase boundaries.

Field width is 256 μm.

3.5.3 Two-point correlation functions

e resulting data les from MicroChar are mostly self-explanatory, as shown in Figs. 26–28.

However, correlation functions may be less familiar to the reader, so they will now be explained

in more detail.

e two-point correlation function for a phase, S2 (r), is the probability that two randomly

selected points in the image, separated by a distance r, both lie within that phase. In random

composite materials, S2 has some generic features that can be used to characterize the spatial

distribution of the phase. For example, when r = 0, the points coincide; the probability of a

randomly selected point lying within a phase in a 2D image is equal to the area fraction, φ, of

that phase, on a total area basis. erefore, S2 (0) = φ, that is, the y-intercept of the plot is equal

27
3 APPLICATION TO SRM 2686A

(a) (b)

Figure 28: (a) Contents of the alite two-point correlation function data le for SRM 2686a image

in Fig. 20 with distance in the rst column and probability density in the second column; (b) plot

of the two-point correlation function S2 for alite and for belite in the SRM 2686a image (right).

to the phase’s area fraction. A second property of S2 (r), which will not be proven here, is that its

limiting slope as r → 0 is directly proportional to the phase’s boundary length per unit area, Σ.

In particular, for granular media the exact relation can be shown to be [3]

lim

r→0

+

dS2

dr = −

Σ

4

(1)

A third property of S2 is that, at suciently large distances between the two points in a

random microstructure, the phase identity of each point should be independent of the other. In

other words, the phase distribution should be uncorrelated at large distances. When two events

are statistically independent, the joint probability is equal to the product of the (unconditional)

individual event probabilities. As stated earlier, the probability of a single point lying within a

phase is the phase’s area fraction, φ, so S2 at suciently large distances should converge to φ

2

.

Examination of the plot in Fig. 27 shows that S2 (0) = 0.4457 for alite, which therefore is alite’s

area fraction on a total image area basis. Furthermore, the slope atr = 0 for alite is approximately

−0.0158 μm−1, which provides a semiquantitative indication of the specic boundary area of alite

in the image. Finally, the plot indicates that S2 for neither alite nor belite appears to be converging

to φ

2

, which is 0.198 for alite and 0.1024 for belite. is can be interpreted to mean that the

spatial distribution of each phase is still somewhat correlated at distances as large as 200 μm. A

cursory inspection of Fig. 21 suggests that both alite and belite are clustered over a fairly large

length scale, so it is not surprising that the phase distributions are still correlated at 200 μm. is

hypothesis can be tested by examining a larger segmented image eld of SRM 2686a, as discussed

in Section 3.6.

Fig. 29 shows the two-point correlation functions for the minority phases arcanite and periclase. is plot highlights two other features of two-point correlation functions. First, the cor28
3 APPLICATION TO SRM 2686A

relation functions appear to be converging to values of about 10−3for periclase and about 10−4

for arcanite at distances no greater than about 100 μm. At r = 0, the values of φ for periclase

and arcanite are 0.0337 and 0.0046, respectively. Squaring these values gives 0.0011 for periclase

and 0.00002 for arcanite, both values of which are reasonably close to the long-distance values.

erefore, both periclase and arcanite are concluded to lose any spatial correlation at distances

of about 100 μm.

Figure 29: Two-point correlation functions for periclase and arcanite in SRM 2686a. Note the

two dierent scales on the y-axes.

e second feature of the two-point correlation functions in Fig. 29 is the appearance of a

local minimum at short distances followed by a small-amplitude oscillation in the curves. e

periodic local minima are characteristic of granular media and of non-overlapping domains in

composites, such as periclase and arcanite in SRM 2686a (Fig. 21). In such materials, the distance

at which the rst local minimum occurs corresponds to about half the average equivalent radius

of the domain, and the distance at which the following local maximum happens corresponds to

the average domain spacing in the microstructure [3]. Inspection of Fig. 29 suggests that the

average equivalent radius of periclase and arcanite domains is about 31 μm and about 18 μm,

respectively. Likewise, the average domain spacings of periclase and arcanite are about 40 μm

and about 27 μm, respectively. Note that these same kinds of features cannot be extracted from

correlation functions when the phases are signicantly correlated over long distances, as appears

to be the case for alite and belite in Figs. 27 and 28.

29
3 APPLICATION TO SRM 2686A

3.6 Larger Image Field

Image analysis statistics may improve by using larger elds of view that can capture more of the

intrinsic heterogeneity within the microstructure. For example, Fig. 30 shows a larger image eld

for SRM 2686a, which was acquired at 200× magnication, with a resolution of 1.5 μm/pixel.

(a) (b)

Figure 30: Larger image eld of SRM 2686a. Magnication 300×, 1024×768 pixels. (a) False color

PNG image, and (b) the phase mass, volume, and area fractions. Field width is 1.536 mm.

Notice that the abundances of alite and belite, especially, (Fig. 30(b)) are signicantly dierent

than those measured on the smaller image in earlier sections (Fig. 20). For example, the volume

fraction of alite is 32 % greater, and that of belite 47 % less, than in the smaller image. One would

ordinarily place more condence in phase fractions measured on larger elds of view, which

should reduce errors due to heterogeneity in the microstructure. As already discussed, alite and

belite in the SRM 2686a clinker microstructure are still signicantly correlated over the eld of

view in the smaller image, as shown by the fact that their two-point correlation functions do

not converge to the square of their volume fraction at longer distances (Fig. 28). erefore, one

should be skeptical that the volume fractions of these phases measured on the smaller image are

representative of the material.

Fig.  shows the two-point correlation functions for alite and belite measured on the larger

eld of view in Fig. 30(a). In contrast to the corresponding functions for the smaller image, the correlation function for both alite and belite do converge, well within the eld of view, to the square

of their respective volume fractions. is result suggests that the larger image in Fig. 30(a) adequately captures all the heterogeneity within the material and, therefore, that phase abundance

estimates based on this image are probably more representative of their true values. In fact, as

shown in Table 3, the clinker phase mass fractions estimated from Fig. 30 are quite close to the

30
3 APPLICATION TO SRM 2686A

certied values6for this material.

Figure 31: Two-point correlation functions calculated for alite and belite using the large image

eld of SRM 2686a shown in Fig. 30.

Table 3: Comparison of clinker phase mass fractions, expressed as a percentage of total solids,

estimated for SRM 2686a from the larger eld of view in Fig. 30 to the certied values [15]. Measurements are based on one image eld, so no uncertainty estimate is given.

Phase Measured (%) Certied (%) Dierence (%)

Alite 61.13 63.35 -2.22

Belite 18.59 18.68 -0.09

Aluminate 2.52 2.46 0.06

Ferrite 13.95 10.76 3.19

Periclase 3.55 3.40 0.15

Arcanite 0.27 0.27 0.00

6Certied values of SRM characteristics, although obviously having some measurement uncertainty, are conventionally treated as “ground truth” for comparison purposes.

31
3 APPLICATION TO SRM 2686A

3.7 Multiple Image Statistics for SRM 2686a

When larger image elds like the one in Fig. 30 are unavailable or are dicult to capture due to

resolution limits of the microscope, intrinsic heterogeneity can also be accommodated by analyzing multiple image elds at higher magnication. Running averages and standard deviations

are tracked and checked for consistency with the SRM certicate and XRD values for the mass

fraction of each clinker phase. 7 e results, by imaging eld, are presented graphically in Fig. 32;

results for each image eld are also presented in Table 4 along with the certied values and 2-σ

uncertainties (95 % coverage). Consideration must be made that we are working with a mediumgrained, heterogeneous material with this clinker and that heterogeneity will be expressed as

a eld-to-eld variability in the measured results. Materials that are heterogeneous at coarser

length scales will have larger eld-to-eld variability which, as already shown in the previous

section, can be reduced by using larger imaging areas. However, image eld enlargement cannot

continue indenitely without compromising the resolution and losing grain boundary denition.

In fact, multi-phase microstructures will always have eld-to-eld variability that generally increases as the image area is reduced. erefore, uncertainties measured on multiple image elds

provide an estimate of the heterogeneity length scale relative to the image area: larger uncertainties correspond to greater heterogeneity at the chosen length scale.

Another means of visualizing the results is with the box plot (Fig. 33), where the median

value is used for its resistance to outliers, and is identied by the “X”, the interquartile range (i.e.,

“middle half”) of the data are represented by the body of the box, and the extremes (minimum and

maximum) are represented by the ends of the straight lines projecting out of the box. e stars

on the plot correspond to the certied values, also tabulated in Table 3. e constant sum eect is

also apparent in Fig. 33 with the values varying inversely with the major phase estimates. is is

one possible source of bias with normalized estimate methodologies. A cumulative mean estimate

(running average) is presented in Table 4 and in graphical form in Fig. 32. For this clinker and

sampling conditions, the values appear to stabilize aer four elds of view and the nal estimates

are close in value to the certied values.

7Consistency with a certied value is taken to mean that the mean value of a measured characteristic lies within

one standard deviation of the certied value.

32
3 APPLICATION TO SRM 2686A

Table 4: SRM 2686a phase fractions expressed as a mass percent by eld and a running average.

Mass Percent by Field

Field Alite Belite Aluminate Ferrite Periclase Alkali sulfate

1 66.4 18.8 1.8 9.4 3.1 0.6

2 68.1 16.4 1.9 9.4 3.1 1.0

3 64.1 18.7 2.2 11.0 3.0 1.0

4 64.4 17.6 2.5 11.7 3.0 0.8

5 65.1 17.0 2.4 11.5 3.0 0.8

6 65.3 16.3 2.6 11.9 3.1 0.8

7 64.3 17.0 2.7 12.2 3.0 0.8

8 63.5 18.0 2.8 12.0 3.0 0.7

9 64.3 17.3 2.7 12.0 3.0 0.6

10 64.7 17.5 2.6 11.6 3.0 0.6

11 64.3 17.6 2.6 11.9 3.0 0.6

12 64.8 17.3 2.5 11.7 3.0 0.6

13 64.8 17.0 2.6 11.9 3.0 0.6

14 64.9 17.1 2.5 11.8 3.1 0.7

Mass Percent: Running Average

1 66.4 18.8 1.8 9.4 3.1 0.6

2 67.2 17.6 1.9 9.4 3.1 0.8

3 66.2 18.0 2.0 9.9 3.1 0.8

4 65.8 17.9 2.1 10.4 3.0 0.8

5 65.6 17.7 2.2 10.6 3.0 0.8

6 65.6 17.5 2.2 10.8 3.1 0.8

7 65.4 17.4 2.3 11.0 3.0 0.8

8 65.2 17.5 2.4 11.1 3.0 0.8

9 65.1 17.5 2.4 11.2 3.0 0.8

10 65.0 17.5 2.4 11.3 3.0 0.8

11 65.0 17.5 2.5 11.3 3.0 0.8

12 64.9 17.5 2.5 11.4 3.0 0.7

13 64.9 17.4 2.5 11.4 3.0 0.7

14 64.9 17.4 2.5 11.4 3.0 0.7

2-σ 2.2 1.4 0.6 1.8 0.2 0.2

Certicate 63.53 18.80 2.46 10.80 3.40 0.86

2Uc 1.04 1.10 0.39 0.84 0.23 0.17

33
A STEP-BY STEP SUMMARY FOR CLINKER ANALYSIS

Figure 32: Running average plots for SEM/IA analysis of SRM2686a shows values beginning to

stabilize aer about ve image elds.

Appendix A Step-by Step Summary for Clinker Analysis

ImageJ is used to preview the SEM BE and XR images and to interactively make some adjustments

to the background and noise by improving the image quality. We always retain the original images

and use the File→Save As command to save a modied image. Appending the lename with an

M to indicate the operation performed on the image indicates that it has been modied. By

appending the lename with an operation code, the images will appear together in the directory,

allowing for easy viewing of the available images. Always save as a TIFF le format to eliminate

any image quality loss with subsequent processing operations.

1. Read images individually or by drag-and-drop with ImageJ.

34
A STEP-BY STEP SUMMARY FOR CLINKER ANALYSIS

Figure 33: Box plots by phase for SEM analysis of SRM 2686a. Stars indicate certied values.

35
A STEP-BY STEP SUMMARY FOR CLINKER ANALYSIS

2. Tile images to arrange all on the monitor (Window→Tile)

3. Assign a gray LUT to each image (Image→Lookup Table→Grays)

4. Adjust each image’s brightness and contrast to assess the noise associated with low element

concentration or the background (image noise): Image→Adjust→Brightness Contrast

5. Subtract background noise with operation (Process→Subtract) using the preview check box

in the Subtract operation window to interactively select the noise level. When acceptable,

select OK.

6. Clean up any noise by ltering. e Despeckle and the Remove Outliers commands under Process→Noise to remove isolated pixels or small pixel groups. e median lter (Process→Filters→Median) or the resholded Blur will also work well to reduce image noise

while retaining edge details. Apply any ltering lightly as edge blur may result in problems during the segmentation phase. Save the modied image using File→Save As→TIFF,

appending the original le name with an “M” to indicate this le has been modied, so that

the original is not over-wrien.

7. Use MultiSpec to generate a segmented, indexed image where each pixel is uniquely assigned to a class to which it most likely belongs based upon a user-developed training set

of data. An abbreviated set of steps, which will be detailed below include (a) reading image

les, (b) linking les, (c) designating phase classes, (c) establishing training elds for each

class, (d) selecting a classication algorithm, (e) evaluating the classication, (f) saving the

project, and (g) saving the classied image as an indexed, segmented TIFF le. is image

le will then be used for the image analysis for determining phase area fraction, phase mass

fraction, and phase perimeter fraction. e resulting image and phase-related data are useful characteristics of the microstructure and for introduction into simulation models like

VCCTL [4, 21] or THAMES [5].

8. Read initial image as a multispectral type (File→Open Image), and check “multispectral

type” box.

Under “Multispectral Display Specications” accept the defaults with the exception of the

“Enhancement” section, where you can opt to not change the image (the usual), clip it on

the low and high end, or user selected clipping.

9. Link the set of X-ray images deemed useful in making a segmentation: File→Open Image.

Note: is operation requires selecting “Link to active image window” from the drop-down

box.

Select the X-ray images by clicking on them individually. Remember to select the modied

ones (M). I use a standard image sequence so I can recall the exact image source in the virtual

linked les (BEI, Al, Mg, K, Na, S) or write it down. ey can be named later if necessary

and excluded in the classication if desired. Click on the Cancel buon to complete the

linking.

36
A STEP-BY STEP SUMMARY FOR CLINKER ANALYSIS

10. Save the image set (Processor→Reformat→Rectify Image)

Use the window defaults (you could aempt to correct any displaced images here if needed),

make sure the path is correct and provide a name, retaining the .LAN extension. Use the

original le name to keep them together in a listing. If you do not include the extension,

the le will be wrien without one and will not be recognized for future Read operations.

11. Open the .LAN image set for review, assigning classes and training elds, and classication

(File→Open Image)

You will then see a window labeled “Multispectral Display Specications,” where you can

assign the channels to a color (R G B) and enhance the images by clipping or transformation. We did the enhancements using ImageJ since Multispec seems to apply a single set of

operations to each image. View the image using dierent channel/color combinations. e

most useful are:

• BE, Mg, Al — alite, belite, aluminate, ferrite, periclase, lime, voids

• BE, K, S — alkali sulfates, calcium sulfates

12. Establish training classes (Processor→Statistics)

Accept the defaults in the “Set Project Options” window by clicking OK.

In the new Project window you will see a Class Designation at the top. Enter the list of

classes here by boxing in a region in the image for each, selecting “Add to List” and providing a name. e classes are listed based on this sequence and the ultimate index value

used for VCCTL or THAMES, following this convention:

1 = alite 2 = belite 3 = aluminate 4 = ferrite

5 = periclase 6 = arcanite 7 = aphthitalite 8 = void

Record the index designations for later use as some phases will be absent from some materials.

Subclasses can be created—for example, alkali sulfate could contain sub-classes arcanite

(potassium sulfate) and aphthitalite (potassium, sodium sulfate).

Add training classes then by selecting the class, boxing in more regions, and specifying

“Add to List”. e more training pixels, the beer the class characteristics are dened.

13. Classify (Processor→Classify)

Specify the “Minimum Euclidian Distance” algorithm as it works consistently well. Check

the box “Write Classication Results” to Disk File and check the path, retaining the .GIS

extension for this image.

37
REFERENCES

Initial quality check is the percent re-assignment of the training set pixels and the quality

of the class assignment seen in the .GIS image. e .GIS image is read into Multispec as a

Multispectral type image and you will see the class assignments with distinct colors. e

colors may be changed by double clinking on the class list followed by seing the color

using the color picker.

14. Save the work as a Project

Even if not successful, you can go back and add training elds, delete elds, run alternate

classication algorithms, and even exclude specic images to see if it can be improved.

15. Save the .GIS le as a TIFF le.

is preserves the color assignments for the phases and the index values. Do not save these

in the .JPG format as the index values and colors will not be saved. You also do not want

to overwrite any les so append the le name with a MS sux.

16. Clean up any mis-assigned pixels using a surface blur lter that retains edge details.

e plugin “resholded Blur” works well for this by ignoring pixels that deviate from the

mean and by allowing you to specify the kernel size (radius, start with 2 or 3), threshold (1

or 2 as this is an indexed image and the values dier by 1), soness (aects the edge blur)

and strength (number of times to apply).

17. Read the TIFF le into ImageJ and write the ASCII Indexed le (File→Save As→Text Image).

You now have an ASCII le with the phase index values. If needed, you can recall the .GIS le

in Multispec or the *MS.TIF le in ImageJ to see which index value corresponds to each phase.

References

[1] “ASTM C 183: Standard Practice for Sampling and the Amount of Testing of Hydraulic

Cement”. In: Annual Book of ASTM Standards. Vol. 4.01. West Conshohocken, PA: ASTM

International, 2014.

[2] P. H. Bates and A. A. Klein. Properties of the Calcium Silicates and Calcium Aluminates

Occurring in Normal Portland Cement. Technological Papers of the Bureau of Standards 78.

Washington, D.C.: U.S. Bureau of Standards, 1917.

[3] J. G. Berryman and S. C. Blair. “Use of Digital Image Analysis to Estimate Fluid Permeability

of Porous Materials: Application of Two-Point Correlation Functions”. In: J. Appl. Phys. 60.6

(1986), pp. 1930–1938.

[4] J. W. Bullard and P. E. Stutzman. “Analysis of CCRL Portland Cement Prociency Samples

Number 151 and Number 152 Using the Virtual Cement and Concrete Testing Laboratory”.

In: Cem. Concr. Res. 36 (2006), pp. 1548–1555.

38
REFERENCES

[5] J. W. Bullard et al. “Coupling ermodynamics and Digital Image Models to Simulate Hydration and Microstructure Development of Portland Cement Pastes”. In: J. Mater. Res. 26

(2011), pp. 609–622.

[6] D. H. Campbell. Microscopical Examination and Interpretation of Portland Cement and Clinker.

2nd. Skokie, IL: Portland Cement Association, 1999.

[7] M. A. Delesse. “Proced ́ e Mecanique pour Determiner la Composition des Roches (Mechan-  ́

ical Methods for Determining the Composition of Rocks)”. In: Annales des Mines 13 (4th

series) (1866), pp. 379–388.

[8] J. I. Goldstein et al. Scanning Electron Microscopy and X-Ray Microanalysis. 2nd. New York:

Plenum Press, 1981.

[9] F. Hofmanner.  ̈ Microstructure of Portland Cement Clinker. Holderbank, Switzerland: Holderbank Management and Consulting, Ltd., 1975.

[10] ImageJ: Image Processing and Analysis in Java. Accessed: 2015-03-16. url: http://imagej.

nih.gov/ij/.

[11] H. Insley and V. D. Frechee.  ́ Microscopy of Ceramics and Cements. New York: Academic

Press, 1955.

[12] D. A. Landgrebe. Signal eory Methods in Multispectral Remote Sensing. New York: Wiley

Interscience, 2003.

[13] Henri Le Chatelier. ˆ Experimental Researches on the Constitution of Hydraulic Mortars (English Translation). New York: McGraw-Hill, 1905.

[14] MultiSpec: A Freeware Multispectral Image Data Analysis System. Accessed: 2015-03-16. url:

https://engineering.purdue.edu/∼biehl/MultiSpec/.

[15] NIST Standard Reference Material 2686a—Portland Cement Clinker. Accessed: 2015-04-21.

url: https://www-s.nist.gov/srmors/view detail.cfm?srm=2686a.

[16] P. E. Stutzman. “Microscopy of Clinker and Hydraulic Cements”. In: Rev. Mineralogy Geochem.

74.1 (2012), pp. 101–146.

[17] P. E. Stutzman. “Scanning Electron Microscopy Imaging of Hydraulic Cement Microstructure”. In: Cem. Concr. Composites 26.8 (2004), pp. 957–966.

[18] P. E. Stutzman and J. R. Clion. “Specimen Preparation for Scanning Electron Microscopy”.

In: Proceedings of the 21st International Conference on Cement Microscopy. Ed. by L. Jany

and A. Nisperos. Las Vega, NV: International Cement Microscopy Association, Apr. 1999,

pp. 10–22.

[19] P. E. Stutzman, G. Lespinasse, and S. Leigh. Compositional Analysis and Certication of NIST

Reference Material Clinker 2686a. Technical Note. Washington, D.C.: National Institute of

Standards and Technology, 2008.

[20] resholded Blur: Another Edge-Preserving Averaging (Smoothing) Algorithm, also Known as

Surface Blur. Accessed: 2015-03-16. url: http://imagejdocu.tudor.lu/doku.php?id=

plugin:filter:thresholded%5C blur:start%5C&do=export%5C pdf.

39
REFERENCES

[21] VCCTL Soware. Accessed: 2015-03-16. url: http : / / www . nist . gov / el / building

materials/evcctl.cfm.

40