[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Newbie tries to read in data file



>>>>> "Kreft" == Jan Kreft <Kreft@Cardiff.ac.uk> writes:

    Kreft> I'm probably missing something simple, but being a biologist I have no
    Kreft> experience with what the data formats supported by gstat actually look
    Kreft> like. 

See gstat manual, section 2.2 "Data formats'. And look at examples in
the directory cmd of gstat distribution.

    Kreft> I have matrices that simply contain, row by row, and column by
    Kreft> column, the data values. A small example is this:
<snip>     
    Kreft> As you might notice, I'm trying to use the ArcInfo asciigrid format, which
    Kreft> is why I have added the header lines to this file.

    Kreft> How can I read this into gstat?

Ypu really have raster data, don't you? Otherwise it makes little
sense to format it to any raster format.

You can use your data in data-statement or in map statement. The
minimal example is:

data(test): 't.dat';

Save it to t.cmd and use with 

gstat t.cmd

t.dat would be the name of file with your data. The above example
starts gstat in interactive mode, giving a possibility to calculate
and fit variogram/variogram model.

Hope it helps.

-- 
Konstantin Malakhanov, research fellow
---
Institut for Hydraulic Engineering & Groundwater Mgmt., RWTH Aachen, Germany
E-Mail: malakhanov@iww.rwth-aachen.de          Tel. +49/(0)241/807343