<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
<html>
<body text="#000000" bgcolor="#FFFFFF" link="#0000EE" vlink="#551A8B" alink="#FF0000">
&nbsp;
<br>Hello all
<p>being unfamiliar with cartographic projections, I hope I won't bother
you with too naive questions.
<p>I have to process together data of different nature and sources:
<br>&nbsp;&nbsp;&nbsp; * satellite (landsat ETM) data, geocorrected (level
1G) , delivered in UTM projection and unspecified ellipsoid -in which case
I assume WGS84
<br>&nbsp;&nbsp;&nbsp; * data given in latitude/longitude for WGS84
<br>&nbsp;&nbsp;&nbsp; * data from ground based stations, delivered in
the local co-ordinate system, in my case (former East Germany):
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Projection: Gauss Krueger in 6&deg; stripes (or Transverse Mercator)
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
projection parameters: central meridian 15E false northing 3,500,000m
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Ellipsoid: Krassowsky with Pulkovo Datum (dX=28, dY=-130 dZ=-95)
<p>then come my questions:
<br>-up to version 4.3, proj was restricted to projections only, within
the same ellipsoid/datum system. It seems that the newer version -and program
cs2cs- is able to do the job. Am I right, or is cs2cs also restricted to
one ellipsoid at a time?
<p>-to check my results I was provided with some checked remarkable points,
which co-ordinates are known in lat/long and in the local system.
<br>I was never able to find these results with proj/cs2cs
<br>Example:
<br>"ground truth": geographical, WGS84: 52.4172632N 13.6254409E becomes
y=3406605m x=5810804m on GK/Krassowsky etc
<p>The command
<br>cs2cs&nbsp; +ellps=WGS84 +proj=latlong +to +ellps=krass +proj=tmerc
+lon_0=15 +x_0=3500000
<br>gives me approximately the same figures, yet slighty different:
<br>y=3406478.97&nbsp;&nbsp; x=5810764.06 (difference is 40m in x and 126m
in y)
<p>-although this difference might be sufficiently small for my application
I would be glad to know if it is normal or much too large, in comparison
to cartography standards
<br>-the difference in x and y looks very much the same as the Pulkovo
datum. I'm not sure of the meaning of this datum: is it a shift from Krassovsky
ellipsoid or from WGS84? In the first case, how can I tell it to cs2cs,
and in the second case, why the same cs2cs command with datum specified
as +towgs84=28,-130,-95
<br>gives me absurd results (e.g. x=3500000.00&nbsp; y=10002137.50)
<p>-last but not least, wanting to check my command before bothering you,
I've asked for more info with option -v, and look at what is said about
ellipsoid:
<br>(with above cs2cs command, option -v added)
<br># ---- From Coordinate System ----
<br>#Lat/long (Geodetic)
<br>#
<br># +ellps=WGS84 +proj=latlong
<br># ---- To Coordinate System ----
<br>#Transverse Mercator
<br>#&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Cyl, Sph&amp;Ell
<br># +ellps=krass +proj=tmerc +lon_0=15 +x_0=3500000
<br>#--- following specified but NOT used
<br># +ellps=WGS84
<p>how do I have to understand the last comment?
<p>Many thanks for your help
<p>Jean-Paul Berroir
<p>PS: I'm using the proj package version 4.4.4 on a lunix i686 machine
<pre>--&nbsp;
----------------------
AIR project - INRIA - BP 105 - F78153 Le Chesnay Cedex
Tel: +33 1 3963 5703 Fax: +33 1 3963 5880
jean-paul.berroir@inria.fr</pre>
&nbsp;
</body>
</html>