[GRASS-dev] i.rectify still chokes on @mapset

Markus Neteler neteler at osgeo.org
Tue Apr 14 17:46:20 EDT 2009


On Tue, Apr 14, 2009 at 8:56 PM, Markus Neteler <neteler at osgeo.org> wrote:
> On Tue, Apr 14, 2009 at 6:39 PM, Michael Barton <michael.barton at asu.edu> wrote:
>> One of my students is testing the GRASSRC4 windows release. Apparently,
>> i.rectify still raises and error if the maps in the group file have the
>> mapset specification (map at mapset). I thought this was fixed. We can try to
>> strip the mapset specifier out in the GUI, but I thought that all modules
>> were supposed to respect it.
>> Michael
>
> The related discussion is in this ticket:
> http://trac.osgeo.org/grass/ticket/70

I have downloaded the "imagery" package for Spearfish,
- run i.target to point the "demo" group to the "spearfish60"
  location,
- run i.points to set 4 GCPs
- run i.rectify:

GRASS 6.5.svn (imagery60):~ > i.rectify demo in=nhap.1 ext=rect order=1
Region N=4928327.747796 S=4922940.933568 E=602038.752218 W=588501.065840
Resolution EW=6.474264 NS=6.293007
Rectified input raster map <nhap.1> will be saved as <nhap.1rect>
 100%
Rectify <nhap.1 at PERMANENT> (location <imagery60>)
into <nhap.1rect at neteler> (location <spearfish60>) ... complete
-----------------------------------------------
i.rectify complete.

and it works... (result is ok in the target UTM location).

head -20 grassdata/imagery60/user1/group/demo/*
==> grassdata/imagery60/user1/group/demo/POINTS <==
#                   image                          target           status
#            east           north            east           north   (1=ok)
#
      2042.513344     -712.821839   600117.032967  4925562.236842    1
      1097.131868     -765.296935   594306.527473  4925394.868421    1
       650.701727     -499.641762   591547.582418  4927194.078947    1
       448.056089     -557.957883   590370.589286  4926762.321514    1
       599.690518    -1082.381863   591988.872607  4923848.733105    1

==> grassdata/imagery60/user1/group/demo/REF <==
nhap.1 PERMANENT b
nhap.2 PERMANENT g
nhap.3 PERMANENT r

==> grassdata/imagery60/user1/group/demo/TARGET <==
spearfish60
neteler

looks all good to me.

And then, just to be sure, created a new group:
GRASS 6.5.svn (imagery60):~ > i.group group=mydemo input=nhap.1,nhap.2,nhap.3
Adding raster map <nhap.1 at PERMANENT> to group
Adding raster map <nhap.2 at PERMANENT> to group
Adding raster map <nhap.3 at PERMANENT> to group
i.group complete.
GRASS 6.5.svn (imagery60):~ > cat grassdata/imagery60/user1/group/mydemo/REF
nhap.1 PERMANENT
nhap.2 PERMANENT
nhap.3 PERMANENT

also ok.

So might it be a Windows specific issue?

Markus


More information about the grass-dev mailing list