<html><head></head><body><div class="ydpa5ec09d1yahoo-style-wrap" style="font-family:verdana, helvetica, sans-serif;font-size:13px;"><div></div>
<div dir="ltr" data-setdir="false">Hi Markus, Having split my input data into smaller sized shp files, these now process fine using v.select. So clearly the problem that I had lies with some file size limitation with the v.select 'disjoint' algorithm - also failing with one input file prior using the v.extract algorithm. The original (large) input files input to grass ok using v.in.ogr, but then failed when using v.select. Might you be able to indicate the file size limitation for the v.select algorithm and the reasons for this? 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 id="yahoo_quoted_1833985760" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
On Wednesday, 19 May 2021, 14:29:56 BST, Christopher Lloyd via grass-user <grass-user@lists.osgeo.org> wrote:
</div>
<div><br></div>
<div><br></div>
<div><div id="yiv5771353777"><div><div class="yiv5771353777ydp3995c95cyahoo-style-wrap" style="font-family:verdana, helvetica, sans-serif;font-size:13px;"><div></div>
<div dir="ltr">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"><br clear="none"></div><div dir="ltr">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"><br clear="none"></div><div dir="ltr">Best wishes, Chris</div>
</div><div class="yiv5771353777yqt1461114019" id="yiv5771353777yqt41360"><div class="yiv5771353777yahoo_quoted" id="yiv5771353777yahoo_quoted_1695427563">
<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 clear="none"></div>
<div><br clear="none"></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 rel="nofollow noopener noreferrer" shape="rect" target="_blank" href="https://grass.osgeo.org/download/data/">https://grass.osgeo.org/download/data/</a>)?<br clear="none"><br clear="none">thanks,<br clear="none">Markus<br clear="none"><div class="yiv5771353777yqt4699886086" id="yiv5771353777yqtfd41237"><br clear="none">On Wed, May 12, 2021 at 1:14 PM Christopher Lloyd via grass-user<br clear="none"><<a rel="nofollow noopener noreferrer" shape="rect" ymailto="mailto:grass-user@lists.osgeo.org" target="_blank" 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 rel="nofollow noopener noreferrer" shape="rect" ymailto="mailto:grass-user@lists.osgeo.org" target="_blank" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br clear="none">> <a rel="nofollow noopener noreferrer" shape="rect" target="_blank" href="https://lists.osgeo.org/mailman/listinfo/grass-user">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 rel="nofollow noopener noreferrer" shape="rect" target="_blank" href="https://www.mundialis.de ">https://www.mundialis.de </a>- free data with free software<br clear="none"><a rel="nofollow noopener noreferrer" shape="rect" target="_blank" href="https://grass.osgeo.org">https://grass.osgeo.org</a><br clear="none"><a rel="nofollow noopener noreferrer" shape="rect" target="_blank" href="https://courses.neteler.org/blog">https://courses.neteler.org/blog</a><div class="yiv5771353777yqt4699886086" id="yiv5771353777yqtfd21069"><br clear="none"></div></div></div>
</div>
</div></div></div></div><div class="yqt1461114019" id="yqt39276">_______________________________________________<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><br clear="none"></div></div>
</div>
</div></body></html>