Active zoom for polyline selection.
Have you ever wanted to carefully select points from an image but could not get close enough to see where you wanted to select? Have you ever wanted to create a smooth, evenly spaced set of parametric contour points from a minimum of user-selected points? Then this function is for you -- it allows the user to select points from an image at a user specified zoom level that moves along with the points selected in a centered frame of reference. The resulting contour can then be interpolated to constant arc-length resolution, smoothed, and even closed into a loop.
Im1 is the input image (matrix) from which points will be selected.
The optional parameter 'interp' creates evenly spaced X Y points with the resolution 'res', where res = 1, will space the X,Y points one pixel apart in arc length. The secondary option 'method' selects between a 'linear' and 'spline' interpolation, the latter being a differentially smooth curve. The default method is 'linear'.
The optional parameter 'loop_path' will loop the path so that the last point is the same as the first point. If the method is 'spline', then the looped path will also be smoothed such that the derivatives match at the beginning and ending points.
The optional parameter 'stationary' will keep the zoom and image position fixed, but all other options will be active.
Hold 'shift' or 'control' and then click to end polyline selection.
This script requires the 'ginput' function.
This script will generate its own figure handle.
SEE ALSO: getpts, getline, getrect, zoom
Credit to John D'Errico for writing the useful 'interparc' function included within this function.
[X0,Y0] = getline_zoom(Im1,'plot');
[X1,Y1] = getline_zoom(Im1,'plot','interp',0.5,'method','spline');
[X2,Y2] = getline_zoom(Im1,'plot','interp',0.5,'method','spline','loop_path');
Really helped me also. I would suggest to change line #141 for : Im_temp=Im1(max([yb,1]):min([yt,size(Im1,1)]),max([xl,1]):min([xr,size(Im1,2)]));
That would avoid the program to crash when the temporary image reaches one of the original image boundary.
Thank you for posting this code -- it was very helpful for my project!
I just have a couple of suggestions to improve the accuracy of your code:
1. In line 198 (calculating the length of the curve), summing across the raw selected points will always underestimate the length of the 'true' curve, since you are subsampling the curve. Instead, I would use interparc.m to generate a smoother curve (roughly at the resolution that is desired by the user), and estimate the length of the curve using those interpolated points. This takes more computation but on certain applications the extra accuracy might be important.
2. I would change line 205 (using interparc.m) to the following:
O1 = interparc(round(Ltot/res) + 1,X,Y,mt1);
to account for the two edge points.
Otherwise, great work!
added smoothing, looping, fixed resolution, and other improvements.
Inspired by: interparc