Krovak Projection

Prof. Ing. Bohuslav Veverka, DrSc.

Czech Technical University - Department of Mapping and Cartography

Thakurova 7, 166 29 Prague 6,

Mobil 607 860 440 Fax:+420 2 71962433 E-mail: VEVERKA.BOHUSLAV@SEZNAM.CZ

MOTTO : there are many information barriers among the cartographically oriented people. One of them is a disability to use various coordinate systems and map projections of other countries.

This text is to geting off on this site

a) Krovak.pdf (72 kB)

b) Krovak.zip (Microsoft Word document packed by WinZip, 166 KB)

1. CZECH COORDINATE SYSTEM S-JTSK

S-JTSK : plane geodetic grid coordinate system called in Czech "System jednotne trigonometricke site katastralni - The System of the Unified Czech/Slovak Trigonometric Cadastral Net".

l.l History

Scale, location and orientation of the S-JTSK on the surface of the Bessel ellipsoid was derived from the results of historical Austro/Hungarian military surveys in the years 1862-98. There are 42 identical points within the Czech territory used for transformation computations. Astronomical orientation was measured only on the Hermannskogel, a trigonometric point in Austria, scale factor was derived from the length geodetic basis in Josefov.

1.2 Use, Applications

Krovak projection and national grid S-JTSK were adopted on the territory of the Czech and Slovak republics (former Czechoslovakia) in 1927. The use of this system in all geodetic surveying and cartographic activities (state mapping) is in accordance with the Act No. 200/94 Dig. of the Czech State Law.

State map series of cadastral and civil topographic maps 1:500 - 1 :200 000 are based on the JTSK System (S-JTSK). There is also a former project of the State Information System of the Czech Republic, where S-JTSK is used for territorial localisation of geographical information subjects.

2. KROVAK PROJECTION AND COORDINATE SYSTEM

While terms projection and coordinate system are often used interchangeably, they do not have the same meaning.

The map projection is a systematic representation of the surface of the Earth, to plane. To the identify the location of points on the Earth, a graticule or network of longitude λ (in the U.S., the abbreviations E - East, or W - West are used) and latitude φ (see N - North or S - South) lines have been superimposed on the surface. They are commonly referred to as meridians and parallels respectively. A projection is defined by an equation or mostly a set of equations containing several parameters where to each of those parameters has been assigned a specific value. The result is a coordinate system.

The Conformal Oblique Conic Projection of Czechoslovakia was prepared by survey engineer Josef Krovak, in the year 1922 for construction of cadastral maps and topographical maps of medium scales for the civil geodetic service of Czechoslovakia.

The projection is based on the Bessel ellipsoid of 1841, which is widely adopted in Central Europe. Longitude λ is assumed from prime Ferro meridian (an island of the Canaries), not from Greenwich meridian! The projection is conformal, so that in the projection plane meridians and parallels intersect at right angles. Angular distance between Ferro and Greenwich is given as

λ greenwich = λferro - 17º 39' 59.7354"

but round value 17º 40' was used for topographic mapping in Czechia, Slovakia and Austria.

The origin of the S-JTSK coordinate system is located in Finnish Basin, X-axis normally coincides with the meridian 42º 30' from Ferro with increasing to the South. Y-axis is perpendicular to the X-axis and increasing to the West.

2.1 The Conformal Oblique Conic Projection of former Czechoslovakia

This projection consists of four steps.

a) Conformal projection of Bessel ellipsoid to Gauss sphere

Projection ellipsoid - sphere is given by relation Tl

Tl : (φ, λferro) Þ  (U,V)

where input values are:

φ: latitude on the Bessel ellipsoid

λferro : longitude on the Bessel ellipsoid (from Ferro)

and output values

U : latitude on the Gauss spheroid

V : longitude on the Gauss spheroid

these general formulas may be written as follows :

U = f (φ)

V = g (λferro)

but because we use a conic projection, this relation may be inverted to

tg (0.5 U + 45° ) = K {tga (0.5 φ + 45° ) * ((1-e sin φ)/(1+e sin φ )) a e/2}

and for value V

V = a * λferro

with constants

K = 1.003419164

a = 1.000597498372

e = 0.0816968303;

b) Change from geographic to cartographic coordinates on the Gauss sphere

Purpose : change from north pole (Sp) to cartographic pole (Q) with better respecting of the borderline of the former Czechoslovak Republic. There is no change of Gauss sphere surface. This transformation is given as relation T2

T2 : (U,V) Þ  (S,D)

where spherical triangle should be solved with the values

S : cartographic latitude on the Gauss sphere (from pole Q)

D : cartographic longitude on the Gauss sphere (from pole Q)

where so called cartographic pole Q (φQ , λQ) on the Bessel ellipsoid has coordinates

φQ = 59° 45' 27"

λQ = 42° 30' 00"

and on the Gauss sphere as the pole Q (UQ ,VQ), with values

U = 59° 42' 42.69689"

V = 42° 31' 31.41725"

Transformation relation T2 is determined by formulas

sin(S)=cos(90°-UQ)*sin(U)+sin(90°-UQ)*cos(U)*cos(VQ -V)

sin(D) = { cos(U) / cos(S) } * sin (VQ - V)

c) Conformal projection of Gauss sphere to oblique tangent cone

The aim is wrapping an oblique cone around the sphere, so that it touches the Gauss sphere surface along a standard cartographic parallel S = 78° 30' where no distance distortion is assumed. This projection is given as relation T3

T3 : (S,D) Þ (r , e )

where

r = r o {tg (0.5 So + 45°) / tg (0.5 S + 45°)}n

and

e = n * D

with constants

r o = 1 298 039.0046 m

So = 78° 30'

n = sin ( So) = 0.9799247046

where polar coordinate e is an angular distance and second polar coordinate r is a radiusvector from the cartographic pole.

d) Projection of oblique cone to grid plane S-JTSK

Polar conic coordinates (r ,e ) are transformed into rectangular coordinates (X,Y) using the relation T4. As it was mentioned in previous text, the origin of the grid X,Y system called S-JTSK is a projection of the cartographic pole Q on plane. This point lies near Estonian Tallin in the Finnish Basin.

The X axis coinciding with the meridian of longitude λferro = 42° 30´ increasing to the South. The Y axis is perpendicular increasing to the West. The whole Czechoslovak territory is situated in the first quadrant with only positive coordinates.

Final relation T4 between polar and rectangular (grid S-JTSK) is given as

T4 = (r ,e ) Þ (X,Y)

where

X = r * cos e

Y = r * sin e

Fig 1 Krovak projection of the Czechoslovak Republic - spheroid situation

3. TEST DATA

FI = 48° 07' 46.2973"

LAMBDAF= 35° 42' 35.2147"

 

U = 48° 05' 29.7061 "

V = 35° 43' 52.0262"

 

S = 77° 43' 29.8385"

D = 21° 49' 09.0823"

 

RO = 1 384 345.283 m

Y = 504 691.675 m

EPSILON = 21° 22' 52.1863"

 

Y = 504 691.675 m

X = 1 289 068.724 m

 

REREFENCES

[1] SRNKA, E. (1986): Mathematical Cartography. Textbook. Military Academy, Brno.

[2]CIMBÁLNÍK, M. - VEVERKA, B. (1992): Transformations between national coordinate systems in the Czech state map series. Geodeticky a kartograficky obzor. 38/80,1992,cislo 10, s.204-208.

[3] BUCHAR, P. - HOJOVEC, V. (1996): Mathematical Cartography. Textbook. Czech Technical University. Praha.

[4] VEVERKA, B. (2001): Topographic and Thematic Cartography. Textbook. Czech Technical University. Praha.[5] Datums, Ellipsoids, Grids and Grid Reference System. Defence Mapping Agency. Technical Manual 8358.1. Fairfax, Virginia, USA.

[6] Geodeticke referencni systemy v Ceske republice. VUGTK Zdiby a VZU Praha. Praha 1998.

[7] Veverka, B. (2000): MATKART 2000. User's manual. GeoSoft. Praha.

 

You can use follow reference:

http://homel.vsb.cz/~hom50/GEOINFRM/KROVAK/KROVAK.HTM