[GRASS-user] Re: Signature made by i.class not accessible?

[Glynn Clements]

Markus Neteler wrote:

> On Thu, Mar 6, 2008 at 12:10 PM,  <Paula.Sankelo at fimr.fi> wrote:
> > The commands I'm using are just:
> >
> > i.class map=varikuva group=testi subgroup=alatesti outsig=signature
> >
> > -Then there's the whole choosing and analyzing the test regions in the
> > graphic window. When that's finished, I say:
> >
> >  i.maxlik group=testi subgroup=alatesti sigfile=signature class=resultmap
> >
> > And I get the error message:
> >
> > ERROR: Can't read signature file [signature].
> >
> > Still, the signature file IS created, and I can see it sitting and lurking
> > there. I just can't use it.
> 
> I have tried to replicate your problem and see that i.group is now
> broken:
> 
> # NC dataset
> i.group mylsat7_2000 sub=mylsat7_2000
> in=lsat7_2000_10,lsat7_2000_20,lsat7_2000_30,lsat7_2000_40,lsat7_2000_50,lsat7_2000_61,lsat7_2000_70,lsat7_2000_80
> Adding raster map <lsat7_2000_10 at landsat> to group
> Adding raster map <lsat7_2000_20 at landsat> to group
> Adding raster map <lsat7_2000_30 at landsat> to group
> Adding raster map <lsat7_2000_40 at landsat> to group
> Adding raster map <lsat7_2000_50 at landsat> to group
> Adding raster map <lsat7_2000_61 at landsat> to group
> Adding raster map <lsat7_2000_70 at landsat> to group
> Adding raster map <lsat7_2000_80 at landsat> to group
> Adding raster map <lsat7_2000_10 at landsat> to subgroup
> Adding raster map <lsat7_2000_20 at landsat> to subgroup
> Adding raster map <lsat7_2000_30 at landsat> to subgroup
> Adding raster map <lsat7_2000_40 at landsat> to subgroup
> Adding raster map <lsat7_2000_50 at landsat> to subgroup
> Adding raster map <lsat7_2000_61 at landsat> to subgroup
> Adding raster map <lsat7_2000_70 at landsat> to subgroup
> Adding raster map <lsat7_2000_80 at landsat> to subgroup
> WARNING: Unable to create file [REF] for subgroup [mylsat7_2000] of group
>          [mylsat7_2000 in landsat]
> i.group complete.
> 
> It only generates
> GRASS 6.4.svn (nc_spm_07):~/grassdata/nc_spm_07/landsat > tree
> group/mylsat7_2000/
> group/mylsat7_2000/
> `-- REF
> 
> 0 directories, 1 file
> 
> which is not sufficient. I suspect that this change broke i.group:
> 
> lib/imagery/fopen.c
> 2008-04-12 03:04  glynn
>         * : Fix handling of group elements to cope with @mapset
> 
> http://trac.osgeo.org/grass/changeset/30948
> 
> Apparently the subdirectory is no longer generated.

Ah; useful information at last.

Yep, it looks like that change broke it (it works if the directory
already exists).

It should hopefully be fixed in r31416 (trunk) and r31418 (6.4
branch); please test.

The original change (r30948) wasn't backported to 6.3, so that isn't
affected.

-- 
Glynn Clements <glynn at gclements.plus.com>