Bonjour,

Voici mon problème, j'ai un raster fait à partir d'une matrice de dimension: nrow=4950, ncol=5890

Quand je fais plot(a) ou levelplot(a), alors ça fonctionne parfaitement.

Mais maintenant, j'ai besoin de la librairie sp (pour superposer plusieurs couches, ouvrir des shapes avec mon raster, etc.)

Mais quand je fais par exemple:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
 
> spplot(a)
Alors je me retrouve avec une carte complètement pixélisée, ce n’est pas net....

Je suis alors allé voir dans la doc et j'ai trouvé l'argument maxpixels donc la logique voudrait que je fasse ceci:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
> spplot(q, col.regions=myPal(100), maxpixels=4950*5890)
Je mets les dimensions de ma matrice pour maxpixels, mais ça charge pendant un certain temps puis ça me fait une erreur d'allocation de la mémoire (variable trop grosse). Si je mets un million de pixels, ça fonctionne, mais je suppose que tous mes pixels ne s'affichent pas, et cela ne convient pas au travail que je fais.

Ce qui est bizarre c'est que le chargement est long pour spplot, mais pas pour plot ou levelplot.

Quelqu'un a une idée ?

Merci