Changes between Version 3 and Version 4 of Ticket #545


Ignore:
Timestamp:
Dec 4, 2014, 5:40:04 PM (10 years ago)
Author:
dac
Comment:

Scripts now exist to create a DEM from ASTER, NextMap or SRTM for use in APL (create_apl_dem.py) and to create a DEM from lidar data, optionally patching with ASTER or NextMap data and performing horizontal and vertical reprojection as required (create_lidar_dem.py)

The paths for DEMs, separation files, etc., are read in from a config file (arsf_dem.cfg) installed alongside the library, these can be overridden by a copy in the home directory or processing folder. The justification for using a config file rather than defining in the Python library is the ultimate aim of releasing the software. A config fill will make it easier for users to set paths for their system.

Scripts are currently being tested.

Legend:

Unmodified
Added
Removed
Modified

  • Ticket #545 – Description

    v3 v4  
    2020Create DEM subset to navigation data for use in APL. Locations of existing DEMs (ASTER and NMB) will be hardcoded but will allow any DEM mosaic (real or virtual raster) to be provided.
    2121
    22 Inputs:   
    23    --aster [bool]
    24    --nmb   [bool]
    25    --demmosaic [File, if not ASTER or NMB]
    26    --nav [Navigation file]
    27    --projectpath [Path to project]
    28 
    29 Outputs:
    30    --outdem [File]
     22{{{
     23  -o Out DEM, --outdem Out DEM
     24                        Output name for DEM
     25  -n Nav file, --nav Nav file
     26                        Navigation data (.sol / .sbet file)
     27  -p Main project directory, --project Main project directory
     28                        Main project directory (default=".")
     29  --aster               Use ASTER data
     30                        (/users/rsg/arsf/aster/aster_15m_dem_mosaic.vrt)
     31  --nextmap             Use Nextmap data (/users/rsg/arsf/nextmap/neodc/nextma
     32                        p_dsm_mosaic_bng.vrt)
     33  --srtm                Use SRTM data (/local1/data/basedata/srtm/srtm_global_
     34                        mosaic_90m.tif)
     35  --demmosaic Input DEM mosaic
     36                        Input DEM mosaic. For non-standard DEM. Use "--aster"
     37                        or "--nextmap" for standard DEMs.
     38  --separation_file Seperation file
     39                        File with Height offset to add if "--demmosaic" is
     40                        used and DEM heights are not relative to WGS-84
     41                        elepsoid. Not required if using "--aster", "--nextmap"
     42                        or "--srtm" for standard DEMs.
     43  -b BIL Navigation Files, --bil_navigation BIL Navigation Files
     44                        Directory containing post-processed navigation files
     45                        in BIL format. By default raw navigation data will be
     46                        used for "--project". If this is not available (e.g.,
     47                        for ARSF delivered data use this option and point to
     48                        "flightlines/navigation" within delivery directory
     49  --keepgrassdb         Keep GRASS database (default=False)
     50}}}
    3151
    3252''create_dem_from_lidar.py''
    3353Create DEM from lidar data
    34 Inputs:
    35    lidar files(s)
    36    --inproj [shouldn't be needed if projection is stored properly with lidar]
    37    --outproj [default same as input]
    38    --aster [bool - patch with aster]
    39    --nmb [bool - patch with nmb]
    40    --mosaic [mosaic if multiple input files are provided]
    41 Outputs:
    42    --outdem [File]
    43    --screenshots [File(s) optional]
     54{{{
     55  -o Out DEM, --outdem Out DEM
     56                        Output name for DEM
     57  -s Out Screenshot File or Directory, --screenshot Out Screenshot File or Directory
     58                        Output directory for screenshots or single file for
     59                        screenshot of mosaic, in JPEG format.
     60  --las                 Input LiDAR data are in LAS format (default=True)
     61  --ascii               Input LiDAR data are in ASCII format (default=False)
     62  -r Resolution, --resolution Resolution
     63                        Resolution for output DEM (default=2)
     64  --in_projection In Projection
     65                        Input projection (e.g., UTM30N; default=UKBNG)
     66  --out_projection Out Projection
     67                        Out projection. Default is same as input
     68  -n Nav file, --nav Nav file
     69                        Navigation data (.sbet / .sol file) used if patching
     70                        with another DEM
     71  -p Main project directory, --project Main project directory
     72                        Main project directory, used if patching with another
     73                        DEM
     74  --demmosaic Input DEM mosaic
     75                        Input DEM mosaic to patch with in GDAL compatible
     76                        format. Vertical datum needs to be the same as output
     77                        projection. Only required for non-standard DEM. Use "
     78                        --aster" or "--nextmap" for standard DEMs.
     79  --aster               Patch with ASTER data
     80                        (/users/rsg/arsf/aster/aster_15m_dem_mosaic.vrt)
     81  --nextmap             Patch with Nextmap data (/users/rsg/arsf/nextmap/neodc
     82                        /nextmap_dsm_mosaic_bng.vrt)
     83  --hyperspectral_bounds
     84                        If patching with another DEM, get extent from
     85                        hyperspectral navigation data, recommended if DEM is
     86                        to be used with APL and navigation data are available.
     87  --lidar_bounds        If patching with another DEM, get extent from lidar
     88                        data plus default buffer of 2000 m. If DEM is not
     89                        required to be used with APL this option is
     90                        recommended.
     91}}}
    4492
    4593'''New DEM Library'''
    … …  
    4896{{{
    4997|-- arsf_dem
     98|   |-- arsf_dem.cfg
     99|   |-- common_functions.py
    50100|   |-- dem_common.py
    51101|   |-- dem_lidar
     102|   |   |-- ascii_lidar.py
    52103|   |   |-- grass_lidar.py
     104|   |   |-- __init__.py
    53105|   |   |-- laspy_lidar.py
    54106|   |   |-- lastools_lidar.py
    55107|   |   |-- lidar_utilities.py
    56108|   |   |-- spdlib_lidar.py
    57 |   |   |-- ascii_lidar.py
    58109|   |-- dem_nav_utilities.py
    59110|   |-- dem_utilities.py
    60 |-- scripts
    61 |   |-- create_apl_dem.py
    62 |   |-- create_dem_from_lidar.py
    63 
     111|   |-- grass_library.py
    64112}}}