<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi,<br>
<br>
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)<br>
<br>
#Code<br>
<i>r.in.gdal in="./00_REFERENCE/reference_map.tif" out=region_ref
--o</i><i><br>
</i><i>g.region rast=region_ref</i><i><br>
</i><i><br>
</i><i>v.neighbors in=pt out=pt_dens300 size=600 meth=count
--overwrite </i><i><br>
</i><br>
<br>
#extract of "v.info
pt": <br>
| Location: CH_1903+ <br>
| Title: <br>
| Map scale: 1:1 <br>
|--------------------------------------------------<br>
| Map format: native <br>
|--------------------------------------------------<br>
| Type of map: vector (level: 2) <br>
| <br>
| Number of points: 10278 Number of centroids:
0 <br>
| Map is 3D:
No <br>
| Number of dblinks:
1 <br>
|
<br>
| Projection: Swiss. Obl. Mercator <br>
|
<br>
| N: 1309490.94577888 S:
1060967.94534242 <br>
| E: 2906024.04080627 W:
2424801.03992381 <br>
|
<br>
| Digitization threshold: 0 <br>
<br>
<br>
#extract of "r.info region_ref": <br>
| Location:
CH_1903+ <br>
|----------------------------------------------------------------------------<br>
|
<br>
| Type of Map: raster Number of Categories:
0 <br>
| Data Type:
CELL <br>
| Rows:
8912 <br>
| Columns:
14016 <br>
| Total Cells:
124910592 <br>
| Projection: Swiss. Obl.
Mercator <br>
| N: 1297000 S: 1074200 Res:
25 <br>
| E: 2834800 W: 2484400 Res:
25 <br>
| Range of data: min = 1 max =
999 <br>
<br>
</body>
</html>