Hello<br><br>The third definition is ok. Since you are specifying both semi-axis the flattening is redundant. The first is wrong because +rf should be infinite in the case of the sphere. Within PROJ syntax you could also use:<br>
<br>proj +proj=laea +lat_0=48 +lon_0=9 +x_0=0 +y_0=0 +R=6378388 +units=m +no_defs<br><br>Regards<br><br>José A. Gonçalves<br><br><br><br><div class="gmail_quote">2010/6/24 Agustin Lobo <span dir="ltr"><<a href="mailto:alobolistas@gmail.com">alobolistas@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi!<br>
<br>
The European Forest Map has the following projection (according to a<br>
file distributed along with the file):<br>
<br>
Projection: Lambert_Azimuthal_Equal_Area<br>
Parameters:<br>
False_Easting: 0.000000<br>
False_Northing: 0.000000<br>
Central_Meridian: 9.000000<br>
Latitude_Of_Origin: 48.000000<br>
Linear Unit: Meter (1.000000)<br>
Coordinate system<br>
Geographic Coordinate System:<br>
Name: GCS_User_Defined<br>
Angular Unit: Degree (0.017453292519943299)<br>
Prime Meridian: Greenwich (0.000000000000000000)<br>
Datum: D_User_Defined<br>
Spheroid: User_Defined_Spheroid<br>
Semimajor Axis: 6378388.000000000000000000<br>
Semiminor Axis: 6378388.000000000000000000<br>
Inverse Flattening: 0.000000000000000000<br>
<br>
I've made several trials to set up a custom projection (in QGIS):<br>
1.<br>
+proj=laea +lat_0=48 +lon_0=9 +x_0=0.000 +y_0=0.000 +a=6378388.0000<br>
+b=6378388.000 +rf=0.000 +units=m +no_defs<br>
Rejected by the custom projection tool as erroneous<br>
<br>
2.<br>
+proj=laea +lat_0=48 +lon_0=9 +x_0=0.000 +y_0=0.000 +a=6378388.0000<br>
+b=6378388.000 +f=0.000 +units=m +no_defs<br>
Accepted by the custom projection tool but the Properties dialog does<br>
not accept it<br>
<br>
3.<br>
+proj=laea +lat_0=48 +lon_0=9 +x_0=0.000 +y_0=0.000 +a=6378388.0000<br>
+b=6378388.000 +units=m +no_defs<br>
Seems to work fine.<br>
<br>
I'd like if someone could actually confirm that the third proj4<br>
definition is actually<br>
in accordance with the definition provided in the server and that I<br>
have reproduced at the beginning of this<br>
message. I'm not sure if option 1 and 2 are rejected because they are wrong or<br>
this is a QGIS problem (message sent to the QGIS list remains unanswered).<br>
<br>
I also would like answer if there is a document instructing how to<br>
translate PROJCS specifications<br>
into proj4 format.<br>
<br>
Thanks!<br>
<br>
Agus<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" target="_blank">http://lists.maptools.org/mailman/listinfo/proj</a><br>
</blockquote></div><br>