1 / 38

A Multi-Stage Expert System for Aerosol Classification

A Multi-Stage Expert System for Aerosol Classification. Statisticians : Raymond Mugno and Wei Zhu Computer Scientists : Peter Imrich and Klaus Mueller Environmental Chemists : Dan Imre and Alla Zelenyuck. Our Project. Working with Environmental Chemists at Brookhaven National Laboratory

zalman
Télécharger la présentation

A Multi-Stage Expert System for Aerosol Classification

An Image/Link below is provided (as is) to download presentation 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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 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

  2. Our Project • 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

  3. Our Project: An Expert System • 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”

  4. Our Project: The Chemistry • 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

  5. 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

  6. 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

  7. The Two-Level Classification Scheme • 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).

  8. Clustering Analysis • Statistical tool used to classify multi-dimensional entities • Hierarchical clustering: • Non-Hierarchical clustering

  9. K-Means Clustering Analysis • 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

  10. K-Means Clustering Analysis • 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

  11. First Level Classification:K-Means Clustering • 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

  12. Distance Function • 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)

  13. Notation X39=20 X41=5 X53=2 Y26=1 Y28=25 Y68=2

  14. Why use Correlation Coefficient? • 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

  15. First Level Results • 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

  16. Cluster Example 1, based on an Organic seed (23) 19 major peaks from 27 to 97

  17. Cluster Example 2, based on an Fe seed(13) 4 major peaks at 54, 55, 56 and 57

  18. New Cluster --- Example 1 (464) 6 Major Peaks at: 23, 24, 28, 30, 36, and 39

  19. New Cluster --- Example 2 (640) 6 Major Peaks at: 53, 54, 55, 56, 57, and 58

  20. New Cluster --- Example 3 (574) 19 Major Peaks from 23 to 131

  21. 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

  22. Comments on the First Level Classifier • 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!

  23. 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

  24. 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

  25. 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)

  26. 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

  27. 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

  28. 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

  29. Second Level Classification Circular Dendrogram Seeds are located around the circumference More similar seeds are merged closer to the outer edge of the circle

  30. Second Level Classification

  31. 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

  32. Second Level Classification

  33. 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

  34. Future Work - Molecule Library

  35. Future Work - Time Series Analysis • Time series distribution of average (or median) aerosol size • Time series distribution of atmospheric composition • Time series models for oxidation rate comparison

  36. 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?

  37. Future Work - Spatial Temporal Analysis • 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 ...

  38. In the News... • Technology Takes on Terrorism • By Earl LaneWASHINGTON BUREAUNewsday: February 26, 2002

More Related