[GRASS5] s.proj site header

[Andreas Lange]

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