-
Notifications
You must be signed in to change notification settings - Fork 0
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
Jbacon ims latlong #53
Conversation
0f2f3b4
to
f78c34a
Compare
Delete icedyno/preprocess/preprocess_netcdf.ipynb Incorrect formatting Add files via upload Delete icedyno/preprocess/preprocess_netcdf.ipynb Improper formatting Add files via upload
…ordinates from x,y values and equations
…st, more to come. Had to include init for tests to see the icedyno project to import
cropped pipeline: Parallelized workflow to crop original IMS netcdf files to centered window at variable window size with luigi
…reprocess/geolocation. Returning x, y in meters and pinning constant values to IMS specific ones. Updated test to reflect change
…ten tests suggests.
…nted out print statements
f78c34a
to
9e5147c
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Are the cells in ims_latlong.ipynb
that's not run still WIP?
…ule, cleanup latlon notebook, addfind closest index in grid function to geolocation module.
…ion for lat/lon crop windows
It's cleaned up now! And ready for review |
… found by manual work in colab. Corrects out of memory error with netcdf rotations with temporary solution to an initial crop before grid correction.
…ression settings in original netcdf file
9d8b147
to
66c6baf
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks great. I have aa question
…ould look like with google maps, added descriptions
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This looks perfect! Thank you!
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Great call on the netcdf crop, this looks good! Thank you!
Major features:
Details on Grid Rotation Solution
We can't rotate the entire netcdf image due to its size, but we also can't select the right window without correcting grid with a rotation first. Temporary solution is to crop to a large, but much smaller than original image first, then rotate and crop to correct window. Throws an assert error if the first crop excludes the selected window.
Details on lat/lon conversion:
Removed northern hemisphere considerations and removes a set of lines in the original
polar_xy_to_lonlat
that causes the longitude to become invalid (>180) for certain lonlat pairs.https://github.com/nsidc/polarstereo-lonlat-convert-py/blob/366069fc65e781c277adf6d2fb0dac60c521d77d/polar_convert/polar_convert.py#L60C1-L61C38
Previously, the tests I wrote to check that the conversions should undo each other doesn't pass with that re-wrapping of the longitude value. Now it does.
The original functions from latlon to xy actually don't have any tests: https://github.com/nsidc/polarstereo-lonlat-convert-py/blob/main/polar_convert/test.py