Geoscience Australia’s Historical Aerial Photography Program currently involves scanning and georeferencing old flight diagrams to enable the digitising and positioning of historical aerial photographs for easy discovery and download. Accurate digital mapping of GA’s aerial photography collection will make catalogue searches easier and the collection more accessible to the public. This story map presents an interactive history of aerial photography, a background of aerial photography in Australia, historical aerial photography use cases and scenarios, and a background on Geoscience Australia's program to digitise flight diagrams and create a catalogue of aerial photographs.

A `weighted geometric median’ approach has been used to estimate the median surface reflectance of the barest state (i.e., least vegetation) observed through Landsat-8 OLI observations from 2013 to September 2018 to generate a six-band Landsat-8 Barest Earth pixel composite mosaic over the Australian continent. The bands include BLUE (0.452 - 0.512), GREEN (0.533 - 0.590), RED, (0.636 - 0.673) NIR (0.851 - 0.879), SWIR1 (1.566 - 1.651) and SWIR2 (2.107 - 2.294) wavelength regions. The weighted median approach is robust to outliers (such as cloud, shadows, saturation, corrupted pixels) and also maintains the relationship between all the spectral wavelengths in the spectra observed through time. The product reduces the influence of vegetation and allows for more direct mapping of soil and rock mineralogy. Reference: Dale Roberts, John Wilford, and Omar Ghattas (2018). Revealing the Australian Continent at its Barest, submitted.

An estimate of the spectra of the barest state (i.e., least vegetation) observed from imagery of the Australian continent collected by the Landsat 5, 7, and 8 satellites over a period of more than 30 years (1983 – 2018). The bands include BLUE (0.452 - 0.512), GREEN (0.533 - 0.590), RED, (0.636 - 0.673) NIR (0.851 - 0.879), SWIR1 (1.566 - 1.651) and SWIR2 (2.107 - 2.294) wavelength regions. The approach is robust to outliers (such as cloud, shadows, saturation, corrupted pixels) and also maintains the relationship between all the spectral wavelengths in the spectra observed through time. The product reduces the influence of vegetation and allows for more direct mapping of soil and rock mineralogy. This product complements the Landsat-8 Barest Earth which is based on the same algorithm but just uses Landsat8 satellite imagery from 2013-2108. Landsat-8’s OLI sensor provides improved signal-to-noise radiometric (SNR) performance quantised over a 12-bit dynamic range compared to the 8-bit dynamic range of Landsat-5 and Landsat-7 data. However the Landsat 30+ Barest Earth has a greater capacity to find the barest ground due to the greater temporal depth. Reference: Exposed Soil and Mineral Map of the Australian Continent Revealing the Land at its Barest - Dale Roberts, John Wilford and Omar Ghattas Ghattas (2019). Nature Communications, DOI: 10.1038/s41467-019-13276-1. https://www.nature.com/articles/s41467-019-13276-1

The Sentinel-2 Bare Earth thematic product provides the first national scale mosaic of the Australian continent to support improved mapping of soil and geology. The bare earth algorithm using all available Sentinel-2 A and Sentinel-2 B observations up to September 2020 preferentially weights bare pixels through time to significantly reduce the effect of seasonal vegetation in the imagery. The result are image pixels that are more likely to reflect the mineralogy and/or geochemistry of soil and bedrock. The algorithm uses a high-dimensional weighted geometric median approach that maintains the spectral relationships across all Sentinel-2 bands. A similar bare earth algorithm has been applied to Geoscience Australia’s deeper Landsat time series archive (please search for "Landsat barest Earth". Both bare earth products have spectral bands in the visible near infrared and shortwave infrared region of the electromagnetic spectrum. However, the main visible and near-infrared Sentinel-2 bands have a spatial resolution of 10 meters compared to 30m for the Landsat TM equivalents. The weighted median approach is robust to outliers (such as cloud, shadows, saturation, corrupted pixels) and also maintains the relationship between all the spectral wavelengths in the spectra observed through time. Not all the sentinel-2 bands have been processed - we have excluded atmospheric bands including 1, 9 and 10. The remaining bands have been re-number 1-10 and these bands correlate to the original bands in brackets below: 1 = blue (2) , 2 = green (3) , 3 = red (4), 4 = vegetation red edge (5), 5 = vegetation red edge (6), 6= vegetation red edge (7), 7 = NIR(8), 8 = Narrow NIR (8a), 9 = SWIR1 (11) and 10 = SWIR2(12). All 10 bands have been resampled to 10 meters to facilitate band integration and use in machine learning.