[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