Ich habe einen Xarray-Datensatz namens „dens“ und möchte ihn grafisch darstellen.

Das ist der Datensatz:

<xarray.dataset>
dimensions:  (time: 641, lat: 30, lon: 30)
coordinates:
  * time     (time) datetime64[ns] 2013-07-01t12:00:00 ... 2013-08-02t12:00:00
  * lon      (lon) float64 32.73 32.83 32.94 33.05 ... 35.53 35.64 35.75 35.85
  * lat      (lat) float64 31.08 31.27 31.47 31.66 ... 36.06 36.25 36.44 36.63
data variables:
    density  (time, lat, lon) float64 2e+03 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

Ich verwende den Befehl

plt.contourf(dens.density.values[-1,:,:]);

hat es eingezeichnet und es hat funktioniert, aber da ich wollte, dass auch die Küstenlinie auf dem Grundstück eingezeichnet wird, habe ich es auch mit

m = basemap(llcrnrlon=data['lon'].min(), llcrnrlat=data['lat'].min(),
                urcrnrlon=data['lon'].max(), urcrnrlat=data['lat'].max(), resolution='i', suppress_ticks=1)
m.drawcoastlines();
m.fillcontinents(color='gray',lake_color='gray');

Aber wenn ich alle Befehle ausführe und dann plt.show() verschwindet das Konturdiagramm und alles, was es mir zeigt, ist die Küstenlinie.

Wie kann dieses Problem gelöst werden, um eine Höhenlinienkarte und eine Küstenlinienkarte in derselben Abbildung zu erhalten?

Tut mir leid, wenn das eine dumme Frage ist, aber ich bin sehr neu in Python

Vielen Dank für Ihre Hilfe,

Yutam

Bearbeiten: Mir ist gerade erst klar geworden, dass ich versuche, zwei verschiedene „Toolkits“ zu kombinieren, und dass es möglich ist, all dies nur mit dem Grundkarten-Toolkit zu tun, aber einfach zu schreiben

m.contourf(dens.density.values[-1,:,:]);

gibt mir diesen Fehler:

---------------------------------------------------------------------------
typeerror                                 traceback (most recent call last)
cell in[21], line 1
----> 1 m.contourf(dens.density.values[-1,:,:])

typeerror: basemap.contourf() missing 2 required positional arguments: 'y' and 'data'

Noch eine Änderung: Ich entdecke immer mehr und nachdem ich die Dokumentation für basemap gelesen hatte, wurde mir klar, dass die Syntax des Befehls so aussehen sollte

m.contourf(dens.lon.values,dens.lat.values,dens.密度.values[-1,:,:]);

Aber jetzt bekomme ich diesen Fehler:

IndexError: too many indices for array: array is 1-dimensional, but 2 were indexed

Ich schätze, das liegt daran, dass mein Dichtearray zweidimensional ist, aber wie extrahiere ich den Dichtewert daraus? Ich verwende [-1] in der Zeitdimension, weil ich eigentlich nur den letzten Zeitschritt benötige

Nochmals vielen Dank im Voraus, Yotam

Letzte Bearbeitung

Dies ist die letzte Handlung. Wie kann ich das umliegende Land grau statt lila machen? Gibt es außerdem eine Möglichkeit, ein größeres, etwas größeres geografisches Gebiet zu beschreiben, ohne die Daten zu verfälschen?

Hier sind die neuen Zahlen für meine tatsächlichen Daten

Richtige Antwort

Das Plotten mit Grundkarte und Xarray wurde hier besprochen.

m = basemap(llcrnrlon=data['lon'].min(), llcrnrlat=data['lat'].min(),
            urcrnrlon=data['lon'].max(), urcrnrlat=data['lat'].max(), 
  resolution='i', suppress_ticks=1)
m.drawcoastlines();
m.fillcontinents(color='gray',lake_color='gray')
dens.density[-1,:,:].plot.contourf()
plt.show()

import xarray as xr
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cf

ds = xr.open_dataset('filename.nc')
fig = plt.figure(figsize=(8,8))
crs=ccrs.platecarree()
ax = fig.add_subplot(1,1,1, projection=crs)
gl = ax.gridlines(crs=crs, draw_labels=true,
                linewidth=0.01, color='gray', alpha=0.5, linestyle='-.')

ax.add_feature(cf.coastline.with_scale("50m"), lw=0.5)
ax.add_feature(cf.borders.with_scale("50m"), lw=0.3)

ds.density[-1,:,:].plot.contourf()
plt.show()

from matplotlib.colors import LinearSegmentedColormap
cm = LinearSegmentedColormap.from_list('', ['white', *plt.cm.Blues(np.arange(255))])
ds.density[-1,:,:].plot.contourf(cmap=cm)