[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