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. <[email protected]> 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
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
