[GRASS-user] r.resamp.rst - ERROR: Input window changed while maps are open for read.

GRASS GRASS user
1

Dear all,

I'm trying to use r.resamp.rst in GRASS7, I would like to interpolate
a map from a resolution of 1m to a resolution of 0.25.
I'm working with the following region:

{{{
GRASS 7.0.svn (edi):~ > g.region -p
projection: 99 (Transverse Mercator)
zone: 0
datum: osgb36
ellipsoid: airy
north: 674050
south: 673911
west: 325025
east: 325233
nsres: 1
ewres: 1
rows: 139
cols: 208
cells: 28912
}}}

Now if I try to run the following command I got:

{{{
GRASS 7.0.svn (edi):~ > r.resamp.rst input=diff ew_res=0.25 nsres=0.25
elev=diff_overl3_ten40 overlap=3 tension=40.
ERROR: Input window changed while maps are open for read. Map name <diff>
}}}

If I change the region to the resolution 0.25 I still got the same error:

{{{
GRASS 7.0.svn (edi):~ > g.region res=0.25 -p
projection: 99 (Transverse Mercator)
zone: 0
datum: osgb36
ellipsoid: airy
north: 674050
south: 673911
west: 325025
east: 325233
nsres: 0.25
ewres: 0.25
rows: 556
cols: 832
cells: 462592
GRASS 7.0.svn (edi):~ > r.resamp.rst input=diff ew_res=0.25 nsres=0.25
elev=diff_overl3_ten40 overlap=3 tension=40.
ERROR: Input map resolution differs from current region resolution!
}}}

the map diff is a FCELL map

{{{
GRASS 7.0.svn (edi):~ > r.info diff -g
north=680000
south=650000
east=340000
west=310000
nsres=1
ewres=1
rows=30000
cols=30000
cells=900000000
datatype=FCELL
}}}

What did I do wrong?
Any hints?

Pietro

2

Pietro wrote:

I'm trying to use r.resamp.rst in GRASS7,

...

Now if I try to run the following command I got:

...

ERROR: Input window changed while maps are open for read.
Map name <diff>

...

What did I do wrong?
Any hints?

you didn't do anything wrong, the module has not been fully ported to grass 7 yet. try it in GRASS 6.

Hamish

3

On Mon, May 6, 2013 at 12:49 PM, Pietro <peter.zamb@gmail.com> wrote:

Dear all,

I'm trying to use r.resamp.rst in GRASS7, I would like to interpolate
a map from a resolution of 1m to a resolution of 0.25.

Note also the module mentioned below.

Re r.resamp.rst usage:

On Thu, May 22, 2008 at 4:02 PM, Paul Kelly wrote:

On Thu, 22 May 2008, Hamish wrote:

Thanks for the answer - I am in the process of trying it out. But:
When I try to use r.resamp.rst, I get the error message "Input map
resolution differs from current region resolution!". I thought that
r.resamp.rst is doing exactly that - changing the resolution?

...

As far as I know it works "the wrong way round", i.e. different from the
other resampling modules - you actually run it with the resolution set to
that of the original map (g.region rast=mapname), and set the desired new
resolution using command-line options.

I'm afraid that this is still the case.
However, like this it should work for you.

Note also this module:
http://grass.osgeo.org/grass70/manuals/r.resamp.interp.html
r.resamp.interp - Resamples raster map layers to a finer grid using
interpolation.

@devs: the best might be to add the RST method to r.resamp.interp...

Markus

4

Dear all,

the issue is still there. I get the same error message. Using the stable
version of GRASS (6.4) is not an option for me since I am working with very
large files (>2GB) for which this module returns the error described here
<http://trac.osgeo.org/grass/ticket/1786&gt; .

Any suggestions?

Many thanks,
Bert

--
View this message in context: http://osgeo-org.1560.x6.nabble.com/r-resamp-rst-ERROR-Input-window-changed-while-maps-are-open-for-read-tp5051545p5129311.html
Sent from the Grass - Users mailing list archive at Nabble.com.

5

Bert,

On Sun, Mar 16, 2014 at 6:39 PM, bwus <b3rt.wu@gmail.com> wrote:

Dear all,

the issue is still there. I get the same error message. Using the stable
version of GRASS (6.4) is not an option for me since I am working with very
large files (>2GB) for which this module returns the error described here
<http://trac.osgeo.org/grass/ticket/1786&gt; .

The ticket reports on r.resamp.rst in GRASS 6. We cannot fix it there
easily since many libraries are involved (see ticket).

Does the problem also exist in GRASS 7 for you?
You can install it in parallel to GRASS 6, see
http://grass.osgeo.org/download/software/linux/#grass70

Markus

6

Markus Neteler wrote:

> the issue is still there. I get the same error message. Using the stable
> version of GRASS (6.4) is not an option for me since I am working with very
> large files (>2GB) for which this module returns the error described here
> <http://trac.osgeo.org/grass/ticket/1786&gt; .

The ticket reports on r.resamp.rst in GRASS 6. We cannot fix it there
easily since many libraries are involved (see ticket).

Does the problem also exist in GRASS 7 for you?
You can install it in parallel to GRASS 6, see
http://grass.osgeo.org/download/software/linux/#grass70

In 6.x, the tempfile code in the RST library is limited to 2 GiB
(ticket #1786).

In 7.x, the RST library isn't working because it hasn't been updated
for the new raster window handling (r40772 etc) (ticket #1718).

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