dividing coastline into equal interval

2025 2 6
1 个回答
更新时间：2025 2 7
11 次查看（30 天）

0 个投票

Hello, I needed to divide the coastline into equal interval (as shown in the image) and identify latitude and longitude of each segments of these coastline.
Basically, for each of these segment, I wish to count the frequency of storms hitting/year.
For this, I did the following:
coastlines = load('C:\Program Files\MATLAB\R2024a\toolbox\map\mapdata\coastlines.mat');
LatLimit = [15 44]; LonLimit = [60 104];
X = coastlines.coastlat((coastlines.coastlat >= LatLimit(1,1) & coastlines.coastlat <= LatLimit(1,2)),:);
Y = coastlines.coastlon((coastlines.coastlon >= LonLimit(1,1) & coastlines.coastlon <= LonLimit(1,2)),:);
However, I am unable to divide it each segmentwise, considering linear distance interval of say, 50 km and identify latitude-longitudes within those segments. Any help in this regard is deeply appreciable.

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

回答（1 个）

Walter Roberson
Walter Roberson 2025-2-6

0 个投票

See John D'Errico's File Exchange contribution interparc https://www.mathworks.com/matlabcentral/fileexchange/34874-interparc which is designed for interpolating curves at fixed distances traveled.

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

Poulomi
Poulomi 2025-2-6
but how to implement this function for the coastline bounded by the fixed latitude and longitude bounds?
Walter Roberson
Walter Roberson 2025-2-6
t = linspace(0,1,250);
pt = interparc(t, X, Y);
Poulomi
Poulomi 2025-2-7
Thanks, but still it is not clear how can I divide the coastline bounded by certain latitude, longitude box, with this function at a linear increment of 50 km.
Walter Roberson
Walter Roberson 2025-2-7

Assuming that the distances are such that Euclidean approximation is good enough, and assuming that the coast is continuous:

You can calculate the total distance along the coast using

 D = sum(sqrt(diff(X).^2 + diff(Y).^2))

This allows you to calculate the fractional increment such that the adjacent distance is 50 km.

Walter Roberson
Walter Roberson 2025-2-7

I said "assuming that the coast is continuous". It appears that works for your situation. It would be a nuisance if it did not, partly because it would not be defined what distance should be measured over any jump in the coast.

Walter Roberson
Walter Roberson 2025-2-7

Looking at the map, I am not certain that using Euclidean approximation of the total distance would be good enough for your purposes. On the other hand, interparc expects Euclidean coordinates I expect.

You might need to convert map coordinates to some kind of linear coordinates before using interparc

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Descriptive Statistics and Visualization 的更多信息

产品

版本

R2024a

标签

Community Treasure Hunt

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

Start Hunting!

Translated by