def pollyxt_cge_display_saturation(tmpFile, saveFolder):
"""
Description
-----------
Display the housekeeping data from laserlogbook file.
Parameters
----------
tmpFile: str
the .mat file which stores the housekeeping data.
saveFolder: str
Usage
-----
pollyxt_cge_display_saturation(tmpFile)
History
-------
2019-01-10. First edition by Zhenping
"""
if not os.path.exists(tmpFile):
print('{filename} does not exists.'.format(filename=tmpFile))
return
# read matlab .mat data
try:
mat = spio.loadmat(tmpFile, struct_as_record=True)
figDPI = mat['figDPI'][0][0]
mTime = mat['time'][0][:]
height = mat['height'][0][:]
SAT_FR_355 = mat['SAT_FR_355'][:]
SAT_FR_532 = mat['SAT_FR_532'][:]
SAT_FR_1064 = mat['SAT_FR_1064'][:]
yLim_FR_RCS = mat['yLim_FR_RCS'][:][0]
pollyVersion = mat['campaignInfo']['name'][0][0][0]
location = mat['campaignInfo']['location'][0][0][0]
version = mat['processInfo']['programVersion'][0][0][0]
fontname = mat['processInfo']['fontname'][0][0][0]
dataFilename = mat['taskInfo']['dataFilename'][0][0][0]
xtick = mat['xtick'][0][:]
xticklabel = mat['xtickstr']
imgFormat = mat['imgFormat'][:][0]
except Exception as e:
print(e)
print('Failed reading %s' % (tmpFile))
return
# set the default font
matplotlib.rcParams['font.sans-serif'] = fontname
matplotlib.rcParams['font.family'] = "sans-serif"
# meshgrid
Time, Height = np.meshgrid(mTime, height)
# load colormap
dirname = os.path.dirname(os.path.dirname(os.path.abspath(__file__)))
sys.path.append(dirname)
try:
from python_colormap import signal_status_colormap
except Exception as e:
raise ImportError('python_colormap module is necessary.')
# display status of 355 FR
fig = plt.figure(figsize=[10, 5])
ax = fig.add_axes([0.11, 0.15, 0.74, 0.75])
pcmesh = ax.pcolormesh(Time,
Height,
SAT_FR_355,
vmin=-0.5,
vmax=2.5,
cmap=signal_status_colormap(),
rasterized=True)
ax.set_xlabel('UTC', fontsize=15)
ax.set_ylabel('Height (m)', fontsize=15)
ax.set_ylim(yLim_FR_RCS.tolist())
ax.yaxis.set_major_locator(MultipleLocator(2500))
ax.yaxis.set_minor_locator(MultipleLocator(500))
ax.set_xticks(xtick.tolist())
ax.set_xticklabels(celltolist(xticklabel))
ax.tick_params(axis='both',
which='major',
labelsize=15,
right=True,
top=True,
width=2,
length=5)
ax.tick_params(axis='both',
which='minor',
width=1.5,
length=3.5,
right=True,
top=True)
ax.set_title('Signal Status at ' +
'{wave}nm Far-Range from {instrument} at {location}'.format(
wave=355, instrument=pollyVersion, location=location),
fontsize=15)
cb_ax = fig.add_axes([0.865, 0.15, 0.02, 0.75])
cbar = fig.colorbar(pcmesh,
cax=cb_ax,
ticks=[0, 1, 2],
orientation='vertical')
cbar.ax.tick_params(direction='in', pad=5)
cbar.ax.set_title('', fontsize=9)
cbar.ax.set_yticklabels(['Good Signal', 'Saturated', 'Low SNR'])
cbar.ax.tick_params(axis='both',
which='major',
labelsize=12,
right=True,
top=True,
width=2,
length=5)
cbar.ax.tick_params(axis='both',
which='minor',
width=1.5,
length=3.5,
right=True,
top=True)
fig.text(0.05,
0.04,
datenum_to_datetime(mTime[0]).strftime("%Y-%m-%d"),
fontsize=15)
fig.text(0.8,
0.04,
'Version: {version}'.format(version=version),
fontsize=14)
fig.savefig(os.path.join(
saveFolder, '{dataFilename}_SAT_FR_355.{imgFmt}'.format(
dataFilename=rmext(dataFilename), imgFmt=imgFormat)),
dpi=figDPI)
plt.close()
# display status of 532 FR
fig = plt.figure(figsize=[10, 5])
ax = fig.add_axes([0.11, 0.15, 0.74, 0.75])
pcmesh = ax.pcolormesh(Time,
Height,
SAT_FR_532,
vmin=-0.5,
vmax=2.5,
cmap=signal_status_colormap(),
rasterized=True)
ax.set_xlabel('UTC', fontsize=15)
ax.set_ylabel('Height (m)', fontsize=15)
ax.set_ylim(yLim_FR_RCS.tolist())
ax.yaxis.set_major_locator(MultipleLocator(2500))
ax.yaxis.set_minor_locator(MultipleLocator(500))
ax.set_xticks(xtick.tolist())
ax.set_xticklabels(celltolist(xticklabel))
ax.tick_params(axis='both',
which='major',
labelsize=15,
right=True,
top=True,
width=2,
length=5)
ax.tick_params(axis='both',
which='minor',
width=1.5,
length=3.5,
right=True,
top=True)
ax.set_title('Signal Status at ' +
'{wave}nm Far-Range from {instrument} at {location}'.format(
wave=532, instrument=pollyVersion, location=location),
fontsize=15)
cb_ax = fig.add_axes([0.865, 0.15, 0.02, 0.75])
cbar = fig.colorbar(pcmesh,
cax=cb_ax,
ticks=[0, 1, 2],
orientation='vertical')
cbar.ax.tick_params(direction='in', pad=5)
cbar.ax.set_title('', fontsize=9)
cbar.ax.set_yticklabels(['Good Signal', 'Saturated', 'Low SNR'])
cbar.ax.tick_params(axis='both',
which='major',
labelsize=12,
right=True,
top=True,
width=2,
length=5)
cbar.ax.tick_params(axis='both',
which='minor',
width=1.5,
length=3.5,
right=True,
top=True)
fig.text(0.05,
0.04,
datenum_to_datetime(mTime[0]).strftime("%Y-%m-%d"),
fontsize=15)
fig.text(0.8,
0.04,
'Version: {version}'.format(version=version),
fontsize=14)
fig.savefig(os.path.join(
saveFolder, '{dataFilename}_SAT_FR_532.{imgFmt}'.format(
dataFilename=rmext(dataFilename), imgFmt=imgFormat)),
dpi=figDPI)
plt.close()
# display status of 1064 FR
fig = plt.figure(figsize=[10, 5])
ax = fig.add_axes([0.11, 0.15, 0.74, 0.75])
pcmesh = ax.pcolormesh(Time,
Height,
SAT_FR_1064,
vmin=-0.5,
vmax=2.5,
cmap=signal_status_colormap(),
rasterized=True)
ax.set_xlabel('UTC', fontsize=15)
ax.set_ylabel('Height (m)', fontsize=15)
ax.set_ylim(yLim_FR_RCS.tolist())
ax.yaxis.set_major_locator(MultipleLocator(2500))
ax.yaxis.set_minor_locator(MultipleLocator(500))
ax.set_xticks(xtick.tolist())
ax.set_xticklabels(celltolist(xticklabel))
ax.tick_params(axis='both',
which='major',
labelsize=15,
right=True,
top=True,
width=2,
length=5)
ax.tick_params(axis='both',
which='minor',
width=1.5,
length=3.5,
right=True,
top=True)
ax.set_title('Signal Status at ' +
'{wave}nm Far-Range from {instrument} at {location}'.format(
wave=1064, instrument=pollyVersion, location=location),
fontsize=15)
cb_ax = fig.add_axes([0.865, 0.15, 0.02, 0.75])
cbar = fig.colorbar(pcmesh,
cax=cb_ax,
ticks=[0, 1, 2],
orientation='vertical')
cbar.ax.tick_params(direction='in', pad=5)
cbar.ax.set_title('', fontsize=9)
cbar.ax.set_yticklabels(['Good Signal', 'Saturated', 'Low SNR'])
cbar.ax.tick_params(axis='both',
which='major',
labelsize=12,
right=True,
top=True,
width=2,
length=5)
cbar.ax.tick_params(axis='both',
which='minor',
width=1.5,
length=3.5,
right=True,
top=True)
fig.text(0.05,
0.04,
datenum_to_datetime(mTime[0]).strftime("%Y-%m-%d"),
fontsize=15)
fig.text(0.8,
0.04,
'Version: {version}'.format(version=version),
fontsize=14)
fig.savefig(os.path.join(
saveFolder, '{dataFilename}_SAT_FR_1064.{imgFmt}'.format(
dataFilename=rmext(dataFilename), imgFmt=imgFormat)),
dpi=figDPI)
plt.close()
def pollyDisplaySigStatus387NR(tmpFile, saveFolder):
"""
Description
-----------
Display the housekeeping data from laserlogbook file.
Parameters
----------
tmpFile: str
the .mat file which stores the housekeeping data.
saveFolder: str
Usage
-----
pollyDisplaySigStatus387NR(tmpFile)
History
-------
2019-01-10. First edition by Zhenping
"""
if not os.path.exists(tmpFile):
print('{filename} does not exists.'.format(filename=tmpFile))
return
# read matlab .mat data
try:
mat = spio.loadmat(tmpFile, struct_as_record=True)
figDPI = mat['figDPI'][0][0]
flagWatermarkOn = mat['flagWatermarkOn'][0][0]
if mat['partnerLabel'].size:
partnerLabel = mat['partnerLabel'][0]
else:
partnerLabel = ''
mTime = mat['time'][0][:]
height = mat['height'][0][:]
SAT_NR_387 = mat['SAT_NR_387'][:]
yLim_NR_RCS = mat['yLim_NR_RCS'][:][0]
pollyVersion = mat['CampaignConfig']['name'][0][0][0]
location = mat['CampaignConfig']['location'][0][0][0]
version = mat['PicassoConfig']['PicassoVersion'][0][0][0]
fontname = mat['PicassoConfig']['fontname'][0][0][0]
dataFilename = mat['PollyDataInfo']['pollyDataFile'][0][0][0]
xtick = mat['xtick'][0][:]
xticklabel = mat['xtickstr']
imgFormat = mat['imgFormat'][:][0]
except Exception as e:
print(e)
print('Failed reading %s' % (tmpFile))
return
# set the default font
matplotlib.rcParams['font.sans-serif'] = fontname
matplotlib.rcParams['font.family'] = "sans-serif"
# meshgrid
Time, Height = np.meshgrid(mTime, height)
# load colormap
dirname = os.path.dirname(os.path.dirname(os.path.abspath(__file__)))
sys.path.append(dirname)
try:
from python_colormap import signal_status_colormap
except Exception as e:
raise ImportError('python_colormap module is necessary.')
# display status of 387 NR
fig = plt.figure(figsize=[10, 5])
ax = fig.add_axes([0.11, 0.15, 0.74, 0.75])
pcmesh = ax.pcolormesh(Time,
Height,
SAT_NR_387,
vmin=-0.5,
vmax=2.5,
cmap=signal_status_colormap(),
rasterized=True,
shading='nearest')
ax.set_xlabel('UTC', fontsize=15)
ax.set_ylabel('Height (m)', fontsize=15)
ax.set_ylim(yLim_NR_RCS.tolist())
ax.yaxis.set_major_locator(MultipleLocator(1000))
ax.yaxis.set_minor_locator(MultipleLocator(200))
ax.set_xticks(xtick.tolist())
ax.set_xticklabels(celltolist(xticklabel))
ax.tick_params(axis='both',
which='major',
labelsize=15,
right=True,
top=True,
width=2,
length=5)
ax.tick_params(axis='both',
which='minor',
width=1.5,
length=3.5,
right=True,
top=True)
ax.set_title('Signal Status at ' +
'{wave}nm Near-Range from {instrument} at {location}'.format(
wave=387, instrument=pollyVersion, location=location),
fontsize=15)
cb_ax = fig.add_axes([0.865, 0.15, 0.02, 0.75])
cbar = fig.colorbar(pcmesh,
cax=cb_ax,
ticks=[0, 1, 2],
orientation='vertical')
cbar.ax.tick_params(direction='in', pad=5)
cbar.ax.set_title('', fontsize=9)
cbar.ax.set_yticklabels(['Good Signal', 'Saturated', 'Low SNR'])
cbar.ax.tick_params(axis='both',
which='major',
labelsize=12,
right=True,
top=True,
width=2,
length=5)
cbar.ax.tick_params(axis='both',
which='minor',
width=1.5,
length=3.5,
right=True,
top=True)
# add watermark
if flagWatermarkOn:
rootDir = os.path.dirname(
os.path.dirname(os.path.dirname(os.path.abspath(__file__))))
im_license = matplotlib.image.imread(
os.path.join(rootDir, 'img', 'by-sa.png'))
newax_license = fig.add_axes([0.58, 0.006, 0.14, 0.07], zorder=10)
newax_license.imshow(im_license, alpha=0.8, aspect='equal')
newax_license.axis('off')
fig.text(0.72,
0.003,
'Preliminary\nResults.',
fontweight='bold',
fontsize=12,
color='red',
ha='left',
va='bottom',
alpha=0.8,
zorder=10)
fig.text(0.84,
0.003,
u"\u00A9 {1} {0}.\nCC BY SA 4.0 License.".format(
datetime.now().strftime('%Y'), partnerLabel),
fontweight='bold',
fontsize=7,
color='black',
ha='left',
va='bottom',
alpha=1,
zorder=10)
fig.text(0.05,
0.04,
datenum_to_datetime(mTime[0]).strftime("%Y-%m-%d"),
fontsize=15)
fig.text(0.2,
0.04,
'Version: {version}'.format(version=version),
fontsize=14)
fig.savefig(os.path.join(
saveFolder,
'{dataFilename}_SAT_NR_387.{imgFmt}'.format(dataFilename=rmext(
os.path.basename(dataFilename)),
imgFmt=imgFormat)),
dpi=figDPI)
plt.close()