inpolygon usage with big matrices or faster function

44 次查看（过去 30 天）
Ahmet Hakan UYANIK
Ahmet Hakan UYANIK 2022-11-9
评论： Ahmet Hakan UYANIK 2022-11-10
Hello,
I have a matrix for latitudes with the size of 12200x6700 and longitudes with same size as well. I have a polygon of a river as well to check if the points are inside the polygon. However when I do
in = inpolygon(longitude,latitude,Poly.Lon,Poly.Lat);
it never ends.
If I check it by each column with for loop(6700 times), it again takes so much time.
Is there any way to fasten this single line. Like maybe using parfor by seperating the longitude and latitude?
Thanks in advance, Cheers

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

采纳的回答

Bruno Luong
Bruno Luong 2022-11-9
Try this https://fr.mathworks.com/matlabcentral/fileexchange/27840-2d-polygon-interior-detection?s_tid=srchtitle
At some point it was 10 time faster than MATLAB stock function.
But you can start to discard points outside the conves hull first as Matt suggests then refine with this non convex polygonal function.
  1 个评论
Ahmet Hakan UYANIK
Ahmet Hakan UYANIK 2022-11-10
Regarding big size of matrices, this tool works way more faster than MATLAB stock function(even much more faster than 10 times :) Thanks in advance.

请先登录，再进行评论。

更多回答（1 个）

Matt J
Matt J 2022-11-9
编辑：Matt J 2022-11-9
Is it a convex polygon? There are much faster alternatives to inpolygon for convex polygons,e.g.,
https://www.mathworks.com/matlabcentral/fileexchange/10226-inhull?s_tid=srchtitle
Even if it is not a convex polygon, perhaps you can frst discard points that are not in the convex hull of the polygon, thereby reducing the amount of work that inpolygon needs to do.
  2 个评论
Bruno Luong
Bruno Luong 2022-11-9
"I have a polygon of a river"
Personally I never seen a convex river. ;-)
Matt J
Matt J 2022-11-10
编辑：Matt J 2022-11-10
Some sections of the Amazon boundary some fairly convex-looking bits of land:

请先登录，再进行评论。

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by