GDAL and PROJ Libraries Integrated with GRASS GIS for Terrain Modelling of the Georeferenced Raster Image

Libraries with pre-written codes optimize the workflow in cartography and reduce labour intensive data processing by iteratively applying scripts to implementing mapping tasks. Most existing Geographic Information System (GIS) approaches are based on traditional software with a graphical user's interface which significantly limits their performance. Although plugins are proposed to improve the functionality of many GIS programs, they are usually ad hoc in finding specific mapping solutions, e.g., cartographic projections and data conversion. We address this limitation by applying the principled approach of Geospatial Data Abstraction Library (GDAL), library for conversions between cartographic projections (PROJ) and Geographic Resources Analysis Support System (GRASS) GIS for geospatial data processing and morphometric analysis. This research presents topographic analysis of the dataset using scripting methods which include several tools: (1) GDAL, a translator library for raster and vector geospatial data formats used for converting Earth Global Relief Model (ETOPO1) GeoTIFF in XY Cartesian coordinates into World Geodetic System 1984 (WGS84) by the 'gdalwarp' utility; (2) PROJ projection transformation library used for converting ETOPO1 WGS84 grid to cartographic projections (Cassini-Soldner equirectangular, Equal Area Cylindrical, Two-Point Equidistant Azimuthal, and Oblique Mercator); and (3) GRASS GIS by sequential use of the following modules: r.info, d.mon, d.rast, r.colors, d.rast.leg, d.legend, d.northarrow, d.grid, d.text, g.region, and r.contour. The depth frequency was analysed by the module 'd.histogram'. The proposed approach provided a systematic way for morphometric measuring of topographic data and combine the advantages of the GDAL, PROJ, and GRASS GIS tools that include the informativeness, effectiveness, and representativeness in spatial data processing. The morphometric analysis included the computed slope, aspect, profile, and tangential curvature of the study area. The data analysis revealed the distribution pattern in topographic data: 24% of data with elevations below 400 m, 13% of data with depths −5000 to −6000 m, 4% of depths have values −3000 to −4000 m, the least frequent data (−6000 to 7000 m)