Artifacts from lat/lon interpolation across -180/180
bnb32 opened this issue · 1 comments
bnb32 commented
Bug Description
When the spatial domain extends across -180/180 longitude plots show an artifact at this boundary. The lat/lon grid could be incorrectly mapped to the fields or it could be that the RBFInterpolate call doesn't work as well across the -180/180 boundary (RBFInterpolate uses n_neighbors and probably doesn't know to use neighbors across the -180/180 boundary, even though they are "close"). It's also possible that the output array isn't getting filled at this boundary for some reason.
Code
plt.scatter(lons.flatten(), lats.flatten(), c=windspeed.flatten())
plt.contourf(lons.flatten(), lats.flatten(), windspeed.flatten())
plt.imshow(windspeed)
All of the above show the artifact.