Tools to measure glacier lengths statistically
Project description
glacier_lengths — Statistical glacier length calculations
Often when glacier lengths are calculated, only the glacier centerline is considered.
This is arguably not a statistically representative measure for the entire front, as it just considers one point on the glacier outline.
The glacier_lengths
package aims to simplify length calculations along an arbitrary amount of lines buffered around the glacier centerline.
Installation
pip install glacier_lengths
pip install git+https://github.com/erikmannerfelt/glacier_lengths.git
Documentation
See the documentation at https://glacier-lengths.readthedocs.io
Example
Calculate the length change of Rhonegletscher:
import geopandas as gpd
import glacier_lengths
# Read the example data
outlines = gpd.read_file("examples/rhone/data/outlines.shp").sort_values("year")
old_outline = outlines.iloc[0]
new_outline = outlines.iloc[1]
centerline = gpd.read_file("examples/rhone/data/centerline.shp").iloc[0]
# Generate ~40 buffered lines around the glacier centerline
old_buffered_lines = glacier_lengths.buffer_centerline(centerline.geometry, old_outline.geometry)
# Cut the newly generated lines to the new_outline
new_buffered_lines = glacier_lengths.cut_centerlines(old_buffered_lines, new_outline.geometry)
# Measure the lengths of the old and new glacier centerlines.
old_lengths = glacier_lengths.measure_lengths(old_buffered_lines)
new_lengths = glacier_lengths.measure_lengths(new_buffered_lines)
# Print the results.
print(f"""
{old_outline['year']}: {old_lengths.mean():.1f}±{old_lengths.std():.1f} m
{new_outline['year']}: {new_lengths.mean():.1f}±{new_lengths.std():.1f} m
""")
prints:
1928: 10783.6±38.8 m
2020: 9699.9±7.6 m
Plot a figure
python examples/rhone/plot_rhone.py
Testing
Run python -m pytest
in the cloned repo base directory.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for glacier_lengths-0.1.2-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | c254b044abae1283a585c0930577cf9402225eaefadf315ad5cd03e858356fce |
|
MD5 | 09ebf1d6f7131266e1eb0b1f82ac0598 |
|
BLAKE2b-256 | 530a3050e6c487d44a8b0004e3139a22ce6d15044605b80d74b6dc5e2a6e7a4f |