I try to interpolate scattered points with a breakline. When i run gstat, i always get errors. When i take out the 'edges:...' command, everything works fine, otherwise, i get the errors.
My command file looks like this: data(zuestoll): '/home/marco/gelaendedaten/zuestoll.txt', x=2, y=3, v=4, d=1; edges: '/home/marco/gelaendedaten/zuestolledges.txt'; variogram(zuestoll): 25000 Sph(400); mask: '/home/marco/gelaendedaten/zuestollmask.agr'; blocksize: dx=40, dy=40; predictions(zuestoll): '/home/marco/gelaendedaten/zuestoll_strat_block_40_40.agr'; and the gstat messages look like this: gstat: Linux version 2.4.3 (04 July 2003) Copyright (C) 1992, 2003 Edzer J. Pebesma using Marsaglia's random number generator data(zuestoll): /home/marco/gelaendedaten/zuestoll.txt (Ascii table file) attribute: col[4] [x:] x_2 : [ 739687, 741308] n: 99 [y:] y_3 : [ 223848, 225310] sample mean: 1933.52 sample std.: 159.452 base functions: intercept, x, y reading plot file `/home/marco/gelaendedaten/zuestolledges.txt', 1 polygons read edges: /home/marco/gelaendedaten/zuestolledges.txt (1 edge(s)) [using universal kriging] initializing maps . "chfactor.c", line 153: sizes of objects don't match in function LDLsolve() gstat caught an error that occurred in the matrix library, the reason for it was: unknown, error code 1 gstat: matrix library error: unknown, error code 1 the file with the edge looks like this: 16 739687 224111 739815 224122 739890 224169 739960 224227 740040 224230 740100 224225 740201 224235 740320 224290 740395 224262 740520 224260 740680 224300 740800 224280 740860 224260 740930 224227 741070 224248 741308 224384 Any idea what i am doing wrong? Thanks in advance, Marco
