Hi,
I understand that you want to distinguish land and ocean in different color. The 'streets' basemap in MATLAB automatically distinguishes land from ocean by using different colors.
For example, you can use the GSHHS dataset or another high-resolution shoreline dataset to add a coastline to your map. You can then use the 'geoplot' function to display the coastline in a contrasting color.
% Define the latitude and longitude limits of the area of interest
latlim = [33.70 33.78]; % limits for Port of Los Angeles
lonlim = [-118.28 -118.22]; % limits for Port of Los Angeles
% Create a geographic axes and set the basemap to 'streets'
ax = geoaxes;
geobasemap(ax, 'streets');
geolimits(latlim, lonlim)
% Load the GSHHS coastline data and display it on the map
load coastlines
% Add some data to the map (replace this with your own data)
geoplot(ax, [33.73 33.75], [-118.25 -118.24], 'k')
In this code, 'geoaxes' is used to create a geographic axes, `geobasemap` is used to set the basemap to 'streets', and `geolimits` is used to set the latitude and longitude limits of the map. The `load` command is used to load the GSHHS coastline data, and `geoplot` is used to display the coastline on the map.
Please note that the `coastlines` data loaded by MATLAB might not be as high resolution as you need. You would need to replace that with your own high-resolution coastline data.
If you want to distinguish land from ocean more clearly, you could consider using a different basemap that provides more contrast, such as 'satellite' or 'topographic'. You can set the basemap using the `geobasemap` function, like this: `geobasemap(ax, 'satellite');`
% Define the latitude and longitude limits of the area of interest
latlim = [33.70 33.78]; % limits for Port of Los Angeles
lonlim = [-118.28 -118.22]; % limits for Port of Los Angeles
% Create a geographic axes and set the basemap to 'streets'
ax = geoaxes;
geobasemap(ax, 'satellite');
geolimits(latlim, lonlim)
% Load the GSHHS coastline data and display it on the map
load coastlines
% Add some data to the map (replace this with your own data)
geoplot(ax, [33.73 33.75], [-118.25 -118.24], 'k')
Please take a look at this documentation for better understanding:
I hope this helps.
