<html><head></head><body><div class="yahoo-style-wrap" style="font-family:verdana, helvetica, sans-serif;font-size:13px;"><div dir="ltr" data-setdir="false">Hi, I am running Grass 7.8.2 on a HPC (linux) to utilise memory not available to me via my local machine.</div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false">It is running using proj 6.3, gdal 3.0.4, geos 3.8, python 3</div><div dir="ltr" data-setdir="false"><span><br></span></div><div dir="ltr" data-setdir="false"><span>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:</span></div><div dir="ltr" data-setdir="false"><span><br></span></div><div dir="ltr" data-setdir="false"><div>g.region n=5.4715700149536133 e=31.4501590728759766 s=-13.5662746429443359 w=12.1540479660034180 -p</div><div><br></div><div>## Input and clean topology</div><div>v.in.ogr -e -o input=file.json output=filev2 snap=-1 --verbose --overwrite</div><div>v.in.ogr -e -o input=file2.gpkg output=file2 snap=-1 --verbose --overwrite encoding=UTF-8</div><div><br></div><div>## Extract all categories for each dataset so as to maintain polygon (building) contiguity where polygons overlap other polygons within each dataset</div><div>v.extract -d input=filev2 layer=1 type=centroid,area output=buildingsdiss --verbose --overwrite</div><div dir="ltr" data-setdir="false">v.extract -d input=<span><span style="color: rgb(0, 0, 0); font-family: verdana, helvetica, sans-serif;">file</span></span>2 layer=1 type=centroid,area output=file2diss1 where="cat < 6000001" --verbose --overwrite</div><div><br></div><div>## Perform disjoint operation to find buildings that do not intersect buildings in primary dataset</div><div>v.select ainput=file2diss1 alayer=file2diss1 binput=buildingsdiss blayer=buildingsdiss output=DISJOINT1 operator=disjoint --verbose --overwrite</div><div><br></div><div>## Output disjoint building layer and cleaned primary dataset</div><div>v.out.ogr -m input=DISJOINT1 output=DISJOINT1.gpkg --verbose --overwrite</div><div>v.out.ogr -m input=filev2 output=filev2.gpkg --verbose --overwrite</div><div><br></div><div>Having input the topologies and extracted the categories I persistently get the following error using various data when trying to run the disjoint operation:<br></div></div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false"><div><div>projection: 3 (Latitude-Longitude)</div><div>zone: 0</div><div>datum: wgs84</div><div>ellipsoid: wgs84</div><div>north: 5:28:17.649311N</div><div>south: 13:33:58.588333S</div><div>west: 8:09:27.521013W</div><div>east: 51:45:42.670764E</div><div>nsres: 1:00:07.170402</div><div>ewres: 1:00:56.104945</div><div>rows: 19</div><div>cols: 59</div><div>cells: 1121</div><div>NS and EW resolutions are different</div><div>Processing features...</div><div> 0%..........100%</div><div>Processing areas...</div><div> 0%ERROR: Unable to seek: Invalid argument</div><div>NS and EW resolutions are different</div><div><br></div></div><div dir="ltr" data-setdir="false">Any advice appreciated on getting v.select to run. Thanks</div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false">Best wishes, Chris</div><div><br></div></div></div></body></html>