<div dir="ltr">Hi André, thank you for your response. <div><br></div><div>I am not using <span style="font-family:arial,sans-serif;font-size:13px">EPSG:900913 anywhere</span></div><div><span style="font-family:arial,sans-serif;font-size:13px"><br>
</span></div><div><span style="font-family:arial,sans-serif;font-size:13px">I got some img files someone else processed in arcgis and I am processing them this way:</span></div><div><span style="font-family:arial,sans-serif;font-size:13px"><br>
</span></div><div><font face="arial, sans-serif">gdal_translate -of GTIFF -co "COMPRESS=LZW" -scale 0 255 0 255 -ot Byte input.img  input.tif</font></div><div><font face="arial, sans-serif"><div>gdalwarp  -t_srs "EPSG:3857" input.tif output_p.tif</div>
<div>gdal2tiles.py output_p.tif  </div><div><br></div><div>I tried forcing the source projection during warp this way</div><div><br></div><div>gdalwarp -s_srs "EPSG:22184" -t_srs "EPSG:3857"  input.tif output_p.tif<br>
</div><div><br></div><div>I am not sure about source projection, I found it out doing:</div><div><br></div><div>gdalinfo -proj4 input.img </div><div><br></div><div>and searching for the proj4 string in the internet. The image is actually from the zone of projection EPSG:22184 so I believe is ok.</div>
<div><br></div><div>In what point of the process you say I should be forcing with a_srs?</div><div><br></div><div>Thanks!</div><div>Marcela</div></font></div><div><span style="font-family:arial,sans-serif;font-size:13px"> </span></div>
</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, May 27, 2014 at 11:53 AM, Andre Joost <span dir="ltr"><<a href="mailto:andre+joost@nurfuerspam.de" target="_blank">andre+joost@nurfuerspam.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Am 27.05.2014 05:31, schrieb Marcela Campo:<div class=""><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
   I am trying to use gdal2tiles with a set of geotiff and I get the<br>
following error<br>
<br>
ERROR 6: EPSG PCS/GCS code 900913 not found in EPSG support files.  Is this<br>
a valid<br>
EPSG coordinate system?<br>
ERROR 6: No translation for an empty SRS to PROJ.4 format is known.<br>
<br>
<br>
I have already transformed to projection EPSG:3857, here is the gdalinfo<br>
output for one of the files<br>
</blockquote>
<br>
</div><div class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
the command generates the tiles, but as the images appear a bit shifted I<br>
wonder if it somewhat misinterpreted the projection.<br>
<br>
</blockquote>
<br></div>
Instead of reprojecting to EPSG:900913 (which will place the points on the same spot on the globe), try to assign the projection with -a_srs. This will effectively shift the points.<br>
<br>
If the shift goes in the wrong direction, try EPSG:3395. This is the projection on a ellipsoid rather than a sphere.<br>
<br>
HTH,<br>
André Joost<br>
<br>
______________________________<u></u>_________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/gdal-dev</a><br>
</blockquote></div><br></div>