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

iber projection lost in the North Pacific #1240

Closed
gerritholl opened this issue Jun 24, 2020 · 2 comments
Closed

iber projection lost in the North Pacific #1240

gerritholl opened this issue Jun 24, 2020 · 2 comments

Comments

@gerritholl
Copy link
Member

Describe the bug

According to its description the satpy built-in iber projection should cover "North half of the Iberian Peninsula and the Gulf of Biscay image 0 degrees". However, as defined it is somewhere in the remote North Pacific.

To Reproduce

import pyresample
ar = pyresample.load_area("/home/gholl/checkouts/satpy/satpy/etc/areas.yaml",
        "iber")
(lons, lats) = ar.get_lonlats()
print(lons.min(), lons.mean(), lons.max(), lats.min(), lats.mean(), lats.max())

Expected behavior

I expect the range of longitudes and latitudes to correspond to the described area, somewhere near the meridian, perhaps a little west of it.

Actual results

In reality the area covers longitude 166 to 180 degrees, seemingly shifted by 180 degrees compared to where it should be expected:

166.3308164211224 173.2332707794299 179.85396117722175 39.629863255964146 42.60892251305624 45.41767241192362

Environment Info:

  • Satpy Version: 0.22.1
  • Pyresample Version: 1.16.0
@djhoese
Copy link
Member

djhoese commented Jun 29, 2020

My guess is the description was a copy/paste error from the sve area right above iber:

sve:
description: North half of the Iberian Peninsula and the Gulf of Biscay image 0
degrees
projection:
init: epsg:3006
shape:
height: 2000
width: 2000
area_extent:
lower_left_xy: [-342379.698, 6032580.06]
upper_right_xy: [1423701.52, 8029648.75]
iber:
description: North half of the Iberian Peninsula and the Gulf of Biscay image 0
degrees
projection:
proj: utm
shape:
height: 1000
width: 2000
area_extent:
lower_left_xy: [-342379.698, 4432580.06]
upper_right_xy: [723701.52, 5029648.75]
units: m

Actually...nevermind that description seems more accurate for iber than sve. I'm a little confused by the projection for iber. It looks like it should be a UTM zone but is missing the actual UTM zone number. PROJ doesn't even define a default: https://proj.org/operations/projections/utm.html

@gerritholl
Copy link
Member Author

This projection was deleted by commit 4126039.

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

2 participants