"""
Count pixel used
======================
"""
from matplotlib import pyplot as plt
from matplotlib.colors import LogNorm
from py_eddy_tracker.observations.tracking import TrackEddiesObservations
import py_eddy_tracker_sample
# %%
# Load an experimental med atlas over a period of 26 years (1993-2019)
a = TrackEddiesObservations.load_file(
py_eddy_tracker_sample.get_path(
"eddies_med_adt_allsat_dt2018/Anticyclonic.zarr"))
c = TrackEddiesObservations.load_file(
py_eddy_tracker_sample.get_path(
"eddies_med_adt_allsat_dt2018/Cyclonic.zarr"))
t0, t1 = a.period
step = 0.1
bins = ((-10, 37, step), (30, 46, step))
kwargs_pcolormesh = dict(cmap="terrain_r",
vmin=0,
vmax=0.75,
factor=1 / (t1 - t0),
name="count")
# %%
# Plot
fig = plt.figure(figsize=(12, 18.5))
ax_a = fig.add_axes([0.03, 0.75, 0.90, 0.25])
if event.button == "down":
self.settings["now"] -= 1
self.update()
def adjust(self, event=None):
self.map._pan_start = None
self.map.end_pan()
def show(self):
self.update()
plt.show()
if __name__ == "__main__":
# a_ = A.load_file(
# "/home/toto/dev/work/pet/20200611_example_dataset/tracking/Anticyclonic_track_too_short.nc"
# )
# c_ = A.load_file(
# "/home/toto/dev/work/pet/20200611_example_dataset/tracking/Cyclonic_track_too_short.nc"
# )
a = A.load_file(
sample.get_path("eddies_med_adt_allsat_dt2018/Anticyclonic.zarr"))
# c = A.load_file(sample.get_path("eddies_med_adt_allsat_dt2018/Cyclonic.zarr"))
# g = GUI(Acyc=a, Cyc=c, Acyc_short=a_, Cyc_short=c_)
g = GUI(Acyc=a)
# g = GUI(Acyc_short=a_)
# g = GUI(Acyc_short=a_, Cyc_short=c_)
g.med()
g.show()
