Note: Suggested to avoid "i" as a vaiable name, as it may conflict with inbuilt "i"
>> i
ans =
0.0000 + 1.0000i
Lets consider image_data as input data, in 2D case, you can directly used the find function, as
[r,c]=find(image_data==max(image_data));
In 3D or higher multidimentional case (Add more variable in LHS, as per dimension of the array), you may use ind2sub to get the desired result
idx=find(image_data==max(image_data(:)));
[rows,colm,sl_num]=ind2sub(size(image_data),idx)
It gives the indices of rows, columns and slice numbers of all maximum pixels elements.
