How I can count the number lines falling inside the green region
1 次查看(过去 30 天)
显示 更早的评论
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
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1609416/image.png)
3 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Map Display 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!