Hi @yichen, this issue occurs because when 'Origin',[0 180 0] is set in axesm, the 180° meridian becomes the center of the map. The landareas.shp file defines longitude values in the range of −180° to 180°, so when features in the shapefile cross the dateline (±180°), they appear split due to the way the coordinates are stored and displayed.
To display these features correctly—without setting 'EdgeColor','none'—the longitude values must be adjusted to match the new map origin.
1. Shift longitude values relative to the new origin (e.g., 180°):
function shifted = shiftLonToCenter(shapes, centerLon)
shifted = shapes;
for k = 1:numel(shapes)
lon = shapes(k).Lon;
if ~isempty(lon)
shifted(k).Lon = wrapTo180(lon - centerLon);
end
end
end
2. Apply the shift to the land data:
land = shaperead('landareas.shp', 'UseGeoCoords', true);
landShifted = shiftLonToCenter(land, 180);
3. Display the adjusted land areas with the correct map origin:
figure;
ax = axesm('eqdcylin', 'Origin', [0 180 0], ...
'MapLatLimit', [-90 90], 'MapLonLimit', [-180 180], ...
'Frame', 'on', 'Grid', 'off');
setm(ax, 'FLineWidth', 1);
geoshow(landShifted, 'FaceColor', [0.8 0.8 0.8], 'EdgeColor', 'black');
This approach ensures that polygons crossing the dateline are properly displayed with the new origin, and that the edges remain visible and correctly connected.
