File Exchange

image thumbnail


version (42.4 KB) by Dan Chavas
Determine if input points (lat/lon) are over land or ocean


Updated 27 Jan 2014

View License

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:
As far as I can tell, this code has not been put into a function file for public use.

Cite As

Dan Chavas (2021). land_or_ocean.m (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (2)


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.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.

MATLAB Release Compatibility
Created with R2013b
Compatible with any release
Platform Compatibility
Windows macOS Linux

Inspired: borders

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!