How to use arrayfun to create an array whose value are relating to index?
显示 更早的评论
I'm now doing a project of image processing, but I have some troubles when using the function arrayfun(). Let me show the core part of my code first. I create an array named resizedImage.
resizedImage = zeros(dr, dc, 'uint8');
Then, I attempt to change the value of each entry of this array, so I use for-loop to achieve this.
for i = 1 : dr
for j = 1 : dc
resizedImage(i, j) = pixel_replication(img, i * dist_r, j * dist_c);
end
end
function val = pixel_replication(img, x, y)
% do something
end
However, for-loop is very time-consuming, so I try to use arrayfun() to make my code more efficient. Anyone can help me?
3 个评论
Dyuman Joshi
2023-9-20
arrayfun is not going to be faster than the for loop, at best it might match it.
Your best bet would be to see if the function pixel_replication() can be optimized.
Can you attach your code? Use the paperclip button to attach.
Duen Chian
2023-9-20
dpb
2023-9-20
for i = 1 : dr
for j = 1 : dc
resizedImage(i, j) = pixel_replication(img, i * dist_r, j * dist_c);
The order of the loops above is processing the array in row major, not column major order; rearranging those as
for j = 1 : dc
for i = 1 : dr
resizedImage(i, j) = pixel_replication(img, i * dist_r, j * dist_c);
could help some, but as @Dyuman Joshi notes, the time consumed will be in the function so whatever can be done there would be the place to look.
采纳的回答
function val = pixel_replication(img, x, y)
val = img(nearest(x), nearest(y));
end
function u = nearest(v)
l = floor(v);
r = ceil(v);
if v - l < r - v
u = l;
else
u = r;
end
end
nearest(v) returns the nearest integer to the input v, and if v is exactly halfway between the two nearest integers, then nearest(v) returns ceil(v). For a non-negative scalar number v, this is the same behavior as the built-in round function, but round also accepts array inputs, so if you use round instead of nearest, then you can easily vectorize your entire code:
function resizedImage = resizeImage_replication(originalImage, scalingFactor)
% load image
img = imread(originalImage);
% get number of rows and columns
[sr, sc] = size(img);
% adjust row and column
dr = ceil(sr * scalingFactor);
dc = ceil(sc * scalingFactor);
% % create an empty array - no longer necessary (and now resizedImage will
% % be the same class as originalImage, be it uint8 or whatever else)
% resizedImage = zeros(dr, dc, 'uint8');
% get dist
dist_r = (sr - 1) / (dr - 1);
dist_c = (sc - 1) / (dc - 1);
% pixel replication: resizedImage = imresize(img, scalingFactor, "nearest");
% for i = 1 : dr
% x = 1 + (i-1) * dist_r;
% for j = 1 : dc
% y = 1 + (j-1) * dist_c;
% resizedImage(i, j) = pixel_replication(img, x, y);
% end
% end
% vectors x and y:
x = 1 + (0:dr-1) * dist_r;
y = 1 + (0:dc-1) * dist_c;
% use round instead of nearest:
resizedImage = img(round(x),round(y));
end
2 个评论
Duen Chian
2023-9-21
Voss
2023-9-21
You're welcome!
The two 1-D arrays are indexing into the 2-D array img. round(x) is the row index and round(y) is the column index.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Image Arithmetic 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
