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

Invalid geometry error #520

Open
shawnpatti opened this issue Nov 12, 2019 · 5 comments
Open

Invalid geometry error #520

shawnpatti opened this issue Nov 12, 2019 · 5 comments

Comments

@shawnpatti
Copy link

The following

gips_process landsat -t 019032 -d 2018-03-22 --chunksize 512.0 --format GTiff --numprocs 1 --verbose 1 -p ndti-toa-coreg crc-toa-coreg ndvi-toa-coreg lswi-toa-coreg cloudmask-coreg ndsi-toa-coreg

Returns error

No handlers could be found for logger "shapely.geos"
Fatal: 1 error(s) occurred:
Problem with running AROP: The operation 'GEOSIntersection_r' could not be performed. Likely cause is invalidity of the geometry <shapely.geometry.polygon.Polygon object at 0x7fadd4aa8490>
[1mGIPS Data Processing (v0.16.0)[0m
utm_zone is 17
@ra-tolson
Copy link
Collaborator

What's the specific commit or branch you're working with? Is it the top of the dev branch? 0.16.0 seems to be unreleased in any case.

@shawnpatti
Copy link
Author

My branch is nwf-ccc-dagan which is based on dev branch.

@ra-tolson
Copy link
Collaborator

I get a different error; I'm not sure how we can be on the same setup. My config is:

  • ran gips_process with -v99 --fetch
  • my config set landsat fetch to usgs and s2 fetch to esa
  • I'm on python 3 & a recent build of the gips 'full test' container
  • my dev branch commit is fb1825e Merge branch 'add-fisk-orcid' into 'dev'
Fatal: 1 error(s) occurred:
Problem with running AROP:
Traceback (most recent call last):
  File "/gips/gips/data/landsat/landsat.py", line 1206, in process
    s2_export = self.sentinel2_coreg_export(tmpdir_fp)
  File "/gips/gips/data/landsat/landsat.py", line 1908, in sentinel2_coreg_export
    landsat_footprint = wkt_loads(self.assets[next(iter(self.assets))].get_geometry())
  File "/gips/gips/data/core.py", line 536, in get_geometry
    return self.get_geofeature().WKT()
  File "/usr/local/lib/python3.6/dist-packages/gippy-1.0.3-py3.6-linux-x86_64.egg/gippy/gippy.py", line 4852, in <lambda>
    __getattr__ = lambda self, name: _swig_getattr(self, GeoFeature, name)
  File "/usr/local/lib/python3.6/dist-packages/gippy-1.0.3-py3.6-linux-x86_64.egg/gippy/gippy.py", line 59, in _swig_getattr
    raise AttributeError(name)
AttributeError: WKT

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/gips/gips/utils.py", line 647, in cli_error_handler
    yield
  File "/gips/gips/data/landsat/landsat.py", line 1208, in process
    except NoSentinelError:
NameError: name 'NoSentinelError' is not defined

@ircwaves
Copy link
Collaborator

@ags-tolson -- your WKT error seems like a portion of the code that has not been migrated to gippy1==1.0.

@ircwaves
Copy link
Collaborator

gips==0.16.0 is tagged, but wasn't promoted to a release.

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