清空筛选条件
清空筛选条件

problem with dealing images with formulas

2 次查看（过去 30 天）
joanna
joanna 2012-11-19
i want to display the results in an axes. i1 i2 i3 i4 i5 i6 are images displayed on individual axes and i want to display a final result which will only display 1 image "diff".
theta = 0.5 * arctan ((i5 - i3) / (i4 - i6));
delta = arctan( 2 * (i5 - i4) / (i1 - i2)*(sin2(theta) - cos2(theta)));
N = delta / (2*3.14);
fsigma = ((8 * 200) / (3.14 * 22.5 * 7));
diff = (N * fsigma) / 22.5;
axes(handles.axes14)
imagesc(abs(diff)),title('Stress distribution in the bone'),xlabel('Pixel'),ylabel('Pixel'),colorbar;
but it gives me these errors:
??? Error using ==> mrdivide
Linear algebra is not supported for integer data types.
Error in ==> VPPP>calButton_Callback at 208
theta = 0.5 * arctan ((i5 - i3) / (i4 - i6));
Error in ==> gui_mainfcn at 96
feval(varargin{:});
Error in ==> VPPP at 42
gui_mainfcn(gui_State, varargin{:});
Error in ==>
@(hObject,eventdata)VPPP('calButton_Callback',hObject,eventdata,guidata(hObject))
??? Error while evaluating uicontrol Callback
I'm not sure what has gone wrong.

请先登录，再回答此问题。

回答（1 个）

Walter Roberson
Walter Roberson 2012-11-19
Are you sure you want matrix division and not element-by-element division which is the ./ operator ?
Anyhow, you should probably be using double() of the images when you are doing numeric calculations.
Note: it is not recommended that you call any variable "diff" as that will interfere with the use of the diff() function.
  2 个评论
Walter Roberson
Walter Roberson 2012-11-19
For example, change
theta = 0.5 * arctan ((i5 - i3) / (i4 - i6));
to
theta = 0.5 * arctan ((i5 - i3) ./ (i4 - i6));

请先登录，再进行评论。

类别

在 Help Center 和 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