[cwhite911]

Description

Adds a raster profile tool to grass.jupyter.

gj.Profile samples one or more rasters along a line with r.profile and returns the data. Its figure has a map panel with the line and numbered vertices, and a chart with matching markers. Nulls become gaps in the line.

Motivation and context

The only profile UI in GRASS is in the wxGUI. d.profile draws a single raster transect as an image and returns no data. In a notebook there was no way to draw a line on the map and get the values back, or to produce a figure showing where along the line each value comes from.

Examples

profile = gj.Profile(["elevation", "elev_lid792_1m"], [(638000, 220000), (640000, 222000)])
profile.data["elevation"][:2]
# [{'easting': 638000.0, 'northing': 220000.0, 'distance': 0, 'value': 109.4}, ...]
profile.vertices        # [{'easting': ..., 'northing': ..., 'distance': ...}, ...]
profile.to_dataframe()  # optional, needs pandas
profile.show()
profile.save("profile.png")

`InteractiveMap` gets a profile button (ipyleaflet backend). Draw a line, pick the rasters, and the figure appears below the map:

m = gj.InteractiveMap()
m.add_raster("elevation")
m.show()
# after drawing a line with the profile tool:
m.profile.data
m.profile.save("profile.png")

Behavior Change

ipyleaflet, InteractiveMap.show() returns a VBox with the map and an output area below it, instead of the bare ipyleaflet.Map. The map widget is still InteractiveMap.map. Folium is unchanged.

Parts of profile.py were written with AI assistance and reviewed by me.

AI Assisted