Skip to content

Problem in eddy identification with py-eddy-tracker #274

@shivashankar11

Description

@shivashankar11

I have installed pyeddytracker-3.6.1 in Ubuntu. When I tried to identify eddies from SLA data using the below code, I could detect only 2 eddies (shown in the attached figure). However, with an independent okubo-weiss parameter-based eddy detection and tracking algorithm, I could found reasonably good number of eddies which are in consistent with the background SLA. Kindly help me out where I am committing mistake in eddy identification also guide me for eddy tracking once identification is successful.

import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.io.shapereader as shpreader
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter
from py_eddy_tracker.dataset.grid import RegularGridDataset
from datetime import datetime

NETCDF_FILE = "sla_20190101_meters.nc"
LON_NAME = "LONGITUDE50_257"
LAT_NAME = "LATITUDE85_190"
DATA_VAR = "SLA_m" # in m

g = RegularGridDataset(NETCDF_FILE, LON_NAME, LAT_NAME)
g.add_uv(DATA_VAR)

EXTENT = [47, 99, 6, 27.5] # whole study area

g.bessel_high_filter(DATA_VAR, 1000, order=1)
a, c = g.eddy_identification(
DATA_VAR, "u", "v",
date,
step=0.002, # 2 mm step
pixel_limit=(4, 8000),
shape_error=70, # META3.2 standard
nb_step_min=2,
nb_step_to_be_mle=0,
sampling=50,
sampling_method="uniform"
)

g.contours.display(ax, step=4)
update_axes(fig, ax)
plt.show()

fig.savefig("figure.png", dpi=300, bbox_inches="tight")

Image

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions