[Proj] Re: Krovak test points

[Martin Landa]

Hi all,

just one example (from [1])

1) WGS-84 -> S-JTSK
12d48'25.15992	49d27'8.14571 -> 868208.52 1095793.96

echo "12d48'25.15992  49d27'8.14571" | cs2cs +proj=latlong +datum=WGS84 \
+to +proj=krovak +ellps=bessel +no_defs \
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56

-868208.53      -1095793.57 -46.96

dy=1cm, dx=39cm (~ bessel -> wgs-84, OK)

2) Bessel (Ferro) -> S-JTSK
32d05'39.77584 50d05'11.32822 -> 742455.8799 1043259.9773

echo "32d05'39.77584 50d05'11.32822" | cs2cs +proj=latlong
+ellps=bessel +pm=ferro +to +proj=krovak +ellps=bessel +no_defs
-742455.88      -1043259.98 0.00

-> OK!

Hope it helps...

Best regards, Martin

[1] http://grass.fsv.cvut.cz/wiki/index.php/S-JTSK#Praktick.C3.A1_uk.C3.A1zka

2006/10/19, Melita Kennedy <mkennedy2 at earthlink.net>:
>
> >----------------------------------------------------------------------
> >
> >Date: Tue, 17 Oct 2006 16:41:49 -0400
> >From: Frank Warmerdam <warmerdam at pobox.com>
> >Subject: Re: [Proj] PROJ 4.5.0 beta (and Krovak!)
> >To: gerald.evenden at verizon.net,        "PROJ.4 and general Projections
> >       Discussions" <proj at lists.maptools.org>
> >
> [contents snipped]
> >
> >I don't really know any test points for Krovak or even what options
> >are required in the two implementations.
>
> While we were implementing Krovak for ESRI software, we were sent
> 7 test points. I don't know how official these are, as the source is
> unclear. Two of the points show up in another document sent by the
> ESRI distributor in  the Czech Republic. The longitudes are Greenwich-based,
> Bessel 1841 and the XY axes are south-west oriented, rather than east-north.
> The sample you posted is EN. When using the SW version, I got
>
> x = 1546506.6922  y = 1164802.1478
>
> longitude,latitude,x,y
> 16 50 59.17900, 50 12 32.44160, 1050538.6433, 568990.9967
> 12 50 23.78945, 49 50 35.78965, 1053301.0476, 859084.7000
> 22 19 38.87658, 49 50 35.78965, 1117926.6680, 180241.3620
> 22 19 38.87658, 48 38 23.46592, 1251684.1815, 184607.1409
> 12 50 23.78945, 48 38 23.46592, 1185471.9460, 880118.0569
> 16 50 59.17900, 50 12 32.44140, 1050538.6494, 568990.9974
> 17 25 51.17920, 50 15 30.70080, 1049221.1899, 527189.1948
>
> Melita Kennedy
> ESRI Product Specialist
>
> _______________________________________________
> Proj mailing list
> Proj at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/proj
>


-- 
Martin Landa <landa.martin at gmail.com> * http://gama.fsv.cvut.cz/~landa *