[GRASS-user] v.select disjoint operation. Grass 7.8.2

Christopher Lloyd chrislloyd2 at yahoo.co.uk
Wed May 12 04:14:29 PDT 2021


Hi, I am running Grass 7.8.2 on a HPC (linux) to utilise memory not available to me via my local machine.
It is running using proj 6.3, gdal 3.0.4, geos 3.8, python 3
I am trying to perform the v.select disjoint operation, which I have used successfully before on a similar setup and similar data using the same version of Grass (I suspect only the supporting modules and input data have changed). My code is:
g.region n=5.4715700149536133 e=31.4501590728759766 s=-13.5662746429443359 w=12.1540479660034180 -p
## Input and clean topologyv.in.ogr -e -o input=file.json output=filev2 snap=-1 --verbose --overwritev.in.ogr -e -o input=file2.gpkg output=file2 snap=-1 --verbose --overwrite encoding=UTF-8
## Extract all categories for each dataset so as to maintain polygon (building) contiguity where polygons overlap other polygons within each datasetv.extract -d input=filev2 layer=1 type=centroid,area output=buildingsdiss --verbose --overwritev.extract -d input=file2 layer=1 type=centroid,area output=file2diss1 where="cat < 6000001" --verbose --overwrite
## Perform disjoint operation to find buildings that do not intersect buildings in primary datasetv.select ainput=file2diss1 alayer=file2diss1 binput=buildingsdiss blayer=buildingsdiss output=DISJOINT1 operator=disjoint --verbose --overwrite
## Output disjoint building layer and cleaned primary datasetv.out.ogr -m input=DISJOINT1 output=DISJOINT1.gpkg --verbose --overwritev.out.ogr -m input=filev2 output=filev2.gpkg --verbose --overwrite
Having input the topologies and extracted the categories I persistently get the following error using various data when trying to run the disjoint operation:

projection: 3 (Latitude-Longitude)zone:       0datum:      wgs84ellipsoid:  wgs84north:      5:28:17.649311Nsouth:      13:33:58.588333Swest:       8:09:27.521013Weast:       51:45:42.670764Ensres:      1:00:07.170402ewres:      1:00:56.104945rows:       19cols:       59cells:      1121NS and EW resolutions are differentProcessing features...   0%..........100%Processing areas...   0%ERROR: Unable to seek: Invalid argumentNS and EW resolutions are different
Any advice appreciated on getting v.select to run. Thanks
Best wishes, Chris
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20210512/e9eecb17/attachment.html>


More information about the grass-user mailing list