Brian Jackson

Associate Professor of Physics at Boise State University

  • About
  • Research
    • CV
    • Joining the Boise State Planetary Science Research Group
    • Ultra-short-period planet database
    • The Short Period Planets Group — S(u)PerP(i)G
    • Google Scholar Page
    • Code
  • Teaching
  • Press
  • Extracurricular Activities
    • Field Trips
  • Public Outreach
    • Boise State’s Astronomical Observatory

Automated pitted cone identification experiment

Posted by admin on August 24, 2014
Posted in: Data Science.
A field of pitted cones in Utopia Planitia on Mars, as observed by the HiRISE instrument onboard the Mars Reconnaissance Orbiter. Taken from http://beautifulmars.tumblr.com/post/82817530060/field-of-cones-in-utopia-planitia.

A field of pitted cones in Utopia Planitia on Mars, as observed by the HiRISE instrument onboard the Mars Reconnaissance Orbiter. Taken from http://beautifulmars.tumblr.com/post/82817530060/field-of-cones-in-utopia-planitia.

A visit from a planetary sciences colleague this weekend got me excited again about automated feature identification for remote-sensing imagery. This colleague and I talked about her work on Martian pitted cones (examples shown at left). These small hills form when a basalt flow passes over volatile (i.e., water) deposit, heating it and producing a steam explosion that uplifts the flow into a cone shape with a crater at the apex.

As for the vast majority of geomorphological studies, pitted cones on Mars are identified by groups of dedicated researchers, sifting by hand through hundreds or thousands of high-resolution images. If, instead, identification could be automated, it would help realize dramatic savings in person-hours and probably significantly increase the number of known features. It could also mitigate potential observational biases introduced by human image processing.

As an experiment, I applied algorithms from the scikit-image python package to find pitted cones in the example image shown at left. Fortunately, the documentation already provides a good example of circular feature detection.

So modifying the example code, I applied it to a small portion of the example image, and the figure below shows the results. Here I’m just showing the top ten most strongly detected circular features.

Unfortunately, the algorithm is not perfect — some of the obvious cones were not picked out, and others were highlighted more than once. But overall, not terrible, considering it took me about an hour and a half to implement (which included upgrading to Python 3 via the Anaconda package because scikit-image wouldn’t work with the Enthought Python version 2.7 — this upgrade will probably adversely affect my ability to use astroml, by the way).

Possible ways to improve things:

  1. The example code only returns the top two most strongly detected features for each circle radius it tries; this restriction would be simple to remove.
  2. The fact that the pitted cones are cone-shaped means you could require that any putative crater be framed by an appropriately shaped shadow. Regardless of the cone’s actual height (probably unknown anyway), the shadow on the cone should darken as the solar incident angle approaches 180 degrees, modulo any nearby morphological features (as evident in the figure below).
My attempt to automatically identify pitted cones using scikit-image.

My attempt to automatically identify pitted cones using scikit-image.

Posts navigation

← Valsecchi+ (2014) — “From Hot Jupiters to Super-Earths Via Roche Lobe Overflow”
Vanderburg & Johnson — “A Technique for Extracting Highly Precise Photometry for the Two-Wheeled Kepler Mission” →
  • Twitter: decaelus

    Brian Jackson
    • RT @ChristinaLords: Hi friends. I’m beyond excited to say I’m starting a new nonprofit online journalism project focused on Idaho state pol… about 27 minutes ago
    • @MarkTurnerNC @ThePlanetaryGuy https://t.co/VojCLI4TGa about 44 minutes ago in reply to MarkTurnerNC
    • @SenJohnThune Exactly! What are people complaining about? Senators in the 1860s made $3k per year and got by just f… https://t.co/yvq8ZNKI4d about 45 minutes ago in reply to SenJohnThune
    @decaelus
  • Recent Posts

    • Third Thursday Planetarium Show – 2021 Feb 18
    • First Friday Astronomy – Psyche: Journey to a Metallic World – 2021 Mar 5
    • Borah High School Career Fair – 2021 Jan 26
    • Third Thursday Planetarium Show – 2021 Jan 21
    • First Friday Astronomy – How To Talk So Science Journalists Will Listen – 2021 Feb 5
  • Archives

    • February 2021
    • January 2021
    • December 2020
    • November 2020
    • October 2020
    • September 2020
    • August 2020
    • July 2020
    • June 2020
    • May 2020
    • April 2020
    • March 2020
    • February 2020
    • January 2020
    • December 2019
    • November 2019
    • October 2019
    • September 2019
    • August 2019
    • July 2019
    • June 2019
    • May 2019
    • April 2019
    • March 2019
    • February 2019
    • January 2019
    • December 2018
    • November 2018
    • October 2018
    • September 2018
    • August 2018
    • July 2018
    • June 2018
    • May 2018
    • April 2018
    • March 2018
    • February 2018
    • January 2018
    • December 2017
    • November 2017
    • October 2017
    • September 2017
    • August 2017
    • July 2017
    • June 2017
    • May 2017
    • April 2017
    • March 2017
    • February 2017
    • January 2017
    • December 2016
    • November 2016
    • October 2016
    • September 2016
    • August 2016
    • July 2016
    • June 2016
    • May 2016
    • April 2016
    • March 2016
    • February 2016
    • January 2016
    • December 2015
    • November 2015
    • October 2015
    • September 2015
    • August 2015
    • July 2015
    • June 2015
    • May 2015
    • April 2015
    • March 2015
    • February 2015
    • January 2015
    • December 2014
    • November 2014
    • October 2014
    • September 2014
    • August 2014
    • July 2014
    • June 2014
    • May 2014
    • April 2014
    • March 2014
    • February 2014
    • January 2014
    • December 2013
    • November 2013
    • October 2013
    • September 2013
    • August 2013
    • July 2013
Proudly powered by WordPress Theme: Parament by Automattic.