How I can count the number lines falling inside the green region

1 view (last 30 days)
figure(2);
rows = 1;
columns = 1;
for subplotIndex = 1 : rows * columns
subaxis(rows, columns, subplotIndex, 'spacing', 0.01, 'margin', 0.01);
m_proj('stereographic', 'lat', 90, 'long', 30, 'radius', 30.1);
m_coast('patch', [.7 .7 .7], 'edgecolor', 'k');
m_grid('xtick', -180:90:180, 'tickdir', 'out', 'ytick', [70 80], 'linest', '-');
end
hold on;
a = 'F:\hysplit_output\ny_alesund\ny_18_19_djf\*.';
files = dir(a);
for jj = 1:numel(files)
fullFilePath = fullfile(files(jj).folder, files(jj).name);
if ~files(jj).isdir
headerLines = 7;
data = readtable(fullFilePath, 'HeaderLines', headerLines);
X = data.Var11;
Y = data.Var10;
[X_map, Y_map] = m_ll2xy(X, Y);
plot(X_map, Y_map, '-b', 'LineWidth', 0.1);
hold on;
M = m_shaperead('first_r');
for k = 1:length(M.ncst)
m_line(M.ncst{k}(:,1), M.ncst{k}(:,2), 'linewidth', 4, 'color', 'g');
end
[xi, yi] = polyxpoly([X_map, Y_map].', [M.ncst{1}(:, 1), M.ncst{1}(:, 2)].');
end
end
  3 Comments

Sign in to comment.

Answers (1)

Matt J
Matt J on 6 Feb 2024
Edited: Matt J on 6 Feb 2024
green=polyshape(vertcat(M.ncst));
Lines=findobj(gca,'Type','line');
clear insideGreen
for i=numel(Lines):-1:1
points=[Lines(i).XData;Lines(i).YData]';
insideGreen(i)= any(isinterior(green,points));
end
total=sum(insideGreen)

Products


Release

R2021a

Community Treasure Hunt

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

Start Hunting!