Mně to přijde taky nekompletní, asi by to chtělo chtělo sehnat ten CD-ROM, o
kterém autor ve své práci píše, že je přílohou. :-) Proto jsem chtěl zkusit
tu cestu přes GDAL.
Udělal jsem bodove_pole.txt -
https://drive.google.com/file/d/0B5J67rBdf34NbVhtWTJHN1hiWUE/edit?usp=sharing
- souřadnice bodů TB, ZhB a CZEPOS v S-JTSK a WGS84
Chybu transformačního klíče lze vizualizovat tak, že se vezmou ty S-JTSK
souřadnice, natransformují se v PostGISu pomocí
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 do WGS84 jako jeden bod
úsečky a ty WGS84 souřadnice jako druhý konec úsečky.
Nebo si lze vyhrát s GRASSem:
http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Chyba_p%C5%99i_transformaci_z_WGS84_do_S-JTSK
MK
"Petr Vejsada" píše v diskusním příspěvku news:1703919.yriTVijfGK@mrnous...
Dobrá, tak jsem se mírně znemožnil ;-). Počítat potřebujeme ten vlastní
grid.
Vstupem je databáze bodů a ta se bude interpolovat.
Kromě Postgisu tedy potřebujeme ještě R a PL/R, k tomu možná něco z CRANu.
Zacházet s R celkem zvládám, i když v jiné oblasti než jsou geo-úlohy. Zdá
se
mi, že k pochopení, dostatečnému ke zprovoznění úlohy, mi chybí porozumění
všem parametrům funkce dataz_tps_raster.
Máme tam:
CREATE OR REPLACE FUNCTION dataz_tps_raster
(
tab character varying,
col character varying,
numx integer,
numy integer,
srid integer,
zoom numeric
)
RETURNS raster AS
$$
...
parametry:
tab - jméno tabulky s body
col - jméno sloupce, jakého?
numx - nevím - používá se jako parametr st_makeemptyraster
numy - nevím - dtto
srid - čeho SRID?
nemohu nalézt tyto funkce:
- dataz_create_vector
- dataz_return_value
struktura tabulky tab:
- the_geom - asi geometrie bodu, kterého? Křovákova nebo WGS84?
- 'col' - ???
Je to někomu jasnější?
--
Petr
Dne Út 13. května 2014 22:49:17, Petr Vejsada napsal(a):
Zdravím,
tato problematika mi není úplně jasná, tak se třeba zeptám blbě - co je na
tom k počítání? IMO jde o to, tu databázi, co máte k dispozici, převést do
zdrojové formy pro nad2bin, pustit na to nad2bin a pak už jen upravit
definici +proj ... a od té doby to Postgis bude umět, ne? Netuším formát
té
zdrojové formy ani nevím, jak vypadá ta databáze bodů. Je někde k
dispozici
k nahlédnutí či sosnutí?
--
Petr
Dne Út 13. května 2014 21:56:26, Martin Kokes napsal(a):
> Mám celou databázi polohového bodového pole TB, ZhB a CZEPOS v
> ETRS89(ETRF2000) <> S-JTSK, nechtěl by mi někdo, kdo se víc orientuje v
> GDAL (Jáchym Čepický?) pomoci vytvoři nový nadgrid ntv2?
> http://lists.maptools.org/pipermail/proj/2013-January/006539.html
> Docela by nám to všem pomohlo v souvislosti s propojením s RÚIAN,
> protože
> přiznejme si, sedmiprvková transformace není v pohraničí to pravé
> ořechové.
>
> viz. http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid
> http://www.kma.zcu.cz/main.php?KMAfile=./STRUCTURE/05_ebooks/04_Zaverecne_
> pr
> ace/zav_prace.php&DRC=./STRUCTURE/05_ebooks/04_Zaverecne_prace/&DRL=CZ&DR
> OF= 0&osCislo=52920
>
> Psal jsem si na to téma s p. Ježkem, zkoušel jsem kontaktovat p. Chlupa,
> ale bez výsledku. Jsem na to ochotný dát volnou výpočetní kapacitu na
> postgisu na Xeon E5 3.6 GHz s SSD.
>
> MK
>
>
>
> _______________________________________________
> Talk-cz mailing list
> [email protected]
> https://lists.openstreetmap.org/listinfo/talk-cz
_______________________________________________
Talk-cz mailing list
[email protected]
https://lists.openstreetmap.org/listinfo/talk-cz
_______________________________________________
Talk-cz mailing list
[email protected]
https://lists.openstreetmap.org/listinfo/talk-cz
