[GRASS5] s.proj site header
Andreas Lange
Andreas.Lange at Rhein-Main.de
Fri Sep 28 05:29:44 EDT 2001
Hi,
Helena wrote:
>
> This appears to be a small but annoing bug - s.proj puts a comment at the
> beginning of a projected
> file, however other programs such as s.info do not like it
>
> s.info lidar.utm
> ERROR: G_site_new_struct: invalid # dims or fields
>
> # s.proj input=lidar99 location=coastll mapset=helena dbase=/data/grassdata
> output=lidar99.utm
> 443873.48614701|3978547.65291465|#1 %5.054
>
> So you have to remove it manually.
>
> To keep the useful info on where the data came from I believe that instead of
> #s.proj.......
> this info should be passed as
> description|s.proj.......
>
> It is in main.c here:
> if(outfile){ /* also write G_recreate_command as comment */
> Site_head sh;
> char buf[80];
>
> fprintf(outfile,"# %s\n", G_recreate_command());
>
> G_site_get_head (infile, &sh);
> G_site_put_head (outfile, &sh);
> }
> Can somebody please fix it, in case that it has not been done already.
I think the correct way should be to read out the sites header data with
G_site_get_head and then add the G_recreate_command output to the desc
part of the header structure and write this back with G_site_put_head to
the new output file. I'll look next week.
cu,
Andreas
--
Andreas Lange, 65187 Wiesbaden, Germany, Tel. +49 611 807850
url: http://mitglied.tripod.de/AndreasLange
mail: Andreas.Lange_at_Rhein-Main.de - A.C.Lange_at_GMX.net
More information about the grass-dev
mailing list