Fernando, The GCTP.FOR is a Fortran77 source file (free from the "USGS.gov" websites) for all of the map projections used by the U.S. Geological Survey. There are two data files associated with the source code that are included also. The math was documented by John P. Snyder (now deceased) in "Map Projections Used by the U.S. Geological Survey" Bulletin 1532, and later revised as "Map Projections - A Working Manual" Bulletin 1535. Bulletin 1535 is better because it has more projections and Mr. Snyder referenced me (ha, ha)! :-) Anyway, Dr. Atef Elassal (now retired), then translated 1532 into Fortran for the USGS. That is GCTP - the General Cartographic Transformation Package which is specifically for cartographic applications within the United States. That's what the data files are for. GCTP is absolutely perfect for what it was intended for - INSIDE THE UNITED STATES OF AMERICA ONLY ! ! ! ! ! Many, many commercial software packages worldwide use this as the basic foundation for their coordinate transformation engine. I give it away to my students as an example of "how not to do it." This is essentially useless for geodetic applications outside of the United States. It can oftentimes be used for cartographic applications outside of the U.S. IF AND ONLY IF the computational accuracy (and precision) is not needed for mapping at scales larger than 1:24,000!!!!!!! If you are going to use this for a NON-geodetic application, this will do just fine. If you are doing geodesy, do not touch this code! ------------------------------------------------------ The ellipsoidal case of the Transverse Mercator was cooked up by Heinrich Lambert in the middle 1700's. It was a mathematical curiousity that was useless for practical applications until the City of Hannover asked Professor Carl Freiderich Gauss to do a geodetic survey of the city in preparation for a new set of accurate tax maps. There are two things you cannot avoid in life; those are death and taxes. Most all geodetic research has been funded (since the late 1700's) for either tax mapping purposes or military purposes looking for better and more efficient ways of killing people ... Anyway, Gauss worked up an expansion of Lambert's formulae that his Ph.D. students could follow in doing the grunt work of adjusting the Hannover Triangulation Net on the Gauss-Conformal Transverse Mercator. Years later, a Prussian Artillery Office named Schreiber used a simplified form of the Gauss-Conformal Transverse Mercator that was a specific truncation called the Gauss-Schreiber Transverse Mercator. Another Prussian Artillery Officer named Krüger came up with a more elaborate expansion of the infinite series. Yup, it is called the Gauss-Krüger Transverse Mercator. In the 1920's or 1930's an Italian Professor in Italy came up with a local version for the Instituto Geografico Militare, and his name was Boaga. Yup, the Italians use the Gauss-Boaga Transverse Mercator. And so on and so forth for ALL the ellipsoidal projections used for Grids on topographic maps. When looking at geodetic accuracy and computational precision at the sub-millimeter level TO THE MULTI-METER LEVEL for coordinates many degrees east or west of the central meridian, the specific truncation of a Transverse Mercator makes a big difference. Doing foreign work for bazillion-dollar exploration, drilling, and production for oil wells in specific countries? Pay attention to your math. If you are doing UTM or DHG (Deutches Herres Gitter) within a plus or minus 3 degree longitude distance from the Central Meridian, it will do fine. Diddling with some X,Y coordinates for a Ph. D. dissertation? Unless your Major Professor is a geodesist or mathematical cartographer, they won't even know the difference. Cliff -- Clifford J. Mugnier ([EMAIL PROTECTED]) The Topographic Engineering Laboratory Department of Civil and Environmental Engineering UNIVERSITY OF NEW ORLEANS New Orleans, Louisiana 70148 Voice and Facsimile: (504) 280-7095 ----------------------------------------------------------------------- Fernando wrote: > > Hello Cliff > > I am doctorate student that needs to program a convert from Gauss > Krüger to lat. long and back. The ideal solution would be to get a > Fortran code for this, but maybe that is too much luck, so I would be > happy with any hint you can give me (if I do not have to buy any > software or module, even better). > > I have already spent a lot of time looking for it in Internet and I have > found nothing, except your name in 'users.netonecom.net' (1998). > > And I have another question that confuses me (I am a beginner on this): > Is Gauss-Krüger the same Transverse Mercator, or there is an important > difference? > > Thanks in advance for any help, > > Fernando ---------------------------------------------------------------------- To unsubscribe from this list, send e-mail to [EMAIL PROTECTED] and put "unsubscribe MAPINFO-L" in the message body, or contact [EMAIL PROTECTED]
MI Flavors of Transverse Mercator
Cliff Mugnier - University of New Orleans Thu, 22 Jul 1999 13:21:15 -0700
- MI RE: Flavours of Transverse Me... Cliff Mugnier - University of New Orleans
- MI RE: Flavours of Transver... Malcolm Jones