Vector plots draw an unintended constant longitude line about dateline #757
Comments
|
The same behavior is observed with geoviews import geoviews as gv
projection = ccrs.Orthographic(central_longitude=180)
crs = ccrs.PlateCarree(central_longitude=180)
gv.Path(gdf_180_lon, crs=crs).opts(projection=projection)
|
|
Can you see if matplotlib does the same? |
Same behavior. import geoviews as gv
import holoviews as hv
hv.extension("matplotlib")
projection = ccrs.Orthographic(central_longitude=180)
crs = ccrs.PlateCarree(central_longitude=180)
gv.Path(gdf_180_lon, crs=crs).opts(projection=projection)
|
|
Oh I meant pure matplotlib :P |
Oops! Here's a Matplotlib plot using a PolyCollection. There's some UXarray logic here, though the GeoDataFrames I generated above were also from UXarray. projection = ccrs.Orthographic(central_longitude=180)
crs = ccrs.PlateCarree(central_longitude=180)
# shift central coordinates by applying a PlateCarree at 180 degrees
pc = uxds['bottomDepth'].to_polycollection(periodic_elements="exclude", projection=crs)
pc.set_antialiased(True)
pc.set_cmap("Blues")
fig, ax = plt.subplots(
1,
1,
figsize=(10, 5),
facecolor="w",
constrained_layout=True,
subplot_kw=dict(projection=projection), # project using matplotlib, with source PlateCarree at 180 lon
)
ax.add_feature(cfeature.COASTLINE)
ax.add_feature(cfeature.BORDERS)
ax.add_collection(pc)
ax.set_global()
Switching |
|
So it seems like there's something wrong with GeoViews projecting it. Moving to geoviews. |
|
Maybe you can poke around here and see if you can spot the issue here: |
Will take a look! |
|
Initially I was thinking that the line may be a reside of the frame surround the plot. After setting the following, the line is still there. backend_opts={
"axes.frame_on": True,
} |
|
It's likely related to cartopy / matplotlib likely does something internally to fill it(?) |
ALL software version info
Software Version Info
Description of expected behavior and the observed behavior
I am working with a mesh of polygons that I want to project and visualize over the dateline (180 degrees longitude). I have two datasets:
I've put a few code examples below. For vector plots, there is line of longitude being drawn at 180 degrees, which is not observed in the rasterized plot or the un-projected one.
With the elements
Complete, minimal, self-contained example code that reproduces the issue
Shapefile of
gdf_0_lon: https://drive.google.com/file/d/1i35d7tCvXOTm1EcmHsKHWn3iTHIHUNm2/view?usp=drive_linkShapefile of
gdf_180_lon: https://drive.google.com/file/d/12hYPfCbDpc0gv7ohzX_yKdxufOC1dWos/view?usp=drive_linkStack traceback and/or browser JavaScript console output
N/A
The text was updated successfully, but these errors were encountered: