File Exchange

image thumbnail

2D Line Curvature and Normals

version 1.3.0.0 (8.97 KB) by Dirk-Jan Kroon
Accurate Curvature and Normals of a line/contour consisting of 2D points

89 Downloads

Updated 17 Oct 2011

View License

Editor's Note: This file was selected as MATLAB Central Pick of the Week

LineCurvature2D, This function calculates the curvature of a 2D line. It first fits polygons to the points. Then calculates the analytical curvature from the polygons.

K = LineCurvature2D(Vertices,Lines)

inputs,
Vertices : A M x 2 list of line points.
(optional)
Lines : A N x 2 list of line pieces, by indices of the vertices
(if not set assume Lines=[1 2; 3 4 ; ... ; M-1 M])

outputs,
K : M x 1 Curvature values

.
.
LineNormals2D, This function calculates the normals, of the line. Using the neighbouring points of each line/contour point, and forward an backward differences on the end points

N=LineNormals2D(V,L)

inputs,
V : List of points/vertices 2 x M
(optional)
Lines : A N x 2 list of line pieces, by indices of the vertices
(if not set assume Lines=[1 2; 3 4 ; ... ; M-1 M])

outputs,
N : The normals of the Vertices 2 x M

Comments and Ratings (32)

Luxiao Cui

Bruno

Dear all, the curvature is calculated simply by the formula in the following link, which is a result of planar differential geometry. https://math.stackexchange.com/questions/1250554/curvature-of-plane-parametric-curves

Bruno

Dear Dirk,
Could you refer as to the math behind the last step of the calculation?Is there any literature? That is, how you get the curvature from the fitted polynomial?
bnscocozza@gmail.com

jian fu

Dear Dirk-Jan Kroon,
thanks for this great function! Is there any chance for me to get the related papers explaining this method?
Thanks a lot!
1197349243@qq.com

Aurea94

Aitor Arana

It's a very nice code but I would like to know the mathematical foundation. Does anybody know it??

ZQ W

you23cef

Anwer khan

Ben King

Bojan

Hi . Can this code be used to break the curved line into small line segments of appropriate size with their ends overlapped to some small extent and then take hough transform?

can the code be used to find singular points in fingerprint? if yes how plz? thx. Ajit

From what i understand the algorithm only uses the direct neighbors of the point to be evaluated correct? is it possible to include multiple nearby points to get a more accurate result?

Khalid Ali

What is meant by lines?

rabeeah

very useful code. i would also like to know the mathematical model used. does it use the freeman chain code Dirk?

Michael

Thank you, good code to find in my case finger tips on a hand scan

Yunfeng

Simon

Is there any literature to back this up??
Also, the grammar in the comments needs serious work....

Willy

Thank you! Very usefull code! I would like to understand the maths behind the code better. Is this code based on some literature?

I would like to point out a small typo:-
*****(correct)****
if(nargin<2)
L=[(1:(size(V,1)-1))' (2:size(V,1))'];
end
*****
In the code I saw,the above was :-
*****(typo)*****
Lines=[(1:(size(Vertices,1)-1))' (2:size(Vertices,1))'];
*****

Updates

1.3.0.0

Corrected typo reported by Ravikiran Janardhana

MATLAB Release Compatibility
Created with R2011a
Compatible with any release
Platform Compatibility
Windows macOS Linux
Acknowledgements

Inspired: thick_line(x,y,thickness)