memls_legacy

Table of Contents

In [1]:
import os

# general import for smrt
from smrt import make_snowpack, make_model, sensor

# import for memls
from smrt.utils import memls_legacy


import matplotlib.pyplot as plt
%matplotlib inline

Prepare the snowpack and the sensor

In [2]:
# prepare snowpack
pc=0.2e-3
snowpack = make_snowpack([10], "exponential",
                        density=[300], temperature=[265], corr_length=pc)

# create the sensor
theta = range(10, 80, 5)
radiometer = sensor.passive(37e9, theta)

Run SMRT with IBA + DORT

In [3]:
# create the EM Model
m = make_model("iba", "dort")

# run the model
sresult = m.run(radiometer, snowpack)

Run MEMLS matlab code

In [4]:
# set the path where memls matlab code is
mresult = memls_legacy.run(radiometer, snowpack)
In [5]:
# outputs
# outputs
plt.plot(theta, sresult.TbV(), 'r-', label='SMRT V')
plt.plot(theta, sresult.TbH(), 'r--', label='SMRT H')
plt.plot(theta, mresult.TbV(), 'b-', label='MEMLS V')
plt.plot(theta, mresult.TbH(), 'b--', label='MEMLS H')
plt.legend(loc='best')
Out[5]:
<matplotlib.legend.Legend at 0x7ff5c40fba58>