<div dir="ltr"><div><div><div>Hi<br><br></div>I created a small test case here:<br><br><a href="http://www.catais.org/tmp/pdal_tindex.zip">http://www.catais.org/tmp/pdal_tindex.zip</a><br><br></div>It contains a small las file and the tileindex file. This is what I've done:<br><br></div><div>1) Create tileindex: <br><br>find . -iname "*.las" | pdal tindex tileindex.shp --a_srs EPSG:21781 --t_srs EPSG:21781  --fast_boundary --stdin<br><br></div><div>2) Commented out line 161 to allow "--t_srs". Pdal now does *not* reproject data anymore to EPSG:4326. But when trying to write the las file it stops with the following error:<br><br>PDAL: Unable to convert scaled value (2.04e+10) to int32 for dimension 'X' when writing LAS/LAZ file tmp.las.<br><br></div><div>3) Commented out writerOptions.add() methods for scales and offsets on line 445 - 450 and added setCommonOptions(writerOptions).  <br><br></div><div>Now I get the desired result by executing:<br><br>pdal tindex --scale "1 1 1" --merge tileindex.shp --lyr_name tileindex --a_srs EPSG:21781 --t_srs EPSG:21781 tmp.las<br><br></div><div>So, this works for me but definitely introduces some evil stuff :) See also attached patch to exactly see what I've done.<br><br></div><div>regards<br></div><div>Stefan<br></div><div><div><div><div><br></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 14, 2015 at 6:16 AM, 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>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<span class="HOEnZb"><font color="#888888"><br></font></span></div><span class="HOEnZb"><font color="#888888"><div>Stefan<br></div><div><div><div><br><br><div><div><br></div></div></div></div></div></font></span></div><div class="HOEnZb"><div class="h5"><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>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><font color="#888888"><div> </div></font></span></div><span><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>
</div></div></blockquote></div><br></div>