Console-Based Mapping of Mongolia Using GMT Cartographic Scripting Toolset for Processing TerraClimate Data

Satellite Image Processing by Python and R Using Landsat 9 OLI/TIRS and SRTM DEM Data on Côte d'Ivoire, West Africa

In this paper, we propose an advanced scripting approach using Python and R for satellite image processing and modelling terrain in Côte d'Ivoire, West Africa. Data include Landsat 9 OLI/TIRS C2 L1 and the SRTM digital elevation model (DEM). The EarthPy library of Python and 'raster' and 'terra' packages of R are used as tools for data processing. The methodology includes computing vegetation indices to derive information on vegetation coverage and terrain modelling. Four vegetation indices were computed and visualised using R: the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index 2 (EVI2), Soil-Adjusted Vegetation Index (SAVI) and Atmospherically Resistant Vegetation Index 2 (ARVI2). The SAVI index is demonstrated to be more suitable and better adjusted to the vegetation analysis, which is beneficial for agricultural monitoring in Côte d'Ivoire. The terrain analysis is performed using Python and includes slope, aspect, hillshade and relief modelling with changed parameters for the sun azimuth and angle. The vegetation pattern in Côte d'Ivoire is heterogeneous, which reflects the complexity of the terrain structure. Therefore, the terrain and vegetation data modelling is aimed at the analysis of the relationship between the regional topography and environmental setting in the study area. The upscaled mapping is performed as regional environmental analysis of the Yamoussoukro surroundings and local topographic modelling of the Kossou Lake. The algorithms of the data processing include image resampling, band composition, statistical analysis and map algebra used for calculation of the vegetation indices in Côte d'Ivoire. This study demonstrates the effective application of the advanced programming algorithms in Python and R for satellite image processing.

Handling Dataset with Geophysical and Geological Variables on the Bolivian Andes by the GMT Scripts

In this paper, an integrated mapping of the georeferenced data is presented using the QGIS and GMT scripting tool set. The study area encompasses the Bolivian Andes, South America, notable for complex geophysical and geological parameters and high seismicity. A data integration was performed for a detailed analysis of the geophysical and geological setting. The data included the raster and vector datasets captured from the open sources: the IRIS seismic data (2015 to 2021), geophysical data from satellite-derived gravity grids based on CryoSat, topographic GEBCO data, geoid undulation data from EGM-2008, and geological georeferences’ vector data from the USGS. The techniques of data processing included quantitative and qualitative evaluation of the seismicity and geophysical setting in Bolivia. The result includes a series of thematic maps on the Bolivian Andes. Based on the data analysis, the western region was identified as the most seismically endangered area in Bolivia with a high risk of earthquake hazards in Cordillera Occidental, followed by Altiplano and Cordillera Real. The earthquake magnitude here ranges from 1.8 to 7.6. The data analysis shows a tight correlation between the gravity, geophysics, and topography in the Bolivian Andes. The cartographic scripts used for processing data in GMT are available in the author’s public GitHub repository in open-access with the provided link. The utility of scripting cartographic techniques for geophysical and topographic data processing combined with GIS spatial evaluation of the geological data supported automated mapping, which has applicability for risk assessment and geological hazard mapping of the Bolivian Andes, South America.