[Proj] Projection conversion problem

Damian Rakowski damian.rakowski at yahoo.com
Thu Jun 30 01:20:48 PDT 2016


That  looks very interesting, thank you for your reply.
One more question:
6371000 - this is an Earth radius0.0093202511 - how did you got this value? I just want to understand :)
Best Regards,
Damian
 

    19:18 środa, 2016-6-29, Andre Joost <andre+joost at nurfuerspam.de> napisał(a):
 

 Am 28.06.2016 um 19:40 schrieb Hermann Peifer:

>
> With the given sample coordinates, this will give:
>
> $ echo 2.2945 48.858222 | cs2cs -v +proj=longlat +datum=WGS84 +no_defs
> +to +proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84
> +units=m +no_defs
> # ---- From Coordinate System ----
> #Lat/long (Geodetic alias)
> #
> # +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
> # ---- To Coordinate System ----
> #Equidistant Cylindrical (Plate Caree)
> #    Cyl, Sph
> #    lat_ts=[, lat_0=0]
> # +proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m
> # +no_defs +ellps=WGS84 +towgs84=0,0,0
> 255422.57    5438872.39 0.00
>
> The expected result was: 27374451 582901293
>

I suggest to take a sphere, and scale with the units:

+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +R=6371000 
+to_meter=0.0093202511 +no_defs

which results to 27374451.23    582901292.36 0.00

HTH,
André Joost


_______________________________________________
Proj mailing list
Proj at lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj

  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/proj/attachments/20160630/380c003f/attachment.html>


More information about the Proj mailing list