Skip to content

Commit

Permalink
Change vertical CRS below to avoid failure
Browse files Browse the repository at this point in the history
  • Loading branch information
rhugonnet committed Dec 6, 2024
1 parent a08b73e commit aedb3ae
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions doc/source/cheatsheet.md
Original file line number Diff line number Diff line change
Expand Up @@ -25,7 +25,7 @@ identify on **a map of elevation differences with another elevation dataset (loo
## Cheatsheet

The patterns of errors categories listed in this spreadsheet **are linked to visual examples further below** that
you use to compare to your own elevation differences.
you can use to compare to your own elevation differences.

```{list-table}
:widths: 1 2 2 2
Expand Down Expand Up @@ -153,7 +153,7 @@ dh.plot(cmap='RdYlBu', vmin=-3, vmax=3, cbar_title="Elevation differences of\nho
(smooth-large-field)=
### Smooth-field offset

Example of smooth large offset field created by a wrong vertical CRS. We here show the difference due to a local
Example of smooth large offset field created by a wrong vertical CRS. We here show the difference due to the EGM96
geoid added on top of the ellipsoid.

```{code-cell} ipython3
Expand All @@ -162,10 +162,10 @@ geoid added on top of the ellipsoid.
: code_prompt_show: "Show code to simulate vertical referencing errors"
: code_prompt_hide: "Hide code to simulate vertical referencing errors"
# Set current vertical CRS
# Set current vertical CRS as ellipsoid
dem.set_vcrs("EGM96")
# Transform to a local reference system from https://cdn.proj.org/
trans_dem = dem.to_vcrs("no_kv_arcgp-2006-sk.tif")
# Transform vertical reference to geoid
trans_dem = dem.to_vcrs("Ellipsoid")
# Plot the elevation differences of the vertical transformation
dh = dem - trans_dem
Expand Down

0 comments on commit aedb3ae

Please sign in to comment.