Geostatistical R Package - Board

Hello,

Je crée une db à partir d'un raster. J'utilise db.sel pour indiquer quelles sont les mailles actives puis j'affiche la carte de la db. Celle-ci se trouve inversée Nord-Sud par rapport au raster d'origine. Ai-je fait une erreur quelque part ?

Question subsidiaire : comment procéder pour obtenir un ratio égal entre les 2 axes (comme avec asp=1) ?

Voici un extrait de mon code :

v1=getValues(grid.in)
cond=(!is.na(v1))
target.db = db.create(flag.grid = TRUE, nx=c(grid.in@ncols,grid.in@nrows),
x0=c(grid.in@extent@xmin,grid.in@extent@ymin), dx=res(grid.in))
target.db = db.locate(target.db,2:3,"x")
target.db = db.sel(target.db,cond,flag.replace=TRUE)

et les images du raster et de la db :

Salut Dimitri

Peux-tu m'envoyer (par mail) toin fichier raster initial STP.
La lecture est tres dependente du contenu du fichier.
Je suppose que la procedure incriminee est db.read.format.

Par ailleurs j'ai corrige db.write.format qui ne prenait pas en compte proprement les selections (la correction sera disponible dans la version 11.1.2).

Didier

Je viens de t'envoyer le raster par mail.

Merci pour ton aide,

Dimitri