[GRASS-dev] Test if raster map and current region match?

[Markus Neteler]

(carried here from grass-stats, discussing improvements of the
 R-GRASS interface)

On Tue, Jun 24, 2008 at 1:48 PM, Roger Bivand <Roger.Bivand at nhh.no> wrote:
> On Tue, 24 Jun 2008, Markus Neteler wrote:
>> On Tue, Jun 24, 2008 at 12:23 PM, Roger Bivand <Roger.Bivand at nhh.no> wrote:
>>> Is there a shell-level utility in GRASS to check whether the region of a
>>> raster is equal to the current working region - the function could call
>>> that and create a copy of the required raster in current resolution if they
>>> differ?
>>
>> The easiest way it to compare the output of
>>
>> r.info -g map127
>> north=230000
>> south=214000
>> east=646000
>> west=628000
>>
>> and
>> g.region -g | grep '^n=\|^s=\|^e=\|^w='
>> n=228513
>> s=214975.5
>> w=629992.5
>> e=645012
>
> But the resolution may differ too - I need a one shot comparison. Isn't
> anything exposed - I've looked in g.copy and r.resample without seeing
> anything. Maybe:
>
> r.info -gs elevation.dem
> g.region -g | grep '^n=\|^s=\|^e=\|^w=\|^ewres=\|^nsres='
>
> but an internal match (where r.resample would be a no-op or not) would be
> helpful.

Any hints for Roger?

Markus