How to change land data to NaN?

5 次查看（过去 30 天）

eevee 2014-10-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/160188-how-to-change-land-data-to-nan

回答： Chad Greene 2014-12-5

Hi,

I have plotted a bathymetry and now I'm given coastline and some islands files where I need to use them to define where the land is and then turn the land data to NaN.

How should I do that?

Thanks, Eve

4 个评论
显示 2更早的评论隐藏 2更早的评论

Geoff Hayes 2014-10-29

Eve - how do you (currently) calculate the volume? What information do you use for that?

eevee 2014-10-30

Hi Geoff,

I use:

dx*dy*zi(~isnan(zi))

However, my zi includes both land and sea data. So now, I need to remove the land data before I could calculate the volume of bathymetry. I don't know how to tell matlab that my land data are those that falls below the coastline (please see image below)

Thanks,

Eve

请先登录，再进行评论。

请先登录，再回答此问题。

回答（2 个）

Kelly Kearney 2014-10-27

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/160188-how-to-change-land-data-to-nan#answer_156742

The inpolygon function will do most of the work. Depending on the scale of your region and the format of the coastlines, you may need to close off polygons, reverse order of vertices, and upsample polygons in order to get the best results; the Mapping Toolbox has several functions ( flatearthpoly, interpm, ispolycw, poly2cw) to help with that.

6 个评论
显示 4更早的评论隐藏 4更早的评论

eevee 2014-10-31

sub14n11test.m

Hi Kelly,

sorry I just realized I forgot to attach my script. :)

Kelly Kearney 2014-10-31

在 MATLAB Online 中打开

Can you upload the coastline file? The answer will depend on whether it contains just the visible line segments you plotted above, or fully-closed polygons. Assuming it's the latter, your code to mask the land points will look something like:

isin = inpolygon(xi, yi, xcoast, ycoast);
zi_diff(isin) = NaN;

But you may need to do some processing of xcoast and ycoast to convert them to closed, clockwise polygons.

请先登录，再进行评论。

Chad Greene 2014-12-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/160188-how-to-change-land-data-to-nan#answer_161437

Another solution is the landmask function on File Exchange.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

类别

Mathematics and Optimization Mapping Toolbox Data Analysis Vector Data

Radar Mapping Toolbox Data Analysis Vector Data

在 Help Center 和 File Exchange 中查找有关 Vector Data 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by