aymeric-spiga / planetoplot Goto Github PK

Le protocole OpeNDAP, m'a l'air tout a fait étudier pour.

http://docs.opendap.org/index.php/QuickStart

Avec la librairie PyDAP on peut faire tourner un serveur OpeNDAP sous le serveur Apache (1) et côté client on peut l'interroger en python gards ou ferret et notamment faire de la lecture on-the-fly (2)

1. https://pydap.readthedocs.io/en/latest/server.html#apache

2. https://publicwiki.deltares.nl/display/OET/Reading+data+from+OpenDAP+using+python

En couplant ça à PlanetoPlot y'a moyen de faire un truc pas mal.

• Pourquoi doit-on avoir un axe des temps qui croit de manière monotone
  ? J'ai des donnés instrumentales avec des gaps temporels qui ne sont pas
  présents explicitement dans le fichier netCDF. Forcer la monotonie de
  l'axe des temps empêche de visualiser correctement ces données avec
  planetoplot.

Quel est l'axe des temps que tu utilises ?

L'astuce monotone c'est juste une facçon de gagner du temps pour lire les fichiers avec beaucoup de sorties temporelles (sinon il est toujours un peu lent). Mais il n'y a pas d'autres raisons. Je peux changer le comportement par défaut. Par exemple pour Ls j'avais ajouté un correctls pour changer ça.

Ha oui l'option 'correctls' correspond à ce que je veux faire. Dans ce cas-ci j'utilise un axe Ls, mais avec des jours on s'en sort aussi avec 'correctlsnoadd', donc tout fonctionne. Peut-être on peut faire une option du genre 'notmonotonous' ou 'readtaxis' qui est un clone de 'correctlsnoadd' pour plus de clarté dans le cas où on utilise pas des Ls.

pour lire une variable 1D dans un fichier NetCDF, il faut préciser x=0, y=0.

suggestion par @jbmadeleine
http://eoimages.gsfc.nasa.gov/ve/2430/land_ocean_ice_2048.jpg

(question par S. Fromang)

Ma question est la suivante: je travaille dans ipython. Lorsque
j'appelle le script, je me recupère donc un objet (mettons "dpsi") en
retour de l'appel, en plus du plot. Si j'appelle la méthode plot() sur
cette objet, il me refait le plot. Mais si j'essaie de visualiser ce qui
est stocké dans dpsi.f, ca ne correspond pas. Voila la sequence:

dpsi=anomaly.anomaly(var='psi',tmin=0.,tmax=0.,tdiff=15.,zlevel=0)

dpsi.plot()

imshow(dspi.f)

C'est l'appel a imshow(dspi.f) qui pose problème. Pourquoi?

Du coup, je ne comprends pas ou est stockée l'info sur le champ que je
veux visualiser. Il s'en sort lors de l'appel a la méthode plot, mais
manifestement l'info n'est pas dans dpsi.f.

Some of the changes brought in matplotlib 2.0+ (https://matplotlib.org/api/api_changes.html) since last year, lead to crashes and incompatibility of planetoplot...

Especially "The axisbg and axis_bgcolor properties on Axes have been deprecated in favor of facecolor." is quite annoying ...

To be short, we cannot plot anything with latest versions of mpl as long as this one (at least) is not fixed ...

for background maps

replace basemap by cartopy proposed by @jbmadeleine

planetoplot could use xarray, which would make a lot more versatile and less complicated library

a fix is e.g. to set proj="cyl"

Quand j'affiche un champ dont les dimensions x et y ont la même
taille, il arrive que le champ affiché soit "transposé" par rapport à
ses axes. En effet, dans ppplot on ne cherche à transposer un champ que
lorsque les axes sont de tailles différentes (cf. ppplot.py, L 774), et
le cas particulier où shape[0] == shape[1] n'est pas pris en compte. Je
ne vois pas de solution élégante à ce problème, si ce n'est d'identifier
les noms des axes et le nom du champ à afficher, ce qui nécessite de
remanier un peu planetoplot.

When trying to use --modx option I get another UnboundLocalError

/path/to/planetoplot/modules/ppplot.py", line 971, in make
self.ax = labelmodulo(ax,self.modx)
UnboundLocalError: local variable 'ax' referenced before assignment

I tried to fix it with self.ax = labelmodulo(self.ax,self.modx) but then x axis is blank ...