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//
//
//v.neighbors in=pt out=pt_dens300 size=600 meth=count --overwrite //
/
#extract of "v.info pt":
| Location: CH_1903+
| Title:
| Map scale: 1:1
|--------------------------------------------------
| Map format: native
|--------------------------------------------------
| Type of map: vector (level: 2)
|
| Number of points: 10278 Number of centroids: 0
| Map is 3D: No
| Number of dblinks: 1
|
| Projection: Swiss. Obl. Mercator
|
| N: 1309490.94577888 S: 1060967.94534242
| E: 2906024.04080627 W: 2424801.03992381
|
| Digitization threshold: 0
#extract of "r.info region_ref":
| Location: CH_1903+
|----------------------------------------------------------------------------
|
| Type of Map: raster Number of Categories: 0
| Data Type: CELL
| Rows: 8912
| Columns: 14016
| Total Cells: 124910592
| Projection: Swiss. Obl. Mercator
| N: 1297000 S: 1074200 Res: 25
| E: 2834800 W: 2484400 Res: 25
| Range of data: min = 1 max = 999
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
