<p>A very neat approach for optimal interpolation, which consider not only breaklines along the coast but also anisotropies caused by currents and horizontal diffusion, was given by <a href="http://dx.doi.org/10.1016/S0278-4343(01)00007-3">Lynch &amp; McGillicuddy (2001)</a>. Their approach is quite elegant and takes advantage of the finite element methodology to avoid the transfer of information across boundaries.</p>

CC BY-SA 3.0
2014-07-17T15:54:55.790Z
Oh, that is clever. It's possibly not fully appropriate for the current task, and there's no way that I have time to try to implement it, but I'll file the paper away for future reference. Thanks!
Semidiurnal Simon
CC BY-SA 3.0
2014-07-17T15:32:59.377Z
2014-07-17T15:32:59.377Z
arkaia
<p>Depending on the nature of the data, the approach of </p>

<blockquote>
  <p>Hess, K., R. Schmalz, C. Zervas, and W. Collier, 2004. 'Tidal Constituent and  Residual Interpolation (TCARI): A New Method for the Tidal Correction of Bathymetric Data'. NOAA Technical Report NOS CS 4, Silver Spring, MD, 112 pp.</p>
</blockquote>

<p>(<a href="http://www.nauticalcharts.noaa.gov/staff/docs/TCARI_CS4.pdf">http://www.nauticalcharts.noaa.gov/staff/docs/TCARI_CS4.pdf</a>)</p>

<p>might be appropriate.  This boils down to solving a diffusion equation on an unstructured grid of the domain (TCARI also adds in some additional terms to help extrapolate trends in the data rather than a zero-gradient condition at the shorelines).  With just the diffusion problem, ~100 lines of python was enough to code this up for a San Francisco Bay domain with 500k cells used in the diffusion problem.</p>

<p>Simon Wood and his then PhD student David Lawrence have also done a fair bit of work in this area, though my impression of their methods is that they are fairly expensive computationally for dense ocean observations.  (I'm going to feel silly if Simon W the asker is actually Simon Wood)</p>

CC BY-SA 3.0
2014-07-27T07:01:35.117Z
*grin* No, I'm not Simon Wood. Thanks, that looks intersting - I don't have time to read it in detail now, but will add it to the to-read pile :-)
Semidiurnal Simon
CC BY-SA 3.0
2014-07-26T14:13:21.660Z
2014-07-26T14:13:21.660Z
rusty
<p>Partial answer: after some further research: It looks as though the terminology that I am looking for is "Breakline". I need to define my coastlines as being breaklines. Now I need to figure out how to implement them in the languages that I am using, but that's probably a question for GIS.SE or Stackoverflow.</p>

CC BY-SA 3.0
2014-07-17T14:35:24.700Z
2014-07-17T14:35:24.700Z
Semidiurnal Simon
<p>I have a dataset consisting of irregularly spaced points in the ocean. For each point I have coordinates and a data value.
I need to interpolate to a regular grid. So far, so easy.</p>

<p>There are areas of land within the area covered by the dataset, and in the original dataset there are simply no points there. I have polygons of the land boundaries.</p>

<p>I need to exclude the land from the output. This is fairly straightforwardly accomplished by checking which of the points on the target grid are within the polygons and removing them.
However, I also need to make sure that whatever interpolation technique is used does not consider points that are (for example) the other side of a peninsula - potentially a short distance away in a straight line, but a long way (or possibly not connected at all) by sea.</p>

<p>I'm sure that this must be a solved problem - I imagine any GIS package will have a way to deal with it - but I am at a loss as to what search terms to use in looking for the solution. </p>

<p>I have been debating whether to ask this question here, or on stats.SE, or on GIS.SE, or on Stackoverflow... but let's see how it goes.</p>

<ocean><geospatial>
Interpolating spatial data in contiguous sea areas only
Semidiurnal Simon
# Interpolating spatial data in contiguous sea areas only

(Asked by: username_0 on 2014-07-17)

I have a dataset consisting of irregularly spaced points in the ocean. For each point I have coordinates and a data value.
I need to interpolate to a regular grid. So far, so easy.

There are areas of land within the area covered by the dataset, and in the original dataset there are simply no points there. I have polygons of the land boundaries.

I need to exclude the land from the output. This is fairly straightforwardly accomplished by checking which of the points on the target grid are within the polygons and removing them.
However, I also need to make sure that whatever interpolation technique is used does not consider points that are (for example) the other side of a peninsula - potentially a short distance away in a straight line, but a long way (or possibly not connected at all) by sea.

I'm sure that this must be a solved problem - I imagine any GIS package will have a way to deal with it - but I am at a loss as to what search terms to use in looking for the solution.

I have been debating whether to ask this question here, or on stats.SE, or on GIS.SE, or on Stackoverflow... but let's see how it goes.

**Tags:** ocean, geospatial


---

## 3 answers

### Answer by username_1 on 2014-07-17. Score: 7

A very neat approach for optimal interpolation, which consider not only breaklines along the coast but also anisotropies caused by currents and horizontal diffusion, was given by [Lynch & McGillicuddy (2001)](http://dx.doi.org/10.1016/S0278-4343(01)00007-3). Their approach is quite elegant and takes advantage of the finite element methodology to avoid the transfer of information across boundaries.

> **Comments:**
> **username_0**: Oh, that is clever. It's possibly not fully appropriate for the current task, and there's no way that I have time to try to implement it, but I'll file the paper away for future reference. Thanks!  

---

### Answer by username_2 on 2014-07-26. Score: 7

Depending on the nature of the data, the approach of

> Hess, K., R. Schmalz, C. Zervas, and W. Collier, 2004. 'Tidal Constituent and Residual Interpolation (TCARI): A New Method for the Tidal Correction of Bathymetric Data'. NOAA Technical Report NOS CS 4, Silver Spring, MD, 112 pp.

(<http://www.nauticalcharts.noaa.gov/staff/docs/TCARI_CS4.pdf>)

might be appropriate. This boils down to solving a diffusion equation on an unstructured grid of the domain (TCARI also adds in some additional terms to help extrapolate trends in the data rather than a zero-gradient condition at the shorelines). With just the diffusion problem, ~100 lines of python was enough to code this up for a San Francisco Bay domain with 500k cells used in the diffusion problem.

username_0 Wood and his then PhD student David Lawrence have also done a fair bit of work in this area, though my impression of their methods is that they are fairly expensive computationally for dense ocean observations. (I'm going to feel silly if username_0 W the asker is actually username_0 Wood)

> **Comments:**
> **username_0**: *grin* No, I'm not username_0 Wood. Thanks, that looks intersting - I don't have time to read it in detail now, but will add it to the to-read pile :-)  

---

### Answer by username_0 on 2014-07-17. Score: 2

Partial answer: after some further research: It looks as though the terminology that I am looking for is "Breakline". I need to define my coastlines as being breaklines. Now I need to figure out how to implement them in the languages that I am using, but that's probably a question for GIS.SE or Stackoverflow.


---

