Thanks Oscar, all was good now.<br><br>Thanks a lot for all your responses.<br><br><div><span class="gmail_quote">2006/8/8, Oscar van Vlijmen <<a href="mailto:ovv@hetnet.nl">ovv@hetnet.nl</a>>:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
From: "Gheysen Thibaut"<br><br>... This is my first attempt to use proj4 and I don't know what I do wrong.<br>Here is the definition I use :<br>WGS 84 (epsg : 4326) : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
<br>no_defs <><br>Lambert 72 Belge (epsg : 31370) : +proj=lcc +lat_1=51.16666723333333<br>+lat_2=49.8333339 +lat_0=90 +lon_0=4.367486666666666 +x_0=150000.013<br>+y_0=5400088.438 +ellps=intl +units=m +no_defs no_defs <>
<br><br>Reply:<br><br>I don't do PROJ command lines, but it seems to me that you first have to<br>perform a datum transform from WGS84 to International for Belgium. This<br>can't be evoked with the epsg:4326 init since it has no transform
<br>parameters. Find the complete transform parameters, check rotation model and<br>direction (from/to), supply the command line with a proper datum transform,<br>then project.<br><br>Belgium 1972 to WGS84:<br>translation in m, rotation in arcsec, scale in ppm.
<br>Position vector rotation (PROJ)<br>dx=-99.059; dy=53.322; dz=-112.486; rx=0.419; ry=-0.83; rz=1.885; ds=-1.00;<br>Derived from ESRI ArcGIS parameters.<br><br>Going in the other direction, from wgs to intl, just negate all parameters.
<br><br><br>_______________________________________________<br>Proj mailing list<br><a href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br><a href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj
</a><br></blockquote></div><br>