Re: [LAStools] Coordinate projection from NAD_1983_CSRS_UTM_Zone_21N to NAD_1983_CSRS_MTM_3
Martin Isenburg
The command you are looking for is:
las2las -v -i lidar.laz ^
-ellipsoid 11 ^
-tm 500000 0 meter 0 -57 0.9996 ^
-target_tm 304800 0 meter 0 -58.5 0.9999 ^
-odix "_mtm3" -olaz
Sometimes it helps to use the GUI to figure out the correct parameters (see attached screen-shot and output).
las2las -v -i "C:\lastools\bin\eh.laz" -target_tm 304800 0 meter 0 -58.5 0.9999899864 -odir "reprojected/" -odix "_mtm3" -olaz -tm 500000 0 meter 0 -57 0.9995999932 -ellipsoid 11
command was manually edited. new command:
las2las -v -i "C:\lastools\bin\eh.laz" -target_tm 304800 0 meter 0 -58.5 0.9999 -odir "reprojected/" -odix "_mtm3" -olaz -tm 500000 0 meter 0 -57 0.9996 -ellipsoid 11
using TM target projection 'false east/north: 304800/0 [m], origin lat/meridian long: 0/-58.5, scale: 0.9999'
using TM projection 'false east/north: 500000/0 [m], origin lat/meridian long: 0/-57, scale: 0.9996'
using ellipsoid '11 - GRS 1980 (NAD-83) (6.37814e+006 0.00669438)'
reading 1360219 and writing all surviving points ...
total time: 41.39 sec. written 1360219 surviving points.
Cheers,
Martin @rapidlasso
--
http://rapidlasso.com - fast tools to catch reality
Hi All,
I am struggling with the conversion of las data in UTM 21N to MTM zone3. I don't know what order the parameters are supposed to be in. The readme file does not specify what the values in the examples are. For example in the las2las readme we have:
-target_utm 12T : output is UTM zone 12T
-target_sp83 CO_S : output is state plane NAD83 Colorado South
-target_sp27 SC_N : output is state plane NAD27 South Carolina North
-target_longlat : output is geometric coordinates in longitude/latitude
-target_latlong : output is geometric coordinates in latitude/longitude
-target_survey_feet : output uses survey feet
-target_feet : output uses feet
-target_meter : output uses meter
-target_elevation_surveyfeet : output uses survey feet for elevation
-target_elevation_feet : output uses feet for elevation
-target_elevation_meter : output uses meter for elevation
-tm 609601.22 0.0 meter 33.75 -79 0.99996 : specifies a transverse mercator projection
-tm 1804461.942257 0.0 feet 0.8203047 -2.1089395 0.99996
-lcc 609601.22 0.0 meter 33.75 -79 34.33333 36.16666 : specifies a lambertian conic confomal projection
-lcc 1640416.666667 0.0 surveyfeet 47.000000 -120.833333 47.50 48.733333
-ellipsoid 23 : use WGS-84 (specify '-ellipsoid -1' for a list)
There is no -target_MTM and I am not sure if -tm means output coordinate system or what order the values are supposed to be in.
I want to go from:
NAD_1983_CSRS_UTM_Zone_21N
WKID: 2962 Authority: EPSG
Projection: Transverse_Mercator
False_Easting: 500000.0
False_Northing: 0.0
Central_Meridian: -57.0
Scale_Factor: 0.9996
Latitude_Of_Origin: 0.0
Linear Unit: Meter (1.0)
Geographic Coordinate System: GCS_North_American_1983_CSRS
Angular Unit: Degree (0.0174532925199433)
Prime Meridian: Greenwich (0.0)
Datum: D_North_American_1983_CSRS
Spheroid: GRS_1980
Semimajor Axis: 6378137.0
Semiminor Axis: 6356752.314140356
Inverse Flattening: 298.257222101
TO
NAD_1983_CSRS_MTM_3
WKID: 2945 Authority: EPSG
Projection: Transverse_Mercator
False_Easting: 304800.0
False_Northing: 0.0
Central_Meridian: -58.5
Scale_Factor: 0.9999
Latitude_Of_Origin: 0.0
Linear Unit: Meter (1.0)
Geographic Coordinate System: GCS_North_American_1983_CSRS
Angular Unit: Degree (0.0174532925199433)
Prime Meridian: Greenwich (0.0)
Datum: D_North_American_1983_CSRS
Spheroid: GRS_1980
Semimajor Axis: 6378137.0
Semiminor Axis: 6356752.314140356
Inverse Flattening: 298.257222101
Thanks.
