Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
6 messages Options
jpr
Reply | Threaded
Open this post in threaded view
|

Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

jpr
Hi All,

I am trying to convert some points from rpoly to utm. I have tried typing Cartesian coordinates and cartographic coordinates and both result in the error: 
<cs2cs>: while processing file: <stdin>, line 1
pj_transformation(): Unknown error

For some reason it works correctly if I go from utm to rpoly, but that is not what I am trying to do.
Here is my code:
cs2cs +proj=rpoly +lat_ts=25.066667 +lon_0=-80.636667 +ellps=bessel +to +proj=utm +zone 17 +ellps=WGS84
I tried inputting: 
-80.4666667 +24.983333333
and got:   *                  *       0.00
17151.59  2799896.42
and got:    *                 *       0.00

Can anyone tell me what I am doing wrong?

Thanks!

_______________________________________________
Proj mailing list
[hidden email]
http://lists.maptools.org/mailman/listinfo/proj
Reply | Threaded
Open this post in threaded view
|

Re: Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

John Buonagurio
Hi,

PROJ.4 only supports the forward transformation (from lon/lat to projected). The inverse transformation isn't implemented unfortunately for rectangular polyconic:

<https://github.com/OSGeo/proj.4/blob/master/src/PJ_rpoly.c#L12>

("no inv.")

Has anyone looked into this? The inverse is described in Yang/Snyder/Tobler.

John


From: [hidden email] [mailto:[hidden email]] On Behalf Of Resnick, Jessica
Sent: Tuesday, September 27, 2016 2:06 PM
To: [hidden email]
Subject: [Proj] Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

Hi All,

I am trying to convert some points from rpoly to utm. I have tried typing Cartesian coordinates and cartographic coordinates and both result in the error: 
<cs2cs>: while processing file: <stdin>, line 1
pj_transformation(): Unknown error

For some reason it works correctly if I go from utm to rpoly, but that is not what I am trying to do.
Here is my code:
cs2cs +proj=rpoly +lat_ts=25.066667 +lon_0=-80.636667 +ellps=bessel +to +proj=utm +zone 17 +ellps=WGS84
I tried inputting: 
-80.4666667 +24.983333333
and got:   *                  *       0.00
17151.59  2799896.42
and got:    *                 *       0.00

Can anyone tell me what I am doing wrong?

Thanks!
_______________________________________________
Proj mailing list
[hidden email]
http://lists.maptools.org/mailman/listinfo/proj
Reply | Threaded
Open this post in threaded view
|

Re: Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

Kristian Evers-2
I had a quick look in Yang, Snyder & Tobler and it should be fairly easy to implement the inverse projection. I've created a ticket for it: https://github.com/OSGeo/proj.4/issues/420

Feel free to submit a pull request!

/Kristian

> -----Oprindelig meddelelse-----
> Fra: [hidden email] [mailto:proj-
> [hidden email]] På vegne af John Buonagurio
> Sendt: 27. september 2016 23:16
> Til: PROJ.4 and general Projections Discussions
> Emne: Re: [Proj] Receiving "pj_transform(): Unknown error" while trying to
> use cs2cs from rpoly to utm
>
> Hi,
>
> PROJ.4 only supports the forward transformation (from lon/lat to projected).
> The inverse transformation isn't implemented unfortunately for rectangular
> polyconic:
>
> <https://github.com/OSGeo/proj.4/blob/master/src/PJ_rpoly.c#L12>
>
> ("no inv.")
>
> Has anyone looked into this? The inverse is described in Yang/Snyder/Tobler.
>
> John
>
>
> From: [hidden email] [mailto:proj-
> [hidden email]] On Behalf Of Resnick, Jessica
> Sent: Tuesday, September 27, 2016 2:06 PM
> To: [hidden email]
> Subject: [Proj] Receiving "pj_transform(): Unknown error" while trying to
> use cs2cs from rpoly to utm
>
> Hi All,
>
> I am trying to convert some points from rpoly to utm. I have tried typing
> Cartesian coordinates and cartographic coordinates and both result in the
> error:
> <cs2cs>: while processing file: <stdin>, line 1
> pj_transformation(): Unknown error
>
> For some reason it works correctly if I go from utm to rpoly, but that is not
> what I am trying to do.
> Here is my code:
> cs2cs +proj=rpoly +lat_ts=25.066667 +lon_0=-80.636667 +ellps=bessel +to
> +proj=utm +zone 17 +ellps=WGS84
> I tried inputting:
> -80.4666667 +24.983333333
> and got:   *                  *       0.00
> 17151.59  2799896.42
> and got:    *                 *       0.00
>
> Can anyone tell me what I am doing wrong?
>
> Thanks!
> _______________________________________________
> Proj mailing list
> [hidden email]
> http://lists.maptools.org/mailman/listinfo/proj
_______________________________________________
Proj mailing list
[hidden email]
http://lists.maptools.org/mailman/listinfo/proj
jpr
Reply | Threaded
Open this post in threaded view
|

Re: Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

jpr
In reply to this post by John Buonagurio
Thanks for your response John.

Is there any way for me to convert rpoly to anything at all that is more modern? I just need to get my coordinates into something that can be used in a GIS software.

Thanks,
Jessica
Reply | Threaded
Open this post in threaded view
|

Re: Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

John Buonagurio
In reply to this post by Kristian Evers-2
Thanks, my thoughts as well. The calculations for this projection are very similar to the regular (American) polyconic.

I'll give it a shot - working on implementing the inverse on the sphere. Will submit via pull request.

John

-----Original Message-----
From: [hidden email] [mailto:[hidden email]] On Behalf Of Kristian Evers
Sent: Wednesday, September 28, 2016 3:57 AM
To: PROJ.4 and general Projections Discussions <[hidden email]>
Subject: Re: [Proj] Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

I had a quick look in Yang, Snyder & Tobler and it should be fairly easy to implement the inverse projection. I've created a ticket for it: https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_OSGeo_proj.4_issues_420&d=DQIGaQ&c=t0wRGL5ICVzH157W8C8Wew&r=Q-Rs-yPWbL0-v_m_Ov9CS1uGraVrKdE4mvUkBfZ3rWE&m=9qi7q5QoQWxg6vZjCV32eGaFCcF3h2BpKvbInYryzts&s=kz8sB_TLq3wDwp1Gvttw4ZxIq8a-L40nxugvAoI-Lps&e= 



Feel free to submit a pull request!



/Kristian



> -----Oprindelig meddelelse-----

> Fra: [hidden email] [mailto:proj-

> [hidden email]] På vegne af John Buonagurio

> Sendt: 27. september 2016 23:16

> Til: PROJ.4 and general Projections Discussions

> Emne: Re: [Proj] Receiving "pj_transform(): Unknown error" while trying to

> use cs2cs from rpoly to utm

>

> Hi,

>

> PROJ.4 only supports the forward transformation (from lon/lat to projected).

> The inverse transformation isn't implemented unfortunately for rectangular

> polyconic:

>

> <https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_OSGeo_proj.4_blob_master_src_PJ-5Frpoly.c-23L12&d=DQIGaQ&c=t0wRGL5ICVzH157W8C8Wew&r=Q-Rs-yPWbL0-v_m_Ov9CS1uGraVrKdE4mvUkBfZ3rWE&m=9qi7q5QoQWxg6vZjCV32eGaFCcF3h2BpKvbInYryzts&s=0CIu0APRS-wSnhnKpZRmHYh2htrzQNE_AQAz-wYulAg&e= >

>

> ("no inv.")

>

> Has anyone looked into this? The inverse is described in Yang/Snyder/Tobler.

>

> John

>

>

> From: [hidden email] [mailto:proj-

> [hidden email]] On Behalf Of Resnick, Jessica

> Sent: Tuesday, September 27, 2016 2:06 PM

> To: [hidden email]

> Subject: [Proj] Receiving "pj_transform(): Unknown error" while trying to

> use cs2cs from rpoly to utm

>

> Hi All,

>

> I am trying to convert some points from rpoly to utm. I have tried typing

> Cartesian coordinates and cartographic coordinates and both result in the

> error:

> <cs2cs>: while processing file: <stdin>, line 1

> pj_transformation(): Unknown error

>

> For some reason it works correctly if I go from utm to rpoly, but that is not

> what I am trying to do.

> Here is my code:

> cs2cs +proj=rpoly +lat_ts=25.066667 +lon_0=-80.636667 +ellps=bessel +to

> +proj=utm +zone 17 +ellps=WGS84

> I tried inputting:

> -80.4666667 +24.983333333

> and got:   *                  *       0.00

> 17151.59  2799896.42

> and got:    *                 *       0.00

>

> Can anyone tell me what I am doing wrong?

>

> Thanks!

> _______________________________________________

> Proj mailing list

> [hidden email]

> https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.maptools.org_mailman_listinfo_proj&d=DQIGaQ&c=t0wRGL5ICVzH157W8C8Wew&r=Q-Rs-yPWbL0-v_m_Ov9CS1uGraVrKdE4mvUkBfZ3rWE&m=9qi7q5QoQWxg6vZjCV32eGaFCcF3h2BpKvbInYryzts&s=5kUXaRL5BcgHzRyfMSdYtVwoAvruAaxQYJiUhoCw6zA&e= 

_______________________________________________
Proj mailing list
[hidden email]
https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.maptools.org_mailman_listinfo_proj&d=DQIGaQ&c=t0wRGL5ICVzH157W8C8Wew&r=Q-Rs-yPWbL0-v_m_Ov9CS1uGraVrKdE4mvUkBfZ3rWE&m=9qi7q5QoQWxg6vZjCV32eGaFCcF3h2BpKvbInYryzts&s=5kUXaRL5BcgHzRyfMSdYtVwoAvruAaxQYJiUhoCw6zA&e= 
_______________________________________________
Proj mailing list
[hidden email]
http://lists.maptools.org/mailman/listinfo/proj
Reply | Threaded
Open this post in threaded view
|

Re: Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

John Buonagurio
In reply to this post by jpr
It's not possible to do that accurately without an inverse transformation. PROJ.4 and GCTP[C] at least don't implement it. I'll work on it over the next few days and keep you updated.

Also, apologies to others on the list for "urldefense.proofpoint.com" links in previous message. Our email server adds these and I forgot to remove.

John

-----Original Message-----
From: [hidden email] [mailto:[hidden email]] On Behalf Of jpr
Sent: Wednesday, September 28, 2016 11:36 AM
To: [hidden email]
Subject: Re: [Proj] Receiving "pj_transform(): Unknown error" while trying to use cs2cs from rpoly to utm

Thanks for your response John.

Is there any way for me to convert rpoly to anything at all that is more modern? I just need to get my coordinates into something that can be used in a GIS software.

Thanks,
Jessica
 
_______________________________________________
Proj mailing list
[hidden email]
http://lists.maptools.org/mailman/listinfo/proj