State Plane conversion to UTM from Oregon Lambert 2992...
Matthew Robbins
las2las -i *.las -sp83 OR_N -feet -elevation_feet -target_utm 10N -target_meter -target_elevation_meter
Evon Silvia
 |
Evon Silvia PLS
Geomatics Specialist 517 SW 2nd Street, Suite 400, Corvallis, OR 97333 P: (541) 752-1204 |
--
Download LAStools at
http://lastools.org
http://rapidlasso.com
Be social with LAStools at
http://facebook.com/LAStools
http://twitter.com/LAStools
http://linkedin.com/groups/LAStools-4408378
Manage your settings at
http://groups.google.com/group/lastools/subscribe
Martin Isenburg
whenever you (or anyone else on this lists) finds an ESPG code that LAStools does not recognize (check the lasinfo report) then please send me an email about it with a link to the file and/or this lasinfo report. Then, as Evon pointed out, check the spatial reference on some online site. For EPSG 2992 google quickly finds this here: http://spatialreference.org/ref/epsg/nad83-oregon-lambert-ft/prettywkt/
PROJCS["NAD83 / Oregon Lambert (ft)",
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4269"]],
UNIT["foot",0.3048,
AUTHORITY["EPSG","9002"]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",43],
PARAMETER["standard_parallel_2",45.5],
PARAMETER["latitude_of_origin",41.75],
PARAMETER["central_meridian",-120.5],
PARAMETER["false_easting",1312335.958],
PARAMETER["false_northing",0],
AUTHORITY["EPSG","2992"],
AXIS["X",EAST],
AXIS["Y",NORTH]]
This tells us we have a lambert conformal conic projection. So either via the GUI (source on the left side, target on the right side) or in the command line we then set the source and target projections (as long as they are on the same ellipsoid). To fill the right number into the right plae of the '-lcc' definition check the README file of las2las as suggested by Evon.
las2las -i points.laz ^
-lcc 1312335.958 0 feet 41.75 -120.5 43 45.5 -elevation_feet ^
-ellipsoid 11 ^
-target_utm 10N -target_elevation_meter ^
-o points_10N.laz
And then check the result with
lasboundary -i points_10N.laz -odix _boundary -okml
las2dem -i points_10N.laz -first_only -hillshade -odix _dsm -opng
The next version of LAStools will know the EPSG code 2992/
Cheers,
Martin @rapidlasso
On Tue, Aug 5, 2014 at 2:09 AM, Evon Silvia <esi...@quantumspatial.com> wrote:
>
> Oregon Statewide Lambert (EPSG 2992 for int'l feet) is not the same as Oregon State Plane North (EPSG 2269 for int'l feet). As far as I know, LAStools does not have an implementation of Oregon Statewide Lambert so you'll need to use user-defined projection parameters as instructed on line 236 of the las2las readme. You can find the Lambert Conformal Conic projection parameters in several places, including the EPSG registry (www.epsg-registry.org) and the oregon.gov website (http://www.oregon.gov/DAS/EISPD/GEO/pages/data/format.aspx).
>
>
> Good luck!
>
>
Heidemann, Hans
H. Karl Heidemann, GISP- Arthur James Balfour
Requested CRS:KarlPhysical Scientist, Lidar ScienceU.S. Geological SurveyMundt Federal Building47914 252nd StreetSioux Falls, SD 57110"Nothing matters very much, and very few things ... matter at all."- Arthur James Balfour
--
Evon Silvia
PROJCS["NAD83 / Oregon Lambert",
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4269"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",43],
PARAMETER["standard_parallel_2",45.5],
PARAMETER["latitude_of_origin",41.75],
PARAMETER["central_meridian",-120.5],
PARAMETER["false_easting",400000],
PARAMETER["false_northing",0],
AUTHORITY["EPSG","2991"],
AXIS["X",EAST],
AXIS["Y",NORTH]]--
Ian Madin
And while you are at it, all new Oregon data will use NAD83 (2011) and ellipsoid 12a
Users can find info about the Oregon projections at
http://www.oregon.gov/DAS/EISPD/GEO/Pages/coordination/projections/projections.aspx
From: last...@googlegroups.com [mailto:last...@googlegroups.com]
On Behalf Of Heidemann, Hans
Sent: Tuesday, August 05, 2014 8:22 AM
To: LAStools; Martin Isenburg
Subject: Re: [LAStools] State Plane conversion to UTM from Oregon Lambert 2992...
Karl
H. Karl Heidemann,
GISP
Physical Scientist, Lidar Science
U.S. Geological Survey
Mundt Federal Building
47914 252nd Street
Sioux Falls, SD 57110
"Nothing matters very much, and very few things ... matter at all."
- Arthur James Balfour
On Tue, Aug 5, 2014 at 10:12 AM, Heidemann, Hans <kheid...@usgs.gov> wrote:
Requested CRS:
Albers Equal Area, CONUS
http://epsg.io/5071
Karl
H. Karl Heidemann,
GISP
david.herries
Thanks Martin
Also be great to have New Zealand Transverse Mercator in the list.
SRID Datum / Projection
---------------------------------------------------------
EPSG:2193 NZGD2000 / New Zealand Transverse Mercator
PROJCS["NZGD2000 / New Zealand Transverse Mercator 2000",
GEOGCS["NZGD2000",
DATUM["New_Zealand_Geodetic_Datum_2000",
SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6167"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4167"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",173],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",1600000],
PARAMETER["false_northing",10000000],
AUTHORITY["EPSG","2193"],
AXIS["Easting",EAST],
Cheers
Dave
Ravi Peters
Op dinsdag 5 augustus 2014 16:45:58 UTC+2 schreef Martin Isenburg:
Chris Inskeep
Heidemann, Hans
H. Karl Heidemann, GISP- Arthur James Balfour
--
Petr.D...@cuzk.cz
Hallo Martin,
Could you add S-JTSK_Easth_North_Projection?
Problem is only about EPSG Code, it was confused by ESRI in 90s, they introduce their own code 1020067. But the right code is 5514 introduced by ČÚZK and authorized by EPSG.
Best Regards,
Petr
PROJCS["S-JTSK_Krovak_East_North",
GEOGCS["GCS_S_JTSK",
DATUM["Jednotne_Trigonometricke_Site_Katastralni",
SPHEROID["Bessel_1841",6377397.155,299.1528128]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Krovak"],
PARAMETER["False_Easting",0],
PARAMETER["False_Northing",0],
PARAMETER["Pseudo_Standard_Parallel_1",78.5],
PARAMETER["Scale_Factor",0.9999],
PARAMETER["Azimuth",30.28813975277778],
PARAMETER["Longitude_Of_Center",24.83333333333333],
PARAMETER["Latitude_Of_Center",49.5],
PARAMETER["X_Scale",-1],
PARAMETER["Y_Scale",1],
PARAMETER["XY_Plane_Rotation",90],
UNIT["Meter",1],
AUTHORITY["EPSG","5514"]]