Re: ordinary block kriging output problem

["Edzer J. Pebesma" <e.pebesma@geog.uu.nl>]

Stef,

Have you tried to remove the radius=1000 in the script?
It seems that the predictions locations are more than 1000
away from the observation locations.
--
Edzer

stef wrote:

> Hello,
> If possible would you give some suggestions?  I have
> created this script and I think it should work, but
> all I get for predictions in -999..I have included
> everything:
>
> Script:
>
> # gstat command file, Linux version 2.4.3 (04 January
> 2004)
> # Sat Apr 17 13:36:35 2004
> #
> data(VarioAllSites): 'AllSites', x=1, y=2, v=1,
> min=20, max=40, radius=1000;
> variogram(VarioAllSites): 197.675 Nug(0) + 91204.4
> Gau(3.76684e+06);
> data(): 'site_lagoon', x=1, y=2;
> #set cutoff = 500000;
> #set fit = 4;
> #set width = 6500;
>
> blocksize: dx=40, dy=40;
> set output = 'lagoon.out';
>
>
>
> result in .out file:
>
>
> 6235844.5|2106311.5|#1 %-9999 %-9999
> 6287419|2095235.9|#2 %-9999 %-9999
> 6157494|2087979.5|#3 %-9999 %-9999
> 6264220|2085700.8|#4 %-9999 %-9999
> 6256589|2084986.6|#5 %-9999 %-9999
> 6257108.5|2084890.1|#6 %-9999 %-9999
> 6257066|2084825|#7 %-9999 %-9999
> 6254950.5|2083291|#8 %-9999 %-9999
> 6260952|2082856.5|#9 %-9999 %-9999
> 6254825.5|2082516.9|#10 %-9999 %-9999
> 6266499.5|2079800.8|#11 %-9999 %-9999
> 6264909|2079132.5|#12 %-9999 %-9999
> 6310347|2077621.8|#13 %-9999 %-9999
>
>
> input file (AllSites):
>
> 6291417.48572|1788184.62821|#1 %90
> 6299808.45022|1785768.56614|#2 %87
> 6304906.65078|1782448.75422|#3 %84
> 6295229.17171|1779921.87627|#4 %78
> 6295289.61946|1778575.56779|#5 %94
> 6307791.2059|1779666.95835|#6 %96
> 6260517.33576|1856289.34721|#7 %99
> 6311148.23853|1781623.14737|#8 %88
> 6314903.94272|1779944.74406|#9 %78
> 6291488.78722|1790034.45407|#10 %98
> 6306458.86845|1779682.93289|#11 %94
> 6304601.17715|1795455.45236|#12 %98
> 6260792.29089|1853948.96805|#13 %99
> 6289427.30344|1862636.65383|#14 %94
> 6293273.52917|1864338.9109|#15 %90
> 6299150.90601|1865563.54416|#16 %80
> 6305298.07969|1874210.28485|#17 %92
>
> Thanks,
> Stef
>
>
>
> 	
> 		
> __________________________________
> Do you Yahoo!?
> Yahoo! Photos: High-quality 4x6 digital prints for 25¢
> http://photos.yahoo.com/ph/print_splash
>