1
function [val_UWattsPerCm2,limit_UWattsPerCm2] = ISO2007MPEComputeType1ContinuousCornealUVUnweightedValue(S,radiance_WattsPerSrM2,stimulusDurationSecs,stimulusAreaDegrees2)
2
% [val_UWattsPerCm2,limit_UWattsPerCm2] = ISO2007MPEComputeType1ContinuousCornealUVUnweightedValue(S,radiance_WattsPerSrM2,stimulusDurationSecs,stimulusAreaDegrees2)
4
% Compute the unweighted UV radiation for Type 1 instruments as given on page 7, Table 2,
7
% Input spectrum is radiance in units of Watts/[sr-m2-wlinterval].
9
% Also return the exposure limit for this quantity.
11
% See page 6 for a definition of a Type 1 instrument. As far as I can tell, the key
12
% criterion is that it doesn't put out more light that exceeds the Type 1 limits.
14
% If the exposure time is longer than 2 hours the specified limits should be reduced by
15
% 1/exposureDuration in hours. This routine implements that adjustment for its returned
16
% limit value. It does not implement a further reduction of of the limit (by a factor of 2)
17
% for microscopes and endoilluminators.
19
% ****************************************************************************
20
% IMPORTANT: Before using the ISO2007MPE routines, please see the notes on usage
21
% and responsibility in PsychISO2007MPE/Contents.m (type "help PsychISO2007MPE"
22
% at the Matlab prompt.
23
% ****************************************************************************
25
% 6/25/13 dhb Wrote it.
27
%% Specify the limit (from table)
29
% Limit in table is in mWatts, we convert to uWatts
30
% for consistency with the weighted routine.
31
exposureDurationHours = stimulusDurationSecs/3600;
32
if (exposureDurationHours <= 2)
33
limit_UWattsPerCm2 = 1*(10^3);
35
limit_UWattsPerCm2 = 1*(10^3)/(exposureDurationHours/2);
38
%% Convert radiance to corneal irradiance
39
cornealIrradiance_WattsPerM2 = RadianceAndDegrees2ToCornIrradiance(radiance_WattsPerSrM2,stimulusAreaDegrees2);
40
cornealIrradiance_UWattsPerM2 = (10^6)*cornealIrradiance_WattsPerM2;
41
cornealIrradiance_UWattsPerCm2 = (10^-4)*cornealIrradiance_UWattsPerM2;
43
%% Get sum betwen 360 and 400.
45
index = find(wls >= 360 & wls <= 400);
47
error('Should not call this routine with no spectral sampling between 360 and 400');
49
val_UWattsPerCm2 = sum(cornealIrradiance_UWattsPerCm2(index));