Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

No areas of sufficient density - no shapes. Decrease density or area size #26

Open
geographika opened this issue Oct 12, 2022 · 5 comments

Comments

@geographika
Copy link

I have built the hexer project on Windows, using GDAL 3.04. The help command works fine, however I've not yet been successful in getting any valid outputs.

curse boundary points.shp --output myboundary.shp produces an empty shapefile with no records.

curse density points.shp --output D:\Temp\myhexagons.shp --edge 100 always ends up returning the following error:

No areas of sufficient density - no shapes. Decrease density or area size

I've played around with the --edge and --count settings, and also run the command on various point shapefiles, and LAS files, including a point shapefile generated with the make_shapes.py script.

Is anyone able to get outputs with recent versions of GDAL, or on Windows? Or maybe I've missed a something basic..

@abellgithub
Copy link
Collaborator

abellgithub commented Oct 12, 2022 via email

@hobu
Copy link
Collaborator

hobu commented Oct 12, 2022

Something must have rotted. This code hasn't been used in this form for quite a while.

@hobu
Copy link
Collaborator

hobu commented Oct 12, 2022

This has really rotted. I took a stab at it and couldn't get it working in 15 minutes of futzing. I would try to use PDAL directly to do this nowadays, as that's where this code has migrated to.

  • Write out a text file
  • Read the text file with readers.text
  • Run pdal boundary or pdal tindex on it.

@geographika
Copy link
Author

Does your input shape file consist of points?

Yes, I tried a shapefile generated from a LAS file, a Python generated one, and a couple more point shapefiles I had to hand all with the same result. I also tried converting to web mercator in case that was an issue (and as used at https://gist.github.com/tmcw/7436261).

This has really rotted. I took a stab at it and couldn't get it working in 15 minutes of futzing. I would try to use PDAL directly to do this nowadays, as that's where this code has migrated to.

* Write out a text file

* Read the text file with `readers.text`

* Run `pdal boundary` or `pdal tindex` on it.

Thanks @hobu for investigating and the reply. As far as I can tell the boundary function has moved to PDAL, but there is no equivalent of the density function for creating the hexagon outputs - is this correct?

A shame this project is now deprecated - the hexer, curse names are inspired!

@geographika
Copy link
Author

Ignore comment above, the hexagons are available as described at https://pdal.io/en/stable/workshop/exercises/analysis/density/density.html

pdal density test.las -o test.shp

image

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants