<div dir="ltr"><div>Hi all,</div><div><br></div><div>Thank Jürgen for your  explanation about how QGIS set the reference system from source.</div><div><br></div><div>I agreement with André Joost about the user should  select your specific reference.</div>
<div><br></div><div>The source of my problem is the creating of shapefile by GDAL/OGR.</div><div><br></div><div>For creating shapefile, the OGR/GDAL not using  all parameters of reference system.</div><div><br></div><div>
I resolve my problem, creating the GML file, in this case, the refefence system is identical of source.</div><div><br></div><div>How suggested Etienne i will raise this issue in gdal-dev list.</div><div><br></div><div>If someone want to see the script example:</div>
<div><a href="http://pastebin.com/PbjwdrHR">http://pastebin.com/PbjwdrHR</a></div><div><br></div><div>Regards,</div><div>Luiz Motta</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014/1/17 Andre Joost <span dir="ltr"><<a href="mailto:andre+joost@nurfuerspam.de" target="_blank">andre+joost@nurfuerspam.de</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Am 17.01.2014 11:56, schrieb Luiz Motta:<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
The question is exact this, why the QGIS read the PRJ and understands the<br>
'towgs84', otherwise, the GDAL don't understands, where , the QGIS use the<br>
GDAL/OGR provider ?<br>
<br>
</blockquote>
<br></div>
This is an intended behaviour.<br>
<br>
Other software (like ESRI) don't write datum shifts in the .prj file, but leave it to the end user to select one of the available shifts. QGIS has to handle these files correctly, so it looks up the EPSG table and selects the best fitting code. If you think it has guessed wrong, you can assign another CRS to the layer.<br>

<br>
GDAL does no guessing. Unless there are not datum shifts provided on source and target CRS (using either +datum, +towgs84 or +nadgrids), it does no datum shift at all. You can specify -s_srs to override the source file's .prj file with an EPSG code that contains a datum shift. This should bring the same result as QGIS does automatically.<div class="HOEnZb">
<div class="h5"><br>
<br>
HTH,<br>
André Joost<br>
______________________________<u></u>_________________<br>
Qgis-developer mailing list<br>
<a href="mailto:Qgis-developer@lists.osgeo.org" target="_blank">Qgis-developer@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/qgis-developer" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/qgis-<u></u>developer</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Luiz Motta<div>Analista Ambiental - IBAMA/COTIG<br></div><div><div>Instituto Brasileiro do Meio Ambiente dos Recursos Naturais - IBAMA </div>
<div>IBAMA - SCEN Trecho 2 - Ed. Sede - Cx. Postal nº 09566 </div><div>CEP 70818-900 - Brasília-DF </div><div>Tel: 61.3316-1807</div><div>Skype: ibama_luiz-motta</div></div><div><br></div><div><br></div></div>
</div>