<html><head></head><body><div class="ydp3995c95cyahoo-style-wrap" style="font-family: verdana, helvetica, sans-serif; font-size: 13px;"><div></div>
<div dir="ltr" data-setdir="false">Hi Markus, Thanks for your response. I note that another user recently had a similar error with a similar (but not the same) module.</div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false">I am currently splitting my input data into smaller sized files, each for input separately to the v.select module, to see if this resolves the issue. I will let you know on progress. If the failure continues then I will supply you with a reproducible example using the NC data as you indicate. Many thanks.</div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false">Best wishes, Chris</div>
</div><div id="yahoo_quoted_1695427563" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
On Tuesday, 18 May 2021, 22:01:58 BST, Markus Neteler <neteler@osgeo.org> wrote:
</div>
<div><br></div>
<div><br></div>
<div><div dir="ltr">Hi,<br clear="none"><br clear="none">Is there a chance to receive the dataset for testing (or, ideally, a<br clear="none">reproducible example with the North Carolina sample dataset from<br clear="none"><a shape="rect" href="https://grass.osgeo.org/download/data/" target="_blank">https://grass.osgeo.org/download/data/</a>)?<br clear="none"><br clear="none">thanks,<br clear="none">Markus<br clear="none"><div class="yqt4699886086" id="yqtfd41237"><br clear="none">On Wed, May 12, 2021 at 1:14 PM Christopher Lloyd via grass-user<br clear="none"><<a shape="rect" ymailto="mailto:grass-user@lists.osgeo.org" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>> wrote:<br clear="none">><br clear="none">> Hi, I am running Grass 7.8.2 on a HPC (linux) to utilise memory not available to me via my local machine.<br clear="none">><br clear="none">> It is running using proj 6.3, gdal 3.0.4, geos 3.8, python 3<br clear="none">><br clear="none">> 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:<br clear="none">><br clear="none">> g.region n=5.4715700149536133 e=31.4501590728759766 s=-13.5662746429443359 w=12.1540479660034180 -p<br clear="none">><br clear="none">> ## Input and clean topology<br clear="none">> v.in.ogr -e -o input=file.json output=filev2 snap=-1 --verbose --overwrite<br clear="none">> v.in.ogr -e -o input=file2.gpkg output=file2 snap=-1 --verbose --overwrite encoding=UTF-8<br clear="none">><br clear="none">> ## Extract all categories for each dataset so as to maintain polygon (building) contiguity where polygons overlap other polygons within each dataset<br clear="none">> v.extract -d input=filev2 layer=1 type=centroid,area output=buildingsdiss --verbose --overwrite<br clear="none">> v.extract -d input=file2 layer=1 type=centroid,area output=file2diss1 where="cat < 6000001" --verbose --overwrite<br clear="none">><br clear="none">> ## Perform disjoint operation to find buildings that do not intersect buildings in primary dataset<br clear="none">> v.select ainput=file2diss1 alayer=file2diss1 binput=buildingsdiss blayer=buildingsdiss output=DISJOINT1 operator=disjoint --verbose --overwrite<br clear="none">><br clear="none">> ## Output disjoint building layer and cleaned primary dataset<br clear="none">> v.out.ogr -m input=DISJOINT1 output=DISJOINT1.gpkg --verbose --overwrite<br clear="none">> v.out.ogr -m input=filev2 output=filev2.gpkg --verbose --overwrite<br clear="none">><br clear="none">> 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 clear="none">><br clear="none">> projection: 3 (Latitude-Longitude)<br clear="none">> zone: 0<br clear="none">> datum: wgs84<br clear="none">> ellipsoid: wgs84<br clear="none">> north: 5:28:17.649311N<br clear="none">> south: 13:33:58.588333S<br clear="none">> west: 8:09:27.521013W<br clear="none">> east: 51:45:42.670764E<br clear="none">> nsres: 1:00:07.170402<br clear="none">> ewres: 1:00:56.104945<br clear="none">> rows: 19<br clear="none">> cols: 59<br clear="none">> cells: 1121<br clear="none">> NS and EW resolutions are different<br clear="none">> Processing features...<br clear="none">> 0%..........100%<br clear="none">> Processing areas...<br clear="none">> 0% ERROR: Unable to seek: Invalid argument<br clear="none">> NS and EW resolutions are different<br clear="none">><br clear="none">> Any advice appreciated on getting v.select to run. Thanks<br clear="none">><br clear="none">> Best wishes, Chris<br clear="none">><br clear="none">> _______________________________________________<br clear="none">> grass-user mailing list<br clear="none">> <a shape="rect" ymailto="mailto:grass-user@lists.osgeo.org" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br clear="none">> <a shape="rect" href="https://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a></div><br clear="none"><br clear="none"><br clear="none"><br clear="none">-- <br clear="none">Markus Neteler, PhD<br clear="none"><a shape="rect" href="https://www.mundialis.de " target="_blank">https://www.mundialis.de </a>- free data with free software<br clear="none"><a shape="rect" href="https://grass.osgeo.org" target="_blank">https://grass.osgeo.org</a><br clear="none"><a shape="rect" href="https://courses.neteler.org/blog" target="_blank">https://courses.neteler.org/blog</a><div class="yqt4699886086" id="yqtfd21069"><br clear="none"></div></div></div>
</div>
</div></body></html>