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
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
