:: 
:: planar patches script for (raw) MLS data in ECEF coordinates
::

::::::::::::::::::::::::
:: define parameters
::::::::::::::::::::::::

set NUM_CORES=8
set UTM_ZONE=34N

:: folder containing input files and format
set INPUT_FOLDER=E:\BKKK\BP006_111_Margit_hid_Budai_MLS
set INPUT_FORMAT=laz
set INPUT_TRAJECTORY=BP006_111_trajectory

:: folder and names for output file and planes
set OUTPUT_FOLDER=D:\lastools\bin\mls_test_output
set OUTPUT_NAME=BP006_111
set POLYGON_NAME=M_006_111_09
set OUTPUT_FORMAT=pef

:: only use points with height from BELOW to ABOVE
set CUTOFF_BELOW=0.5
set CUTOFF_ABOVE=10

:: lasplanes parameters
set CELL_SIZE=1.0
set CELL_POINTS=100
set EIGEN_RATIO_SMALLEST=0.0001
set EIGEN_RATIO_LARGEST=0.9
set PLANE_THICKNESS=0.01
set PLANE_EXCLUSION=5
set PLANE_POINTS=100
set POLYGON_AREA=0.5
set POLYGON_POINTS=8
set POLYGON_STDDEV=0.01

:: convert trajectory from ASCII text to LAZ

txt2las -i %INPUT_FOLDER%\%INPUT_TRAJECTORY%.txt ^
        -parse txyz -skip 1 ^
        -odir %OUTPUT_FOLDER% -o %INPUT_TRAJECTORY%.laz

:: convert trajectory from ECEF to UTM

las2las -i %OUTPUT_FOLDER%\%INPUT_TRAJECTORY%.laz ^
        -ecef ^
        -target_utm %UTM_ZONE% ^
        -target_precision 0.001 ^
        -odir %OUTPUT_FOLDER% -o utm_%INPUT_TRAJECTORY%.laz

:: convert file(s) from ECEF to UTM, preserve high precision, zero classifications

las2las -i %INPUT_FOLDER%\*.%INPUT_FORMAT% ^
        -set_classification 0 ^
        -ecef ^
        -target_utm %UTM_ZONE% ^
        -target_precision 0.001 ^
        -odir %OUTPUT_FOLDER% -odix _utm -olaz ^
        -cores %NUM_CORES%

:: keep only points that are between 0.5 and 10 meters above the track

lastrack -i %OUTPUT_FOLDER%\*_utm.laz ^
         -track %OUTPUT_FOLDER%\utm_%INPUT_TRAJECTORY%.laz ^
         -drop_below 0.5 -drop_above 10 ^
         -odix _cut -olaz ^
         -cores %NUM_CORES%

:: convert file from UTM34 back to ECEF

las2las -i %OUTPUT_FOLDER%\*_utm_cut.laz ^
        -target_ecef ^
        -target_precision 0.001 ^
        -odir %OUTPUT_FOLDER% -odix _ecef -olaz ^
        -cores %NUM_CORES%

:: find planes 

lasplanes -v -i %OUTPUT_FOLDER%\*_utm_cut_ecef.laz -merged ^
          -cell_size %CELL_SIZE% ^
          -cell_points %CELL_POINTS% ^
          -eigen_ratio_smallest %EIGEN_RATIO_SMALLEST% ^
          -eigen_ratio_largest %EIGEN_RATIO_LARGEST% ^
          -plane_thickness %PLANE_THICKNESS% ^
          -plane_exclusion %PLANE_EXCLUSION% ^
          -plane_points %PLANE_POINTS% ^
          -polygon_area %POLYGON_AREA% ^
          -polygon_points %POLYGON_POINTS% ^
          -polygon_stddev %POLYGON_STDDEV% ^
          -polygon_name %POLYGON_NAME% ^
          -odir %OUTPUT_FOLDER% -o %OUTPUT_NAME%.%OUTPUT_FORMAT% -o%OUTPUT_FORMAT%

GOTO:EOF

:: clean-up temp files

del %OUTPUT_FOLDER%\%INPUT_TRAJECTORY%.laz
del %OUTPUT_FOLDER%\utm_%INPUT_TRAJECTORY%.laz
del %OUTPUT_FOLDER%\*_utm.laz
del %OUTPUT_FOLDER%\%OUTPUT_NAME%_thinned.laz
del %OUTPUT_FOLDER%\%OUTPUT_NAME%_classified.laz
del %OUTPUT_FOLDER%\%OUTPUT_NAME%_ground.laz
del %OUTPUT_FOLDER%\*_utm_cut.laz
del %OUTPUT_FOLDER%\*_utm_cut_ecef.laz