hdrvdp_joint_rod_cone_sens_diff.m 1.1 KB

1234567891011121314151617181920212223
  1. function Sd = hdrvdp_joint_rod_cone_sens_diff( la, metric_par )
  2. % Copyright (c) 2011, Rafal Mantiuk <mantiuk@gmail.com>
  3. % Permission to use, copy, modify, and/or distribute this software for any
  4. % purpose with or without fee is hereby granted, provided that the above
  5. % copyright notice and this permission notice appear in all copies.
  6. %
  7. % THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
  8. % WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
  9. % MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
  10. % ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
  11. % WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
  12. % ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
  13. % OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
  14. p1 = metric_par.csf_sa(1);
  15. p2 = metric_par.csf_sa(2); % in the paper - p6
  16. p3 = metric_par.csf_sa(3); % in the paper - p7
  17. p4 = metric_par.csf_sa(4); % in the paper - p8
  18. Sd = (p1*p2*p3*p4*(p2./la).^(p3 - 1))./(la.^2.*((p2./la).^p3 + 1).^(p4 + 1));
  19. end