Add DGT Lidar tutorial in portuguese

Author: joa-quim

_extensions/gmtref/gmtref.lua

@@ -24,7 +24,7 @@ return {
     local utilities = {
       "ablines", "anaglyph", "append2fig", "autocov", "autocor", "blendimg", "cart2pol", "cart2sph",
       "coastlinesproj", "colorzones", "cpt4dcw", "crop", "cube", "cubeplot", "cylinder",
-      "date2doy", "delrows", "df2ds", "dodecahedron", "doy2date", "ecmwf", "ellipse3D",
+      "date2doy", "delrows", "df2ds", "dgt_lidar", "dgt_mosaic", "dodecahedron", "doy2date", "ecmwf", "ellipse3D",
       "era5time", "era5vars", "extrude", "findpeaks", "flatfv", "fv2fv",
       "gadm", "geocoder", "geodetic2enu", "getbyattrib", "gmtread", "gmtwrite", "graticules",
       "grid2tri", "gridit", "gunique", "hampel", "hlines", "icosahedron",

_quarto.yml

@@ -586,7 +586,8 @@ project:
     - tutorials/Landsat8/histogram_stretch.qmd
     - tutorials/longest_sail/longestsail.qmd
     - tutorials/mosaics/mosaics.qmd
-    - tutorials/dgt_lidar/dgt_lidar.qmd
+    - tutorials/dgt_lidar/tut_dgt_lidar.qmd
+    - tutorials/dgt_lidar/tut_dgt_lidar_pt.qmd
     - tutorials/remote_datasets/remote_datasets.qmd
     - tutorials/signal/signal.qmd
     - tutorials/signal/autocorr_chloro/autocorr_chloro.qmd
@@ -792,7 +793,7 @@ website:
           - text: "Mosaics"
             file: tutorials/mosaics/mosaics.html
           - text: "DGT LIDAR Portugal"
-            file: tutorials/dgt_lidar/dgt_lidar.html
+            file: tutorials/dgt_lidar/tut_dgt_lidar_pt.html
           - text: "Insets"
             file: tutorials/insets/insets.html
           - text: "Remote Grids"

documentation/modules/gravprisms.md

@@ -25,7 +25,7 @@ gravity gradient anomalies, or geoid anomalies. Options are available to control
 axes units and direction.
 
 ~~~
-<img src="/assets/GMT_seamount_prisms.png" alt="Simple" width="500" class="center"/>
+<img src="/assets/GMT_seamount_prisms.png" alt="Simple" width="600" class="center"/>
 ~~~
 
    Three density models modeled by prisms for a truncated Gaussian seamount via **C**: (left) Constant

documentation/utilities/dgt_lidar.qmd

@@ -164,7 +164,7 @@ using GMT, HTTP
 dgt_lidar(-7.607, 40.322; neighbors=1, mosaic="estrela.tiff")
 ```
 
-For a complete worked example see the [DGT LIDAR tutorial](../../tutorials/dgt_lidar/dgt_lidar.qmd).
+For a complete worked example see the [DGT LIDAR tutorial](../../tutorials/dgt_lidar/tut_dgt_lidar_pt.qmd).
 
 ```{julia}
 #| echo: false

documentation/utilities/dgt_mosaic.qmd

@@ -101,4 +101,4 @@ See Also
 
 {{< gmtref dgt_lidar >}}, {{< gmtref mosaic >}}
 
-For a complete worked example see the [DGT LIDAR tutorial](../../tutorials/dgt_lidar/dgt_lidar.qmd).
+For a complete worked example see the [DGT LIDAR tutorial](../../tutorials/dgt_lidar/tut_dgt_lidar_pt.qmd).

documentation/utilities/mat2img.qmd

@@ -27,14 +27,14 @@ If `stretch` is a scalar, scale the values > `stretch` to [0 255]
 The `kw...` kwargs search for [:layout :mem_layout], [:names] and [:metadata]
 
 
-See Also
---------
-
-{{< gmtref mat2ds >}}, {{< gmtref mat2grid >}}
-
 ```{julia}
 #| echo: false
 #| output: asis
 using GMT
 GMT.doc_source_links("mat2img")
 ```
+
+See Also
+--------
+
+{{< gmtref mat2ds >}}, {{< gmtref mat2grid >}}

tutorials/dgt_lidar/dgt_lidar.qmd tutorials/dgt_lidar/tut_dgt_lidar.qmdtutorials/dgt_lidar/dgt_lidar.qmd renamed to tutorials/dgt_lidar/tut_dgt_lidar.qmd

@@ -5,8 +5,8 @@ provides free high-resolution LIDAR-derived elevation data covering the entire P
 through its [Collaborative Data Distribution (CDD)](https://cdd.dgterritorio.gov.pt/) portal.
 GMT.jl provides two functions for working with this dataset:
 
-- **`dgt_lidar`** — authenticates against the DGT portal and downloads all tiles intersecting a bounding box.
-- **`dgt_mosaic`** — assembles downloaded tiles into a single, clipped GeoTIFF.
+- {{< gmtref dgt_lidar >}} — authenticates against the DGT portal and downloads all tiles intersecting a bounding box.
+- {{< gmtref dgt_mosaic >}} — assembles downloaded tiles into a single, clipped GeoTIFF.
 
 Five collections are available:
 
@@ -53,7 +53,7 @@ You can also create or edit the file directly with any text editor.
 
 ### Load the extension
 
-`dgt_lidar` lives in the `GMTDGTLidarExt` package extension, activated automatically when `HTTP` is loaded:
+{{< gmtref dgt_lidar >}} lives in the `GMTDGTLidarExt` package extension, activated automatically when `HTTP` is loaded:
 
 ```{julia}
 using GMT, HTTP
@@ -408,6 +408,22 @@ dgt_lidar([-8.995, -8.968, 38.473, 38.486]; collection="LAZ", output_dir="_arrab
 pc = lasread(joinpath(homedir(), ".gmt/DGT/arrabida_laz/LAZ/LO-148174-06-2024.laz"))
 ```
 
+And we can make a little detour and drape a Bing aereal image fetch with mosaic on top of the DEM:
+
+```{julia}
+#| eval: false
+# Get the DEM as a grid in memory, reprojected to geographic coordinates
+G = dgt_lidar([-8.995, -8.968, 38.473, 38.486]; output_dir = "_arrabida", mosaic="grid", proj="geog");
+
+# Use its limits to fetch a Bing image of the same area. We use a zoom=16 to get a approximate  resolution to G
+I = mosaic(G, zoom=16);
+
+# Visualize the draped result
+viz(G, drape=I, view=(155,35), zsize=3)
+```
+
+![](/assets/portinho_draped.png){width="700" .center}
+
 
 ## Large Areas — Automatic Subdivision