'improfile' is a missing function in octave, what other options could I use to make my loop works?
2 views (last 30 days)
Show older comments
BA on 28 Jul 2022
Commented: Jan on 29 Jul 2022
[x, y] = ginput(2);
% preallocate matrices
cog = zeros(size(matrix,3),1);
% loop start
for i = 1:size(matrix,3)
I = matrix(:,:,i);
%what to do to make this code works in octave
test = improfile(I,[x(1) x(2)],[y(1) y(2)]);
cog(i) = sum((1:length(test)).*test')/sum(test);
% loop end
scog = (cog - min(cog)) / (max(cog) - min(cog));
Jan on 29 Jul 2022
@Image Analyst: I've re-opened the question. Although Octave is off-topic, the answer works for Matlab users without the Image Processing Toolbox.
Jan on 28 Jul 2022
The command interp2(I, 2) does not consider the selected line. Better:
n = ceil(norm([diff(x), diff(y)])); % A rough estimation of number of points
test = interp2(I, 2, linspace(x(1), x(2), n), linspace(y(1), y(2), n));
Walter Roberson on 28 Jul 2022
improfile figures out whether the row difference is greater than the column difference, and does the interp2 along the longer dimension.
Find more on Octave in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!