How to return answer in 8 columns of 2D matrices through division

Question

Yun Inn 2012-7-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/44285-how-to-return-answer-in-8-columns-of-2d-matrices-through-division

I need to create 8 columns of matrices based on division of 2 separate sets of matrices. My code below:

PhysicalThickness=[0,5,10,15,20,25,30,35];
for i=1:length(PhysicalThickness)
h=num2str(i-1);
fileextension='.dcm';
filenameEN=strcat(XIOM0,h,fileextension);
filenameEX=strcat(XIO,h,fileextension);
[OPEN]=dicomread(filenameEN);
[data]=dicomread(filenameEX);
rTPS=data./OPEN;
end
%______________________________________________________________________%
PhysicalThickness=[0,5,10,15,20,25,30,35];
for j=1:length(PhysicalThickness)
h=num2str(j-1);
fileextension='.dcm';
EPIDfilenameEN=strcat(EPIDM0,h,fileextension);
EPIDfilenameEX=strcat(EPID,h,fileextension);
[OPEN]=dicomread(EPIDfilenameEN);
[data]=dicomread(EPIDfilenameEX);
rEPID=data./OPEN;
end
%____________________________________________________________________%
ratio=rTPS./rEPID;

rTPS and rEPID, both returned answer in 8 columns of PhysicalThickness (which is what I want). However, when I divide the product of the two (rTPS./rEPID), I get only 1 answer (belonging to the last PhysicalThickness).

How do I get the answer for all i/j (PhysicalThickness) in 8 columns instead of just one?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Jan 2012-7-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/44285-how-to-return-answer-in-8-columns-of-2d-matrices-through-division#answer_54294

编辑：Jan 2012-7-23

在 MATLAB Online 中打开

The problem appears, because you overwrite rTPS and rEPID in each iteration. The solution is to read both input and calculate the results inside the loop:

PhysicalThickness=[0,5,10,15,20,25,30,35];
ratio = cell(1, length(PhysicalThickness));
for i=1:length(PhysicalThickness)
  h=num2str(i-1);
  fileextension='.dcm';
  filenameEN=strcat(XIOM0,h,fileextension);
  filenameEX=strcat(XIO,h,fileextension);  
  OPEN=dicomread(filenameEN);
  data=dicomread(filenameEX);
  rTPS=data./OPEN;
  %
  EPIDfilenameEN=strcat(EPIDM0,h,fileextension);
  EPIDfilenameEX=strcat(EPID,h,fileextension);
  OPEN=dicomread(EPIDfilenameEN);
  data=dicomread(EPIDfilenameEX);
  rEPID=data./OPEN;
  %
  ratio{i} = rTPS./rEPID;
end