inpoly mex file

Find points inside or on the border of a polygon

Description

A fast mex version of the modified inpoly function.

please run mexme_inpoly.m to compile mex file on your plateform

please run test_inpoly.m for demo

Acknowledgements

The author wishes to acknowledge the following in the creation of this submission:
Fast points-in-polygon test
This submission has inspired the following:
2D random paths generator integrating leg's contraints

• MATLAB release: MATLAB 7.5 (R2007b)
• Other requirements: A C compiler such as MSVC/Intel CPP/GCC

[Darren Engwirda]

Well, I guess this is what you get when you "mex" someone else's code...

This mex file is based on an old version of my "inpoly.m" function. Unfortunately for the author I've made significant improvements to the algorithm since then and a different algorithm is used in my new "inpoly.m" function.

This mex file appears to be slower than my new "inpoly.m" function by a factor of between 3 and 4 on the examples I've tested, even though I have still only coded it as a *.m file.

There is also no "help inpoly" available unless you directly open the *.c file and it also seems that the inputs are required to be transposed when compared with the original "inpoly.m" function.

As with all versions/copies of "inpoly" the point of this code is that it's faster than the MATLAB "inpolygon.m" function, especially for large inputs.

[Sébastien Paris]

Yes in fact, I did this mex file a long times ago. I realized that you update the algorithm. Maybe if I had 1 hour, I'll recreate the new version of this mex-file. If not, I'll remove this submission

[Sébastien Paris]

sorry for my poor english level skill ...

[Martin Podt]

Very useful. On the examples I have tested, the current mex version of inpoly is about 2 times faster than the *.m file.

[Luigi Giaccari]

Very fast!!! well done.

I thing I found an error in the demo file:
The second demo model is buggy: the polygon is self-interecant!!

[Zhang Yanxiang]

it sounds like inpoltgon.

[Korea Univ Saleh]

some points are on the edges of the polygon and the function can not detect that. what is the reason?

[Sebastien Paris]

Please give me an example of this bahaviour.

[John Kua]

Thanks! This has been very useful on a number of projects!

One problem - the function does not handle an empty input array of [2 x 0] and crashes MATLAB. Adding a check of the input dimensionality with a warning if the input is empty and perhaps returning an empty array would be nice.

[Sebastien Paris]

You are right, I'll upload a version with argument checking. Thanks

[Nicolas Billard]

Very good function used on my project. Thanks.

[Matthew Behrend]

great! this was 50 times faster than inpoly.m for my application