Dear gstat users, I am trying to specify a mask (Surfer ASCII grid) map for kriging prediction in the command file. I only get:
"gstat: Win32/Cygwin version 2.3.6 (24 September 2001) Copyright (C) 1992, 2001 Edzer J. Pebesma gstat: read failed on file `out.grd` " Would you be so kind as to give me some advice at where to look or what to do to get the file working. As I gather from the "gstat user's manual", I can use the grid file in the Surfer ascii format directly. I have joined a copy of my command file and the first few lines of the grid file at the end of this message. Your help is greatly appreciated. Cordialement, Daniel Chamberland-Tremblay Biologiste Ministère des Ressources naturelles Direction de la conservation des forêts Division des analyses et stratégies 880, chemin Sainte-Foy, 6e étage Québec (Québec) Canada G1S 4X4 Téléphone : (418) 627-8642 poste 4060 Télécopieur : (418) 643-2368 Here a copy of my command file : ########################################### # # # Description du fichier mauricie_population.data : # "sta.romer","xmtm","ymtm","ins.log","ins.log.222","ins.log.cmax","ins.log.cm in","ins.log.pmid","ins.log.detrend3","ins.log.detrend2" # # 1. sta_romer : numéro romer # 2. xmtm : longitude en mètre # 3. ymtm : latitude en mètre # 4. ins.log : population d'insecte log(x+1) # 5. ins.log.222 : ins.log valeurs/222 # 6. ins.log.cmax : ins.log valeurs pondérées par cellule (max mean) # 7. ins.log.cmin : ins.log valeurs pondérées par cellule (min mean) # 8. ins.log.pmid : ins.log valeurs pondérées par polygonale # 9. ins.log.detrend3 : ins.log sans tendance, polynomiale 3 termes # 10. ins.log.detrend2 : ins.log sans tendnace, polynomiale 2 termes set logfile='gstat.log'; data(log.ins) : '/home/daniel/work/master/gstat/mauricie_population.data', x=2, y=3, v=4; set cutoff=70596.01; set width=4706.401; variogram(log.ins) : 0.09 Nug() + 0.07 Sph(39800); # ou 37800 mask: 'out.grd'; predictions(log.ins): 'log_okpr'; ########################################### The start of my Surfer Ascii grid data file DSAA 100 73 206249.8 377999.8 5130584 5255834 782.85101010133 37858.1489899 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 782.85101010133 1171.457070707 1560.0631313118 1948.6691919193 2337.2752525252 2725.8813131315 3114.4873737403 3503.0934343435 3891.6994949497 4280.3055555558 4668.9116161616 5057.5176767677 5446.1237373834 5834.7297979797 6223.3358585859 6611.9419192035 7000.5479797979 7389.154040404 7777.7601010092 8166.3661616165 8554.9722222222 8943.5782828284 9332.1843434331 9720.7904040404 10109.396464647 10498.002525253 10886.608585858 11275.214646465 11663.820707083 12052.426767677 12441.032828283 12829.638888874 13218.244949495 13606.851010101 13995.457070709 14384.063131313 14772.669191919 15161.275252525 15549.881313132 15938.487373737 16327.093434343 16715.699494948 17104.305555556 17492.911616162 17881.517675732 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038 1.70141E+038