Along track views of height change over time #99
Merged
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Towards making transect/cross-section plots! Ensure that our ATL11 xarray dataset has the referencegroundtrack number in it, for easier tracking later on! Builds on f30110c. I've updated the Kamb subset for my friend, and it's good to see that hvplot works interactively when zooming in now if datashade is on (but rasterize shows up blank?). Ideally we would be using intake.open_ndzarr for this, but there are too many missing pieces in the puzzle to get it to work seamlessly now. So I've reverse engineered the add_path_to_ds mechanism and used it to extract the reference ground track number from the filename when loading from the Zarr files! The coordinates are stored as int64 which seems overkill (int8 would suffice), but that can be amended properly at a later date, hopefully upstream at intake/xarray.
|
Check out this pull request on Review Jupyter notebook visual diffs & provide feedback on notebooks. Powered by ReviewNB |
Getting closer to detecting subglacial lake drainage/filling events! With the referencegroundtrack variable available with efe93e0, and interactive hvplot working again (somewhat), we can finally do some fun exploratory data analysis on examining active subglacial lake movements in Antarctica. Made a small tweak to the `add_path_to_ds` format_string to exclude version and revision numbers, since they add little value. Showed an example at Whillans Ice Stream (EPSG:3031, x: -444300, y: -537500) where a ~7 metre filling event occured. Probably spent a bit too much time perfecting the PyGMT plot, complete with the ReferenceGroundTrack number, legends for each cycle colour point, mean datetimes for those points, etc. Can't wait to go 3D soon!
Lots of updates across the PyViz stack to support Bokeh 2.1. HvPlot interactivity works great now too!
Made a cool ICESat-2 rate of height change over time (dhdt) dashboard using Panel and HvPlot! Has buttons to select variable to plot (dhdt_slope, referencegroundtrack or h_corr), sliders to filter datapoints by cycle_number and dhdt_range, and checkboxes to turn on/off rasterize/datashade! Actually my first time using Panel really, but quite pleased with the results.
weiji14
added a commit
that referenced
this pull request
Jul 16, 2020
Closes #99 Along track views of height change over time.
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Towards making transect/cross-section plots, along a single reference ground track! For example, being able to make a scatter plot with
x_atcon the x-axis andh_corron the y-axis, and overlay differentcycleson top of one another.Also made an interactive dashboard map view plot to locate interesting reference ground tracks:
TODO:
References:
Notes:
Currently using some 'hacks' to get this work, until upstream libraries like intake-xarray, xarray, zarr, et al. release some features to improve PyData interoperability (see #100).