This paper presents a cartographic framework based on algorithms of GMT codes for mapping seismically active areas in Venezuela. The data included raster grids from GEBCO, EGM-2008, and vector geological layers from the USGS. The data were iteratively processed in the console of GMT, converted by GDAL, formatted, and mapped for geophysical data visualisation; the QGIS was applied for geological mapping. We analyzed 2000 samples of the earthquake events obtained from the IRIS seismic database with a 25-year time span (1997–2021) in order to map the seismicity. The approach to linking geological, topographic, and geophysical data using GMT scripts aimed to map correlations among the geophysical phenomena, tectonic processes, geological setting, seismicity, and earthquakes. The practical application of the GMT scripts consists in automated mapping for the visualization of geological risks and hazards in the mountainous region of the Venezuelan Andes. The proposed method integrates the approach of GMT scripts with state-of-the-art GIS techniques, which demonstrated its effectiveness as a tool for mapping spatial datasets and rapid data processing in an iterative regime. In this context, using GMT and GIS to find similarities between the regional earthquake distribution and the geological and topographic setting is essential for hazard risk assessment. This study can serve as a basis for predictive seismic analysis in geologically vulnerable regions of Venezuela. In addition to a technical demonstration of GMT algorithms, this study also contributes to geological and geophysical mapping and seismic hazard assessments in South America. We present the full scripts used for mapping in a GitHub repository.