[GRASS-user] v.neighors error:region do not overlap

patrick s. patrick_gis at gmx.net
Sun Aug 2 08:03:40 PDT 2015


Hi Markus

thanks for your answer. epsg:2056 is correct for CH_1903+;
And here the print you ask for, comming out of mapset in location 1903+.

Patrick

#g.region raster=region_ref -p
projection: 99 (Swiss. Obl. Mercator)
zone:       0
datum:      ch1903
ellipsoid:  bessel
north:      1297000
south:      1074200
west:       2484400
east:       2834800
nsres:      25
ewres:      25
rows:       8912
cols:       14016
cells:      124910592


On 01.08.2015 23:56, Markus Neteler wrote:
> On Wed, Jul 29, 2015 at 4:40 PM, patrick s. <patrick_gis at gmx.net> wrote:
>> Hi,
>>
>> I need to run v.neighbors on multiple datasets for the identic region
>> (extend and resolution).  One dataset has observations outside of this
>> region, and causes an error: "ERROR: Input vector and computational region
>> do not overlap." However, viewing the extend in v.info and r.info (see
>> below) and overlaying the layers in the GUI, I see overlaps for 90% of my
>> points. Is this expected behavior or do I miss something? (GRASS70 from
>> Repos, OS:Ubuntu 15.04)
>>
>> #Code
>> r.in.gdal in="./00_REFERENCE/reference_map.tif" out=region_ref --o
>> g.region rast=region_ref
> Just to be sure, please post
>
> g.region raster=region_ref -p
>
> I generated a location with epsg:2056 (is that the right one for
> CH_1903+ ?) and there is indeed map overlap.
> But just to be sure that the computational region is right, too.
>
> Markus

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20150802/fcc261d4/attachment.html>


More information about the grass-user mailing list