manipulation de fichier netcdf avec {raster}

[Sophie Dubois]

Bonjour à tous,

Je travaille sur un fichier au format NC classique qui contient une variable exprimée en fonction de la latitude, de la longitude et du temps (variable mesurée tous les jours pendant plusieurs années).
J'importe le fichier .nc avec la fonction raster:

data <- raster("nom_fichier.nc")

Ca c'est ok.
Ensuite je "lis" mon fichier avec la fonction readAll:
class: rasterLayer
band: 1 (of 3546 bands)
dimensions: 50, 130, 6500 (nrow, ncol, ncell)

Mon problème est que je veux pouvoir réaliser une opération sur les valeurs comprises dans les différentes "bands" de mon fichier (la moyenne) et représenter cette moyenne et non uniquement la première "band". Ma question est donc: comment accéder à ces différentes "bands" pour pouvoir faire des opérations dessus ou ne sélectionner que certaines d'entre elles, si par exemple je ne veux m'intéresser qu'à une certaine période de mon étude?

J'en profite pour vous demander si vous connaissez de bons tutoriels de cartographie de ce type de données avec le package raster ou avec les packages ncdf/ncdf4 ?

Merci d'avance de votre aide!

Sophie

[dicko ahmadou]

Salut,

Utilises stack à la place de raster.

Code : Tout sélectionner

data <- stack("file.nc")

[Sophie Dubois]

Génial! Merci beaucoup!!

Une autre petite question: je suppose qu'il est possible de mettre un fond de carte vierge (juste un continent par exemple) sous une "projection" raster mais je n'y arrive pas avec le package raster...
J'ai essayé d'utiliser les cartes de fond que je créé avec maps mais elles sont "écrasées" par ma "projection" de ma variable d'intérêt...
Est-ce qu'avec un shapefile, ça marche dans le package raster? Si non, y-a-t-il une petite fonction ou un argument style add=TRUE qui fonctionnerait pour lui dire de ne pas faire un nouveau graphique mais de mettre les 2 "couches" sur le même?

Merci par avance!

Sophie