diff -r 673843d7f9d8 -r f22be112deba example_script/convert_lilas.py --- a/example_script/convert_lilas.py Thu Feb 16 09:22:05 2017 +0200 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,74 +0,0 @@ -""" Sample script to convert licel files to SCC netcdf format. - -The script assumes the following things: - -1. You have already installed the atmospheric_lidar module (e.g. using pip). -2. You have create a class in a file "lilas" describing your system (described the readme file). - If you want to use it for a different system you need to change the script to use your own class - -Run the script using: python convert_lilas.py - -Examples --------- -# This will read all files starting with "l" from "my_dir". -# It will guess the measurement ID based on the date of the files, and will assume measurement number 00. -# For example, the new measurment id could be 20120101mb00 -python convert_lilas.py my_dir l*. - -# This will use the measurement id you defined. -python convert_lilas.py my_dir l*. -m 20120204mb32 # It will create the file 20120204mb32.nc - -Help string ------------ -# To get this, run: python convert_lilas.py -h -usage: convert_lilas.py [-h] [-m MEASUREMENT_ID] [-n MEASUREMENT_NUMBER] - [directory] [searchstring] - -positional arguments: - directory Directory with licel files. - searchstring Processing system id. - -optional arguments: - -h, --help show this help message and exit - -m MEASUREMENT_ID, --measurement_id MEASUREMENT_ID - The new measurement id - -n MEASUREMENT_NUMBER, --measurement_number MEASUREMENT_NUMBER - The measurement number for the date, if no id is - provided - -""" - -import glob -import os -import argparse - -from atmospheric_lidar import lilas - - -if __name__ == "__main__": - - # Define the command line arguments. - parser = argparse.ArgumentParser() - parser.add_argument("directory", nargs='?', help="Directory with licel files.", default='.') - parser.add_argument("searchstring", nargs='?', help="Processing system id.", default="*.*") - parser.add_argument("-m", "--measurement_id", help="The new measurement id", default=None) - parser.add_argument("-n", "--measurement_number", help="The measurement number for the date, if no id is provided", default="00") - args = parser.parse_args() - - - # Get a list of files to convert - search_str = os.path.join(args.directory, args.searchstring) - files = glob.glob(search_str) - - if files: - # Read the files - print "Reading {0} files from {1}".format(len(files), args.directory) - measurement = lilas.LilasLidarMeasurement(files) - - #Save the netcdf - print "Saving netcdf." - measurement.set_measurement_id(args.measurement_id, args.measurement_number) - measurement.save_as_netcdf() - print "Created file ", measurement.scc_filename - else: - print "No files found when searching for ", search_str \ No newline at end of file