[Talk-cz] Nový S-JTSK grid, kdo pomůže?

[Petr Vejsada]

Ahoj,

ještě se k tomu vracím.

Nejprve převeďme tvé výsledky na počítačovější formu. Jestli počítám správně, 
tak tvé výsledky pro transformaci bez gridu jsou:

18.7939444444  49.5737863889

a pro transformaci s gridem

18.793978333333 49.5737916667


což je rozdíl:

select 
import.distance_meters(st_setsrid(st_makepoint(18.7939444444,49.5737863889),4326),st_setsrid(st_makepoint(18.793978333333,49.5737916667),4326));

 distance_meters
-----------------
     2.520413388
(1 řádka)

to je trochu moc, ne? 2.5 metru?

A teď moje výsledky:

bez gridu:

select st_astext(st_transform(st_setsrid(st_makepoint(-436231.78, 
-1133608.56),5514),4326));

                st_astext
------------------------------------------
 POINT(18.7939443539245 49.5737863847451)
(1 řádka)


s gridem:

select st_astext(st_transform(st_setsrid(st_makepoint(-436231.78, 
-1133608.56),999),4326));
                st_astext
------------------------------------------
 POINT(18.7939504934423 49.5737917877839)
(1 řádka)

Rozdíl:

select 
import.distance_meters(st_setsrid(st_makepoint(18.7939443539245,49.5737863847451),4326),st_setsrid(st_makepoint(18.7939504934423,49.5737917877839),4326));

 distance_meters
-----------------
     0.747197172
(1 řádka)

To už odpovídá deklarovaným 0.7m

Seknul jsem se při převodu těch tvých výsledků? Přepočítával jsem to 2x, 
stupně+(minuty/60)+(vteřiny/3600)

Dejme tomu, že je posun při transformaci s gridem správným směrem a že je tedy 
výsledek přesnější než bez gridu. Když si dáme v JOSM přes sebe průhlednou KM 
a k tomu moje vrstvy budov, proč se tedy lépe kryjí růžové budovy (t.j. bez 
gridu) než budovy modré (t.j. s gridem) ???Růžové se lépe kryjí i s tím, co je 
v OSM. Znamená to, že máme vše v OSM špatně? A hlavně - co s tím dál?

Poznámka: funkce distance_meters vrací:

st_distance(st_transform(geom_a,4326)::geography,st_transform(geom_b,4326)::geography)


Dne Út 20. května 2014 10:23:32, hanoj napsal(a):

> echo "-436231.78 -1133608.56" | cs2cs +init=esri:102067
> +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to +init=epsg:4326
> 18d47'38.2"E    49d34'25.631"N 41.921
> 
> nyni transformujme gridem[4]:
> 
> echo "-436231.78 -1133608.56" | cs2cs +proj=krovak +ellps=bessel
> +nadgrids=jezek_czech08.llb +to +proj=longlat +datum=WGS84
> 18d47'38.222"E    49d34'25.65"N 0.000