[Proj] large difference while projecting utm to utm back and forth
Clifford J Mugnier
cjmce at lsu.edu
Wed Feb 16 12:26:43 EST 2011
For that much of a zone extension, you need to have additional terms in the classical expansion, say, to the fourteenth derivative for ±24°of the Central Meridian; - or - use the closed form with nested iterative loops as now implemented in ESRI software. Most people don't use that because of the extreme convergence angle between Grid North and True North is counter-productive for large-scale mapping. If it's not large-scale (1:50,000 or larger), why use UTM?
Clifford J. Mugnier, C.P., C.M.S.
Chief of Geodesy,
Center for GeoInformatics
Department of Civil Engineering
Patrick F. Taylor Hall 3223A
LOUISIANA STATE UNIVERSITY
Baton Rouge, LA 70803
Voice and Facsimile: (225) 578-8536 [Academic]
Voice and Facsimile: (225) 578-4578 [Research]
Cell: (225) 328-8975 [Academic & Research]
Honorary Life Member of the
Louisiana Society of Professional Surveyors
Fellow Emeritus of the ASPRS
Member of the Americas Petroleum Survey Group
From: proj-bounces at lists.maptools.org on behalf of Mikhail Tchernychev
Sent: Wed 16-Feb-11 11:13
To: PROJ.4 and general Projections Discussions
Subject: Re: [Proj] large difference while projecting utm to utm back and forth
I might be wrong but I'm afraid it is not simply possible. In my
experience with Transverse Mercator
formulas (not with cs2cs) they are ill-posed when you get far from the
In your case you have 5 zones difference, which is 30 degrees. In other
is just not defined when you get so far.
It would work if you have points at the border of the zones, say 30 and
35 is way far.
On 2/16/2011 5:19 AM, Christoph Dohmen wrote:
> Dear list,
> I know, I should not do it but I have to :-(
> Given is a point 600000 east and 8260000 north in utm zone 35.
> I have to project this location to utm zone 30.
> Doing some calculations I have to project this value back from
> zone 30 to zone 35.
> But I get differences of about 3km east and more than 200m north.
> Those are the steps I made using cs2cs:
> cs2cs +proj=utm +zone=35 +datum=WGS84 +to +proj=utm +zone=30 +datum=WGS84
> 600000 8260000
> 1451785.04 8532758.68 0.00
> and the way back:
> cs2cs +proj=utm +zone=30 +datum=WGS84 +to +proj=utm +zone=35 +datum=WGS84
> 1451785.04 8532758.68
> 597050.49 8260246.21 0.00
> Can anyone describe a way to achive the startvalues?
> Best regards
> I use proj 4.7.1 and the latest from svn.
Proj mailing list
Proj at lists.maptools.org
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Proj