<div dir="ltr"><div><div><div><div>find . -iname "*.laz" | pdal tindex tileindex.shp --a_srs EPSG:21781 --t_srs EPSG:21781 --fast_boundary --stdin<br><br></div>ogrinfo -al tileindex.shp:<br><br>Layer name: tileindex<br>Geometry: Polygon<br>Feature Count: 6<br>Extent: (592000.000000, 226000.000000) - (594000.000000, 230000.000000)<br>Layer SRS WKT:<br>PROJCS["CH1903_LV03",<br> GEOGCS["GCS_CH1903",<br> DATUM["CH1903",<br> SPHEROID["Bessel_1841",6377397.155,299.1528128]],<br> PRIMEM["Greenwich",0],<br> UNIT["Degree",0.017453292519943295]],<br> PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],<br> PARAMETER["latitude_of_center",46.95240555555556],<br> PARAMETER["longitude_of_center",7.439583333333333],<br> PARAMETER["azimuth",90],<br> PARAMETER["scale_factor",1],<br> PARAMETER["false_easting",600000],<br> PARAMETER["false_northing",200000],<br> PARAMETER["rectified_grid_angle",90],<br> UNIT["Meter",1]]<br>location: String (254.0)<br>srs: String (254.0)<br>modified: Date (10.0)<br>created: Date (10.0)<br>OGRFeature(tileindex):0<br> location (String) = /home/stefan/tmp/lidar/srs/./LAS_593227.laz<br> srs (String) = EPSG:21781<br> modified (Date) = 2015/08/13<br> created (Date) = 2015/08/13<br> POLYGON ((593000 227000,593000 228000,594000 228000,594000 227000,593000 227000))<br><br></div>pdal info LAS_593227.laz --metadata:<br><br>{<br> "filename": "LAS_593227.laz",<br> "metadata":<br> {<br> "comp_spatialreference": "PROJCS[\"CH1903 / LV03\",GEOGCS[\"CH1903\",DATUM[\"CH1903\",SPHEROID[\"Bessel 1841\",6377397.155,299.1528128,AUTHORITY[\"EPSG\",\"7004\"]],TOWGS84[674.4,15.1,405.3,0,0,0,0],AUTHORITY[\"EPSG\",\"6149\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4149\"]],PROJECTION[\"Hotine_Oblique_Mercator_Azimuth_Center\"],PARAMETER[\"latitude_of_center\",46.95240555555556],PARAMETER[\"longitude_of_center\",7.439583333333333],PARAMETER[\"azimuth\",90],PARAMETER[\"rectified_grid_angle\",90],PARAMETER[\"scale_factor\",1],PARAMETER[\"false_easting\",600000],PARAMETER[\"false_northing\",200000],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Y\",EAST],AXIS[\"X\",NORTH],AUTHORITY[\"EPSG\",\"21781\"]]",<br> "compressed": true,<br> "count": 14724002,<br> "creation_doy": 224,<br> "creation_year": 2015,<br> "dataformat_id": 3,<br> "dataoffset": 2181,<br> "filesource_id": 0,<br> "global_encoding": "AAA=",<br> "header_size": 227,<br> "major_version": 1,<br> "maxx": 593999.98999999999,<br> "maxy": 227999.98999999999,<br> "maxz": 1667.6600000000001,<br> "minor_version": 2,<br> "minx": 593000,<br> "miny": 227000,<br> "minz": 622.80000000000007,<br> "offset_x": 0,<br> "offset_y": 0,<br> "offset_z": 0,<br> "project_id": "00000000-0000-0000-0000-000000000000",<br> "scale_x": 0.01,<br> "scale_y": 0.01,<br> "scale_z": 0.01,<br> "software_id": "PDAL 1.0.0.b1 (e412bd)",<br> "spatialreference": "PROJCS[\"CH1903 / LV03\",GEOGCS[\"CH1903\",DATUM[\"CH1903\",SPHEROID[\"Bessel 1841\",6377397.155,299.1528128,AUTHORITY[\"EPSG\",\"7004\"]],TOWGS84[674.4,15.1,405.3,0,0,0,0],AUTHORITY[\"EPSG\",\"6149\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4149\"]],PROJECTION[\"Hotine_Oblique_Mercator_Azimuth_Center\"],PARAMETER[\"latitude_of_center\",46.95240555555556],PARAMETER[\"longitude_of_center\",7.439583333333333],PARAMETER[\"azimuth\",90],PARAMETER[\"rectified_grid_angle\",90],PARAMETER[\"scale_factor\",1],PARAMETER[\"false_easting\",600000],PARAMETER[\"false_northing\",200000],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Y\",EAST],AXIS[\"X\",NORTH],AUTHORITY[\"EPSG\",\"21781\"]]",<br> "system_id": "PDAL"<br> },<br> "pdal_version": "PDAL 1.0.0.b1 (e412bd)"<br>}<br><br></div>In the mergeFile() method on line 413 you assign the output srs for the reprojection stage. But this cannot be set by the user and will therefore always be EPSG:4326. As far as I understand.<br><br></div><div>regards<br></div><div>Stefan<br></div><div><div><div><br><br><div><div><br></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 14, 2015 at 12:10 AM, Andrew Bell <span dir="ltr"><<a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="">On Thu, Aug 13, 2015 at 2:46 PM, Stefan Ziegler <span dir="ltr"><<a href="mailto:stefan.ziegler.de@gmail.com" target="_blank">stefan.ziegler.de@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div>Hi<br><br></div>I'm trying to merge some las files with the following command:<br><br>pdal tindex --merge /home/stefan/tmp/lidar/srs/tileindex.gpkg --lyr_name tileindex tmp.las<br><br></div>The input las files have an srs assigned (EPSG:21781). But the merged one is EPSG:4326.<br></div></div></div></div></div></div></blockquote><div><br></div></span><div>Can you use ogrinfo to verify the SRS of the layer? Can you provide the command used to create the index file?</div><div><br></div><div>Thanks,</div><span class="HOEnZb"><font color="#888888"><div> </div></font></span></div><span class="HOEnZb"><font color="#888888">-- <br><div>Andrew Bell<br><a href="mailto:andrew.bell.ia@gmail.com" target="_blank">andrew.bell.ia@gmail.com</a></div>
</font></span></div></div>
</blockquote></div><br></div>