EPSG:1042

Note: These formulas have been transcribed from EPSG Guidance Note #7-2. Users are encouraged to use that document rather than the text which follows as reference because limitations in the transcription will be avoided.

From the defining parameters the following constants for the projection may be calculated :

A = a(1 - e^2)^0.5 / [1 - e^2 sin^2.(latC)]
B = {1 + [e^2 * cos^4(latC) / (1 - e^2)]}^0.5
gammao = asin[sin(latC) / B]
to = tan(pi/4 + gammao/2).[(1 + e sin(latC)) / (1 - e sin(latC))]^(e.B/2) / [tan(pi/4 + latC/2)]^B
n = sin(latp)
ro = kp.A / tan(latp)

To derive the projected Southing and Westing coordinates of a point with geographical coordinates (lat, lon) the formulas for the Krovak are:

U = 2(atan{to.tan^B(lat/2 + pi/4) / [(1 + e sin(lat)) / (1 - e sin(lat))]^[e.B/2]} - pi/4)
V = B(lonO - lon) where lonO and lon must both be referenced to the same prime meridian.
T = asin[cos(alphaC).sin(U) + sin(alphaC).cos(U). cos(V)]
D = asin[cos(U).sin(V)/cos(T)]
theta = n.D
r = ro.tan^n(pi/4 + latp/2) / tan^n(T/2 + pi/4)
Xp = r.cos(theta)
Yp = r.sin(theta)
Xr = Xp – Xo
Yr = Yp – Yo
dX = C1 + C3.Xr – C4.Yr – 2.C6.Xr.Yr + C5.(Xr^2 – Yr^2) + C7.Xr.(Xr^2 – 3.Yr^2) – C8.Yr.(3.Xr^2 – Yr^2) + 4.C9.Xr.Yr.(Xr^2 – Yr^2) + C10.(Xr^4 + Yr^4 – 6.Xr^2.Yr^2)
dY = C2 + C3.Yr + C4.Xr + 2.C5.Xr.Yr + C6.(Xr^2 – Yr^2) + C8.Xr.(Xr^2 – 3.Yr^2)+ C7.Yr.(3.Xr^2 – Yr^2) – 4.C10.Xr.Yr.(Xr^2 – Yr^2) + C9.(Xr^4 + Yr^4 – 6.Xr^2.Yr^2)
Southing   X = FN + Xp – dX
Westing    Y = FE + Yp – dY

The reverse formulas to derive the latitude and longitude of a point from its Southing and Westing values are:

Xr' = (Southing – FN) – Xo
Yr' = (Westing – FE) – Yo
dX' = C1 + C3.Xr' – C4.Yr' – 2.C6.Xr'.Yr' + C5.(Xr'^2 – Yr'^2) + C7.Xr'.(Xr'^2 – 3.Yr'^2) – C8.Yr'.(3.Xr'^2 – Yr'^2) + 4.C9.Xr'.Yr'.(Xr'^2 – Yr'^2) + C10.(Xr'^4 + Yr'^4 – 6.Xr'^2.Yr'^2)
dY' = C2 + C3.Yr' + C4.Xr' + 2.C5.Xr'.Yr' + C6.(Xr'^2 – Yr'^2) + C8.Xr'.(Xr'^2 – 3.Yr'^2) 
+ C7.Yr'.(3.Xr'^2 – Yr'^2) - 4.C10.Xr'.Yr'.(Xr'^2 – Yr'^2) + C9.(Xr'^4 + Yr'^4 – 6.Xr'^2.Yr'^2)
Xp' = (Southing – FN) + dX'
Yp' = (Westing – FE) + dY'
r' = [(Yp')^2 + (Xp')^2]^(1/2)  
theta' = atan2[Yp' , Xp']   (see GN7-2 implementation notes in preface for atan2 convention)
D' = theta' / sin(latp)
T' = 2{atan[((ro / r')^(1/n)).tan(pi/4 + latp/2)] - pi/4}
U' = asin[cos(alphaC).sin(T') - sin(alphaC).cos(T').cos(D')]
V' = asin(cos(T').sin(D') / cos(U'))

Then latitude lat is found by iteration using U' as the value for lat(j-1) in the first iteration:
lat(j) = 2*(atan{tO^(-1/B) tan^(1/B).(U'/2 + pi/4).[(1 + e sin(lat(j-1)) / (1 - e sin(lat(j-1))]^(e/2)} - pi/4) 

Then
lon = lonO - V' / B where lon is referenced to the same prime meridian as lonO.