File takes as input vector of points (lat/lon) and coastal resolution and returns vector of indicators for whether these points are over land or ocean. Makes geographic plot too, if desired.
Primary code was given in a response to a question by Brett Shoelson in Feb 2011: http://www.mathworks.com/matlabcentral/answers/1065-determining-whether-a-point-on-earth-given-latitude-and-longitude-is-on-land-or-ocean
As far as I can tell, this code has not been put into a function file for public use.
Dan Chavas (2021). land_or_ocean.m (https://www.mathworks.com/matlabcentral/fileexchange/45268-land_or_ocean-m), MATLAB Central File Exchange. Retrieved .
Inspired: borders
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Create scripts with code, output, and formatted text in a single executable document.
Great function!
I found some of locations that produce the following error message:
---------------------------------------
Error using delaunayTriangulation/pointLocation
The triangulation is empty.
Error in inpolyfast (line 20)
si = dt.pointLocation(x(:),y(:));
Error in vec2mtx>cellCentersInPolygon (line 157)
in = inpolyfast(glon, glat, lon, lat);
Error in vec2mtx (line 109)
Z(~cellCentersInPolygon(lat, lon, sizeZ, R1)) = 2;
Error in land_or_ocean (line 65)
[Z, R] = vec2mtx(coast.lat, coast.long, ...
---------------------------------------
Some instances are the (lat,lon) couples (-19,41), (-68,163), (47,-127), (74,-72).
It looks like these locations might be at the same (or very close to same) distance from two different shores.