Skip to content

a script that makes error ellipse plots of ensemble atcfunix files#16

Open
SamuelTrahanNOAA wants to merge 4 commits into
dtcenter:developfrom
SamuelTrahanNOAA:ellipses
Open

a script that makes error ellipse plots of ensemble atcfunix files#16
SamuelTrahanNOAA wants to merge 4 commits into
dtcenter:developfrom
SamuelTrahanNOAA:ellipses

Conversation

@SamuelTrahanNOAA

Copy link
Copy Markdown
Contributor

An ensemble error ellipse script with no proper home yet. The script will probably end up somewhere else, but at least this PR gives us a forum for code review and discussion.

Given an ensemble hurricane forecast, this script plots error ellipses. The ellipses are plotted at one standard deviation, tilted to match the spread of the dataset. It uses the eigendecomposition of the covariance matrix of lat & lon differences from the mean location. (See references.) Ellipses are plotted every 24 hours by default, but this can be changed.

I expect that, once this is put inside some larger package, the atcfunix reader and data structures will be replaced with whatever the rest of that package uses.

References:

Idea came from here:

How to do an error ellipse plot:

ellipse-revised

@SamuelTrahanNOAA

Copy link
Copy Markdown
Contributor Author

@KathrynNewman @CPKalb @mkavulich @willmayfield - This is the draft PR I mentioned with the error ellipse plotting script.

@mkavulich

Copy link
Copy Markdown
Collaborator

@SamuelTrahanNOAA This looks good to me, at least as a temporary home, as we discussed last meeting. If you mark this PR ready for review I'm happy to approve it. If you'd like a more thorough review of the new plotting script I can also test it out.

@SamuelTrahanNOAA

Copy link
Copy Markdown
Contributor Author

Someone other than me should test the script.

@SamuelTrahanNOAA
SamuelTrahanNOAA marked this pull request as ready for review May 28, 2026 21:10
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants