example_script/convert_lilas.py

changeset 49
f22be112deba
parent 48
673843d7f9d8
child 50
510d7ac14280
--- 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 <your options>
-
-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

mercurial