1  %'ima_ratio': take the ratio of two input images with same size


2  %


3  % the two fields are subtstracted when of the same nature (scalar or


4  % vector), if the coordinates do not coincide, the second field is


5  % interpolated on the cooridintes of the first one


6  %


7  % when scalar and vectors are combined, the fields are just merged in a single matlab structure for common visualisation


8  %


9  % function SubData=sub_field(Field,XmlData,Field_1)


10  %


11  % OUPUT:


12  % SubData: structure representing the resulting field


13  %


14  % INPUT:


15  % Field: matlab structure representing the first field


16  % XmlData: input calibration parameter, not used


17  % Field_1:matlab structure representing the second field


18 


19  %=======================================================================


20  % Copyright 20082015, LEGI UMR 5519 / CNRS UJF GINP, Grenoble, France


21  % http://www.legi.grenobleinp.fr


22  % Joel.Sommeria  Joel.Sommeria (A) legi.cnrs.fr


23  %


24  % This file is part of the toolbox UVMAT.


25  %


26  % UVMAT is free software; you can redistribute it and/or modify


27  % it under the terms of the GNU General Public License as published


28  % by the Free Software Foundation; either version 2 of the license,


29  % or (at your option) any later version.


30  %


31  % UVMAT is distributed in the hope that it will be useful,


32  % but WITHOUT ANY WARRANTY; without even the implied warranty of


33  % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the


34  % GNU General Public License (see LICENSE.txt) for more details.


35  %=======================================================================


36 


37  function DataOut=ima_ratio(DataIn,XmlData,DataIn_1)


38 


39  %% option to introduce input parameters when function is selected, skipped here


40  DataOut=DataIn;


41  if ~isstruct(DataOut)


42  % if ~isfield(DataIn,'A')


43  % msgbox_uvmat('ERROR','ima_ratio requires two images of the same size as input')


44  return


45  end


46 


47 


48  %% actual calculation


49  if exist('DataIn_1','var')


50  DataIn_1.A(DataIn_1.A==0)=1;


51  DataOut.A=double(DataIn.A)./double(DataIn_1.A);


52  end

