# 'improfile' is a missing function in octave, what other options could I use to make my loop works?

2 Ansichten (letzte 30 Tage)
BA am 28 Jul. 2022
Kommentiert: Jan am 29 Jul. 2022
imshow(matrix(:,:,1))
%identify axes
[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
end
scog = (cog - min(cog)) / (max(cog) - min(cog));
##### 4 Kommentare2 ältere Kommentare anzeigen2 ältere Kommentare ausblenden
Jan am 28 Jul. 2022
@Bilal Alsharif: Do you see the set of icons above the field for typing the text in the forum? They help you to apply a nice formatting. I've formatted the code for you in the question already. Please use these formatting tools to distinguish text, Matlab code and programs.
Jan am 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.

Melden Sie sich an, um zu kommentieren.

### Antworten (1)

Jan am 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));
##### 1 Kommentar-1 ältere Kommentare anzeigen-1 ältere Kommentare ausblenden
Walter Roberson am 28 Jul. 2022
improfile figures out whether the row difference is greater than the column difference, and does the interp2 along the longer dimension.

Melden Sie sich an, um zu kommentieren.

### Kategorien

Mehr zu Octave finden Sie in Help Center und File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by