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

[Markus Neteler]

Hi,

Is there a chance to receive the dataset for testing (or, ideally, a
reproducible example with the North Carolina sample dataset from
https://grass.osgeo.org/download/data/)?

thanks,
Markus

On Wed, May 12, 2021 at 1:14 PM Christopher Lloyd via grass-user
<grass-user at lists.osgeo.org> wrote:
>
> 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 topology
> v.in.ogr -e -o input=file.json output=filev2 snap=-1 --verbose --overwrite
> v.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 dataset
> v.extract -d input=filev2 layer=1 type=centroid,area output=buildingsdiss --verbose --overwrite
> v.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 dataset
> v.select ainput=file2diss1 alayer=file2diss1 binput=buildingsdiss blayer=buildingsdiss output=DISJOINT1 operator=disjoint --verbose --overwrite
>
> ## Output disjoint building layer and cleaned primary dataset
> v.out.ogr -m input=DISJOINT1 output=DISJOINT1.gpkg --verbose --overwrite
> v.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:       0
> datum:      wgs84
> ellipsoid:  wgs84
> north:      5:28:17.649311N
> south:      13:33:58.588333S
> west:       8:09:27.521013W
> east:       51:45:42.670764E
> nsres:      1:00:07.170402
> ewres:      1:00:56.104945
> rows:       19
> cols:       59
> cells:      1121
> NS and EW resolutions are different
> Processing features...
>    0%..........100%
> Processing areas...
>    0% ERROR: Unable to seek: Invalid argument
> NS and EW resolutions are different
>
> Any advice appreciated on getting v.select to run. Thanks
>
> Best wishes, Chris
>
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user



-- 
Markus Neteler, PhD
https://www.mundialis.de - free data with free software
https://grass.osgeo.org
https://courses.neteler.org/blog