Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in
  • arosics arosics
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 14
    • Issues 14
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 3
    • Merge requests 3
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Packages and registries
    • Packages and registries
    • Container Registry
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Repository
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • Daniel Scheffler
  • arosicsarosics
  • Issues
  • #27
Closed
Open
Issue created Feb 13, 2019 by Daniel Scheffler@danschefOwner

KeyError in case input datasets do not have WGS84 as geographic datum

This was raised in case the input datasets have EPSG:3035 as projection:

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "C:\conda\arosics\lib\site-packages\arosics\CoReg.py", line 1330, in calculate_spatial_shifts
    self._get_updated_map_info()
  File "C:\conda\arosics\lib\site-packages\arosics\CoReg.py", line 1341, in _get_updated_map_info
    original_map_info = geotransform2mapinfo(self.shift.gt, self.shift.prj)
  File "C:\conda\arosics\lib\site-packages\py_tools_ds\geo\map_info.py", line 216, in geotransform2mapinfo
    return Geocoding(gt=gt, prj=prj).to_mapinfo()
  File "C:\conda\arosics\lib\site-packages\py_tools_ds\geo\map_info.py", line 49, in __init__
    self.from_geotransform_projection(gt, prj)
  File "C:\conda\arosics\lib\site-packages\py_tools_ds\geo\map_info.py", line 111, in from_geotransform_projection
    self.datum = 'WGS-84' if proj4['datum'] == 'WGS84' else proj4['datum']  # proj4['ellps']?
KeyError: 'datum'
Assignee
Assign to
Time tracking