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

5 visualizzazioni (ultimi 30 giorni)
BA il 28 Lug 2022
Commentato: Jan il 29 Lug 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 CommentiMostra 2 commenti meno recentiNascondi 2 commenti meno recenti
Jan il 28 Lug 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 il 29 Lug 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.

Accedi per commentare.

### Risposte (1)

Jan il 28 Lug 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 CommentoMostra -1 commenti meno recentiNascondi -1 commenti meno recenti
Walter Roberson il 28 Lug 2022
improfile figures out whether the row difference is greater than the column difference, and does the interp2 along the longer dimension.

Accedi per commentare.

### Categorie

Scopri di più su Octave in Help Center e File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by