[GRASS-user] v.surf.idw G_Malloc
Paul Kelly
paul-grass at stjohnspoint.co.uk
Mon Jan 14 13:11:36 EST 2008
On Tue, 15 Jan 2008, Matthew James wrote:
> Hi Paul - thanks for the fast reply.
>
> Ok - my import command is -
>
> v.in.ascii input=/home/Matthew/XYGrp45.csv output=vectorOut format=point
> fs=, skip=0 'columns=x double precision, y double precision, yield double
> precision' x=2 y=1 z=3 cat=0
>
> output of that command is -
>
> Maximum input row length: 28
> Maximum number of columns: 3
> Minimum number of columns: 3
> Building topology ...
> Registering lines: 1000 2000 3000 4000
> 5000 6000 7000 8000 9000 10000
> 11000 12000 13000 14000 15000 16000
> 17000 18000 19000 20000 21000 22000
> 23000 24000 25000 26000 27000 28000
> 29000 30000 31000 32000 33000 34000
> 35000 36000 37000 38000 39000 40000
> 41000 42000 43000 44000 45000 46000
> 47000 48000 49000 50000 51000 52000
> 53000 54000 55000 56000 57000 58000
> 59000 60000 61000 62000 63000 64000
> 65000 66000 67000 68000 69000 70000
> 71000 72000 73000 74000 75000 76000
> 77000 78000 79000 80000 81000 82000
> 83000 84000 85000 86000 87000 88000
> 89000 90000 91000 92000 93000 94000
> 95000 96000 97000 98000 99000 100000
> 101000 102000 103000 104000 105000 106000
> 107000 108000 109000 110000 111000 112000
> 113000 114000 115000 116000 117000 118000
> 119000 120000 121000 122000 123000 124000
> 125000 126000 127000 128000 129000 130000
> 131000 132000 133000 134000 135000 136000
> 137000 138000 139000 140000 141000 142000
> 143000 144000 145000 146000 147000
> 147909 primitives registered
> Building areas:
>
>
> 0 areas built
> 0 isles built
> Attaching islands:
> Attaching centroids:
>
> Topology was built.
> Number of nodes : 147872
> Number of primitives: 147909
> Number of points : 147909
> Number of lines : 0
> Number of boundaries: 0
> Number of centroids : 0
> Number of areas : 0
> Number of isles : 0
>
> MY SURF.IDW COMMAND -
>
> v.surf.idw input=vectorOut output=rasterOut npoints=12 column=3
I think you need column=yield rather than column=3. 150000 points isn't a
huge amount and should be easy enough to work with. But for much large
datasets it would be easier to work with if you create a 3-D vector by
using the -z flag with v.in.ascii (and leave out the columns option). Then
you could use layer=0 for v.surf.idw to interpolate using the z
co-ordinate.
No idea where the G_malloc error is coming from but hopefully it will go
away with the correct column parameter...
Paul
More information about the grass-user
mailing list