Atmospheric Lidar Package: comparison licel.py

-:870fc8f65eeb
+:28d7b0974fe6
 import numpy as np
 import datetime
-from generic import BaseLidarMeasurement, Lidar_channel
+from generic import BaseLidarMeasurement, LidarChannel
 import musa_netcdf_parameters
 import musa_2009_netcdf_parameters
 licel_file_header_format = ['Filename',
 'Site StartDate StartTime EndDate EndTime Altitude Longtitude Latitude ZenithAngle',
 'LS1 Rate1 LS2 Rate2 DataSets', ]
 licel_file_channel_format = 'Active AnalogPhoton LaserUsed DataPoints 1 HV BinW Wavelength d1 d2 d3 d4 ADCbits NShots Discriminator ID'
 class LicelFile:
 def __init__(self, filename):
 self.start_time = None
 self.stop_time = None
 self.import_file(filename)
 self.calculate_physical()
 a = np.fromfile(f, 'b', 1)
 b = np.fromfile(f, 'b', 1)
 if (a[0] != 13) | (b[0] != 10):
 print "Warning: No end of line found after record. File could be corrupt"
-channel = LicelFileChannel(current_channel_info, raw_data)
+channel = LicelFileChannel(current_channel_info, raw_data, self.duration())
 channel_name = channel.channel_name
 if channel_name in channels.keys():
 # If the analog/photon naming scheme is not enough, find a new one!
 raise IOError('Trying to import two channels with the same name')
 channels[channel_name] = channel
 f.close()
 self.raw_info = raw_info
 self.channels = channels
+def duration(self):
+""" Return the duration of the file. """
+dt = self.stop_time - self.start_time
+return dt.seconds
 class LicelFileChannel:
-def __init__(self, raw_info = None, raw_data = None):
+def __init__(self, raw_info = None, raw_data = None, duration = None):
 self.raw_info = raw_info
 self.raw_data = raw_data
+self.duration = duration
 @property
 def wavelength(self):
 if self.raw_info is not None:
 wave_str = self.raw_info['Wavelength']
 wavelength = wave_str.split('.')[0]
 def calculate_physical(self):
 data = self.raw_data
 number_of_shots = float(self.raw_info['NShots'])
-norm = data/number_of_shots
+norm = data / number_of_shots
 dz = float(self.raw_info['BinW'])
 if self.raw_info['AnalogPhoton']=='0':
 # If the channel is in analog mode
 ADCb = int(self.raw_info['ADCbits'])
 # c = 300 # The result will be in MHZ
 # SR  = c/(2*dz) # To account for pulse folding
 # channel_data = norm*SR
 # CHANGE:
 # For the SCC the data are needed in photons
-channel_data = norm*number_of_shots
+channel_data = norm *number_of_shots
 #print res,c,cdata,norm
 #Calculate Z
 number_of_bins  = int(self.raw_info['DataPoints'])
 self.z = np.array([dz*bin_number + dz/2.0 for bin_number in range(number_of_bins)])
 class LicelLidarMeasurement(BaseLidarMeasurement):
 '''
 '''
 extra_netcdf_parameters = musa_netcdf_parameters
-raw_info = {} # Keep the raw info from the files
+raw_info = {}   # Keep the raw info from the files
+durations = {}  # Keep the duration of the files
 def import_file(self, filename):
 if filename in self.files:
 print "File has been imported already:" + filename
 else:
 current_file = LicelFile(filename)
 self.raw_info[current_file.filename] = current_file.raw_info
+self.durations[current_file.filename] = current_file.duration()
 for channel_name, channel in current_file.channels.items():
 if channel_name not in self.channels:
-self.channels[channel_name] = Licel_channel(channel)
+self.channels[channel_name] = LicelChannel(channel)
 self.channels[channel_name].data[current_file.start_time] = channel.data
 self.files.append(current_file.filename)
 def append(self, other):
 def _get_duration(self, raw_start_in_seconds):
 ''' Return the duration for a given time scale. If only a single
 file is imported, then this cannot be guessed from the time difference
 and the raw_info of the file are checked.
 '''
 if len(raw_start_in_seconds) == 1: # If only one file imported
-raw_info = self.raw_info.itervalues().next() # Get the first (and only) raw_info
+duration = self.durations.itervalues().next() # Get the first (and only) raw_info
-start_str = "%s %s" % (raw_info['StartDate'], raw_info['StartTime'])
-end_str = "%s %s" % (raw_info['EndDate'], raw_info['EndTime'])
-start_time = datetime.datetime.strptime(start_str, "%d/%m/%Y %H:%M:%S")
-end_time = datetime.datetime.strptime(end_str, "%d/%m/%Y %H:%M:%S")
-duration = end_time - start_time
 duration_sec = duration.seconds
 else:
 duration_sec = np.diff(raw_start_in_seconds)[0]
 return duration_sec
-class Licel_channel(Lidar_channel):
+class LicelChannel(LidarChannel):
 def __init__(self, channel_file):
-c = 299792458.0 #Speed of light
+c = 299792458.0 # Speed of light
 self.wavelength = channel_file.wavelength
 self.name = channel_file.channel_name
-self.binwidth = channel_file.dz * 2 / c # in microseconds
+self.binwidth = channel_file.dz * 2 / c  # in seconds
 self.data = {}
 self.resolution = channel_file.dz
-self.z = np.arange(channel_file.number_of_bins) *self.resolution + self.resolution/2.0 #Change: add half bin in the z
+self.z = np.arange(channel_file.number_of_bins) * self.resolution + self.resolution/2.0 #Change: add half bin in the z
 self.points = channel_file.number_of_bins
 self.rc = []
-self.duration  = 60
+self.duration  = channel_file.duration
 def append(self, other):
 if self.info != other.info:
 raise ValueError('Channel info are different. Data can not be combined.')

Mercurial > public > atmospheric_lidar / file comparison

comparison: licel.py

licel.py