Rodolfo
Cucchiella
a,
Giuseppe
Falini
b,
Massimo
Ferri
a,
Milena
Stracquadanio
c and
Claudio
Trombini
*b
aDipartimento di Matematica - Università di Bologna, Piazza di Porta S. Donato 5, 40126, Bologna, Italy
bDipartimento di Chimica “G. Ciamician”, Università di Bologna,
via Selmi 2, 40126, Bologna, Italy. E-mail: claudio.trombini@unibo.it; Fax: +39 051 2099456; Tel: +39 051 2099513
cDipartimento di Chimica Industriale e dei Materiali, Università di Bologna, Viale del Risorgimento 4, 40136, Bologna, Italy
First published on 7th November 2008
Mathematical morphology is a tool for extracting image components that are useful for representation and description. The technique consists of a set-theoretic method of image analysis providing a quantitative description of geometrical structures. A simple application of mathematical morphology to a bi-dimensional processing of TEM images of airborne particles allows us to distinguish between particles grown and/or transported in atmosphere under dry conditions or in rainy days by a simple comparison of the corresponding image form factors. The form factors range in the 0.385–0.031 interval in the case of particles sampled in rainy days, and in the 0.103–0.006 interval in the case of non-rainy conditions. The same classification criterion was applied to filters collected under dry conditions and plunged in water. The results demonstrate that a morphological change may be artificially induced to the particle structure. The artificially wet particles, indeed, display an apparent contraction of their structures evidenced by a two-fold increase of the average values of their form factors. The last experiment roughly simulates the impact of particles on membranes of the respiratory tract.
A major contribution to the total budget of atmospheric particulate matter in urban and industrial areas, is given by gasoline and diesel powered engines exhausts, which deliver soot, mainly consisting of carbonaceous primary nanoparticles.11–13Combustion-generated nanoparticles are composed of carbon spherules typically sized between 20 and 50 nm,14–16 and held together by non-directional cohesive forces (van der Waals interactions, electrostatic interactions, π-stacking interactions) to form branched aggregates and clusters of variable size that contain up to hundreds of individual spherules.15,17,18 Once emitted to the atmosphere, the condensation of primary particles or the formation of new volatile particles through nucleation gives rise to secondary particles.
In addition to composition and size distribution, particle morphology is an important factor affecting their toxicity, optical properties,19,20 and lifetime in the atmosphere.21
After the first morphological studies on airborne particulate collected in urban areas using transmission electron microscopy (TEM),22 efforts have been focused on determining the effect on morphology of primary sources, fuels and combustion conditions. For instance, the engine load affects particles morphology and fractal dimensions of particles.23 Particles from exhausts of a direct-injected heavy-duty diesel engine working under different operating conditions,20,24 as well as particles emitted from spark ignition engines of light-duty vehicles have been examined by TEM.25 Moreover, morphological investigations by electron microscopy has been extended to particles deriving from combustion of various fuels such as residual oil,25 or coal and tires in power plants.26
Besides to its primary source, particle morphology depends on general weather conditions such as temperature and relative humidity.27
The research question in this study was to verify if the structure of airborne particles is affected by water and to describe this effect in a quantitative way.
The experimental design involved collecting particles by sampling PM2.5 in quartz fiber filters (QFF), subjecting particles to an extensive TEM analysis, and processing TEM images with mathematical morphology (MM) techniques.
Mathematical morphology (MM) provides a powerful tool for TEM image processing. In this study MM gives us a solution to quantitatively distinguish between particles grown and transported under dry conditions and particles exposed to the action of liquid water during transport (rain) or by impact on an aqueous phase.
A grey-tone image is a matrix M, whose entries are integer numbers in the interval [0, 255], where 0 represents black, and 255 white. A binary image is a matrix N whose entries are either 0 (black) or 1 (white). Thresholding consists in deciding a threshold value L and substituting a grey-tone image M with a binary image N of the same size, where each N(i,j) = 0 if M(i,j) < L, and N(i,j) = 1 if M(i,j) ≥ L.
Erosion and dilation are substitutions of a given binary image N with a new one. They both depend on the choice of a structuring element, i.e. an auxiliary matrix B (usually small with respect to N) whose indices may be null or negative. Erosion substitutes the matrix N with a matrix N′ where N′(i,j) = 1 if and only if N(i + h,j + k) = 1 for all h,k such that B(h,k) = 1, and N′(i,j) = 0 otherwise. Dilation substitutes N with a matrix N″ where N″(i,j) = 1 if and only if there exists at least one h and one k such that N(i + h,j + k) = 1 and B(h,k) = 1, and N″(i,j) = 0 otherwise. Erosion and dilation are the basic manipulations in mathematical morphology.28
On all images of this research, both of training and test sets, the threshold was set at two different values: T1 = 193 and T2 = 200, in order to catch different detail levels. The choice of T1 and T2 had been done by first applying four different thresholds to some randomly chosen images. The lowest and highest thresholds turned out to be unfit, since they destroyed the relevant information of some of the images. The two middle thresholds were then progressively tuned until they gave satisfactory results on all of them. We decided to keep two different (although close) values in order to test the form factor method, both with greater and rougher detail.
The morphological processing was then performed by applying three dilations in sequence, followed by three erosions. The structuring element was always a square of side three whose entries were all equal to 1.
As a parameter for assessing the compactness of the image, we choose the form factor (ff) expressed as:
ff = 4πA/p2 |
In order to test the hypothesis that “wet” particles are more compact, we have computed the ff of all images after thresholding both with T1 and T2, and carrying out the sequence of dilations and erosions. The whole set of black pixels was then considered to be the region of interest of each image. Manual post-processing turned out to be necessary, for getting rid of the grid support from the images. Nothing could be done for correcting the images in which the region of interest is not totally contained in the frame. This is of course a drawback that affects the computation of the ff.
What we obtain is a two-class classifier. It is based on a cutoff value F of the ff. Explicitly, the program classifies an image as positive (i.e. “wet” particle) if ff ≥ F, negative (i.e. “dry” particle) otherwise. Of course, even without the aforementioned problems, one should always expect some “dry” particle images to have ff ≥ F (false positive) and some wet ones to have ff < F (false negative). Two parameters are used for assessing the quality of a classifier: sensitivity (Se) and specificity (Sp). They are based on an a priori knowledge of the real status of the sample as dry or wet (ground truth) which will be considered as absolute.
Se is defined as the ratio
Se = (# of wet samples classified as such)/(total # of wet samples) |
Sp = (# of dry samples classified as such)/(total # of dry samples) |
So, a high sensitivity means that the classifier yields few false negatives, and a high Sp that there are few false positives.
Of course, it would be desirable to have both parameters equal to 1, but this is very rarely possible in practice. A very low cutoff F may yield a high Se, even equal to 1, but at the cost of a low Sp. By gradually raising F, we shall be able to increase Sp, but eventually Se will lessen. Actually, by using F as a parameter, we can trace a curve, called Receiver Operating Characteristic (ROC curve), traditionally plotted with 1-Sp as the abscissa x and Se as the ordinate y, which depicts the global behavior of the classifier. The farther from the diagonal y = x it is, the better the classifier.
Given the ROC curve of a classifier, it can be used in some different ways for choosing a “good” cutoff value F. In the ideal case of a perfect classifier, the ROC curve would simply follow the left and upper sides of the square [0,1] × [0,1] and the optimal cutoff, actually splitting the two classes, would correspond to the upper left corner. With a real ROC curve, one must be content with an approximation. At least two cutoff values are available: the value F′ corresponding to the intersection of the ROC curve with the diagonal y = 1 − x, and the value F″ corresponding to the point of the ROC curve nearest to the upper left corner.
PM2.5 samples were collected on quartz fiber filters (QFF) (QF 20 Schleicher & Schuell) by means of a Bravo PM HV TCR Tecora instrument equipped with a EN LVS PM2.5 size selective inlet, which complies with the European norm EN 12341 (air flux = 2.3 m3 h−1). Sampling details are given in ref. 30.
Henceforth, samples and particles collected in non-rainy days are operationally defined as “dry” samples or particles, those collected in rainy days are defined as “naturally wet”. No attempts were done to correlate particles to the mean relative humidity of the sampling day. To test whether impact with liquid water brings about morphological changes to particles, a few QFF's sampled in non-rainy conditions were cut in two portions. A half filter was dipped either in milliQ®water or in a physiological solution for 12 hours in a Petri dish, then dried in a dessicator for 48 hours at room temperature (T ∼ 20–25 °C). These last samples were referred as “artificially wet”.
To collect particles on a suitable TEM support, two approaches were used: (i) ten 35 mm copper grids, leaned down on the QFF surface, captured impacting particles during sampling; (ii) particles were mechanically transferred from a QFF to grids by softly rubbing grids on the QFF surface. Grids were stored at room temperature in a TEM grid storage box until the analysis. The digitized TEM images were collected by means of a Philips CM100 TEM instrument at a voltage of 80 KV, using a high resolution closed-coupled device (CCD) camera. The TEM was operated in bright field mode. The images were divided in two sets. The first set consisted of 31 images of “dry” particles and 32 images of “wet” particles and constituted the training set. The second one consisted of 32 images randomly extracted from the entire pool of images from dry and wet samples. It constituted the test set. Of course, the images used in the test set were not present in the training set to avoid any bias effect.
All the TEM images confirmed that particles collected in this urban site mainly derive from combustion processes, since they consist of branched chains aggregations of elementary units having the shape of spherules with diameter ranging from 20 to 50 nm, as widely documented in the literature.14–18 In the samples examined, most of submicronic particles were carbonaceous particles. Fig. 1a shows the typical TEM image of branching ultrafine soot aggregate, which is formed from individual primary spherules of 20–50 nm size. Generally, carbon speciation of diesel exhaust and urban particulate matter differentiates between two principal fractions: elemental carbon (EC) and organic carbon (OC). EC normally has almost no hydrogen content and is usually considered to be a product of combustion of fossil fuels. OC has a higher hydrogen content and lower oxidation resistance. Morphology of OC particles is that of big spherical particles (above 50 nm) either isolated or attached to the soot aggregates. In Fig. 1b an OC-type carbonaceous particle attached to soot aggregate is indicated by the arrow.31,32 The TEM images collected from wet samples (Fig. 1c, d) do not show the presence of OC particles. Only spherules from 20 to 50 nm are observed, with the same size of the dry samples, but with a different shape, aggregation and branching.
![]() | ||
Fig. 1 Bright-field TEM images of particulate matter. (a) and (b) show two magnifications of “dry” particles; (c) “naturally wet” particles; (d) “artificially wet” particles. The OC-type carbonaceous particle attached to the soot aggregate is indicated by the arrow. Scale bars: (a) 500 nm; (b), (c) and (d) 100 nm. |
(i) to prove that airborne nano-structured particles show different morphology depending on they are sampled in rainy or non-rainy days; more compact structures are observed in rainy days with respect to those of particles sampled under dry conditions.
(ii) To propose a simple bi-dimensional image processing providing simple numerical descriptors for the ffs that adequately characterize and classify the nanostructured particles on the basis of a contact/non-contact with water.
(iii) To verify that dry particles after contact with water by plunging them into pure water or a pH 7 phosphate-buffered solution, undergo an appreciable and measurable morphological change, which consists, from a qualitative point of view, of a particle shrinking. The result is a substantial decrease of the empty volume between particle chains.
The experimental design was as follows. Two sets of images were provided; the ground truth (dry/wet) of the training set was explicit; for the test set it was hidden. All images of both sets were processed as described, with the two thresholds T1 = 193 and T2 = 200. For each of these thresholds the ROC curve was traced, and the cutoff values were derived: F′ (the intersection of the ROC curve with the diagonal y = 1 − x) is denoted as F1 for T1 thresholding and F3 for T2; F″ (the point of the ROC curve nearest to the upper left corner) is denoted respectively as F2 and F4, as shown in Fig. 2.
![]() | ||
Fig. 2 ROC curves for T1 (a) and T2 (b) thresholding. Se and Sp indicate sensitivity and specificity, respectively. |
Because of closeness of the two pairs of cutoff values, only one of the two points was chosen for our computations. Thus, the images of the test set have been classified by using F2 and F4 as cutoff values for the ff. Table 1 collects the ffs of the training set of images, while Table 2 shows ffs and classification for each image of the test set, against the ground truth; false negative and false positive are evidenced. The use of the lower threshold for binarization of images seems more efficient in keeping the relevant features sufficiently unaltered. However, with both thresholding values, there are three images which escape classification (Table 2, #18, #21 and #23).
Dry set | Wet set | ||
---|---|---|---|
Form factor | Form factor | Wetting modea | |
a Nat. indicates sample collected under “naturally wet” conditions; art. indicates “artificially wet” samples. | |||
1 | 0.0393 | 0.0700 | nat. |
2 | 0.0431 | 0.1809 | nat. |
3 | 0.0374 | 0.3395 | nat. |
4 | 0.0541 | 0.3848 | nat. |
5 | 0.0575 | 0.0498 | nat. |
6 | 0.0394 | 0.1949 | nat. |
7 | 0.0133 | 0.1419 | nat. |
8 | 0.0064 | 0.1962 | nat. |
9 | 0.0133 | 0.0689 | nat. |
10 | 0.0195 | 0.0581 | nat. |
11 | 0.0144 | 0.1266 | nat. |
12 | 0.0092 | 0.1667 | nat. |
13 | 0.0380 | 0.3113 | nat. |
14 | 0.0351 | 0.0602 | nat. |
15 | 0.0502 | 0.0308 | nat. |
16 | 0.1034 | 0.1429 | nat. |
17 | 0.0806 | 0.1113 | nat. |
18 | 0.0611 | 0.0366 | nat. |
19 | 0.0347 | 0.0228 | nat. |
20 | 0.0605 | 0.0294 | nat. |
21 | 0.0162 | 0.0765 | nat. |
22 | 0.0161 | 0.0319 | nat. |
23 | 0.0110 | 0.1041 | nat. |
24 | 0.0059 | 0.1211 | nat. |
25 | 0.0093 | 0.0525 | nat. |
26 | 0.0200 | 0.0734 | art. |
27 | 0.0336 | 0.0214 | art. |
28 | 0.0116 | 0.0306 | art. |
29 | 0.0174 | 0.0550 | art. |
30 | 0.0687 | 0.1074 | art. |
31 | 0.0310 | 0.0675 | art. |
32 | 0.0889 | art. |
Form factor | Classification F2 = 0.058 | Classification F3 = 0.051 | Ground truth | |
---|---|---|---|---|
a Wet indicates sample collected in wet conditions. b Dry indicates samples collected in dry conditions. | ||||
1 | 0.2521 | weta | wet | wet |
2 | 0.1591 | wet | wet | wet |
3 | 0.1504 | wet | wet | wet |
4 | 0.0515 | dryb | wet | dry |
5 | 0.0399 | dry | dry | dry |
6 | 0.0253 | dry | dry | dry |
7 | 0.0272 | dry | dry | dry |
8 | 0.0548 | dry | wet | dry |
9 | 0.0463 | dry | dry | dry |
10 | 0.0395 | dry | dry | dry |
11 | 0.0847 | wet | wet | wet |
12 | 0.0600 | wet | wet | wet |
13 | 0.0090 | dry | dry | dry |
14 | 0.0243 | dry | dry | dry |
15 | 0.0824 | wet | wet | wet |
16 | 0.0964 | wet | wet | wet |
17 | 0.2309 | wet | wet | wet |
18 | 0.0326 | dry | dry | wet |
19 | 0.0364 | dry | dry | wet |
20 | 0.1081 | wet | wet | wet |
21 | 0.0419 | dry | dry | wet |
22 | 0.0316 | dry | dry | dry |
23 | 0.0734 | wet | wet | dry |
24 | 0.0391 | dry | dry | dry |
25 | 0.0189 | dry | dry | dry |
26 | 0.0086 | dry | dry | dry |
27 | 0.0131 | dry | dry | dry |
28 | 0.0329 | dry | dry | dry |
29 | 0.0083 | dry | dry | dry |
30 | 0.2372 | wet | wet | wet |
31 | 0.1049 | wet | wet | wet |
32 | 0.1902 | wet | wet | wet |
Using this bi-dimensional image processing to determine the ff of airborne carbonaceous particles, a simple method to measure the particle shrinking effect due to exposure to water is proposed. Indeed, the ff of “dry” particles varies around the average value of 0.032 with a standard deviation of 0.023, the highest value observed being 0.103 (Table 3). On the other hand, “wet” particles show an average ff of 0.121 with a standard deviation of 0.087 and a maximum value observed of 0.385.
Average value | St. dev. | Maxa | Mina | |
---|---|---|---|---|
a max and min indicate the maximum and minimum value, respectively; wet indicates sample collected in wet conditions; dry indicates samples collected in dry conditions; wet-nat. indicates sample collected under “naturally wet” conditions; wet-art. indicates “artificially wet” samples. | ||||
dry | 0.032 | 0.023 | 0.103 | 0.003 |
wet | 0.121 | 0.087 | 0.385 | 0.021 |
wet-nat. | 0.137 | 0.092 | 0.385 | 0.023 |
wet-art. | 0.063 | 0.027 | 0.107 | 0.021 |
Even though a partial overlap of the ff ranges exists (borderline images where an unambiguous assignment is not possible), values higher than 0.087 (the 90% of wet samples) clearly are associated to airborne particles grown and/or transported in rainy days, while particles possessing a ff < 0.03 are definitely “dry particles” (Fig. 3).
![]() | ||
Fig. 3 Histogram showing the form factors (ff) calculated from “dry” and “wet” particles. The dashed line indicates the ff threshold below which the particle is classified as dry. |
This simple application of MM to TEM images looks as a straightforward tool to classify a particle originated and/or transported under dry or rainy conditions.
Let's now come to the last question, do morphology changes occur in dry particles after contact with liquid water?
The interaction of fine airborne particles with liquid water has a terrific interest when biological membranes are involved. Recently, light has been shed onto the interaction between lung lining liquid and inhaled ambient particles.33 An aggregation process has been observed to occur in this aqueous micro-layer, and it was claimed to be driven by the proteins present in the lining liquid. The new greater and more dense particles are more easily phagocytized by macrophages. The neat result is a more efficient particle clearance in the lung lining fluid, thus rendering particle aggregation a kind of lung protective process.
To prove that water brings morphological changes to particles one half of few QFF’s sampled in non-rainy conditions has been dipped either in milliQ®water or in a physiological solution and then dried. Images taken from the untreated half QFF and the wetted half have been compared.
The study carried out in a set of 7 filters sampled in dry days revealed that the simple action of liquid water shrinks aggregates and this phenomenon can be measured by comparing the ff of particles collected from the original and water-treated portions of the filters. While the average ff for the untreated particles is 0.032, the average value of particles ff from the water treated QFFs increases to 0.063, a clear indication that the direct contact with water is responsible for structure compression.
This journal is © The Royal Society of Chemistry 2009 |