The ASTER Spectral Library provides over 2,000 high-resolution spectra for numerous classes of materials. [Baldridge2009] The spectra and associated metadata are provided as a large set of ASCII text files. SPy provides the ability to import the ASTER library spectra and a subset of the associated metadata into a relational databased that can be easily accessed from Python.
You will first need to get a copy of the ASTER Spectral Library, which is distributed on CD-ROM and can be ordered by filling out a simple online order form. Once you have the library, you can import the data in to an sqlite database as follows:
In : db = AsterDatabase.create("aster_lib.db", "/CDROM/ASTER2.0/data") Importing /CDROM/ASTER2.0/data/jhu.becknic.manmade.asphalt.paving.solid.0095uuu.spectrum.txt. Importing /CDROM/ASTER2.0/data/jhu.becknic.manmade.asphalt.paving.solid.0096uuu.spectrum.txt. Importing /CDROM/ASTER2.0/data/jhu.becknic.manmade.asphalt.paving.solid.0674uuu.spectrum.txt. Importing /CDROM/ASTER2.0/data/jhu.becknic.manmade.asphalt.tar.solid.0099uuu.spectrum.txt. ---// snip //--- Importing /CDROM/ASTER2.0/data/usgs.perknic.rock.sedimentary.shale.solid.phop011.spectrum.txt. Importing /CDROM/ASTER2.0/data/usgs.perknic.rock.sedimentary.shale.solid.phop02a.spectrum.txt. Processed 2443 files. Ignored the following 2 bad files: /CDROM/ASTER2.0/data/jhu.nicolet.mineral.silicate.tectosilicate.fine.albite1.spectrum.txt /CDROM/ASTER2.0/data/usgs.perknic.rock.igneous.mafic.colid.me3.spectrum.txt
Several files with malformed spectral were are skipped during the import process. If those spectra are important to you, you can always repair or adjust the data in the files as you see fit and repeat the import process.
The AsterDatabase create method was written specifically for version 2.0 of the ASTER Spectral Library and will not work for previous versions of the library.
Once the database has been created, it can be accessed by instantiating an AsterDatabase object for the database file (it can also be accessed by using sqlite external to Python). The current implementation of the SPy database contains two tables: Samples and Spectra. There is a one-to-one relationship between rows in the two tables but they have been separate to support potential future changes to the database. The schemas for the tables are in the schemas attribute of the database object:
In : db = AsterDatabase('aster_lib.db') In : for s in db.schemas: ...: print s ...: CREATE TABLE Samples (SampleID INTEGER PRIMARY KEY, Name TEXT, Type TEXT, Class TEXT, SubClass TEXT, ParticleSize TEXT, SampleNum TEXT, Owner TEXT, Origin TEXT, Phase TEXT, Description TEXT) CREATE TABLE Spectra (SpectrumID INTEGER PRIMARY KEY, SampleID INTEGER, SensorCalibrationID INTEGER, Instrument TEXT, Environment TEXT, Measurement TEXT, XUnit TEXT, YUnit TEXT, MinWavelength FLOAT, MaxWavelength FLOAT, NumValues INTEGER, XData BLOB, YData BLOB)
Descriptions of most of the Sample table columns can be found in the ASTER Spectral Library documentation. The sample spectra and bands are stored in BLOB objects in the database, so you probably don’t want to access them directly. The recommended method is to use either the get_spectrum or get_signature method of AsterDatabase, both of which take a SpectrumID as their argument.
As a simple example, suppose you want to find all samples that have “stone” in their name (this isn’t the best way to find stones in the library but it makes for an easy example). There are three ways you can issue queries through the AsterDatabase object. You can call its print_query method to print query results to the command line. You can call its query method to return query results in tuples. Lastly, you can use the cursor attribute of the AsterDatabase object to issue the query.
In : db.print_query('SELECT COUNT() FROM Samples WHERE Name LIKE "%stone%"') 78 In : db.print_query('SELECT SampleID, Name FROM Samples WHERE Name LIKE "%stone%" limit 10') 43|slate stone shingle 250|fossiliferous limestone 251|dolomitic limestone 252|limestone 253|oolitic limestone 254|lithographic limestone 255|argillaceous limestone 256|oolitic limestone 257|fossiliferous limestone 258|dolomitic limestone
Next, let’s retrieve and plot one of the results (we will take the last one).
In : f = plt.figure() In : s = db.get_signature(2436) In : import pylab as plt In : plt.plot(s.x, s.y) Out: [<matplotlib.lines.Line2D at 0x10d3e310>] In : plt.title(s.sample_name) Out: <matplotlib.text.Text at 0x10d222d0> In : plt.grid(1) In : plt.show()
The sqlite3 module is included with Python 2.5+. Details on sqlite3 connection and cursor objects can be found in the Python sqlite3 documentation.
While the AsterDatabase provides a Python interface to the ASTER Spectral Library, there may be times where you want to repeatedly access a small, fixed subset of the spectra in the library and do not want to repeatedly query the database. The ENVI file format enables storage of spectral libraries (see ENVI Headers). SPy can read these files into a SPy SpectralLibrary object.
To enable easy creation of custom spectral libraries, the AsterDatabase has a create_envi_spectral_library method that generates a spectral library that can easily be saved to ENVI format. Spectra in the ASTER library have varying numbers of samples over varying spectral ranges. To generate the library we want to save to ENVI format, we need to specify a band discretization to which we want all of the desired spectra resampled. Let’s pick the bands from our sample hyperspectral image.
In : bands = aviris.read_aviris_bands('92AV3C.spc') In : print bands.centers, bands.centers[-1] 400.019989 2498.959961
We see from the output above that the bands range from about 400 - 2,500 nm (we’re ignoring the fact that the bands at both ends have a finite width). We would like to find library spectra that cover the entire spectral range for our image, so we’ll check the band limits for the library spectra. But first, let’s check the units of the spectra in the library.
In : db.print_query('SELECT DISTINCT Measurement, XUnit, YUnit FROM Samples, Spectra ' + ....: 'WHERE Samples.SampleID = Spectra.SampleID AND ' + ....: 'Name LIKE "%stone%"') ....: directional (10 degree) hemispherical reflectance|wavelength (micrometers)|reflectance (percent) directional hemispherical reflectance|wavelength (micrometers)|reflectance (percent) hemispherical reflectance|wavelength (micrometers)|reflectance (percent) reflectance|wavelength (micrometers)|reflectance (percent)
We see that all spectra are measures of reflectance but wavelengths are in units of micrometers, whereas our sample image bands are in nanometers. To properly query the spectra, we will need to specify the band limits in micrometers.
In : db.print_query('SELECT COUNT() FROM Samples, Spectra ' + ....: 'WHERE Samples.SampleID = Spectra.SampleID AND ' + ....: 'Name LIKE "%stone%" AND ' + ....: 'MinWavelength <= 0.4 AND MaxWavelength >= 2.5') ....: 59
So it appears that 59 of the 78 “stone” spectra cover the desired band limits. To create a new, resampled spectral library for these spectra, we call the AsterDatabase create_envi_spectral_library method, passing it the list of spectrum IDs and our output band schema (in micrometers). Since our bands are defined in nanometers, we will convert them before calling the method.
In : rows = db.query('SELECT SpectrumID FROM Samples, Spectra ' + ....: 'WHERE Samples.SampleID = Spectra.SampleID AND ' + ....: 'Name LIKE "%stone%" AND ' + ....: 'MinWavelength <= 0.4 AND MaxWavelength >= 2.5') ....: In : ids = [r for r in rows] In : bands.centers = [x / 1000. for x in bands.centers] In : bands.bandwidths = [x / 1000. for x in bands.bandwidths] In : lib = db.create_envi_spectral_library(ids, bands)
Now that we’ve created a resampled library of spectra, let’s plot one. We’ll plot the same limestone spectrum plotted above, which happens to be the last spectrum in the resampled library.
In : import pylab as plt In : s = db.get_signature(2436) In : plt.plot(s.x, s.y, 'k-', label='original'); In : plt.hold(1) In : plt.plot(bands.centers, lib.spectra[-1], 'r-', label='resampled'); In : plt.grid(1) In : plt.gca().legend(loc='upper left'); In : plt.xlim(0, 3); In : plt.title('Resampled %s spectrum' % lib.names[-1]); In : plt.show()
The resampled spectral library can be used with any image that uses the same band calibration to which we resampled the spectra. We can also save the library for future use. But before we save the library, we need to change the band units to the units used in the band calibration (we need to convert from micrometers to nanometers).
In : lib.bands.centers = [1000. * x for x in lib.bands.centers] In : lib.bands.bandwidths = [1000. * x for x in lib.bands.bandwidths] In : lib.bands.band_unit = 'nm' In : lib.save('stones', 'Stone spectra from the ASTER library')
Saving the library with the name “stones” creates two files: “stones.sli” and “stones.hdr”. The first file contains the resampled spectra and the second is the header file that we use to open the library.
In : mylib = envi.open('stones.hdr') In : mylib Out: <spectral.io.envi.SpectralLibrary instance at 0xfd5c638>
|[Baldridge2009]||Baldridge, A. M., Hook, S. J., Grove, C. I. and G. Rivera, 2008(9). The ASTER Spectral Library Version 2.0. In press Remote Sensing of Environment|