In [1]:
import accuread as ar
import matplotlib.pyplot as plt
%matplotlib inline

In [2]:
moddir = '../tests/testdata/'
d = ar.ReadART('radiop', 
               basefolder=moddir, # folder where the Output-folder is located
               radiance=True, 
               iops=True, # read iop file
              )

Plot radiance

At the moment only a rather simple polar contour plot is available.


In [3]:
d.plot_rad_contour(wl=341,z=1e5)


Out[3]:
<matplotlib.axes._subplots.PolarAxesSubplot at 0x7f297352bc88>

If the wavelength is specified with the wl keyword, as above, the closest wavelength is chosen. Similarly, if the vertical position is defined with the z keyword, the detector closest to that position is chosen.

One can alternatively specify the index for wavelength and/or height:


In [4]:
d.plot_rad_contour(wl_index=5,z_index=1)


Out[4]:
<matplotlib.axes._subplots.PolarAxesSubplot at 0x7f2943fbe278>

IOPs

A single figure showing absorption and scattering coefficients, asymmetry factor and single scattering albedo, is shown. For materials with very strong forward scattering, AccuRT applies delta-fit scaling. In other words, the phase function is made a little less forward peaked, and the scattering coefficient is reduced. "Scaled" variables are after the delta-fit scaling.


In [5]:
d.plot_iops()


Out[5]:
((<matplotlib.axes._subplots.AxesSubplot at 0x7f2943e811d0>,
  <matplotlib.axes._subplots.AxesSubplot at 0x7f2943e6ceb8>),
 (<matplotlib.axes._subplots.AxesSubplot at 0x7f2943db5da0>,
  <matplotlib.axes._subplots.AxesSubplot at 0x7f2943deb828>))

As with radiance, one can choose wavelength and layer. But unlike the radiance plot, the z key corresponds to layer depths instead of detector depths.


In [6]:
d.plot_iops(wl=520,z=1e5+0.1)


Out[6]:
((<matplotlib.axes._subplots.AxesSubplot at 0x7f2944139278>,
  <matplotlib.axes._subplots.AxesSubplot at 0x7f2944133198>),
 (<matplotlib.axes._subplots.AxesSubplot at 0x7f29440b2080>,
  <matplotlib.axes._subplots.AxesSubplot at 0x7f2943c30da0>))

In [ ]: