Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Plot tropical cyclone tracks from BUFR data #10

Open
sandorkertesz opened this issue Jun 28, 2024 · 0 comments
Open

Plot tropical cyclone tracks from BUFR data #10

sandorkertesz opened this issue Jun 28, 2024 · 0 comments
Labels
enhancement New feature or request

Comments

@sandorkertesz
Copy link
Collaborator

Is your feature request related to a problem? Please describe.

The BUFR file https://get.ecmwf.int/repository/test-data/earthkit-data/test-data/A_JSXX01ECEP251200_C_ECMP_20240625120000_tropical_cyclone_track_70W_132degE_6p4degN_bufr4.bin contains tropical cyclone track forecast for multiple ensemble members.

The following code extracts the minimum mslp locations for a given member into a pandas dataframe:

import earthkit.data

f_path = "A_JSXX01ECEP251200_C_ECMP_20240625120000_tropical_cyclone_track_70W_132degE_6p4degN_bufr4.bin"
ds = earthkit.data.from_source("file", f_path)
df = ds.to_pandas(columns=("latitude", "longitude", "pressureReducedToMeanSeaLevel"),
    filters={"ensembleMemberNumber": 2, "meteorologicalAttributeSignificance": [1,4]},
    required_columns=True)

df = df.dropna()
   latitude  longitude  pressureReducedToMeanSeaLevel
0       6.7      132.8                       100700.0
1       5.5      131.2                       100600.0
2       6.5      130.0                       100600.0
3       8.7      128.3                       100400.0
6      11.9      122.9                       100900.0
7      13.0      121.8                       100800.0

The task is to generate a plot containing a singe track or multiple tracks.

The data also contains the maximum 10m wind locations, which can be extracted with the following code:

df = ds.to_pandas(columns=("latitude", "longitude", "windSpeedAt10M"),
    filters={"ensembleMemberNumber": 2, "meteorologicalAttributeSignificance": 3},
             required_columns=True)

df = df.dropna()

Describe the solution you'd like

No response

Describe alternatives you've considered

No response

Additional context

No response

Organisation

ECMWF

@sandorkertesz sandorkertesz added the enhancement New feature or request label Jun 28, 2024
@sandorkertesz sandorkertesz changed the title Plot tropical cyclon tracks from BUFR data Plot tropical cyclone tracks from BUFR data Jun 28, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

1 participant