[GRASS-dev] [GRASS GIS] #1964: m.proj on wingrass: protect g.proj's path name to the grid file

[GRASS GIS]

#1964: m.proj on wingrass: protect g.proj's path name to the grid file
---------------------------------------------+------------------------------
 Reporter:  hamish                           |       Owner:  grass-dev@…              
     Type:  defect                           |      Status:  new                      
 Priority:  major                            |   Milestone:  6.4.4                    
Component:  Shell Scripts                    |     Version:  svn-releasebranch64      
 Keywords:  wingrass, m.proj, cs2cs, g.proj  |    Platform:  MSWindows XP             
      Cpu:  x86-32                           |  
---------------------------------------------+------------------------------
 Hi,

 this has been a known thing for a while, but I think it deserves its own
 ticket.  e.g. see #820's comment 28.

 The m.proj shell script uses cs2cs to do the coord sys transform on the
 list of points. Unfortunately due to the way it parses its command line
 args, it's impossible to get cs2cs to work with quoting, so a +nadgrids=
 with spaces in the /path/to/grid/file side gets broken.

 To add to this, msys helpfully decides to translate the string-which-
 looks-like-a-path and breaks it even further.

 e.g. on wingrass 6.4.3svn make a new location, epsg 27200, for the datum
 transform option choose the grid file. (if the loc'n wizard is broken you
 may need to use g.setproj from the text startup with proj=nzmg
 datum=nzgd49)

 that gets you this PROJ_INFO file:
 {{{
 name: New Zealand Map Grid
 proj: nzmg
 lat_0: -41
 lon_0: 173
 x_0: 2510000
 y_0: 6023150
 no_defs: defined
 datum: nzgd49
 ellps: international
 nadgrids: nzgd2kgrid0005.gsb
 }}}

 and g.proj translates to these proj4 terms:
 {{{
 GRASS 6.4> g.proj -j
 +proj=nzmg
 +lat_0=-41
 +lon_0=173
 +x_0=2510000
 +y_0=6023150
 +no_defs
 +a=6378388
 +rf=297
 +nadgrids=c:/Program Files/GRASS GIS 6.4.3svn/etc/nad/nzgd2kgrid0005.gsb
 +to_meter=1
 }}}

 m.proj does this to that:
 {{{
 if [ $GIS_FLAG_I -eq 1 ] ; then
    OUT_PROJ=`g.proj -j | (
        OUT_PROJ=
        while read line ; do
            OUT_PROJ="$OUT_PROJ '$line'"
        done
        echo "$OUT_PROJ"
    )`
 fi
 }}}


 msys sees a "/" and decides to act, e.g.:
 {{{
 GRASS 6.4> g.message "/"
 C:/PROGRA~1/GRASSG~1.3SV/msys
 }}}

 and the broken result of all this is:
 {{{
 GRASS 6.4> echo "167 -45" | m.proj -i
 input parameters=[+proj=longlat +datum=WGS84]
 output parameters=[+proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000
 +y_0=6023150 +no_defs +a=6378388 +rf=297
 +nadgrids=c;C:\PROGRA~1\GRASSG~1.3SV\msys\Program Files\GRASS GIS
 6.4.3svn\etc\nad\nzgd2kgrid0005.gsb +to_meter=1]
 Rel. 4.8.0, 6 March 2012
 <cs2cs.exe>:
 Sys errno: 2: No such file or directory
 Files/GRASS
 Rel. 4.8.0, 6 March 2012
 <cs2cs.exe>:
 Sys errno: 2: No such file or directory
 GIS
 Rel. 4.8.0, 6 March 2012
 <cs2cs.exe>:
 Sys errno: 2: No such file or directory
 6.4.3svn/etc/nad/nzgd2kgrid0005.gsb
 }}}


 perhaps the cygpath util could help?
    http://cygwin-lite.sourceforge.net/html/cygpath.html

 or a couple of .bat file one-liners to make short 8.3 space-free path
 names on Windows:
 {{{
 for /d %I in (*) do @echo %~sI
   or
 for %I in (.) do echo %~sI
 }}}


 Hamish

-- 
Ticket URL: <https://trac.osgeo.org/grass/ticket/1964>
GRASS GIS <http://grass.osgeo.org>