[GRASS-dev] 6.4.0 blocker bugs

Hamish hamish_b at yahoo.com
Thu Jul 29 09:20:59 EDT 2010


> Hamish wrote:
> > next problem: r.in.wms calls wms.request which calls r.tileset.
> > r.tileset calls cs2cs*, with +proj4 terms taken from SRS=`g.proj -j`.
> > 
> > for spearfish those terms include the conus grid file, which is in
> > $GISBASE/etc/nad/ ... and $GISBASE of course contains a space, and
> > then cs2cs dies a horrible death.
> > 
> > so how to quote the right side of the individual +nadgrids= line? and
> > will cs2cs accept that?

Glynn wrote:
> You need to construct the cs2cs command line "the hard way"; simply
> inserting the "g.proj -j" output won't work.
> 
> Something like:
> 
>     opts=`g.proj -j | (
>         opts=
>         while read line ; do
>             opts="$opts '$line'"
>         done
>         echo "$opts"
>     )`

thanks, applied in 6.5svn in r.in.wms/wms.request and r.tileset.

... still no-go. at the 6.4.0svn msys prompt:

GRASS 6.4> g.proj -j
+proj=utm
+zone=13
+a=6378206.4
+rf=294.9786982
+no_defs
+nadgrids=c:/Program
Files/GRASS-64-SVN/etc/nad/conus
+to_meter=1.0


.. so g.proj is broken ..

g.proj/output.c  print_proj4():

    for (i = proj4mod; *i; i++) {
        /* Don't print the first space */
        if (i == proj4mod && *i == ' ')
            continue;

! -->   if (*i == ' ' && !(dontprettify))
            fputc('\n', stdout);
        else
            fputc(*i, stdout);
    }
    fputc('\n', stdout);





> BTW: the Python versions of m.proj and r.in.aster will have problems
> with this, as they use "g.proj -jf ...".

ok



Hamish


      


More information about the grass-dev mailing list