[Proj] Convergence

[Kees Krikke]

Noel,

No i'm not in the south. The line with zone= and south appears because i 
use the -v flag.

I always use the -v option with proj. It gives some additional 
information about parameters that can be used with a projection (like 
zone and south for utm) and a complete list of parameters that are 
actually used. This allows you to see all the parameters including those 
not set by directly by the user.


Regards,

Kees


Noel Zinn wrote:
> Pardon my proj4 syntax density, but are you in the south with TM and the
> north with LCC?  If so, that's the reason.  -Noel
> 
> -----Original Message-----
> From: proj-bounces at lists.maptools.org
> [mailto:proj-bounces at lists.maptools.org] On Behalf Of Kees Krikke
> Sent: Monday, September 15, 2008 10:28 AM
> To: PROJ.4 and general Projections Discussions
> Subject: [Proj] Convergence
> 
> Hello,
> 
> 
> I've made a piece of s/w that uses proj's convergence calculation to 
> find the true north for a given grid north.
> 
> When testing with a utm and lcc projection i found some confusing 
> results regarding the sign (direction) of the convergence.
> 
> To demonstrate this i modified proj.exe so it prints the convergence (in 
> deg) at the end of the line when using the -S flag:
> 
> 
> C:\proj-4.6.1\src>proj -vS +proj=utm +datum=WGS84 +lon_0=3
> #Universal Transverse Mercator (UTM)
> #       Cyl, Sph
> #       zone= south
> # +proj=utm +datum=WGS84 +lon_0=3 +ellps=WGS84 +towgs84=0,0,0
> 2 52
> 431350.30       5761510.32      <0.999658 0.999658 0.999316 1.20783e-006 
> 0.999658 0.999658 -0.788041>
> 4 52
> 568649.70       5761510.32      <0.999658 0.999658 0.999316 8.54066e-007 
> 0.999658 0.999658 0.788041>
> 4 52.1
> 568496.62       5772632.28      <0.999658 0.999658 0.999315 0 0.999658 
> 0.999658 0.789115>
> 
>  From the above results it seems that for utm the convergence is 
> negative when west and positive when east of lon_0.
> If you calculate two utm positions along the same longitude east of 
> lon_0, you see that the easting reduces when going northwards.
> 
>  From this is seems that TrueNorth = GridNorth + convergence
> 
> 
> I used the same points for a lcc projection:
> 
> C:\proj-4.6.1\src>proj -vS +proj=lcc +datum=WGS84 +lat_0=52 +lon_0=3 
> +lat_1=50 +lat_2=54
> #Lambert Conformal Conic
> #       Conic, Sph&Ell
> #       lat_1= and lat_2= or lat_0
> # +proj=lcc +datum=WGS84 +lat_0=52 +lon_0=3 +lat_1=50 +lat_2=54 +ellps=WGS84
> # +towgs84=0,0,0
> 2 52
> -68634.11       472.08  <0.999392 0.999392 0.998785 0.000355057 0.999395 
> 0.999389 0.788173>
> 4 52
> 68634.11        472.08  <0.999392 0.999392 0.998785 0.000355708 0.999395 
> 0.999389 -0.788173>
> 4 52.1
> 68481.15        11591.10        <0.999393 0.999393 0.998787 0.000434814 
> 0.999397 0.99939 -0.788173>
> 
> 
> The convergence sign west and east of the origin are opposite of the 
> signs found with the utm. But what confuses me is that when i calculate 
> two positions along a longitude east of the origin, the easting reduces 
> when going northwards just like with the utm projection.
> 
>  From this is seems that TrueNorth = GridNorth - convergence
> 
> 
> I think that either the convergence for lcc or utm has the wrong sign. 
> But i can't tell which one because i don't know how convergence is defined.
> 
> 
> Any thoughts about this?
> 
> Regards,
> 
> Kees Krikke
> 
> _______________________________________________
> Proj mailing list
> Proj at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/proj
> 
>