Python-based analysis of geospatial data

Replication of results described in ¹.

conda create -c conda-forge -n kossou-hillshade earthpy matplotlib misaka numpy rasterio
pip install kossou-hillshade-cli

kossou-hillshade --colormap cividis Figure_3_Hillshade_Kossou_1609.jpg
kossou-hillshade --azimuth 230 Figure_4_Azimuth_hillshade_Kossou.jpg
kossou-hillshade --altitude 10 --colormap magma Figure_5_Angle_hillshade_Kossou.jpg

Download elevation file

To download the elevation file, first register an account on the USGS EROS registration system to get download access via EarthExplorer and/or GloVis.

EarthExplorer

Under the Datasets tab, under the Digital Elevation category, under the SRTM subcategory, select SRTM Void Filled.

Under the Additional Criteria tab, set Entity ID to SRTM3N07W006V2.

The results tab will then contain the SRTM Void Filled DEM entity with ID SRTM3N07W006V2, as described in ¹.

You can select the Download Options icon and download the GeoTIFF 3 Arc-second file, n07_w006_3arc_v2.tif.

GloVis

Under the Interface Controls, click the "+" next to "Selected Data Set(s) and add the SRTM Void Filled data set.

Add a Dataset Metadata Filter for the SRTM Void Filled dataset with Entity ID SRTM3N07W006V2, and apply the filter.

You may need to scroll the map view to west Africa in order to find the result. The scene will appear in the Scene Navigator, where you can view the download options and download the GeoTIFF 3 Arc-second file n07_w006_3arc_v2.tif.

Download landsat files

Similar to above, use EarthExplorer or GloVis to search the Landsat 8-9 OLI/TIRS C2 L1 dataset for a result with Landsat Product Identifier L1 of LC09_L1TP_197055_20220111_20230502_02_T1.

Generate plots

The results of ¹ may be replicated by running:

kossou-hillshade --geotiff n07_w006_3arc_v2.tif --colormap cividis Figure_3_Hillshade_Kossou_1609.jpg
kossou-hillshade --geotiff n07_w006_3arc_v2.tif --azimuth 230 --colormap plasma Figure_4_Azimuth_hillshade_Kossou.jpg
kossou-hillshade --geotiff n07_w006_3arc_v2.tif --altitude 10 --colormap magma Figure_5_Angle_hillshade_Kossou.jpg

In fact, the geotiff file n07_w006_3arc_v2.tif is included in the kossou-hillshade-base package, so this will also work:

kossou-hillshade --colormap cividis Figure_3_Hillshade_Kossou_1609.jpg
kossou-hillshade --azimuth 230 --colormap plasma Figure_4_Azimuth_hillshade_Kossou.jpg
kossou-hillshade --altitude 10 --colormap magma Figure_5_Angle_hillshade_Kossou.jpg