- 58 Views
- Uploaded on
- Presentation posted in: General

A Multi-Stage Expert System for Aerosol Classification

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

A Multi-Stage Expert System forAerosol Classification

Statisticians:

Raymond Mugno and Wei Zhu

Computer Scientists:

Peter Imrich and Klaus Mueller

Environmental Chemists:

Dan Imre and Alla Zelenyuck

- Working with Environmental Chemists at Brookhaven National Laboratory
- Design an Expert System to classify aerosols
- Existing tools are not accurate enough
- Existing tools are not fast enough

- Use statistical techniques to reduce data
- Use visualization techniques to show reduction
- Give experts a tool to view classification
- Experts can change how data is reduced, by moving particles and making “classification rules”

- Aerosols are collected by a SPLAT mass spectrometer
- 5200 voltage readings are collected for each particle
- The 5200 voltage drops are converted to a 450 dimensional vector of intensities
- These vectors are the mass spectra

- Each spectra can be represented as vector, v = (v1, v2, …, v450)
- The subscript, j=1,2, … 450, represents the atomic weight of an element in the particle
- The value of vj represents the amount of that element

Dataset

- Over 1 million mass spectra collected from Houston, Texas
- Filtered to 238,160 actual particles
- Each particle has a unique time stamp down to millisecond
- Each particle has a time of flight, from which mass can be determined
- Each particle has a mass spectra, an ordered array of 450 integer peaks

- Level 1 (dimension reduction): Classify large numbers of mass spectra into clusters with very similar particles, i.e. similar mass spectra, using K-means clustering.
- Level 2 (class determination): Guided by chemical experts, combine clusters using binary classifiers. Determine the particle/cluster membership (acid, aromatics, or finer classes).

- Statistical tool used to classify multi-dimensional entities
- Hierarchical clustering:
- Non-Hierarchical clustering

- Start with a k seeds (representative entities of the k clusters) and a threshold distance
- For each entity find the distance to each seed
- If the minimum distance is less than the threshold, add the entity to that cluster
- If the minimum distance is greater than the threshold, the entity becomes the seed of a new cluster

- After iterating through all the entities, update the seeds
- Iterate though the particles again
- Continue iterating through the particles until none of the entities change clusters or other criterion is met
- Differs from classification, because end number of classes and the classes themselves are not set in advance

- Start with 25 seeds, average spectra of particle class, from experts
- For each particle, the distance between its spectra and each seed’s spectra is calculated
- If the minimum distance is less than a threshold distance, the particle is put into that corresponding cluster
- If the minimum distance is greater than the threshold distance, the current particle is set as a new seed

- 1 – r
- r is the Pearson Correlation Coefficient
- We label seed spectra as (x1,x2,…,xn) and the particle spectra as (y1,y2,…,yn) -- where Xi or Yi represents the magnitude of the ‘peak’ at location i is mass to charge ratio (i=13,14,…,250)

X39=20

X41=5

X53=2

Y26=1

Y28=25

Y68=2

- Spectra with “peaks” of the same proportion at the same locations will have small distance between them
- Classify similar shaped spectra together for dimension reduction

- Started with 25 seeds
- Threshold to create a new cluster set to 0.3
- Processed 238,160 particles (5 iterations)
- Finished with 2000 clusters/seeds
- Seeds are updated to be average of spectra in cluster
- Dimension reduction of 120 fold

19 major peaks from 27 to 97

4 major peaks at 54, 55, 56 and 57

6 Major Peaks at: 23, 24, 28, 30, 36, and 39

6 Major Peaks at: 53, 54, 55, 56, 57, and 58

19 Major Peaks from 23 to 131

Measuring Within-Cluster Similarity

Calculated the average distance for each particle to its clusters center

Calculated the standard deviation of the distances to the cluster’s center for each cluster

Found the particle furthest from the cluster’s center

- Using a distance of threshold of 0.3 yielded clusters where the within cluster similarity level is very high
- Ideal for dimension reduction – Now instead of working with 238,160 original spectra, we can work with the 2000 seeds for the clusters for a second level classification
- A dimension reduction of 120 fold!

Second Level Classification

- 2000 clusters is too many
- Find clusters that are very similar and combine them
- Find clusters that are very similar and classify them into a general group and have the chemical experts sub divide the general groups

Second Level Classification

Hierarchical Clustering

- Find the pair wise distance between each entity
- Merge the two “closest” entities
- Repeat the procedure until there is only 1 entity, or until merging distance threshold is met

Second Level Classification

Hierarchical Clustering

- Simple Linkage (closest entities)
- Average Linkage (average distance over all entities)
- Centroid Linkage (distance between average entities)
- Complete Linkage (distance between furthest elements)

Second Level Classification

- Find clusters that are very similar and classify them into a general group and have the chemical experts sub divide the general groups
- Using Centroid Linkage
- Each cluster is represented by a seed, the average spectra of that cluster

Second Level Classification

- Use Binary Matching metric to group clusters
- From a particle spectra vector v create a binary vector w.
- If vi > Peak_Threshold of total of vi’s peaks, wi = 1, else wi=0
- Experts gave Peak_Threshold =10 to filter out noise

Second Level Classification

Metric for comparing 2 binary vectors, w and x

Binary score = Number of peaks in common

Max Peaks = Maximum number of peaks between x or w

Distance = 1 – (Binary Score)/Max Peaks

Second Level Classification

Circular Dendrogram

Seeds are located around the circumference

More similar seeds are merged closer to the outer edge of the circle

Second Level Classification

Second Level Classification

User can zoom on different area of the dendrogram

From dendrogram, user can obtain seed spectra

From dendrogram, user can obtain cluster information, such as number of particles

Second Level Classification

Conclusion

Chemists now have a tool to view distribution of the data

Can look at 2000 seeds instead of 238,000 particles

Gives chemists insight to distribution of particle classes

Chemists can give feedback that will improve metrics

- Time series distribution of average (or median) aerosol size
- Time series distribution of atmospheric composition
- Time series models for oxidation rate comparison

Future Work - Comparison of Classifiers

1. Hierarchical clustering (Hinz et al, 1995)

2. Non-hierarchical clustering (Trieger et al, 1995)

3. Fuzzy clustering (Hinz et al, 96/99; Trieger et al, 95)

4. Discriminant analysis (Alsberg et al,1998)

5. Neural network (Song et al, 1999)

6. Classification tree (Harrington et al, 1989)

- Which approach is the best?

- Analyzing the spatial & temporal evolution of air-borne particles
- Random field theory
- Interactive classification & analysis via the graphical user interface
The Human-Machine Interface: examine chemical significance of clusters; adjust classifiers for better classification; explore the spatial/temporal trend & model goodness-of-fit graphically ...

- Technology Takes on Terrorism
- By Earl LaneWASHINGTON BUREAUNewsday: February 26, 2002