[GRASS-dev] i.segment: Invalid region id -1

[Nikos Alexandris]

NA:
> > > --%<--
> > > i.segment msx_hpf out=segments_msx_hpf_seeded_t0.02 threshold=0.02
> > > minsize=4 seed=segments_pan_t0.01 memory=3000 iterations=1000
> > > -->%--
> > > It worked in one case (repeated to be sure) and it failed in another!

MM:
> > Different computational regions?

NA:
> Absolutely.  Two different, independent trials. Same Location, different
> Mapsets, different regions ("physically", if I may say so, and
> computationally).

> > > In the failing case, before the ERROR message, there are multiple
> > > WARNINGS issued:
> > > WARNING: Region consists of only one cell, nothing to update

> > This should not happen, a bug in the region growing algorithm.

Is a Warning message always a bug (in the case of i.segment)?

> > > I went after looking all of the details of the involved maps.  The only
> > > "strange" thing I can see (which I caused) is that the region is 0.6,
> > > the seed (segments_pan_t0.01) is also 0.6 while the group of Pan-
> > > Sharpened images are (each) of 0.60017817 (ns) x 0.60016801 (we)
> > > resolution.  Is this my mistake? The resolution(s) should be identical,
> > > right?

> > Yes. I guess in the process of pansharpening, the region was set to
> > 0.6, then the resolution was adjusted to the extents (for g.region,
> > extents have precedence over resolution).

-- useful --
> > The correct way of adjusting the region would be to either set the region
> > to the pan band or align the region to the band band (g.region align=pan).
-- useful --

> > Note that g.region res=0.6 -a can introduce a pixel shift.
> > Can you reproduce that with sample data?

> I'll try.

I can't reproduce this anymore (with my working images) after being careful 
with extent and resolution. Two examples, both worked fine:

1) Seed map (Pan), HPF-Sharpened MSes and Region: all of 0.59998996
2) Seed map (Pan), group of HPF-Sharpened MSes and Region: all of 0.60017817 x 
0.60016801

However, in a completely different setup (Location, Mapsets), I can see 
multiple Warnings, like above, no error though. The process completes fine and 
I get to see the segments.  Is this ok?

If not, could it be related with sub-sequent processing errors? For example, 
like invalid classification signatures from i.cluster and failure of i.maxlik?  
I use(d) fancy segments-stats (segments used as a cover, group of Multi-
Spectral images used as source for stats) like skewness, kurtosis, etc.

Nikos