GeoTIFF RGB

This example demonstrates how to render an OGC GeoTIFF RGB image as a geolocated mesh.

📋 Summary

Creates a mesh from a GeoTIFF encoded RGB image.

The resulting mesh contains quad cells.

The GeoTIFF RGB image is first pre-processed using rasterio.features.sieve() to remove several unwanted masked regions within the interior of the image, which are due to a lack of dynamic range in the uint8 image data.

Note that, the RGB image pixel data is located on the mesh nodes/points.

The masked RGB pixels are then removed by extracting only mesh cells with no masked points.

🏷 Tags: GeoTIFF Lighting RGB

Attention

Optional package dependency rasterio is required.

---

Static Scene

Interactive Scene

from __future__ import annotations

import geovista as gv
from geovista.pantry import fetch_raster
import geovista.theme


def main() -> None:
    """Plot a GeoTIFF RGB image.

    Notes
    -----
    .. versionadded:: 0.5.0

    """
    # Fetch the example GeoTIFF file.
    fname = fetch_raster("bahamas_rgb.tif")

    # Plot the RGB image.
    plotter = gv.GeoPlotter()

    # Load the GeoTIFF image, which requires the optional package
    # dependency 'rasterio'.
    mesh = gv.Transform.from_tiff(fname, rgb=True, sieve=True, extract=True)

    plotter.add_mesh(mesh, lighting=False, rgb=True)
    plotter.add_axes()
    plotter.add_text(
        "Bahamas, RGB GeoTIFF",
        position="upper_left",
        font_size=10,
        shadow=True,
    )
    plotter.view_xz()
    plotter.camera.zoom(1.3)
    plotter.show()


if __name__ == "__main__":
    main()