Search

International audience; Physically accurate spectral rendering generally brings further computationalcomplexity to RGB rendering pipelines which precludes real-time rates. In thisreport, we show that by building on several physical a-prioris for the spectral phe-nomena in standard astronomical scenes, we achieve to obtain a model greatlyreducing the complexity of the spectral aspect for volumetric rendering - by up to83 factor in our examples. By projecting the total local contribution on a small op-timal basis, thus providing scalability for the method, our rendering scheme thenmakes it possible to render a volumetric scene with any number and spectra oflights in real-time and no significant visual discrepancies.

This paper focuses on the environment of Ethiopia, a country highly sensitive to droughts severely affecting vegetation. Vegetation monitoring of Ethiopian Highlands requires visualization of environmental parameters to assess droughts negatively influencing agricultural sustainable management of crops. Therefore, this study presented mapping of several climate and environmental variables including Palmer Drought Severity Index (PDSI). The data were visualized and interpreted alongside the topographic data to evaluate the environmental conditions for vegetation. The datasets included WorldClim and GEBCO and Digital Chart of the World (DCW). Research has threefold objectives: i) environmental mapping; ii) technical cartographic scripting; iii) data processing. Following variables were visualized on seven new maps: 1) topography; 2) soil moisture; 3) T °C minimum; 4) T °C maximum; 5) Wind speed; 6) Precipitation; 7) Palmer Drought Severity Index (PDSI). New high-resolution thematic environmental maps are presented and the utility of GMT for mapping multi-source datasets is described. With varying degrees of soil moisture (mean value of 15.0), min T°C (−1.8°C to 24°C), max T°C (14.4°C to 40.2°C) and wind speed (0.1 to 6.1 m/s), the maps demonstrate the variability of the PDSI fields over the country area (from −11.7 to 2.3) induced by the complex sum of these variables and intensified by the topographic effects notable over the Ethiopian Highlands which can be used for vegetation analysis. The paper presents seven new maps and contributes to the environmental studies of Ethiopia.

International audience

Summary The study presents a comparative analysis of eight Vegetation Indices (VIs) used to examine vegetation greenness over the northern coasts of Iceland. The geographical extent of the study area is set by the coordinates of the two fjords, Eyjafjörður and Skagafjörður, notable for their agricultural significance. Vegetation in Iceland is fragile due to the harsh climate, climate change, overgrazing and volcanic activity, which increase soil erosion. The study was conducted on a Landsat TM image using SAGA GIS as a technical tool for raster bands calculations. The NDVI dataset shows a range from -0.56 to 0.24, with 0 indicating ‘no vegetation’, and negative values – ‘other surfaces’ (e.g. rocks, open terrain). The DVI, compared to the NDVI, shows statistically non-normalized values ranging from -112 to 0, with extreme negative values while the coastal vegetation areas are badly distinguished from the water areas. The NRVI shows an extent from -0.24 to 0.48 with higher values for vegetation. The NRVI reduces topographic, solar and atmospheric effects and creates a normal data distribution. RVI shows a range in a dataset from 0.2 to 3.2 with vegetation in the river valleys clearly visible and depicted, while the water areas have values 0.8 to 1.0. The CTVI shows corrected TVI, in a data range -0.10 to 1.10, as the dataset of NDVI were negative. The TVI dataset ranges from 0.44 to 0.80 with the ice-covered areas and glaciers distinguishable and water values within a range from 0.60 to 0.64 and the vegetation from 0.60 to 0.44. The TTVI dataset ranges from 0.40 to 0.80 performing similarly to the TVI, but more refined with vegetation values 0.64 to 0.68. SAVI dataset ranges from -0.80 to 0.30 with minimized effects of soil on the vegetation through a constant soil adjustment factor added into the NDVI formula. The paper presents a comparison of eight VIs for Arctic vegetation monitoring. The overall behavior of SAGA GIS in calculation and mapping of the VIs is effective in terms of their use for vegetation mapping of the region.

Geologic and geophysical mapping has been so far limited to the traditional single-method GIS-based mapping. A new approach combining integrated analysis of data on geology, gravity, topography and geomorphology is presented for regional characterization of the geophysical setting in Mongolia: the Gobi Altai Mountains, the Khangai Mountains and Khentii Mountains with surrounding areas. Nine new maps have been produced from the high-resolution datasets: GEBCO, gravity raster, USGS geological data and SRTM-90 DEM geomorphological grid. Methodology includes three tools for cartographic data visualization: i) Generic Mapping Tools (GMT), ii) R programming language (‘raster’ and ‘tmap’ libraries); iii) QGIS. The results demonstrated strong agreement between the estimated values in gravity and topography grids, distribution of geological units and provinces over the country and geomorphological landforms with respect to the mountain ranges: Altai, Khangai and Khentii Mountains. The highest values in the gravity anomalies correspond to the mountain ranges in the Altai Mountains and Khangai Mountains (

International audience; Space has fascinated and inspired people since almost the beginning of the human era. More recently, enormous and complex structures like galaxies, nebulae or gaz clouds became observable by humans, using fairly advanced technology. This has only lead to more attraction and many artistic representations of them were created. Even more recently began the computer era, and with it, computer graphics. Many applications of computer graphics were found, ranging from video games to movies and scientific simulations. However, the complexity of these space structures makes them really hard to model in computer graphics. In this paper, we present a way of generating and visualizing these kind of structures, efficiently enough to be done in real-time.

International audience; This paper examines spatial variations in the geomorphology of the Ninety East Ridge (NER), located in the Indian Ocean. The NER is an extraordinary long linear bathymetric feature with topography reflecting complex geophysical setting and geologic evolution. The research is based on a compilation of high-resolution bathymetric, geological, and gravity datasets clipped for the study area extent (65° - 107°E, 35°S - 21°N): General Bathymetric Chart of the Oceans (GEBCO), Earth Gravitational Model (EGM2008, EGM96). The submarine geomorphology of the NER was modeled by digitized cross-sectional profiles using Generic Mapping Tools (GMT). The availability of the method is explained by 1) the free datasets; 2) the open source GMT toolset; 3) the available tutorials of the GMT and the codes explained in this work. Three segments of the NER were selected, digitized, and modeled: 1) northern 89°E, 7°S to 90°E, 7°N; 2) central 88.4°E, 14.7°S to 88.8°E, 8.2°S; 3) southern 87.9°E, 17°S to 87.5°E, 27°S. Measured depths were visualized in graphs, compared, and statistically analyzed by the histograms. The northern segment has a steepness of 21.3° at the western slopes, and 14.5° at the eastern slope. The slopes on the eastern flank have dominant SE orientation. The central segment has a bell-shaped form, with the highest steepness comparing to the northern and southern segments. The eastern flank has a steepness of 49.5°. A local depression at a distance of 50 km off from the axis (90°E) continues parallel to the NER, with the shape of the narrow minor trench. The western slope has a steepness of 57.6°, decreasing to 15.6°. The southern segment has a dome-like shape form. Compared to the northern and central segments, it has a less pronounced ridge crest, with a steepness of 24.9° on the west. The eastern flank has a steepness of 36.8° until 70 km, gradually becoming steeper at 44.23°. A local minor trench structure can be seen on its eastern flank (100 km off the axis). This corresponds to the very narrow long topographic depressions stretching parallel to this segment of the NER at 90.5°E. The study contributes to regional geographic studies of Indian Ocean geomorphology and cartographic presentation of GMT functionality for marine research and oceanographic studies.

International audience; This paper explores the geomorphological features of the Ninety East Ridge (NER), Indian Ocean. NER presents especially long and linear topographic structure formed as a result of complex regional geophysical and geologic development. The research is based on high precision bathymetric, geological and gravity data. The submarine geomorphology of NER was digitized as three cross-sectional profiles. The profiles were selected in northern, central and southern segments. The depths were visualized in graphs, compared and statistically analyzed by histograms. The study contributes to the geophysical studies of the Indian Ocean.; В данной статье исследуются пространственные вариации геоморфологии хребта Девяностого градуса, также известного как Восточно-Индийский хребет (ВИХ) в Индийском океане в трех сегментах. ВИХ представляет собой необычайно протяженный линейный батиметрический объект с топографией, отражающей сложные геофизические условия и геологическую эволюцию. Исследование основано на компиляции наборов батиметрических, геологических и гравиметрических данных высокого разрешения. Подводная геоморфология была смоделирована с помощью оцифрованных профилей поперечного сечения. Были выбраны три сегмента по геоморфологическим особенностям ВИХ: 1) северный; 2) центральный; 3) южный. Измеренные глубины были визуализированы в виде графиков, проведен сравнительный анализ статистических гистограмм распределения и повторяемости глубин. Работа вносит вклад в региональные исследования геоморфологии дна Индийского океана.

Mapping seafloor geomorphology is a complex task requiring the integration of advanced cartographic technologies and high-resolution spatial data. This paper provides a comparative geomorphological analysis of the Hellenic Trench (HT) and the Pliny Trench (PT) located in the Eastern Mediterranean Sea, southern Greece. These trenches were formed as a result of the tectonic plate subduction in the Eastern Mediterranean Sea: the northward motion of the African and Arabian plates, complicated by the regional geological settings, such as active faults and earthquakes, which resulted in their different geomorphological forms and bathymetric shapes. Technically, this paper presents an example of the advanced scripting mapping by scripting the cartographic toolset of Generic Mapping Tools (GMT), which performs mapping through shell scripts. The maps are based on the high-quality topographic, geological and geophysical data: GEBCO, EGM96, geoid, and marine free-air gravity fields. The GMT builds upon the modules used for data processing. The region was subsetted by ‘grdcut’, analysed by the Geospatial Data Abstraction Library (GDAL) (gdalinfo utility), and visualized by ‘grdimage’. Two segments of the trenches formed in a condition of varying geological and geophysical settings, have been modelled, visualized and compared, as representative cross-sections. As a result of the automated digitizing, spatial interpolation and sequential aggregating of GMT codes, the segments of the cross-sections were represented. The HT (Ionian segment) has an asymmetric one-sided shape; a steepness of 56.8° on the NE side and 7° on the SW flank. The PT has a more symmetric view; a steepness of 42.14° on its NW flank and 26.66° on its SE flank. The PT has a clear peak of the depths at -2600 to -2800 m and the most representative data range at -5000 to -6000 m. The HT has a bimodal data distribution with two peaks. The most frequent data for HT is -3200 m to -3400 m. Compared to PT, the HT is deeper. The histogram shows the peak of data for HT in the interval between -3200 to -3400 m (135 samples) while the PT shows the peak of data in the interval at -2600 to -2800 m (310 samples). Besides, 105 samples of the HT have depths exceeding 4000 m, while only 20 samples were recorded for PT in the same interval. This paper contributes to the geomorphological studies of the general Eastern Mediterranean Sea region, particularly relating to regional seafloor mapping of the Hellenic and Pliny trenches., Mapiranje geomorfologije morskoga dna složen je zadatak koji zahtijeva integraciju naprednih kartografskih tehnologija i prostornih podataka visoke rezolucije. Ovaj rad donosi usporednu geomorfološku analizu Helenskoga i Plinijskoga jarka smještenih u istočnome Sredozemlju, na jugu Grčke. Te su strukture nastale kao rezultat podvlačenja tektonske ploče, tj. kretanja afričke i arapske ploče prema sjeveru, usložnjenoga regionalnim geološkim postavkama poput aktivnih rasjeda i potresa, što je rezultiralo njihovim različitim geomorfološkim oblicima i batimetrijom. Prikazano je njihovo napredno kartiranje primjenom alata Generic Mapping Tools (GMT), koji koriste programske skripte. Karte se temelje na visokokvalitetnim topografskim, geološkim i geofizičkim podatcima, kao što su GEBCO, EGM96, geoid te gravitacijsko polje na razini mora. GMT se nadovezuje na module koji se koriste za obradu podataka. Područje je podijeljeno skriptom grdcut, a na temelju podataka iz zbirke apstraktnih geoprostornih podataka (GDAL) (skripta gdalinfo), te je na kraju vizualizirano (skripta grdimage). Dva dijela jarka oblikovana su u različitim geološkim i geofizičkim postavkama te su uspoređeni njihovi reprezentativni presjeci. Kao rezultat automatizirane digitalizacije, prostorne interpolacije i sekvencijalnoga agregiranja GMT kodova predstavljeni su segmenti presjeka. Helenski jarak (jonski segment) ima asimetričan jednostrani oblik, nagiba 56,8° na sjeveroistočnoj i 7° na jugozapadnoj strani. Plinijski je jarak simetričniji, nagiba 42,14° na SZ i 26,66° na JI boku. Plinijski jarak ima jasno izraženo tjeme od –2600 do –2800 m, a najreprezentativniji su podatci na –5000 do –6000 m. Helenski jarak ima bimodalnu distribuciju podataka s dvama tjemenima, a najčešće su dubine od –3200 do –3400 m. Taj je jarak dublji. Histogram prikazuje okupljanje podataka od –3200 do –3400 m, dok Plinijski (135 uzoraka) prikazuje maksimalne vrijednosti od –2600 do –2800 m (310 uzoraka). Osim toga, 105 uzoraka Helenskoga jarka ima dubine veće od –4000 m, dok je u istome intervalu za Plinijski zabilježeno samo 20 uzoraka. Članak doprinosi geomorfološkim studijama regije istočnoga Sredozemnog mora općenito i regionalnomu mapiranju morskoga dna, posebno Helenskoga i Plinijskoga jarka.

International audience; Landsat-TM of 2001 covering Iceland (15.5°W-21°W, 64.5°N-67°N) was processed using SAGA GIS for testing distance-based Vegetation Indices (VIs): four approaches of Perpendicular Vegetation Index (PVI) and two approaches of Transformed Soil Adjusted Vegetation Index TSAVI. The PVI of vegetation from the soil background line indicated healthiness as a leaf area index (LAI). The results showed that the reflectance for vegetation has a linear relation with soil background line. Four PVI models and two TSAVI shown coefficients of determination with LAI. The dataset demonstrate variations in the calculated coefficients. The mode in the histograms of the PVI based on four different algorithms show the difference:-7.1,-8.36, 2.78 and 7.0. The dataset for the two approaches of TSAVI: first case ranges in 4.4.-80.6 with a bellshape mode of a histogram (8.09 to 23.29) for the first algorithm and an irregular shape for the second algorithm with several modes starting from 0.11 to 0.2 and decreasing to 0.26. SAGA GIS permits the calculation of PVI and TSAVI by computed NDVI based on the intersection of vegetation and soil background. Masking the NIR and R, a linear regression of grids was performed using an equation embedded in SAGA GIS. The advantages of the distance-based PVI and TSAVI consists in the adjusted position of pixels on the soil brightness line which refines it comparing to the slope-based VIs. The paper demonstrates SAGA GIS application in agricultural studies.

International audience; The study focused on the comparative analysis of the submarine geomorphology of three oceanic trenches: Hikurangi Trench (HkT), Puysegur Trench (PT) and Hjort Trench (HjT), New Zealand region, Pacific Ocean. HjT is characterized by an oblique subduction zone. Unique regional tectonic setting consist in two subduction zones: northern (Hikurangi margin) and southern (Puysegur margin), connected by oblique continental collision along the Alpine Fault, South Island. This cause variations in the geomorphic structure of the trenches. PT/HjT subduction is highly oblique (dextral) and directed southwards. Hikurangi subduction is directed northwestwards. South Island is caught in between by the “subduction scissor”. Methodology is based on GMT (The Generic Mapping Tools) for mapping, plotting and modelling. Mapping includes visualized geophysical, tectonic and geological settings of the trenches, based on sequential use of GMT modules. Data include GEBCO, ETOPO1, EGM96. Comparative histogram equalization of topographic grids (equalized, normalized, quadratic) was done by module ’grdhisteq’, automated cross-sectioning – by ’grdtrack’.Results shown that HjT has a symmetric shape form with comparative gradients on both western and eastern slopes. HkT has a trough-like flat wide bottom, steeper gradient slope on the North Island flank. PT has an asymmetric V-form with steep gradient on the eastern slopes and gentler western slope corresponding to the relatively gentle slope of a subducting plate and steeper slope of an upper one. HkT has shallower depths < 2,500 m, PT is 6,000 m for HjT. The surrounding relief of the HjT presents the most uneven terrain with gentle slope oceanward, and a steep slope on the eastern flank for PT, surrounded by complex submarine relief along the Macquarie Arc. Data distribution for the HkT demonstrates almost equal pattern for the depths from −600 m to −2,600 m. PT has a bimodal data distribution with 2 peaks: 1) −4,250 to −4,500 m (18%); 2) −2,250 to −3,000 m, < 7,5%. The second peak corresponds to the Macquarie Arc. Data distribution for HjT is classic bell-shaped with a clear peak at −3,250 to −3,500 m. The asymmetry of the trenches resulted in geomorphic shape of HkT, PT and HjT affected by geologic processes.; Studium poświęcone jest analizie porównawczej rzeźby dna trzech rowów oceanicznych: Hikurangi (HkT), Puysegur (PT) i Hjort (HjT), położonych w pobliżu Nowej Zelandii na południowym Pacyfiku. HjT charakteryzuje się skośną strefą subdukcji. Unikalna sy- tuacja geotektoniczna regionu polega na rozdzieleniu dwóch stref subdukcji: północnej (Hikurangi) i południowej (Puysegur), strefą kolizji kontynentalnej wzdłuż uskoku Alpine Fault na Wyspie Południowej. Subdukcja na południe od Wyspy Południowej zachodzi pod dużym kątem w kierunku południowo-wschodnim (PT i HjT), podczas gdy w strefie północnej (Hikurangi) odbywa się na pół- nocny zachód. W konsekwencji Wyspa Południowa jest ujęta w swego rodzaju „nożyce subdukcyjne”. Metodologia oparta na GMT (The Generic Mapping Tools) posłużyła do skartowania, wykreślenia i modelowania obszaru. Kartowanie obejmuje wizualizację danych geofizycznych oraz pozycji tektonicznej i geologicznej rowów, opartą na sekwencyjnym użyciu modułów GMT. Dane obej- mują GEBCO, ETOPO1, EGM96. Porównawcza korekcja histogramu siatek topograficznych (wyrównana, znormalizowana, kwadra- towa) została wykonana przez moduł „grdhisteq”, zaś zautomatyzowane przekroje – przez moduł „grdtrack”.Analiza wykazała , że rów Hjort ma symetryczną formę z porównywalnymi nachyleniami zarówno na zachodnich, jak i wschod- nich zboczach. Rów Hikurangi ma podobne do koryta płaskie szerokie dno, a stok od strony zachodniej (przylegający do Wyspy Pół- nocnej) jest nachylony pod większym kątem od stoku wschodniego. Rów Puysegur ma asymetryczną V-kształtną formę ze stromo na- chylonym zboczem wschodnim i łagodniejszym zachodnim. Rów HkT jest relatywnie płytki < 2500 m, PT osiąga głębokość 6000 m) stwierdzono dla rowu Hjort. Rzeźba dna w otoczeniu HjT jest najbardziej zróżnicowana, a w przypadku położonego bardziej na północ PT zaznacza się wyraźna dysproporcja pomiędzy łagodnym oceanicznym zboczem na zachodzie i stromym zboczem grzbietu Puysegur (północny odcinek Łuku Macquarie) na wschodniej flance rowu. Rozkład da- nych batymetrycznych dla HkT jest stosunkowo zrównoważony dla głębokości od 600 m do 2600 m. PT ma bimodalny rozkład danych z 2 pikami: 1) 4250 do 4500 m (18%); 2) 2250 do 3000 m, < 7,5%. Druga koncentracja danych odpowiada łukowi Macquarie. Rozkład danych dla HjT ma klasyczny kształt dzwonu z wyraźnym ekstremum odpowiadającym głębokościom 3250 do 3500 m. Asymetria zaprezentowanych rowów oceanicznych jest uwarunkowana przez procesy geotektoniczne.

The study presents an analysis of the sediment thickness compared with bathymetric and geophysical settings in the Bay of Bengal and Andaman Sea, Indian Ocean. It uses a combination of the high-resolution data: topographic GEBCO, satellite and marine gravity anomalies, EGM2008 geoid and GlobSed to visualize the correlation between relief, gravity and trends in continent-ocean sediment transport. The results include thematic maps and 3D model showing increased sediment thickness in the Bengal Fan (8,0 to 8,2 km) in NE direction with maximum in Ganges Fan (16,2 km), and southward decrease in the Andaman Sea from Irrawaddy river mouth (6-7 km) to the Strait of Malacca (1-2 km). All maps and 3D model have been plotted by Generic Mapping Tools (GMT) cartographic scripting toolset version 6.0.0.

International audience; The study presents the insights of the tectonic development and geological settings of the Atlantic Ocean supported by cartographic visualization in Generic Mapping Tools (GMT). The aim is to study geologic situation and trends in the tectonic development of the Mid-Atlantic Ridge and Atlantic Ocean seafloor. The objective is to find out impact of various factors (such as volcanic, tectonic, hydrothermal and sedimentary processes) that sculpt seafloor geomorphology, and correlation between early history of crust formation, geological processes and present submarine landforms. Other assignments in this work refer to mutual comparison of raster grids on sedimentation, topography, geology, seafloor fabric and highlighting similarities among the landforms and sediment thickness. Asymmetry in crustal accretion is explained by the tectonic history of the lithosphere formation. Correlation between plate subduction and development of the submarine landforms is explained by the Earth's crust extension resulting in formation of cracks, elongations, faults, rifts. Ocean seafloor geomorphology is shaped by a variety of factors that impact its form at different scales. These drivers (tectonic evolution, oceanic currents, hydrology, sedimentation) have effects on geomorphic landforms of the seafloor in context of historic geological development and during Quaternary. Technical part of this work was performed by GMT scripting toolset with all maps plotted in American polyconic projection. The results are received by overlay, cartographic analysis and synthesis of the multi-source geodata through mapping and interpreting grids (ETOPO1, EGM96, GlobSed, crustal age). This work contributes to expand the knowledge on geological and tectonic development of the Atlantic Ocean seabed in order to complete the view of its submarine geomorphology

Cartographic visualisation is a central concept in geoinformatics, and Generic Mapping Tools (GMT) functionality provides a variety of modules for effective mapping. However, due to its specific scripting approach, there is not enough reported experience of GMT mapping, comparing to traditional GIS. This contribution introduces steps that can be taken to perform cartographic mapping and modelling using GMT. Geographically, this paper investigates the Izu-Bonin Trench in the Pacific Ocean. The aim was to compare its geomorphology in two segments, and each was modelled by a series of profiles. The comparative analysis shows that the southern segment is deeper and is a more precisely V-shaped form with a steeper gradient slope. The northern part has an asymmetric slope with submarine terraces to the west and a straight shape to the east. The northern profile is based on 407 samples (13.5%) at depths of −5,600 to −5,800 m, followed by 304 samples at −5,800 to −6,000 m (10%). The southern histogram has a bimodal distribution with two peaks: 523 samples (20%) with depths of −5,800 to −6,000 m. The second peak (10%) is on the Bonin Ridge. The GMT proved to be an effective instrument for marine geological mapping, 3D and 2D modelling, statistical analysis and graphical plotting.

The main purpose of this article is to present the use of R programming language in cartographic visualization demonstrating using machine learning methods in geographic education. Current trends in education technologies are largely influenced by the possibilities of distance-learning, e-learning and self- learning. In view of this, the main tendencies in modern geographic education include active use of open source GIS and publicly available free geospatial datasets that can be used by students for cartographic exercises, data visualization and mapping, both at intermediate and advanced levels. This paper contributes to the development of these methods and is fully based on the datasets and tools available for every student: the R programming language and the free open source datasets. The case study demonstrated in this paper show the examples of both physical geographic mapping (geomorphology) and socio-economic geography (regional mapping) which can be used in the classes and in self-learning. The objective of this research includes geomorphological modelling of the terrain relief in Italy and regional mapping. The data include DEM SRTM90 and datasets on regional borders of Italy embedded in R packages ‘maps’ and ‘mapdata’. Modelling references to the characteristics of slope, aspect, hillshade and elevation, their visualization using R packages: ‘raster’ and ‘tmap’. Regional mapping of Italy was made using main package ‘ggmap’ with the ‘ggplot2’ as a wrapper. The results present five thematic maps (slope, aspect, hillshade, elevation and regions of Italy) created in R language. Traditionally used in statistical analysis, R is less known as a perfect tool in geographic education. This paper contributes to the development of methods in geographic education by presenting new technologies of the machine learning methods of mapping.

International audience; U ovom istraživanju razvijen je integrisani okvir za analizu podataka o Zemlji (“Big Earth data”) u kon- tekstu geomorfologije Jordana. Istraživanje se bavi korelacijom između nekoliko tematskih skupova podataka, uključujući mašinsko učenje i multidisciplinarne geoprostorne podatke. GIS kartiranje se široko koristi u geološkom mapiranju kao najadekvatnije tehničko sredstvo za vizualizaciju i anal- izu podataka. GIS aplikacije podstiču geološko prospektivno modeliranje vizualizacijom podataka usmjerenih na prognozu mineralnih resursa. Međutim, automatizacija pomoću mašinskog učenja za obradu velikih podataka o Zemlji pruža brzinu i tačnu obradu masivnih skupova podataka iz više izvora. To je moguće primjenom skriptiranja i programiranja u kartografskim tehnikama. Ova stu- dija predstavlja kombinovane metode mašinskog učenja kartografske analize i modeliranja velikih podataka o Zemlji. Cilj studije je da se analizira povezanost faktora koji utiču na geomorfološki oblik Jordana u odnosu na rasjed Mrtvog mora i geološku evoluciju. Tehnička metodologija uključuje tri nezavisna alata: 1) Generic Mapping Tools (GMT); 2) odabrane biblioteke programskog jezika R; 3) QGIS. Konkretno, GMT skriptni program korišćen je za topografsko, seizmičko i geofizičko mapiranje; QGIS - za geološko mapiranje; Programski jezik R - za geomorfometrijsko modeliranje. Shodno tome, tok rada je logično strukturisan kroz ova tri tehnička alata, predstavljajući različite kartografske pri- stupe za obradu podataka. Podaci i materijali uključuju skupove podataka iz više izvora različite rezo- lucije, prostornog opsega, porijekla i formata. Rezultati su predstavili kartografske rasporede kvalita- tivnih i kvantitativnih karata sa statističkim rezimeom (histogrami). Novost ovog pristupa objašnjava se potrebom da se smanji tehnički jaz između tradicionalnog GIS-a i skripti za kartiranje, koje je bit- no za mapiranje velikih podataka, gdje su presudni faktori brzina i preciznost rukovanja podacima i efikasna vizuelizacija postignuta mašinskom grafikom. U radu se analiziraju osnovni geološki procesi koji utiču na formiranje geomorfoloških oblika terena u Jordanu sa 3D vizualizacijom izabranog frag- menta zone rasjeda Mrtvog mora. Istraživanje predstavlja prošireni opis metodologije, uključujući objašnjenja isječaka koda iz modula GMT i primjere upotrebe R-biblioteka ‚raster‘ i ‚tmap‘. Rezultati su otkrili snažnu korelaciju između geoloških i geofizičkih karakteristika terena i geomorfoloških obrazaca. Integrisano proučavanje geomorfologije Jordana zasnovano je na skupovima podataka koji obrađuju više scenarija. Temeljnom analizom predstavljene su regionalne korelacije između geomorfološkog, geološkog i tektonskog okruženja u Jordanu. Rad je doprinio razvoju kartografskog inženjerstva uvođenjem skriptnih tehnika i regionalnim studijama Jordana, uključujući Mrtvo more kao posebnu regiju Jordana. Rezultati uključuju 12 novih tematskih mapa, uključujući 3D model.; In this research, an integrated framework on the big Earth data analysis has been developed in the context of the geomorphology of Jordan. The research explores the correlation between several thematic datasets, including machine learning and multidisciplinary geospatial data. GIS mapping is widely used in geological mapping as the most adequate technical tool for data visualization and analysis. GIS applications encourage geological prospective modeling by visualizing data aimed at the prognosis of mineral resources. However, automatization using machine learning for big Earth data processing provides the speed and accurate processing of multisource massive datasets. This is enabled by the application of scripting and programming in cartographic techniques. This study presents the combined machine learning methods of cartographic analysis and big Earth data modeling. The objective is to analyze a correlation between the factors affecting the geomorphological shape of Jordan with respects to the Dead Sea Fault and geological evolution. The technical methodology includes the following three independent tools: 1) Generic Mapping Tools (GMT); 2) Selected libraries of R programming language; 3) QGIS. Specifically, the GMT scripting program was used for topographic, seismic and geophysical mapping, while QGIS was used for geologic mapping and R language for geomorphometric modeling. Accordingly, the workflow is logically structured through these three technical tools, representing different cartographic approaches for data processing. Data and materials include multisource datasets of the various resolution, spatial extent, origin and formats. The results presented cartographic layouts of qualitative and quantitative maps with statistical summaries (histograms). The novelty of this approach is explained by the need to close a technical gap between the traditional GIS and scripting mapping, which is wider for big data mapping and where the crucial factors are speed and precision of data handling, as well as effective visualization achieved by the machine graphics. The paper analyzes the underlying geologic processes affecting the formation of geomorphological landforms in Jordan with a 3D visualization of the selected fragment of the Dead Sea Fault zone. The research presents an extended description in methodology, including the explanations of code snippets from the GMT modules and examples of the use of R libraries 'raster' and 'tmap'. The results revealed strong correlation between the geological and geophysical settings which affect geomorphological patterns. Integrated study of the geomorphology of Jordan was based on multisource datasets processed by scripting. A thorough analysis presented regional correlations between the geomorphological, geological and tectonic settings in Jordan. The paper contributed both to the development of cartographic engineering by introducing scripting techniques and to the regional studies of Jordan including the Dead Sea Fault as a special region of Jordan. The results include 12 new thematic maps including a 3D model.

International audience; The paper presents the cartographic processing of the Landsat TM image by the two unsupervised classification methods of SAGA GIS: ISODATA and K-means clustering. The approaches were tested and compared for land cover type mapping. Vegetation areas were detected and separated from other land cover types in the study area of southwestern Iceland. The number of clusters was set to ten classes. The processing of the satellite image by SAGA GIS was achieved using Imagery Classification tools in the Geoprocessing menu of SAGA GIS. Unsupervised classification performed effectively in the unlabeled pixels for the land cover types using machine learning in GIS. Following an iterative approach of clustering, the pixels were grouped in each step of the algorithm and the clusters were reassigned as centroids. The paper contributes to the technical development of the application of machine learning in cartography by demonstrating the effectiveness of SAGA GIS in remote sensing data processing applied for vegetation and environmental mapping.

International audience; Along with the development of the scripting technology in cartography, such as the Generic Mapping Tools (GMT) and libraries of R programming language, geologic and geophysical mapping is being vigorously promoted, where the integration of the thematic data, such as GEBCO/SRTM, EGM-2008 and open geological raster and vector layers is one of the primary datasets that provides the highresolution raw sources for cartographic visualization in the geologically complex regions like Pakistan. This study aims to integrate scripting methods of automated cartography, methods of applied geoinformatics for geomorphometric analysis and technical data processing (formatting, projecting, plotting), to provide a synthesis of the geological, geophysical and geomorphological maps of Pakistan they as new information supporting analysis of the geospatial variations of geology, geomorphology, tectonics and gravity fields with a special focus on the geologically remarkable region of Pakistan: Sulaiman Fold Belt and Muslim Bagh ophiolite complex. This study presents new 12 thematic maps, which are technically made using scripting approaches and open tools. All maps cover the region of Pakistan and they are made using open source tools: GMT, R and QGIS. A GMT and R based scripting mapping is applied for mapping Pakistan, and its algorithm steps are presented stepwise as code snippets. A system complex approach of the data integration and formats reshaping, data conversion and reformatting for a single project of the geology of Pakistan is designed and developed based on the combination of the programming and scripting techniques and with additional QGIS based mapping, which effectively integrates the thematic geospatial multi-origin datasets. Various color palettes and cartographic visualization approaches have been used to achieve the best visualization. The resulting maps are explained and discussed. Correlation between spatial phenomena of Earth's gravity, geologic evolution and tectonic movements were pointed out and commented on. New 12 maps present the regional geologic setting of the country.

International audience; Cartographic data visualization is crucial for geological mapping in seismically active areas of Malawi. Representation of the datasets from IRIS database is beneficial for hazard risk assessment, forecasting and damage mitigation, especially for mapping the distribution and magnitude of earthquakes. Scripting cartography using Generic Mapping Tools (GMT) was used to map seismicity in Malawi Rift Zone based on IRIS catalogue for 1972-2021. The maps show relations among the elevation heights, distribution of seismic events and volcanos in the north of the Malawi Lake. The results show correlation between elevations and seismicity with local topographic depressions. Seismic data show variability of the seismicity level with more earthquakes recorded in the north of the country. This paper contributes to the regional studies of Malawi for risk assessment and geological hazard analysis.

3D restitution is now a well-known tool to validate hypotheses on historical buildings that do not exist anymore. The present project takes the method a step further in order to explore the art historical topic of ornament reuse in 19th century revival architecture, particularly in Cairo. The case study is the Saint-Maurice residence, built 1875-79, for which an extensive collection of documents in varied formats, and from multiple locations and disciplines, has been conducted. The paper presents some preliminary results on the 3D restitution, the remaining open questions and the challenges they raise., Eurographics Workshop on Graphics and Cultural Heritage, Short Papers II, 117, 120, Vincent Baillet, Pascal Mora, Corentin Cou, Sarah Tournon-Valiente, Mercedes Volait, Xavier Granier, Romain Pacanowski, and Ga��l Guennebaud, CCS Concepts: Applied computing --> Arts and humanities; Information systems --> Digital libraries and archives

The functionality of Generic Mapping Tools (GMT) to process and visualize geospatial information is crucial to the development of the advanced cartographic method. This paper presents modelling and spatial analysis of the marine geological data using GMT shell scripting. GMT demonstrated effective cartographic solutions for visualization of the georeferenced data. The particular feature of GMT consists in its scripting modular approach that enables to use machine learning to explore reliable georeferenced data. Here, the study applies a sequential shell scripting to devise GMT modules for depicting marine geological data on the Mariana Trench. The data cover bathymetry, geophysics, tectonics and geology. The first method makes use of the 'nearneighbor' GMT module for grid contour modelling using Nearest Neighbor algorithm. This form of modelling classifies the geospatial data based on a similarity. The second method presents surface modelling from the initial XYZ-ASCII dataset by a combination of the 'blockmean' and 'surface' modules. The third method includes the use of the modules 'grdimage', 'psbasemap' and 'grdcontour' for plotting. Compared to GIS methods in which data are processed in a menu, GMT presents the console-based approach which automates cartographic data processing. The results present seven new maps and explanations of scripts.A combination of visual approaches applied using a color fill and various textures to represent data, which is effective in allowing readers to assess geophysical setting. The study demonstrated the effectiveness of GMT in geodata visualization.

In this paper, we propose to address the problem of tracking a deformable surface, typically the free surface of a liquid. This surface domain sees its geometry and topology evolve over time by displacement of its vertices, so the elements of the mesh (edges and facets) are all potentially contracted or expanded and require a remeshing. To do so, we propose to use a technique based on restricted Voronoi diagrams. Voronoi diagrams offer a space partition and, more particularly, a partition of the surface domain considered that allow us, among other things, to optimize the distribution of a sample set among the domain and to define a specific triangulation : the restricted Delaunay triangulation. This remeshing solution is only effective when certain conditions are met. Therefore, the first work consisted in implementing an analysis of the restricted cell configurations to ensure that the dual object meets the definition of triangulated manifold and that it is homeomorphic to the initial domain. For less favourable configurations, we have developed a method to correct the partition automatically. Based on the previous analysis, the method proposes a new minimal approximation for each of the faulty cells, thus we can limit the number of vertices used, contrary to the classicalDelaunay refinement. A second work proposes to improve the proximity between the initial mesh and the result of the remeshing : new vertex positions are finely adjusted to minimize the approximation error which is here expressed as local volume differences. These tools are combined with a sampling strategy that allows to maintain a constant sampling density throughout the deformation and thus we propose a new method to track free surfaces in incompressible fluid simulation.; Dans ce mémoire, nous proposons une technique de remaillage adaptée au suivi de surfaces déformables, typiquement la surface libre d'un liquide. Ce domaine surfacique voit sa géométrie et sa topologie évoluer au cours du temps par déplacement de ces sommets aussi, les éléments constitutifs du maillage (arêtes et facettes) sont tous potentiellement contractés ou dilatés et requièrent une ré-évaluation.Nous proposons ici d'exploiter une technique de remaillage fondée sur les diagrammes de \voronoi ; ces derniers offrent une une partition de l'espace et, plus particulièrement, du domaine surfacique considéré. Cette décomposition permet entre autre d'optimiser la répartition d'un ensemble d'échantillons sur ce domaine et d'en définir une triangulation remarquable : la triangulation de Delaunay restreinte. Cette solution de remaillage n'est efficace que lorsque certaines conditions sont réunies. Aussi, le premier travail a consisté à mettre en oeuvre une technique d'analyse des configurations de cellules permettant d'assurer que l'objet dual répond à l'appellation de variété triangulée et qu'il est homéomorphe au domaine initial. Pour les configurations moins favorables, nous avons développé une méthode de correction automatique du partitionnement ; elle s'appuie sur la précédente analyse et propose une nouvelle approximation minimale pour chacune des cellules miseen défaut. Un second travail propose d'améliorer la proximité entre le maillage initial et le résultat du remaillage : il s'agit de minimiser l'erreur d'approximation qui est ici exprimée sous forme de différences locales de volume.Ces deux outils sont associés à une stratégie d'échantillonnage qui permet de maintenir une densité d'échantillonnage constante tout au long de la déformation et propose ainsi une nouvelle méthode de suivi d'une surface libre pour la simulation d'écoulement de fluides incompressibles.

International audience; This paper presents the GRASS GIS-based thematic mapping of Antarctica using scripting approach and associated datasets on topography and geophysics. The state-ofthe art in cartographic development points at two important aspects. The first one comprises shell scripting promoted repeatability of the GIS technique, increased automatization in cartographic workflow, and compatibility of GRASS with Python, PROJ and GDAL libraries which enables advanced geospatial data processing: converting formats, re-projecting and spatial analysis. The second aspect is that data visualization greatly influences geologic research through improving the interpretation between the Antarctic glaciation and surface. This includes the machine learning algorithms of image classification enabling to distinguish between glacier and non-glacier surfaces through automatically partitioning data and analysis of various types of surfaces. Presented detailed maps of Antarctic include visualized datasets from the ETOPO1, GlobSed, EGM96 and Bedmap2 projects. The grids include bed and surface elevation, ETOPO1-based bathymetry and topography, bed, ice and sediment thickness, grounded bed uncertainty, subglacial bed elevation, geoid undulations, ice mask grounded and shelves. Data show the distribution of the present-day glacier, geophysical fields and topographic landforms for analysis of processes and correlations between the geophysical and geological phenomena. Advances in scripting cartography are significant contributions to the geological and glaciological research. Processing high-resolution datasets of Southern Ocean retrieved by remote sensing methods present new steps in automatization of the digital mapping, as presented in this research, and promotes comprehensive monitoring of geological, permafrost and glacial processes in Antarctica. All maps have been plotted using GRASS GIS version 7.8. with technical details of scripts described and interpreted.

International audience; Seismicity in Yemen and the Gulf of Aden in a geological context. The study presents geologic investigation of Yemen and the Gulf of Aden with a special focus on geophysical, seismic, tectonic and topographic mapping performed by the integrated approach of QGIS and GMT scripting. Cartographic visualization is crucial in geologic analysis, data processing and prognosis of mineral resource prospects. The region of Yemen and Gulf of Aden was formed as a result of Arabian and African plates movements and still tectonically active. Besides, the Gulf of Aden contains mineral resources of hydrocarbons which makes this region actual for investigation. The IRIS database on earthquakes was used for visualization of the magnitude of submarine earthquakes in the Gulf of Aden for the period of 2007-2020. The paper presents 6 new thematic maps for the region of Yemen and Gulf of Aden. The research presented an analysis of correlation between the geological, topographic and geophysical settings. Through combined approach of cartographic high-resolution data visualization and geologic analysis, this paper contributed to the regional geological studies of Yemen, Gulf of Aden and the Middle East.

International audience; This paper presents using GDAL utilities and GRASS GIS for topographic analysis of the raster grids based on GEBCO DEM as NetCDF file at 15 arc-second intervals. The focus study area encompasses the area around Okinawa Trough, Ryukyu trench-arc system, southern Japan, East China Sea and the Philippine Sea, west Pacific Ocean. Several GDAL utilities were applied for data processing: gdaldem, gdalwarp, gdalinfo, gdal_translate. The data were imported to GRASS GIS via r.in.gdal. Data visualization highlighted high resolution and accuracy of GEBCO grid, enabling topographic modelling at the advanced level. The algorithm of DEM processing, implemented in GDAL utility gdaldem, was used for generating multipurpose topographic models: aspect, hillshade, roughness and topographic indices, such as Topographic Position Index (TPI), Terrain Ruggedness Index (TRI). Thematic maps (topography, geoid, marine free-air gravity) were visualized using GRASS GIS modules for raster (d.rast, r.colors, r.contour) and vector (d.vect, v.in.region, d.legend) data processing. The results demonstrated smoother bathymetry in the East China Sea and rugged relief in the Philippine Sea which corresponds to their different geological and geophysical settings. The presented methodology of the topographic analysis based on DEM demonstrated technical aspects of GDAL and GRASS as scripting approach of advanced cartography.

Uganda is traversed by the East African Rift System (EARS), which controls its topography and induces high seismicity in the southwestern and western parts of the country, around the Rwenzori Mountains and Albertine Graben. This paper explores the complexity of the geophysical, topographic and geomorphological settings of Uganda for linking the distribution and magnitude of earthquakes to the controlling factors: geomorphology, topography and geophysics. The methodological workflow for an automated plotting of the topographic, geophysical and seismic data includes the following important codes: 'gmt grdimage ug_relief.nc-Cgeo.cpt-R29/35/-1.5/4.3-JM6.5i-I+a15+ne0.75-t50-Xc-P-K > $ps' for topographic visualisation, 'gmt psxy-R-J quakes_UG.ngdc-Wfaint-i4,3,6,6s0.1-h3-Scc-Csteps.cpt-O-K >> $ps' for mapping earthquakes, 'gmt img2grd grav_27.1.img-R29/35.5/-1.5/4.3-Ggrav_UG.grd-T1-I1-E-S0.1-V' for mapping gravity. The steps of the scripting workflow using 'raster' and 'tmap' libraries of R and a variety of GMT modules includes the following most important codes: 'alt = getData("alt", country = "Uganda", path = tempdir()); slope = terrain(alt, opt = "slope"); aspect = terrain(alt, opt = "aspect"); hill = hillShade(slope, aspect, angle = 40, direction = 270)'. The study uses high-resolution spatial datasets GEBCO, EGM-2008, IRIS, gravity grids and SRTM DEM aimed to undertake thematic mapping of Uganda. Nine new maps have been generated using shell scripts of GMT and R for environmental monitoring. Depth of the earthquakes ranges from 4 to 40 Km with frequency of seismic events near Rwenzori Mountains and Albertine Graben for destructive earthquake is less than about five decades. Earthquake magnitudes varied from 3.7 Richter Magnitude Scale (Lake Edward region) to 6.2 Richter local magnitude Makerere scale (M L) (Lake Albert region). Lake Kivu region has earthquakes with magnitudes between 3.8-5.3. The region of Lake Albert had the highest magnitudes and frequency of seismic events, compared to Lakes Kivu and Edward. Earthquakes around the Lake Albert reached the highest magnitude at 6.2 M L. Earthquakes around Lake Victoria reach the maximal level of earthquake magnitude at 5.0 M l , near Lake George-5.3 M L. Earthquakes were mostly recorded in the Albertine Graben of EARS in the Albert, Edward and George lakes and occasionally in Lake Victoria.

International audience; The paper presents the use of the Landsat TM image processed by the ArcGIS Spatial Analyst Tool for environmental mapping of southwestern Iceland, region of Reykjavik. Iceland is one of the most special Arctic regions with unique flora and landscapes. Its environment is presented by vulnerable ecosystems of highlands where vegetation is affected by climate, human or geologic factors: overgrazing, volcanism, annual temperature change. Therefore, mapping land cover types in Iceland contribute to the nature conservation, sustainable development and environmental monitoring purposes. This paper starts by review of the current trends in remote sensing, the importance of Landsat TM imagery for environmental mapping in general and Iceland in particular, and the requirements of GIS specifically for satellite image analysis. This is followed by the extended methodological workflow supported by illustrative print screens and technical description of data processing in ArcGIS. The data used in this research include Landsat TM image which was captured using GloVis and processed in ArcGIS. The methodology includes a workflow involving several technical steps of raster da ta processing in ArcGIS: 1) coordinate projecting, 2) panchromatic sharpening, 3) inspection of raster statistics, 4) spectral bands combination, 5) calculations, 6) unsupervised classification, 7) mapping. The classification was done by clustering technique using ISO Cluster algorithm and Maximum Likelihood Classification. This paper finally presents the results of the ISO Cluster application for Landsat TM image processing and concludes final remarks on the perspectives of environmental mapping based on Landsat TM image processing in ArcGIS.The results of the classification present landscapes divided into eight distinct land cover classes: 1) bare soils; 2) shrubs and smaller trees in the river valleys, urban areas including green spaces; 3) water areas; 4) forests including the Reykjanesfólkvangur National reserve; 5) ice-covered areas, glaciers and cloudy regions; 6) ravine valleys with a sparse type of the vegetation: rowan, alder, heathland, wetland; 7) rocks; 8) mixed areas. The final remarks include the discussion on the development of machine learning methods and opportunities of their technical applications in GIS-based analysis and Earth Observation data processing in ArcGIS, including image analysis and classification, mapping and visualization, machine learning and environmental applications for decision making in forestry and sustainable development.

International audience; The paper presents the use of scripting cartographic techniques for the visualization of topographic and geomorphological maps by R and GMT. The thematic maps intend to analyze the region of Peru with specific focus on its geomorphology: slope, aspect, elevation and hillshade performed in R libraries ‘raster’ and ‘tmap’. The integrated free materials include several datasets embedded in R library ‘OpenStreetMap’: Stamen, ESRI World Imagery, Bing Maps, National Park Service. The research shows a particular example of the use of open source datasets and free tools available for distance-based online education and research, which becomes an actual tendency in modern geographic research.; El artículo presenta el uso de técnicas cartográficas de scripting para la visualización de mapas topográficos y geomorfológicos por R y GMT. Los mapas temáticos pretenden analizar la región del Perú con un enfoque específico en su geomorfología: pendiente, aspecto, elevación y sombreado realizado en las bibliotecas R «raster» y «tmap». Los materiales integrados incluyen varios conjuntos de datos integrados en la biblioteca R "OpenStreetMap": Stamen, ESRI World Imagery, Bing Maps, National Park Service. La investigación muestra un ejemplo particular del uso de conjuntos de datos de código abierto y herramientas gratuitas disponibles para la educación e investigación en línea a distancia, que se convierte en una tendencia real en la investigación geográfica moderna.

International audience; This study evaluates the geomorphometric parameters of the topography in Ethiopia using scripting cartographic methods by applying R languages (packages 'tmap' and 'raster') and Generic Mapping Tools (GMT) for 2D and 3D topographic modelling. Data were collected from the open source repositories on geospatial data with high resolution: GEBCO with 15 arc-second and ETOPO1 with 1 arc-minute resolution and embedded dataset of SRTM 90 m in 'raster' library of R. The study demonstrated application of the programming approaches in cartographic data visualization and mapping for geomorphometric analysis. This included modelling of slope steepness, aspect and hillshade visualized using DEM SRTM90 to derive geomorphometric parameters of slope, aspect and hillshade of Ethiopia and demonstrate contrasting topography and variability climate setting of Ethiopia. The topography of the country is mapped, including Great Rift Valley, Afar Depression, Ogaden Desert and the most distinctive features of the Ethiopian Highlands. A variety of topographical zones is demonstrated on the presented maps. The results include 6 new maps made using programming console-based approach which is a novel method of cartographic visualization compared to traditional GIS software. The most important fragments of the codes are presented and technical explanations are provided. The presented series of 6 new maps contributes to the cartographic data on Ethiopia and presents the methodology of scripting mapping techniques.

Zambia recently experienced several environmental threats from climate change such as droughts, temperature rise and occasional flooding and they all affect agricultural sustainability and people wellbeing through negative effects on plants and growing crops. This paper is aimed at showing variations in several climate and environmental parameters in Zambia showing spatial variability and trends in different regions of Zambia's key environmental areas (Zambezi River and tributaries), Livingstone near the Victoria Falls and central region with Muchinga Mountains. A series of 10 maps was plotted using data from TerraClimate dataset: precipitation, soil moisture, Palmer Drought Severity Index (PDSI), downward surface shortwave radiation, vapor pressure deficit and anomalies, potential and actual evapotranspiration and wind speed with relation to the topographic distribution of elevations in Zambia plotted using GEBCO/SRTM data. The data range of the PDSI according to the index values ranged from minimum at -5.7 to the maximum at 16.6 and mean at 7.169, with standard deviation at 4.278. The PDSI is effective in quantifying drought in long-term period. Because PDSI index applies temperature data and water balance model, it indicates the effect of climate warming on drought by correlation with potential evapotranspiration. The maximum values for soil moisture of Zambia show minimum at 1 mm/m, maximum at 413 mm/m, mean at 173 mm/m. This study is technically based on using the Generic Mapping Tools (GMT) as cartographic scripting toolset. The paper contributes to the environmental monitoring of Zambia by presenting a series of climate and environmental maps that are beneficial for agricultural mapping of Zambia.

Current paper presents the results of the experimental analysis on permafrost uppermost soil samples with various physical properties (moisture, porosity) tested with varied external pressure and time. The aim of this work is to test properties of the soil samples intended for the construction of buildings, railways and objects of civil infrastructure by modeled external pressure, data visualization and analysis. Variations in the soil samples were studied by analysis of the equivalent soil cohesion (Ceq) in frozen soil samples. Methods include integrated application of the laboratory experiments, methods of the statistical data analysis and 3D plotting performed by the selected LaTeX packages. Laboratory experiments were performed using KrioLab equipment ‘Sharikovy Stamp PSH-1’. The 15 series of experiments have been tested. Models of the soil strength are graphically presented and statistically analyzed showing the results of the experiment.

The paper concerns GMT application for studies of the geophysical and geomorphological settings of the Weddell Sea. Its western part is occupied by the back-arc basin developed during geologic evolution of the Antarctic. The mapping presents geophysical settings reflecting tectonic formation of the region, glaciomarine sediment distribution and the bathymetry. The GlobSed grid highlighted the abnormally large thickness of sedimentary strata resulted from the long lasting sedimentation and great subsidence ratio. The sediment thickness indicated significant influx (>13,000m) in the southern segment. Values of 6,000-7,000 m along the peninsula indicate stability of the sediments influx. The northern end of the Filchner Trough shows increased sediment supply. The topography shows variability-7,160-4,763 m. The ridges in the northern segment and gravity anomalies (>75 mGal) show parallel lines stretching NW-SE (10°-45°W, 60°-67°S) which points at the effects of regional topography. The basin is dominated by the slightly negative gravity >-30 mGal. The geoid model shows a SW-NE trend with the lowest values 20m in the NE and along the Coats Land. The ripples in the north follow the geometry of the submarine ridges and channels proving correlation with topography and gravitational equipotential surface.

International audience; The South Sandwich Trench located eastward of the Drake Passage in the Scotia Sea between Antarctica and South America is one of the least studied deep-sea trenches. Its geomorphological formation and present shape formed under the strong influence of the tectonic plate movements and various aspects of the geological setting, i.e., sediment thickness, faults, fracture zones and geologic lineaments. The aim of this paper is to link the geological and geophysical setting of the Scotia Sea with individual geomorphological features of the South Sandwich Trench in the context of the phenomena of its formation and evolution. Linking several datasets (GEBCO, ETOPO1, EGM96, GlobSed and marine freeair gravity raster grids, geological vector layers) highlights correlations between various factors affecting deep-sea trench formation and development, using the Generic Mapping Tools (GMT) for cartographic mapping. The paper contributes to the regional studies of the submarine geomorphology in the Antarctic region with a technical application of the GMT cartographic scripting toolset.

International audience; Detailed mapping based on the high-resolution grids, such as GEBCO, ETOPO1, GlobSed, EGM-2008 is crucial for various domains of Earth sciences: geophysics, glaciology, Quaternary, sedimentology, geology, environmental science, geomorphology, etc. The study presented a GMTbased scripting technique of the cartographic data processing aimed at the comparative analysis of the bathymetry, sediment thickness, geologic objects and geophysical settings in the study area based on various datasets. The study area is located in the Ross Sea, Antarctic. The highest values of the sediment thickness over 7,500 m are dominating in the southwest segment of the Ross Sea closer to the Victoria Land, followed by the region over the Ross Ice Shelf with values between 5,500 to 7,000 m (170°-175°W). The increased sediment thickness (2,500 to 3,000 m) was also mapped seen in the region NE off the Sulzberger Bay (70-75°S to 140-155°W), caused by the closeness of the Marie Bird Land ice coasts. A remarkable correlation between the gravity and the topography of the sea-land border in the Marie Bird Land area is well reflected in the coastal line and a set of the higher values in the free-air gravity. On the contrary, negative values (-60 to-80 mGal) are notable along the submarine troughs stretching parallel in the western part of the basin: e.g. the trough stretching in NW-SE direction in the 170°W-175°E, 65°S-68°S, between the 167°W-175°W, 70°S-72°S. Such correlations are clearly visible on the map, indicating geological lineaments and bathymetric depressions correlating with gravity grids. The paper contributes to the regional studies of the Ross Sea, the Antarctic and Polar region, and development of the technical cartographic methodologies by presenting an application of the GMT for thematic mapping.; O mapeamento detalhado com base em grades de alta resolução, como GEBCO, ETOPO1, GlobSed, EGM-2008 é fundamental para vários campos das ciências da terra: geofísica, glaciologia, quaternário, sedimentologia, geologia, ciências ambientais, geomorfologia, etc. O estudo apresenta métodos de processamento de dados de mapas com script baseado em GMT visando a análise comparativa de batimetria, espessura de sedimentos, características geológicas e condições geofísicas na área de estudo com base em vários conjuntos de dados. A Área de Estudo está localizada no Mar de Ross, na Antártica. Os valores mais altos de espessura de sedimento de mais de 7.500 m prevalecem no sudoeste do Mar de Ross próximo a Victoria Land, seguido pela área acima da plataforma de gelo de Ross em valores entre 5.500 e 7.000 m (170 ° -175 ° W). ... O aumento da espessura do sedimento (de 2.500 a 3.000 m) também foi mapeado na área nordeste da Baía Sulzberger (de 70-75 ° S a 140-155 ° W), devido à proximidade das costas de gelo da Terra Marie Bird. A notável correlação entre a gravidade e a topografia da fronteira marítima perto da Terra Marie Bird é bem refletida na linha costeira e em muitos dos maiores valores de gravidade ao ar livre. Em contraste, os valores negativos (-60 a -80 mGal) são perceptíveis ao longo das muralhas submarinas que se estendem paralelamente na parte oeste da bacia: por exemplo, uma trincheira que se estende na direção noroeste-sudeste a 170 ° W-175 ° E, 65 ° S-68 ° S, entre 167 ° W-175 ° W, 70 ° S-72 ° S Essas correlações são claramente visíveis no mapa, indicando lineamentos geológicos e depressões batimétricas correlacionadas com grades gravimétricas. O documento contribui para estudos regionais do Mar de Ross, Antártica e região polar, bem como para o desenvolvimento de metodologias técnicas de mapeamento, introduzindo a aplicação GMT para mapeamento temático.

International audience; Detailed mapping based on the high-resolution grids, such as GEBCO, ETOPO1, GlobSed, EGM-2008 is crucial for various domains of Earth sciences: geophysics, glaciology, Quaternary, sedimentology, geology, environmental science, geomorphology, etc. The study presented a GMTbased scripting technique of the cartographic data processing aimed at the comparative analysis of the bathymetry, sediment thickness, geologic objects and geophysical settings in the study area based on various datasets. The study area is located in the Ross Sea, Antarctic. The highest values of the sediment thickness over 7,500 m are dominating in the southwest segment of the Ross Sea closer to the Victoria Land, followed by the region over the Ross Ice Shelf with values between 5,500 to 7,000 m (170°-175°W). The increased sediment thickness (2,500 to 3,000 m) was also mapped seen in the region NE off the Sulzberger Bay (70-75°S to 140-155°W), caused by the closeness of the Marie Bird Land ice coasts. A remarkable correlation between the gravity and the topography of the sea-land border in the Marie Bird Land area is well reflected in the coastal line and a set of the higher values in the free-air gravity. On the contrary, negative values (-60 to-80 mGal) are notable along the submarine troughs stretching parallel in the western part of the basin: e.g. the trough stretching in NW-SE direction in the 170°W-175°E, 65°S-68°S, between the 167°W-175°W, 70°S-72°S. Such correlations are clearly visible on the map, indicating geological lineaments and bathymetric depressions correlating with gravity grids. The paper contributes to the regional studies of the Ross Sea, the Antarctic and Polar region, and development of the technical cartographic methodologies by presenting an application of the GMT for thematic mapping.; O mapeamento detalhado com base em grades de alta resolução, como GEBCO, ETOPO1, GlobSed, EGM-2008 é fundamental para vários campos das ciências da terra: geofísica, glaciologia, quaternário, sedimentologia, geologia, ciências ambientais, geomorfologia, etc. O estudo apresenta métodos de processamento de dados de mapas com script baseado em GMT visando a análise comparativa de batimetria, espessura de sedimentos, características geológicas e condições geofísicas na área de estudo com base em vários conjuntos de dados. A Área de Estudo está localizada no Mar de Ross, na Antártica. Os valores mais altos de espessura de sedimento de mais de 7.500 m prevalecem no sudoeste do Mar de Ross próximo a Victoria Land, seguido pela área acima da plataforma de gelo de Ross em valores entre 5.500 e 7.000 m (170 ° -175 ° W). ... O aumento da espessura do sedimento (de 2.500 a 3.000 m) também foi mapeado na área nordeste da Baía Sulzberger (de 70-75 ° S a 140-155 ° W), devido à proximidade das costas de gelo da Terra Marie Bird. A notável correlação entre a gravidade e a topografia da fronteira marítima perto da Terra Marie Bird é bem refletida na linha costeira e em muitos dos maiores valores de gravidade ao ar livre. Em contraste, os valores negativos (-60 a -80 mGal) são perceptíveis ao longo das muralhas submarinas que se estendem paralelamente na parte oeste da bacia: por exemplo, uma trincheira que se estende na direção noroeste-sudeste a 170 ° W-175 ° E, 65 ° S-68 ° S, entre 167 ° W-175 ° W, 70 ° S-72 ° S Essas correlações são claramente visíveis no mapa, indicando lineamentos geológicos e depressões batimétricas correlacionadas com grades gravimétricas. O documento contribui para estudos regionais do Mar de Ross, Antártica e região polar, bem como para o desenvolvimento de metodologias técnicas de mapeamento, introduzindo a aplicação GMT para mapeamento temático.

The Cascadia Trench is stretching along the convergent plate boundaries of Pacific Plate, North America Plate and Juan De Fuca Plate. It is an important geomorphological structural feature in the north-east Pacific Ocean. The aim of the paper is to analyse the geomorphology of the Cascadia Trench west of Vancouver Island (Canada and USA) using the GMT cartographic scripting toolset. The unique geomorphological feature of the Cascadia Trench is that the thick sediment layer completely obscures the subduction zone and abyssal hills. This results in the asymmetric profile in the cross-section of the trench. Bathymetric data were extracted from the GEBCO 2019 dataset (15 arc-second grid), sediment thickness by the GlobSed dataset. Due to the dominance of high sedimentary rate and complexity of the tectonic processes and geologic settings, Cascadia Trench develops very specific asymmetric geomorphic shape comparing to the typical V-form. The results of the geomorphic modelling show that eastern side of the trench has a gentle curvature (slope: 35.12°), partially stepped, due to the tectonic movements and faults. The opposite, oceanward side is almost completely leveled. The trench is narrow with maximal depth at the selected segment -3489 m and for the whole dataset -6201 m. The most repetitive depth is in a range -2500 to -2400 m (267 samples) and -2500 to -2600 m (261 samples). The bottom is mostly flat due to the high sedimentation rates indicating the accumulative leveling processes. Marine free-air gravity anomalies along the Cascadia Subduction Zone are characterized by weakly positive values (20 mGal) increasing rapidly in the zone of the continental slope (>200 mGal), which is associated with a decrease in thickness of the Earth’s crust.

International audience; The paper discusses the geomorphology of the Sunda Trench, an oceanic trench located in the eastern Indian Ocean along the Sumatra and Java Islands of the Indonesian archipelago. It analysis difference in depths and variation in slope steepness between the two segments of the trench: southern Java transect (108.8°E 10.10°S - 113.0°E 10.75°S) and northern Sumatra transect (97.5°E 1.1°S - 101.0°E 5.5°S). The maps and geomorphological modelling were plotted using Generic Mapping Tools (GMT). The data include high-resolution grids on topography, geology, geodesy and geophysics: GEBCO, EGM2008 EGM-2008, GlobSed. The results include modelled segments, slope gradients, and cross-section profiles. The geological processes take place in the Indian Ocean at different stages of its evolution and influence the nature of the submarine geomorphology and geomorphology of the trench that differs in two segments. Java segment has a bell- shaped data distribution in contrast to the Sumatra with bimodal pattern. Java segment has the most repetitive depths at -2,500 to -5,200 m. Sumatra transect has two peaks: 1) a classic bell-shaped peak (-4,500 m to -5,500 m); 2) shelf area (0 to -1,750 m). The data at middle depths (-1,750 to -4,500 m) have less than 300 samples. The most frequent bathymetry for the Sumatra segment corresponds to the -4,750 m to -5,000 m. Comparing to the Sumatra segment, the Java segment is deeper. For depths > -6,000 m, there are only 138 samples for Sumatra while 547 samples for Java. Furthermore, Java segment has a more symmetrical geometric shape while Sumatra segment is asymmetric, one-sided. The Sumatra segment has a steepness of 57.86° on its eastern side (facing Sumatra Island) and a contrasting 14.58° on the western part. The Java segment has a steepness of 64.34° on its northern side (facing Java Island) and 24.95° on the southern part (facing the Indian Ocean). The paper contributes to the studies of the submarine geomorphology in Indonesia.; Резиме: У раду су детаљно анализиране геоморфолошке карактеристике Сундског рова, који се пружа дуж острва Суматра и Јава у оквиру индонежанског архипелага, у источном делу Индијског океана. Познато је да је геоморфологија океанских ровова настала као резултат субдукције тектонских плоча током сложене геолошке еволуције. Упркос покушајима да се детаљно проучи његова геодинамика, природа и структура рова још увек нису довољно јасни. Конкретно, геоморфолошке карактеристике рова, формиране кроз сложене интеракције геофизичких, геолошких и тектонских фактора, у фокусу су многих истраживања. Сврха овог рада је да истражи разлике у геоморфолошким карактеристикама два различита сегмента Сундског рова (јужни део дуж острва Јава и северни део дуж острва Суматра) како би се олакшало боље разумевање структуре морског дна у североисточном делу Индијског океана. Методолошки приступ овог рада заснован је на употреби GMT, односно универзалних картографских инструмената, који су примењени у циљу конструисања графичких прилога. У овом квалитативном и квантитативном истраживању израђено је осам карата и обављен је преглед бројне литературе о истраживаном подручју.Јужни, јавански сегмент рова простире се дуж линије са следећим координа- тама (почетна и крајња тачка): 108,8°Е 10,10°S - 113,0°Е 10,75°S. Северни, суматрански сегмент лоциран је дуж линије са следећим координатама (почетна и крајња тачка) 97,5°Е 1,1°S - 101,0°Е 5,5°S. Подаци укључују растерску координатну мрежу високе резолуције о топографији, геологији, геодезији и геофизици: GEBCO, EGM2008EGM-2008, GlobSed. Представљено је и визуализовано неколико тематских карата са циљем демонстрирања просторне варијације у дистрибуцији геофизичких поља на проучаваном подручју. Карте посебно илуструју геолошке процесе који су се одвијали у источном делу Индијског океана у различитим фазама његове еволуције, а које су утицале на изглед подморске геоморфологије рова. Попречни профили су дигитализовани помоћу приступа машинског учења којег нуди GMT. Упоредна анализа оба сегмента извршена je визуелизацијом података у виду хистограма. Резултати истраживања показују да је геоморфологија рова зависна од локалних геолошких, геофизичких и тектонских карактеристика. Конкретно, јавански сегмент има дистрибуцију дубина у "облику звона", за разлику од суматранског са бимодалном дистрибуцијом. У сегменту Јаве најчешће је распрострањен појас дубина између -2.500 и -5.200 m. Сегмент Суматре има два јасно видљива врха: 1) врх у облику звона (-4.500 m до -5.500 m); 2) подручје шелфа (0 до -1.750 m). У поређењу са сегментом Суматре (северни), сегмент Јаве (јужни) је дубљи. На пример, за дубине преко -6.000 m постоји само 138 узорака за сегмент22Суматре, док постоји 547 узорака за сегмент Јаве. Даље, јавански сегмент има симетричан геометријски облик, док је сегмент Суматре углавном асиметричан. Сегмент Суматре има нагиб од 57,86° на источној страни (страна окренута према острву Суматра) и 14,58° на западној страни. Сегмент Јаве има нагиб од 64,34° на северној страни (страна окренута према острву Јава) и 24,95° у јужном делу (према Индијском океану).Значај картографске визуелизације је добро познат, јер тачно представљање података помаже да се изврши корелација и укаже на разлике у дистрибуцији геолошких и географских објеката, процеса и појава. Најчешће, традиционални GIS софтвери преовлађују у картографским истраживањима. Овај рад указује на ефикасну функционалност софтвера GMT за визуализацију геофизичких и геолошких података и за проучавање геоморфолошких карактеристика. Нејасност и неразумевање геолошких и геодинамичких процеса често проистиче из недостатка картографске визуализације која покрива одређено истраживачко подручје у правилној пројекцији и палети боја. Да би се повећала ефикасност у анализи података, представљене карте конструисане су у три различите пројекције: поликонусној, меркаторовој и азимутној пројекцији Ламберта, док су попречни профили урађени са декартовим xy координатама. Значај овог рада састоји се у мултидисциплинарном приступу који комбинује савремена картографска представ- љања података и геолошку анализу индонежанског региона. Резултати ове студије показују да су коришћење картографских метода за обраду геоинформација драгоцено средство за даља геофизичка и геоморфолошка моделовања и помоћ у бољем разумевању морског дна, јединствене скривене структуре Земљине површи- не.

International audience; المستخلصتم في هذه الدراسة اعداد خرائط ثنائية الأبعاد وثلاثية الأبعاد لحوض مابين النهرين في العراق والمناطق المجاورة من. خلال تكامل لبيانات متعددة المصادر بواسطة لغة برمجة TMG و R و STGQ. اعتمد رسم الخرائط الجيومورفولوجية والجيولوجية على مجموعات البيانات التي تحتوي على طبقات نقطية ومتجهة للبيانات الجيولوجية والجيوفيزيائية والطوبوغرافية والصخرية والجيوديسية التي تم جمعها من مستودعات مفتوحة المصدر. تتضمن النتائج 8 خرائط جديدة. تقدم هذه الدراسة طريقة رسم خرائط متكاملة تجمع بين أدوات مختلفة لتصور البيانات: STGQ لرسم الخرائط الجيولوجية، TMG لتخطيط الجيود، شذوذ الجاذبية، التصور ثلاثي الأبعاد ورسم الخرائط الطبوغرافية ثنائية الأبعاد، ولغة البرمجة R لنمذجة التضاريس الجيومورفولوجية لأراضي العراق. بينت الخريطة الطوبوغرافية تباين في الإرتفاع من -69 متر تحت مستوى سطح البحر في الخليج العربي الى 3754 متر فوق مستوى سطح البحر في الأجزاء الشمالية الشرقية من العراق. أظهرت الخريطة الجيولوجية ان ترسبات العصر الرباعي (Q) وترسبات العصر الثلاثي (T) هي الترسبات السائدة والتي تقسم العراق بشكل تقريبي الى جزئين (شمالي شرقي وجنوبي غربي). ترسبات العصر الرباعي (Q) هي الترسبات الغالبة التي تغطي حوض مابين النهرين. قدم النموذج ثلاثي الأبعاد رسم منظوري لحوضمابين النهرين وجبال زاكروس.عند تدوير زاوية السمت بمقدار °30/°165. تزداد تموجات الجيود بالإتجاه الشمالي الشرقي: القيمة السالبة الأكثرانخفاضاً للجيود هي (-17 م) وكانت على امتداد الخليج العربي، وتعكس استمرارية الحد الأدنى للجغرافي الهندية. الحد الاعلى للجيود هو 19 في الاتجاه الجنوبي الغربي والشمالي الشرقي. يتباين ارتفاع الجيود في حوض ما بين النهرين من -15 م الى 5 م متوافقتاً مع المنخفض الطوبوغرافي. نموذج جاذبية الهواء الحر أظهر توافق بين قيم الجاذبية والتضاريس والتي تثبت الإتجاه العام لتناقص الكتلة التي تكون سبب لقيم الجاذبية المنخفضة أو السالبة. تختلف قيم انحرافات الجاذبية لـ Faye بين الحد الادنى -124.945 ملي كال الى الحد الأعلى 327.724 ملي كال. من الواضح ان قيمالجاذبية في حوض مابين النهرين هي قيم سالبة.تساهم مخرجات رسم هذه الخرائط (8 خرائط جديدة) في زيادة المعرفة حول الجيومورفولوجيا والجيولوجيا والشذوذ الجيوفيزيائي في العراق مع التركيز على حوض - ما بين النهرين في العراق. تتعلق دراسة أسئلة البحث الفني بتنسيق البيانات الخرائطية ومعالجتها وتكاملها ونمذجتها بواسطة مجموعات بيانات إنشائية متعددة المصادر تغطي العراق. تساهمالورقة بالإضافة الى اعداد الخرائط في الدراسات التكاملية الإقليمية للعراق والشرق الأوسط وحوض مابين النهرين.; The study presents the 2D and 3D mapping of the Mesopotamia Foredeep Basin, Iraq, and the adjacent areas using integrated cartographic approach of multi-source data visualization by Generic Mapping Tools (GMT), R language and QGIS. The geomorphological and geological mapping is based on multiple datasets containing raster and vector layers of the geological, geophysical, topographic, lithological and geodetic data captured from the open source repositories and resulted in 8 new maps. This study develops and presents an integrated cartographic method combining the use of various tools for different tasks in data visualization: QGIS for geological mapping, GMT for mapping geoid, gravity anomalies, 3D mesh visualization and 2D topographic mapping, R language (libraries 'raster' and 'tmap') for modelling relief (slope and aspect) of Iraq. The topography shows variation in data from minimum at-69 m below sea level in the Arabian Gulf, to maximum at 3754 m asl in the NE of Iraq. The map of geologic units of Iraq shows dominating types are Quaternary (Q) and Tertiary (T) that roughly divide the country into two segments (SW and NE). The Mesopotamia Foredeep Basin is dominated by Quaternary (Q). The 3D model of the Iraqi relief presented perspective view of the Mesopotamia Foredeep Basin and Zagros Mts with azimuth rotation at 165°/30°. The geoid undulations increase in NW direction: the lowest, negative values of geoid (-17 m) are over the Arabian Gulf which reflects the continuum of the Indian geopotential minimum. Maximal values of geoid reach 19 m in the SW and NW. The Mesopotamia Foredeep Basin has geoid variations ranging from-15 to 5 m showing correlation with topographic depressions. Free-air gravity model shows correlation of the values with topography proving general trend of the deficit of landmasses causing the negative or lower values in gravity. The values of Faye's gravity anomalies differ from minimum at-124.945 mGal to maximum at 327.724 mGal. The Mesopotamia Foredeep Basin has clearly negative values. The generated cartographic outputs (8 new maps) contribute to the increase of knowledge on the geomorphology, geology and distribution of geophysical anomalies of Iraq with a special focus on Mesopotamia Foredeep Basin. The technical research questions covered by this study relate to the cartographic data formatting, processing, integration and modelling using thematic multi-source datasets covering Iraq. The paper contributes both to the cartographic aspect of integrated mapping and to the regional studies of Iraq, Middle East and the Mesopotamia Foredeep Basin.

International audience; The study present a case study of the Generic Mapping Tools (GMT) applied for cartographic modelling, mapping and comparative analysis of the deep-sea trenches located in southwest Pacific Ocean: the New Britain Trench (NBT) and the San Cristobal Trench (SCT). The aim was to evaluate their geomorphic variation using scripting cartographic approach of GMT. The data was processed using a sequence of the GMT modules with the main module 'grdtrack' used to visualize crosssection profiles along the trenches for their geomorphological modelling. The main grid used for topographic mapping is the SRTM DEM with 15-arc second resolution. The statistical analysis shown variability in depths of both trenches by samples in two transects. The cartographic analysis demonstrated following results. The SCT is generally deeper reaching-9,000 m, while the median for the NBT less then 7,000 m. The gradient slope of SCT is more symmetric with accurate 'V' form. In a cross-section graph, the NBT landward slope is markedly asymmetric U-shaped form and has a crescent form in the east. The NBT slope dips westwards with 35° eastward, and 41° westward, while the SCT slope has 33° oceanwards and 33,69° landwards. The difference between the geomorphology of the trenches is explained by the effects of the geotectonic evolution and actual sedimentary processes affected their formation and sculptured their structure. The marine free-air gravity anomaly illustrated density anomalies at the bathymetry in the region of NBT and SCT with range

The paper discusses geomorphology of the Sunda Trench, an oceanic trench located in eastern Indian Ocean along the Sumatra and Java Islands of the Indonesian archipelago. In particular, it analysis the difference in depths and variation in slope steepness between the two segments of the trench: the southern Java transect (coordinates 108.8°E 10.10°S to 113.0°E 10.75°S) and the northern Sumatra transect (97.5°E 1.1°S to 101.0°E 5.5°S). The thematic maps and geomorphological modelling were plotted using Generic Mapping Tools (GMT). The materials include high-resolution data on topography, geology and geophysics: GEBCO 15 arc-minute resolution grid, EGM2008 2.5 minute Earth Gravitation Model of 2008, GlobSed global 5‐arc‐minute total sediment thickness and vector geological datasets. In addition to the GEBCO-based bathymetric data, geological, topographic and geophysical maps, the results include enlarged transects for the Java and Sumatra segments, their slope gradients and cross-section profiles, derived from the bathymetric GEBCO dataset. The geomorphology framework of the Sunda Trench is largely controlled by the subduction of the Australian plate underneath the Sunda microplate. The geological processes take place in basin of the Indian Ocean at different stages of its evolution and influence the nature of the submarine geomorphology and geometric shape of the trench. Sunda Trench is seismically active part of the Pacific Ring of Fire. A large number of the catastrophic earthquakes are recorded around the trench. The histograms shows variation in depths along the segments of the Sumatra and Java. The Java segment has a bell-shaped data distribution in contrast to the Sumatra with bimodal pattern. The Java segment has the most repetitive depths at -2,500 to -5,200 m. The Sumatra transect has two peaks: 1) a classic bell-shaped peak at depths -4,500 m to -5,500 m; 2) shelf area with a peak from 0 to -1,750 m. The data at middle depths (-1,750 to -4,500 m) have a frequency -6,000 m, there are only 138 samples for the Sumatra while 547 samples for Java. Furthermore, Java segment has more symmetrical geometric shape while Sumatra segment is asymmetric, one-sided. The Sumatra segment has a steepness of 57.86° on its eastern side (facing Sumatra Island) and a contrasting 14.58° on the western part. The Java segment has a steepness of 64.34° on its northern side (facing Java Island) and 24.95° on the southern part (facing Indian Ocean). The paper contributes to the studies of the submarine geomorphology in Indonesia.

International audience; GRASS GIS Modules for Topographic and Geophysical Analysis of the ETOPO1 DEM and Raster Data: North Fiji Basin, Pacific Ocean The paper presents topographic analysis based on the raster data (NetCDF and grid formats) and visualization of the geophysical datasets that resembles the topographic surface. Geographically, the research focuses on the region of North Fiji Basin, South Pacific Ocean. North Fiji Basin is one of the marginal basins located at the converging boundary between the Pacific and Australian tectonic plates and is notable for complex geological settings. Methodology is based on GRASS GIS and includes scripting and cartographic visualization. Data include ETOPO1, GEBCO, EGM96 gravity and geoid raster grids imported to GRASS GIS via GDAL library ('r.in.gdal' module) from NetCDF and GRD formats, evaluated and visualized. Several GRASS modules were used as a sequential scripting for displaying and modelling data. Geomorphometric analysis (slopes, curvature, aspect, elevation) was performed by module 'r.slope.aspect' based on ETOPO1 grid. Statistical data analysis (histograms, polar diagram) for topographic and geoid gravitational fields were visualized by combination of GRASS modules d.rast, g.region, d.polar, d.histogram, d.grid, d.legend, d.text. Topographic analysis reveal that free-air gravity anomaly mark the continental shelf-slope border showing that generally, geophysical settings are dominated by the topographic effect. These observations are consistent with the inferred prominent role of the geophysical settings. The histograms reveal that topography (trenches, abyssal plains, shelf areas) has values between-5,000 and-1,000 m, and a clear peak with major values between 100 and 350 m above sea level. This distribution indicates that dominating depths of the ocean floor in North Fiji Basin are between-1,500 to-3,000 m. The aspect map of the ETOPO1 initial raw raster grid show that the surface has mostly 50°-110° oriented slope followed by data range in 250°-300° oriented slopes which is well correlating with submarine topography and general relief forms. The steepening of the gravity values between 50 and 55 mGal and 58 to 63 mGal indicates geoid undulations suggesting that bathymetry is controlled by the geological settings in the study area. Consolebased scripts of GRASS GIS are presented together with cartographic outputs showing its effectiveness for the geographic spatial analysis.

International audience; Insights on the Indian Ocean Tectonics and Geophysics Supported by GMT. This paper presented analyzed and summarized data on geological and geophysical settings about the tectonics and geological structure of the seafloor of the Indian Ocean by thematic visualization of the topographic, geophysical and geological data. The seafloor topography of the Indian Ocean is very complex which includes underwater hills, isolated mountains, underwater canyons, abyssal and accumulative plains, trenches. Complex geological settings explain seismic activity, repetitive earthquakes, and tsunami. Understanding and prognosis of the disastrous and catastrophic geological events is strongly based on correct data analysis, modelling and visualization. An important feature of this paper is mapping multi-source high-resolution data by GMT. Data include raster grids in NetCDF and GRD formats: ETOPO1, geologic and marine free-air gravity data, EGM96, age, spreading rates, and spreading asymmetry of the ocean crust by NOAA, total sediment thickness. Data were visualized by GMT modules to compare and analyze geophysical and geological settings of the Indian Ocean. Visualization revealed correlations between high bathymetric variations of the oceanic seafloor, distribution of main geological seafloor fabric: Southwest, Southeast, Mid and Carlsberg ridges. Tectonic maps were plotted to perform comparative analysis of several variables: crust age, spreading half rates (mm/yr), asymmetries in crustal accretion on conjugate ridge flanks (%), variations in the geopotential and gravimetric models. Being the warmest of the world's ocean, Indian Ocean has specific climatic conditions (repetitive monsoons, tsunamis, cyclones and storms), complex geologic seafloor structure with triple junction and unique geographic settings. Presented paper contributed to the regional studies of the Indian Ocean.

This article analyzes topographical and geological settings in the Japan Archipelago for comparative raster data processing using GRASS GIS. Data include bathymetric and geological grids in NetCDF format: GEBCO, EMAG2, GlobSed, marine free‐air gravity anomaly and EGM96. Data were imported to GRASS by gdalwarp utility of GDAL and projected via PROJ library. Method‐ ology includes data processing (projecting and import), mapping and spatial analysis. Visualization was done by shell scripting using a sequence of GRASS modules: ‘d.shade’ for relief mapping, ‘r.slope.aspect’; for modelling based on DEM, ‘r.contour’ for plotting isolines, ‘r.mapcalc’ for classification, ‘r.category’ for associating labels, and auxiliary modules (d.vect, d.rast, d.grid, d.legend). The results of the geophysical visualization show that marine free‐air gravitational anomalies vary in the Sea of Japan (–30 to above 40 mGal) reflecting density inhomogeneities of the tectonic structure, and correlating with the geo‐ logical structure of the seafloor. Dominating values of geoid model are 30–45 m reflecting West Pacific rise, determined by deep density inhomogeneities associated with the mantle convention. Sediment thickness varies across the sea reflecting its geological development with density of 2 km in its deepest part and thinner in central part (Yamato Rise). The aspect map and reclassified map express gradient of the steepest descent.

The paper aims to evaluate the presence and condition of vegetation by SAGA GIS. The study area covers northern coasts of Iceland including two fjords, the Eyjafjörður and the Skagafjörður, prosperous agricultural regions. The vegetation coverage in Iceland experience the impact of harsh climate, land use, livestock grazing, glacial ablation and volcanism. The data include the Landsat TM image. The methodology is based on computing raster bands for simulating Tassel Cap Transformation (wetness, greenness and brightness) and Enhanced Vegetation Index (EVI) sensitive to high biomass. The results include modelled three bands of brightness, greenness and wetness. Greenness variation shows the least values in ice-covered areas (-56.98 to -18.69). High values (-23.48 to 9.12) are in the valleys with dense vegetation, correlating with the geomorphology of the river network, the vegetation-free areas and ocean which corresponds to the peak of 30.87 to 41.19. The bell-shaped data distribution shows frequency 43.19–141.74 for vegetation indicating healthy state and canopy density. Maximal values are in ice-covered regions and glaciers (64°N-65°N). Very low values (0 to -20) show desertification and mountainous rocks. Moderate values (20-40) indicate healthy vegetation. The most frequent data: -28,17 to 11,8. The EVI shows data variations (-0.14 to 0.04). The study contributes both to the regional studies of Arctic Iceland and methodological approach of remote sensing data processing by SAGA GIS.

The study investigated geomorphology of the Japan Trench located east of Japan, Pacific Ocean. A high-resolution GEBCO Gridded Bathymetric Dataset was used for modeling, mapping and visualization. The study aimed to compare and analyse variations in the geomorphic structures of the two parts of the trench and to visualize variations in the geological, geophysical and bathymetric settings. Technically, the cartographic work was performed using scripting based on the Generic Mapping Toolset (GMT). Modelled cross-sectioning orthogonal profiles transecting the trench in a perpendicular direction were automatically digitized and graphed in the two segments. The results of the bathymetric analysis shown that the southern part is shallower: with deeper values in absolute (139 samples between –7000 to –8000 m) and statistical records (the most frequent values are within –5500 to –5800 m) comparing to the northern segment (–5300 to –5500 m). The geomorphological analysis shows a more complicated relief in the northern part of the trench, which has a higher seismic activity. The southern part has a gentler slope on the Honshu island side. The geoid modeling along the trench ranges in 0–20 mGal. The highest values are recorded by the Honshu Island (>40 mGal). The rest of the area has rather moderate undulations (20–40 mGal). The free-air marine gravity of the Sea of Japan is

The research problem is about to generate artificial fractal landscape surfaces from the Digital Elevation Model (DEM) using a stochastic algorithm by Geographic Resources Analysis Support System Geographic Information System (GRASS GIS) software. Fractal surfaces resemble appearance of natural topographic terrain and its structure using random surface modelling. Study area covers Kuril-Kamchatka region, Sea of Okhotsk, North Pacific Ocean. Techniques were included into GRASS GIS modules (r.relief, d.rast, r.slope.aspect, r.mapcalc) for raster calculation, processing and visualization. Module 'r.surf.fractal' was applied for generating synthetic fractal surface from ETOPO1 DEM GeoTIFF using algorithm of fractal analysis. Three tested dimensions of the fractal surfaces were automatically mapped and visualized. Algorithm of the automated fractal DEM modelling visualized variations in steepness and aspect of the artificially generated slopes in the mountains. Controllable topographic variation of the fractal surfaces was applied for three dimensions: dim=2.0001, 2.0050, 2.0100. Auxiliary modules were used for the visualization of DEMs (d.rast, r.colors, d.vect, r.contour, d.redraw, d.mon). Modules 'r.surf.gauss' and 'r.surf.random' were applied for artificial modelling as Gauss and random based mathematical surfaces, respectively. Univariate statistics for fractal surfaces were computed for comparative analysis of maps representing continuous fields by module 'r.univar': number of cells, min/max, range, mean, variance, standard deviation, variation coefficient and sum. The paper includes 9 maps and GRASS GIS codes used for visualization., Araştırma problemi, GRASS GIS yazılımı ile stokastik bir algoritma kullanılarak Sayısal Yükseklik Modeli'nden (SYM) yapay fraktal yüzeylerin üretilmesidir. Fraktal yüzeyler, doğal topografik arazinin görünümüne ve yapısına rastgele yüzey modellemesi kullanarak benzerler. Çalışma alanı Kuril-Kamçatka bölgesini, Okhotsk Denizi'ni, Kuzey Pasifik Okyanusu'nu kapsamaktadır. Raster hesaplama, işleme ve görselleştirme için kullanılan yöntemler GRASS GIS modüllerini (r.relief, d.rast, r.slope.aspect, r.mapcalc) içermektedir. Fraktal analiz algoritması kullanılarak ETOPO1 DEM GeoTIFF'den sentetik fraktal yüzey oluşturmak için 'r.surf.fractal' modülü uygulanmıştır. Fraktal yüzeylerin test edilen üç boyutu otomatik olarak haritalanmış ve görselleştirilmiştir. Otomatik fraktal DEM modellemesinin algoritması kullanılarak dağlık alanlarda yapay olarak üretilen yamaçların dikliği ve yönü bakımından oluşturulan varyasyonlar ile görselleştirmeler yapılmıştır. Fraktal yüzeylerin kontrol edilebilir topografik varyasyonu üç boyut için uygulanmıştır: dim = 2.0001, 2.0050, 2.0100. DEM'lerin görüntülenmesi için yardımcı modüller kullanılmıştır (d.rast, r.colors, d.vect, r.contour, d.redraw, d.mon). Yapay modelleme için 'r.surf.gauss' ve 'r.surf.random' modülleri Gauss ve rasgele tabanlı matematiksel yüzeyler olmak üzere sırasıyla uygulanmıştır. Fraktal yüzeyler için tek değişkenli istatistikler 'r.univar' modülüne göre sürekli alanları temsil eden haritaların karşılaştırmalı analizi için hesaplanmıştır: hücre sayısı, min / maks, aralık, ortalama, varyans, standart sapma, varyasyon katsayısı ve toplam. Makalede 9 harita ve görselleştirme için kullanılan GRASS GIS kodları bulunmaktadır.

Computer-generated imagery (CGI) takes a growing place in our everyday environment. Whether it is in video games or movies, CGI techniques are constantly improving in quality but also require ever more qualitative artistic content which takes a growing time to create. With the emergence of virtual and augmented reality, often comes the need to render or re-render assets that exist in our world. To allow widespread use of CGI in applications such as telepresence or virtual visits, the need for manual artistic replication of assets must be removed from the process. This can be done with the help of Image-Based Rendering (IBR) techniques that allow scenes or objects to be rendered in a free-viewpoint manner from a set of sparse input photographs. While this process requires little to no artistic work, it also does not allow for artistic control or editing of scene content.In this dissertation, we explore Multi-view Image Editing and Rendering. To allow casually captured scenes to be rendered with content alterations such as object removal, lighting edition, or scene compositing, we leverage the use of optimization techniques and modern deep-learning. We design our methods to take advantage of all the information present in multi-view content while handling specific constraints such as multi-view coherency.For object removal, we introduce a new plane-based multi-view inpainting algorithm. Planes are a simple yet effective way to fill geometry and they naturally enforce multi-view coherency as inpainting is computed in a shared rectified texture space, allowing us to correctly respect perspective. We demonstrate instance-based object removal at the scale of a street in scenes composed of several hundreds of images.We next address outdoor relighting with a learning-based algorithm that efficiently allows the illumination in a scene to be changed, while removing and synthesizing cast shadows for any given sun position and accounting for global illumination. An approximate geometric proxy built using multi-view stereo is used to generate illumination and shadow related image buffers that guide a neural network. We train this network on a set of synthetic scenes allowing full supervision of the learning pipeline. Careful data augmentation allows our network to transfer to real scenes and provides state of the art relighting results. We also demonstrate the capacity of this network to be used to compose real scenes captured under different lighting conditions and orientation.We then present contributions to image-based rendering quality. We discuss how our carefully designed depth-map meshing and simplification algorithm improve rendering performance and quality of a new learning-based IBR method.Finally, we present a method that combines relighting, IBR, and material analysis. To enable relightable IBR with accurate glossy effects, we extract both material appearance variations and qualitative texture information from multi-view content in the form of several IBR heuristics. We further combine them with path-traced irradiance images that specify the input and target lighting. This combination allows a neural network to be trained to implicitly extract material properties and produce realistic-looking relit viewpoints. Separating diffuse and specular supervision is crucial in obtaining high-quality output., Les images de synthèse (CGI) prennent une place grandissante dans notre environnement. Que ce soit dans les jeux vidéos ou les films, leur qualité ne cesse de s’accroître nécessitant la création fastidieuse de contenus artistiques. L’émergence de la réalité virtuelle et augmentée, entraine la nécessité de rendre des environnements existants. Pour permettre l’utilisation généralisée des images de synthèse dans des applications telles que la télé-présence ou les visites virtuelles, la digitalisation manuelle des contenus par des artistes se doit d’être évitée. Une des solutions peut provenir des techniques de Rendu à Base d’Images (IBR) qui permettent de rendre des scènes, depuis un point de vue libre, à partir d’un ensemble de photographies parcimonieux. Bien que ces méthodes ne nécessitent que peu de travail artistique, elles n’autorisent cependant pas le contrôle ou l’édition du contenu.Dans cette thèse, nous explorons l’Edition et le Rendu d’Images Multi-vues. Afin de permettre à des scènes, capturées avec le moins de contraintes possibles, d’être rendues avec des altérations telles que la suppression d’objets, l’édition d’éclairage, ou la composition de scènes, nous exploitons les techniques d’optimisation et d’apprentissage profond. Nous concevons nos méthodes afin qu’elles tirent pleinement avantage de l’information présente dans le contenu multi-vues, tout en en respectant ses contraintes spécifiques.Pour la suppression d’objets, nous introduisons un algorithme de remplissage automatique, multi-vues cohérent, utilisant une représentation planaire. Les plans sont des objets simples et efficaces pour combler la géométrie, dont la cohérence multi-vues émerge naturellement lorsque le remplissage est effectué dans un espace texture rectifié et partagé. Ils permettent aussi le respect des effets de perspective. Nous démontrons la capacité d’enlever des objets, à grande l’échelle, dans des scènes contenant plusieurs centaines d’images.Nous traitons ensuite le problème du rééclairage des scènes extérieures par une méthode d’apprentissage profond. Elle permet de modifier l’illumination, en enlevant et synthétisant les ombres portées, pour une position du soleil quelconque, tout en tenant compte des variations d’illumination globale. Une représentation géométrique approximative, reconstruite en utilisant la stéréo multi-vues, est utilisée pour générer des images tampons d’illumination et d’ombres qui guident un réseau de neurones. Nous entrainons ce réseau sur un ensemble de scènes synthétiques, permettant une supervision complète. Une augmentation des données minutieuse permet à notre réseau de généraliser aux scènes réelles et de produire l’état de l’art en terme de résultats. Nous démontrons ensuite, la capacité du réseau à être utilisé pour composer des scènes réelles, capturées dans des conditions d’orientation et d’éclairages différentes. Nous présentons ensuite des contributions à la qualité de l'IBR. Nous introduisons un algorithme de maillage de cartes de profondeur et de leur simplification. Nous démontrons son impact sur la qualité et les performances d’une nouvelle méthode d’IBR utilisant l’apprentissage.Enfin, nous introduisons une méthode qui combine rééclairage, IBR, et analyse de matériaux. Afin de permettre un rendu à base d’images, rééclairable et tenant compte des effets spéculaires, nous extrayons du contenu multi-vues les variations d’apparence des matériaux et l’information de texture haute résolution, sous la forme de plusieurs rendus IBR heuristiques. Nous les combinons ensuite avec des rendus d’irradiance, obtenus par lancer de rayons, qui spécifient les conditions d’éclairage initiales et désirées. Cette combinaison permet d’entrainer un réseau de neurones à extraire implicitement les propriétés des matériaux et à produire des points de vue rééclairés réalistes [...]

International audience; This work aim is to perform modelling and spatial analysis of the marine geological data using combination of the QGIS and Python programming. Selecting proper cartographic software is important part of the geospatial research. QGIS provides organizing data in a GIS project for mapping and spatial visualization through vector and raster layers stored in GIS. Study area is Mariana Trench, west Pacific Ocean. A series of cross-section profiles were digitized in QGIS and used for further data processing in Python. Mariana Trench has complex geomorphic structure and unevenness in profiles stretching south-westwards. The geomorphology is subjected to various phenomena that affect its shape. These include bathymetry, geodesy, gravimetry, tectonics plates and geological settings, studied in this paper. To understand the structure of the trench, a data modelling using bathymetric analysis was performed by combination of QGIS mapping and statistical analysis in Python's library Seaborn. Statistical data modelling aimed at the analysis of the spatial variation of the geomorphology of the trench using following methods: multiple facet grids, area charts for the data frames, regression analysis, letter-value plots, hexagonal and Kernel density estimation. The results of the geospatial data analysis show spatial unevenness of the geomorphic structure, gravimetric, geodetic and bathymetric settings of the Mariana Trench. The study demonstrated effectiveness of Python application in geographic data analysis with Python codes provided for repeatability.

The paper presents a comparison of the two languages Python and R related to the classification tools and demonstrates the differences in their syntax and graphical output. It indicates the functionality of R and Python packages {dendextend} and scipy.cluster as effective tools for the dendrogram modelling by the algorithms of sorting and ranking datasets. R and Python programming languages have been tested on a sample dataset including marine geological measurements. The work aims to detect how bathymetric data change along the 25 bathymetric profiles digitized across the Mariana Trench. The methodology includes performed hierarchical cluster analysis with dendrograms and plotted clustermap with marginal dendrograms. The statistical libraries include Matplotlib, SciPy, NumPy, Pandas by Python and {dendextend}, {pvclust}, {magrittr} by R. The dendrograms were compared by the model-simulated clusters of the bathymetric ranges. The results show three distinct groups of the profiles sorted by the elevation ranges with maximal depths detected in a group of profiles 19-21. The dendrogram visualization in a cluster analysis demonstrates the effective representation of the data sorting, grouping and classifying by the machine learning algorithms. The programming codes presented in this study enable to sort a dataset in a similar research aimed to group data based on the similarity of attributes. Effective visualization by dendrograms is a useful modelling tool for the geospatial management where data ranking is required. Plotting dendrograms by R, comparing to Python, presented functional and sophisticated algorithms, refined design control and fine graphical data output. The interdisciplinary nature of this work consists in application of the coding algorithms for spatial data analysis.