[Proj] Projection conversion problem

Damian Rakowski damian.rakowski at yahoo.com
Mon Jun 27 04:43:18 EST 2016


Hi,

my names is Damian, I'm new to proj and projections topic. 
I need to perform transformation between two projections, and I have a problem defining projection string properly. 

I have coordinates in WGS84 like:

2.2945 48.858222

and I need to transform them to WGS84-like projection with different unit - encoded on integer. 
Basically what I need to do is to transform coordinates in the following way:

floor(coordinate * 2^32/360)
2^32/360 = 11930464,711111111111111111111111


After conversion I expect to get the following coordinates:

27374451 582901293


I tried:

cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs  +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +k_0=11930464,711111111111111111111111 << EOF
> 2.2945 48.858222
> EOF
2d17'40.2"E    48d51'29.599"N 0.000


cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs  +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +k=11930464,711111111111111111111111 << EOF
>2.2945 48.858222
>EOF

2d17'40.2"E    48d51'29.599"N 0.000

cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs  +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +m=11930464,711111111111111111111111 << EOF
2.2945 48.858222
EOF

The same for reverse conversion:

cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +k_0=11930464,711111111111111111111111  +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs << EOF         
27374451 582901293
EOF

27374451dE    582901293dN 0.000


So it seems that k, k_0, m does not bring any effect. I must not understand or missing something. Could you please help me?

Best Regards,

Damian 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/proj/attachments/20160627/a660bd90/attachment.htm 


More information about the Proj mailing list