Tag Archives: EO-1

Big Ben Eruption Update 2017-02-27

Mawson Peak's summit crater glows orange in this false-color infrared image (bands 7-6-5) taken February 27, 2017.  Image credit: Bill Mitchell (CC-BY) using data from USGS LANDSAT 8 (public domain).
Mawson Peak’s summit crater glows orange in this false-color infrared image (bands 7-6-5) taken February 27, 2017. Image credit: Bill Mitchell (CC-BY) using data from USGS LANDSAT 8 (public domain).

It has been three weeks since I reported on an active eruption on Heard Island seen by Landsat 8. Since then, the presence of lava at or near the surface in the summit crater of Mawson Peak has continued, and a thermal anomaly is present both in the February 27 Landsat 8 image shown above and in February 20 imagery. It is difficult to discern in the true-color imagery from February 27 whether there are any new lava/debris flows present. The two MODIS instruments (one on Aqua, one on Terra) have not picked up any thermal anomalies since early February.

Unfortunately, one of the best tools I’ve had at my disposal for keeping an eye on Mawson Peak is no longer available: NASA/USGS satellite EO-1 was decomissioned last week. EO-1 provided 10 m/pixel true-color imagery, which is significantly higher resolution than 15 m/pixel of Landsat. Archival data for both satellites remains available, but no new EO-1 data will be taken. New data from Landsat 8 typically comes in a few times each month (every 7-16 days), and I’ll be keeping an eye on it.

Mawson Peak in true color, February 27, 2017.  Image processing: Bill Mitchell (CC-BY) using USGS Landsat 8 data (public domain).
Mawson Peak in true color, February 27, 2017. Image processing: Bill Mitchell (CC-BY) using USGS Landsat 8 data (public domain).

Cloud-Free Heard Island

Composite, cloud-free satellite imagery of Heard Island, being produced in QGIS.  Image credit: Bill Mitchell (CC-BY), using USGS (Landsat 8, EO-1) data (public domain).
Composite, cloud-free satellite imagery of Heard Island, being produced in QGIS. Image credit: Bill Mitchell (CC-BY), using USGS (Landsat 8, EO-1) data (public domain).

Heard Island is a pretty cloudy place most of the time. However, there are occasional times when the weather clears, particularly on the southeastern (leeward) side of the island. On rare occasions, the northwest and southwest sides of the island come out from the clouds as a satellite passes over.

For the past two years, I have been watching Heard Island using true-color imagery from four satellites: Terra, Aqua, Landsat 8, and EO-1. I have posted previously about satellite imagery from these instruments. Although every image of the Island I have seen has clouds in it covering a portion of the island, I was curious whether or not I had accumulated clear imagery of the entirety of Heard Island.

In part, this question was spurred by a follower on Twitter asking about eruptive activity at Heard. I had to admit I didn’t really know whether the activity was low-level and continuous (like Kilauea) or more intermittent. Given that our knowledge of its eruptive activity is primarily from satellite observations, do the satellite “thermal anomalies” correspond to short eruptive events, or simply a cloud-free view of the volcano?

For high-resolution imagery of Heard Island, the datasets of interest are from EO-1 ALI, and Landsat 8 OLI. The two MODIS instruments (one on Aqua, one on Terra) are moderate-resolution, and while 250-m resolution is sufficient for some purposes, this one needs more. Looking through the archives, I was able to find EO-1 ALI data primarily for Mawson Peak and points southeast, and Landsat 8 OLI covered much of the island, particularly the northwest.

Not only is having cloud-free, high-resolution data important for me, but I want the data to be recent. There has been a retreat of up to 5.5 km for some of the glaciers since 1947, and the Google Maps imagery of that area (Stephenson Lagoon) is horribly outdated. Fortunately, I found most of the island covered in large swaths with images from 2014 onward, and mostly 2016. There was even good imagery from when I was on Heard Island! Our ship, the Braveheart, is visible as a few white pixels in Atlas Roads (just north of Atlas Cove), slightly closer to the Azorella Peninsula than to the Laurens Peninsula. The tents and campsite are too small and darkly colored to be visible on this image.

Braveheart in Atlas Roads and the campsite (non-contrasting) at Atlas Cove, Heard Island.  Satellite image pixels are 15 m across, and the Azorella Peninsula isthmus (along Campsite label) is 1 km wide.  Image credit: Bill Mitchell (CC-BY), using USGS (Landsat 8 OLI) data (public domain).
Braveheart in Atlas Roads and the campsite (non-contrasting) at Atlas Cove, Heard Island. Satellite image pixels are 15 m across, and the Azorella Peninsula isthmus (along Campsite label) is 1 km wide. Image credit: Bill Mitchell (CC-BY), using USGS (Landsat 8 OLI) data (public domain).

A small portion of the island between Atlas Cove and Mawson Peak was the most difficult to find. With the topography of the island, the steady stream of wind, and the humid air, the 2.5 km by 2.5 km region was cloudy pretty much all the time. Eventually, using the EO-1 ALI instrument and going back to early 2010, I found a reasonably clear image of it.

Once I had the images (after combining true-color and panchromatic brightness data in QGIS), I needed to stitch them together. Thanks to the wonderful QGIS training manual, I was able to create vector (polygon) layers which corresponded to the clear region of each image (plus some surrounding ocean). At this point the troublesome mostly-cloudy spot became evident, and the search was on for imagery to fill the void.

Creating polygons for clipping the satellite imagery using QGIS.  Four polygons are shown here, including the small polygon of much cloudiness.  A fifth dataset was subsequently incorporated.
Creating polygons for clipping the satellite imagery using QGIS. Four polygons are shown here, including the small polygon of much cloudiness. A fifth dataset was subsequently incorporated.

Finally, I tried to put them together. This turned out to be more trouble than it was worth for my purposes, having only five images. Several of the images had differing resolutions (10 m/pixel for EO-1 ALI, 15 m/pixel for Landsat 8 OLI). Additionally, since I was handling these in their raw format, color balances/exposures were not consistent across images. I decided it best, then, to leave them separate, and sent them around to the Heard Island Expedition team.

Soon I had an email from the expedition leader: he was very interested in the imagery, but it wasn’t opening in Google Earth. Some searching later, I found that Google Earth works best with a certain map projection (EPSG:4326), and when exporting the GeoTIFF, I needed to select “rendered image” rather than “raw data”. I re-exported the images, zipped them up, and tested it out on another computer: success! This Google Earth friendly imagery is now available here (17 MB zip).

One continuation of this project would be to keep looking through the documentation on GeoTiffs to find out how to make the rendered images use a transparent, not white, border where there is no data. That would likely let me create a virtual raster catalog to load all of them in one go, rather than having to load them separately.

Show Me the Data: Satellite Observations

Heard Island, March 27, 2013.  Elephant Spit imagery (at right) is from March 3, 2013.  Image has been adjusted to increase bring out detail in exposed land.  Image credit: processed by Bill Mitchell (CC-BY) using data from NASA/EO-1 (public domain).
Heard Island, March 27, 2013. Elephant Spit imagery (at right) is from March 3, 2013. Image has been adjusted to increase bring out detail in exposed land. Image credit: processed by Bill Mitchell (CC-BY) using data from NASA/USGS/EO-1 ALI (public domain).

With Heard Island being remote and uninhabited, studying it can be a bit difficult. However, as readers of this blog (and my Twitter followers) are aware, one of the ways I have been preparing for the expedition is by keeping an eye on it using various satellites and their remote sensing capabilities. Sure, there is often cloudcover at Heard, but some days it’s clearer and on a few of those days, the satellites pass over.

Most of my information comes from NASA’s MODIS instruments, aboard the Terra and Aqua satellites. These have at least every-other-day coverage of everywhere on Earth, although with a moderate resolution of 250 m/pixel. In the morning when I’m catching up on email and comics, I’ll check the near-real-time MODIS image page to see whether there are clear images of Heard Island from either instrument. Finding Heard Island can be difficult: I still usually find the Kerguelen Islands first, then look to the south-southeast. Many times there are indications such as vortices or gravity waves (not gravitational waves, those are different).

A related page is MODVOLC, which uses MODIS for volcano monitoring. In addition to visible light, MODIS can detect several wavelengths of infrared, and the signature from those wavelengths can be used to determine whether there is a likely volcanic eruption occurring at a given place.

MODIS is a great instrument in that it has daily or every-other-day coverage. However, the 250 m/pixel resolution can be quite limiting. For higher-resolution imagery, I look to the ALI instrument on NASA’s EO-1 satellite. These images are available (free registration required) from EarthExplorer, a data search portal from the USGS. ALI has a 30 m/pixel resolution on its color imagery, and 10 m/pixel resolution on the panchromatic image (total light intensity). These can be combined using QGIS into 10 m/pixel color images. By exploring the EO-1 page I found that members of the public can make requests for image targets! Imaging requests are subject to a bunch of conditions (availability of satellite, >30-day lead time, recommended >3 month window for imaging), but the request and any data generated from fulfillment of the request are free.

How did I come to know about these great resources? It takes time, searching, and some attention to detail. MODIS I learned about as a graduate student, from friends who used data products (not the true-color imagery) in their doctoral research on atmospheric chemistry. I came across EO-1 ALI from searching for images of Heard Island: I found some higher-than-MODIS resolution images from NASA which were good about indicating the source satellite/instrument. Citing image sources is incredibly useful, and I’m always disappointed when images (at least, non-screenshot images) are given without any sort of source information.

MODVOLC I learned about from the Smithsonian’s Global Volcanism Program, which cites the sources of their eruption reports. Information about the source plus a little searching yielded an interesting and useful data source.

Satellite Image Processing Revisited

Heard Island on Nov. 20, 2015, with image processing underway in QGIS.  Image credit: Bill Mitchell (CC-BY) with satellite imagery from USGS (EO-1 satellite, ALI instrument).
Heard Island on Nov. 20, 2015, with image processing underway in QGIS. Image credit: Bill Mitchell (CC-BY) with satellite imagery from USGS (EO-1 satellite, ALI instrument).

Following up on my earlier post about satellite image processing, I am happy to report that I have made progress in being able to process images myself! Through a fortunate combination of search terms, timing, and luck, I managed to come across two key pieces of information that I needed.

First, I found out how to make RGB images from raster data layers, such as different spectral bands on a satellite, fairly easily with QGIS. That was a big step forward from how I had been doing it previously, which was inelegant, inefficient, and only mostly worked. Stacking three layers (one each for red, green, and blue) into a virtual raster catalog was just a few clicks away (Raster | Miscellaneous | Build Virtual Raster (Catalog)).

Encouraged by the success with that project, I continued clicking around and stumbled across some mention of pan-sharpening (also pan sharpening), where a panchromatic (all-color) detector at high resolution is used to enhance the resolution of a colored image (sharpen). Alternately, you can think of it in the complementary way, where lower-resolution color data is added to a high-resolution greyscale image. So thanks to this blog post, I was able to find out what I needed to do to make that happen in QGIS (and Orfeo Toolbox).

Of course, it would be too easy for that to work. I didn’t have the Orfeo Toolbox installed which that needed, and ended up having to compile that from source code.* When the compiler finished and the program was installed, I went to tell QGIS where it was—but a bug in QGIS prevented me from entering the folder location. First, having just installed and compiled stuff, I attempted to get the latest version of QGIS and many of the tools on which QGIS relies. Being unsuccessful in making all of those and some of the compiler configuration software play nicely with each other, I eventually remembered I could get updated packages through apt-get, which gets pre-compiled binary files put out by the maintainers of Debian Linux. That solution worked, I added the folder location, and now I can have my pan-sharpened images.

Here for your viewing pleasure is my first properly pan-sharpened image: Heard Island on Nov. 20, 2015, seen in “true color” by the Advanced Land Imager (ALI) on the EO-1 satellite.** I’m not convinced it’s right, and I think the contrast needs to be brought down a bit, but I think it’s close.

Heard Island in true color on Nov. 20, 2015.  Image processing: Bill Mitchell (CC-BY) using data from USGS/EO-1.
Heard Island in true color on Nov. 20, 2015. Image processing: Bill Mitchell (CC-BY) using data from USGS (EO-1/ALI).

* Knowing how to compile software from source code is a rather handy skill.
** Emily Lakdawalla has written a great explanation of what “true color” means.