-
Notifications
You must be signed in to change notification settings - Fork 20
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
Coordinate priority with non-spherical CRS #490
Comments
I was able to bypass the error by using your suggestion
but it only worked for a bit. I'm getting the same error again.
I'm attempting to sum the variable by day/month/year for every shape in the shapefile using the following command:
Thanks in advance! |
Bummer! A few quick questions:
I cannot think of any recent changes that would cause this. It's possible (if the shapfile changed) that you do have some empty subsets in there. In that case, you can try setting |
It must be the shapefile. I get the error I posted about with one. It succeeds with a few others. |
It was definitely the shapefile. It appears very small features were giving me the error. |
Thanks for the update. Sorry for the slow reply - was taking some PTO. Do you think this is a problem that needs to be fixed in ocgis? I could take a look at the shapefile if that would help. |
Take your much deserved PTO! That's completely up to you. I ended up iterating through selections of my shapefile. If it worked I selected more and if it didn't I deselected a few I previously added. It may be nice to have ocgis keep running the calculation and report a list of features that didn't process because of whatever exception. |
Done! 😉 Your approach makes sense. When you get a chance, could you pass along the shp and netcdf (or just the spatial resolution of the grid if that's easier)? What you described would be a nice feature. |
@mkmitchell reported an issue with a NARR dataset having a Lambert Conformal Conic projection. The lat/lon coordinate variables were marked with axis x/y attributes. The Cartesian x/y coordinate variables had no axis attributes.
ocgis
errored out during a subset as it chose the LCC projection correctly but picked the wrong coordinates because of axis attribute priority during coordinate selection. This was fixed by updating the dimension map prior to the subset.To fix this, coordinate variable selection should attempt to match coordinate variable units against the units recommended by the coordinate system prior to using the axis attributes.
The text was updated successfully, but these errors were encountered: