<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 TRANSITIONAL//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=UTF-8">
<META NAME="GENERATOR" CONTENT="GtkHTML/3.6.2">
</HEAD>
<BODY>
Hi list,<BR>
<BR>
I am trying to get a raster map interpolating from a vector file with 14 nodes by means of v.surf.rst. I use <B>GRASS 6.0.0</B> (according to g.version). The steps I made were the nexts:<BR>
<BR>
1. Create a Grass 3D vector file from an ASCII file using v.in.ascii. The ASCII file (named kc.dat) is<BR>
<BR>
502130|4111284|0.75|1<BR>
502243|4111212|0.23|2<BR>
502532|4111109|0.54|3<BR>
502334|4111674|0.67|4<BR>
502639|4111461|0.34|5<BR>
502905|4111461|0.87|6<BR>
503062|4111460|0.42|7<BR>
503325|4111154|0.74|8<BR>
503111|4111107|0.67|9<BR>
503539|4112190|0.38|10<BR>
503828|4112280|0.46|11<BR>
503658|4112414|0.98|12<BR>
505529|4114325|0.67|13<BR>
508391|4116254|0.87|14
<PRE>
where the first two columns are the coords, the third is the z-value and the last one a cat column.
Using the command
<B>> cat kc.dat | v.in.ascii -z z=3 cat=4 out=kctest columns='x double, y double, z double, cat int' </B>
it seems there is no problem and the vector is created. Actually with v.info kctest I get the next info
<B>> v.info kctest</B>
+----------------------------------------------------------------------------+
| Layer: kctest Organization: |
| Mapset: joseantonio Source Date: |
| Location: huenejaUTM100 Name of creator: |
| Database: /home/joseantonio/grassdata |
| Title: |
| Map Scale: 1:1 |
| Map format: native |
|----------------------------------------------------------------------------|
| Type of Map: Vector (level: 2) |
| |
| Number of points: 14 Number of areas: 0 |
| Number of lines: 0 Number of islands: 0 |
| Number of boundaries: 0 Number of faces: 0 |
| Number of centroids: 0 Number of kernels: 0 |
| |
| Map is 3D: 1 |
| Number of dblinks: 0 |
| |
| Projection: UTM (zone 30) |
| N: 4116254.000 S: 4111107.000 |
| E: 508391.000 W: 502130.000 |
| B: 0.230 T: 0.980 |
| |
| Digitize threshold: 0.00000 |
| Comments: |
| |
+----------------------------------------------------------------------------+
2. Set the region by means of
<B>> g.region -p res=100</B>
3. Interpolate with v.surf.rst getting the next error message:
<B>> v.surf.rst input=kctest elev=kcsurface</B>
<B>Authors: original version - H.Mitasova, L.Mitas, I. Kosinovsky, D.P. Gerdes</B>
<B>See manual pages for reference and publications T3</B>
<B>Percent complete: Loading data from attribute table ...</B>
<B>ERROR: Cannot get layer info</B>
Someone knows where the problem is?
Thanks in advance,
José A. Ruiz
Departament of Physics
University of Jaen
</PRE>
</BODY>
</HTML>