[FreeGeoCZ] Re: reprojekce do krovaka

["Jachym Cepicky" <jachym.cepicky AT gmail.com>]

gdalwarp -s_srs "+ellipsoid=WGS84 +proj=utm +zone=33" -t_srs
"epsg:2065" QBsubset.TIF QBsubset-jtsk.tif

gdalinfo QBsubset-jtsk.tif
Driver: GTiff/GeoTIFF
Size is 355, 355
Coordinate System is:
PROJCS["S-JTSK (Ferro) / Krovak",
   GEOGCS["S-JTSK (Ferro)",
       DATUM["S_JTSK_Ferro",
           SPHEROID["Bessel 1841",6377397.155,299.1528128000008,
               AUTHORITY["EPSG","7004"]],
           AUTHORITY["EPSG","6818"]],
       PRIMEM["Ferro",-17.66666666666667],
       UNIT["degree",0.0174532925199433],
       AUTHORITY["EPSG","4818"]],
   UNIT["metre",1,
       AUTHORITY["EPSG","9001"]],
   AUTHORITY["EPSG","2065"]]
Origin = (-464078.173189933761023,-1131377.828972819726914)
Pixel Size = (2.398827733050111,-2.398827733050111)
Metadata:
 AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  ( -464078.173,-1131377.829)
Lower Left  ( -464078.173,-1132229.413)
Upper Right ( -463226.589,-1131377.829)
Lower Right ( -463226.589,-1132229.413)
Center      ( -463652.381,-1131803.621)
Band 1 Block=355x11 Type=UInt16, ColorInterp=Gray
Band 2 Block=355x11 Type=UInt16, ColorInterp=Undefined
Band 3 Block=355x11 Type=UInt16, ColorInterp=Undefined
Band 4 Block=355x11 Type=UInt16, ColorInterp=Undefined

to samy dostanes po
gdalwarp -s_srs "+ellipsoid=WGS84 +proj=utm +zone=33" -t_srs
"+proj=krovak +ellipsoid=bessel" QBsubset.TIF QBsubset-jtsk.tif

problem je v tom, jestli to dobře chapu, ze aktualni stabilni gdal
neumi ulozit "towgs84" atribut do rastroveho souboru pro krovaka. bude
to snad opraveny v gdal 1.5.0

gdalwarp -s_srs "+ellipsoid=WGS84 +proj=utm +zone=33" -t_srs
"+proj=krovak +ellipsoid=bessel
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56" QBsubset.TIF
QBsubset-jtsk.tif

snad by melo zabrat, provest transformaci v GRASSu:

# zalozit location utm33, naimportovat snimek
GRASS6.3.cvs utm33/jachym: tmp > g.region -g
n=5496931.2
s=5496170.4
w=746558.4
e=747319.2
nsres=2.4
ewres=2.4
rows=317
cols=317
cells=100489

GRASS6.3.cvs utm33/jachym: tmp >g.proj -w
PROJCS["UTM Zone 33, Northern Hemisphere",
   GEOGCS["wgs84",
       DATUM["unknown",
           SPHEROID["WGS_1984",6378137,298.257223563]],
       PRIMEM["Greenwich",0],
       UNIT["degree",0.0174532925199433]],
   PROJECTION["Transverse_Mercator"],
   PARAMETER["latitude_of_origin",0],
   PARAMETER["central_meridian",15],
   PARAMETER["scale_factor",0.9996],
   PARAMETER["false_easting",500000],
   PARAMETER["false_northing",0],
   UNIT["meter",1]]

# zjisteni hranicnich souradnic pro location jtsk:
GRASS6.3.cvs utm33/jachym: tmp > echo "746558.4 5496170.4"|cs2cs
+proj=utm +ellipsoid=WGS84 +zone=33 +to +proj=krovak +ellipsoid=bessel
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56
-463951.58      -1132011.59 -42.36

GRASS6.3.cvs utm33/jachym: tmp > echo "747319.2 5496931.2"|cs2cs
+proj=utm +ellipsoid=WGS84 +zone=33 +to +proj=krovak +ellipsoid=bessel
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56
-463099.38      -1131355.44 -42.34


# zalozit location jtsk:
# coord system: krovak
# specify datum: yes
# datum: hermannskogel
# transformation parameters: 3
GRASS6.3.cvs jtskutm/jachym: tmp >  g.region n=-1131355.44
s=-1132011.59 e=-463099.38 w=-463951.58

# doporučuju odzoomovat, rastry jsou pootočený a nevejdou se

GRASS6.3.cvs jtskutm/jachym: tmp > g.proj -w
PROJCS["Krovak",
   GEOGCS["bessel",
       DATUM["Militar_Geographische_Institut",
           SPHEROID["Bessel_1841",6377397.155,299.1528128],
           TOWGS84[570.8,85.7,462.8,4.998,1.587,5.261,3.56]],
       PRIMEM["Greenwich",0],
       UNIT["degree",0.0174532925199433]],
   PROJECTION["Krovak"],
   PARAMETER["latitude_of_center",0],
   PARAMETER["longitude_of_center",0],
   PARAMETER["azimuth",0],
   PARAMETER["pseudo_standard_parallel_1",0],
   PARAMETER["scale_factor",1],
   PARAMETER["false_easting",0],
   PARAMETER["false_northing",0],
   UNIT["meter",1]]

GRASS6.3.cvs jtskutm/jachym: tmp > for i in `seq 4`; do r.proj
in=qb.$i out=qb.$i location=utm33 --o; done


mno a závěrečný test:

GRASS6.3.cvs jtskutm/jachym: tmp > r.out.gdal in=qb.4 out=qb4.tif type=Int16
ERROR 6: SetColorTable() only supported for Byte or UInt16 bands in
TIFF format. 100%
r.out.gdal complete.
GRASS6.3.cvs jtskutm/jachym: tmp > gdalinfo qb4.tif
Driver: GTiff/GeoTIFF
Size is 1987, 1531
Coordinate System is `'
Origin = (-464519.016431930009276,-1130918.535061219939962)
Pixel Size = (1.000234741781585,-1.000228658536903)
Metadata:
 AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  ( -464519.016,-1130918.535)
Lower Left  ( -464519.016,-1132449.885)
Upper Right ( -462531.550,-1130918.535)
Lower Right ( -462531.550,-1132449.885)
Center      ( -463525.283,-1131684.210)
Band 1 Block=1987x2 Type=Int16, ColorInterp=Gray
 NoData Value=-32768
 Metadata:
   COLOR_TABLE_RULES_COUNT=1
   COLOR_TABLE_RULE_RGB_0=2.250000e+02 2.007000e+03 0 0 0 255 255 255
GRASS6.3.cvs jtskutm/jachym: tmp >

qgis qb4.tif QBsubset.TIF

a skutečně - je to posunytý

jdu psát bugreport

jestli to dobře chápu, tak GRASS transformaci "umí" i s
transformačními parametry, kdežto gdalwarp ne

jachym


12.4.07, HÁJEK Filip Ing. <Hajek.Filip AT uhul.cz>:
> Jasne, my to muzem overit (aspon do urcite miry) na vektorech LHP - ackoli 
> tyvsou stejne asi zdrble.
> Snimek je podle me ve WGS 84 UTM (33N).prikladam textak, co byl u origos QB.
>
> Dekuju. F.
>
> ::::::::::::::::::::::::::::::::::::::::::::::::::
> Ing. Filip Hájek
> zpracování dat DPZ
> ÚHÚL Brandýs nad Labem, pobočka Olomouc
> +420724935916
> -----Original Message-----
> From: Jachym Cepicky [mailto:jachym.cepicky AT gmail.com]
> Sent: Wednesday, April 11, 2007 5:20 PM
> To: HÁJEK Filip Ing.
> Subject: Re: reprojekce do krovaka
>
> copak já
>
> šlo mi spíš o to, jesli máš někde v mapě změřený souřadnice v alespoň
> jednom ze systémů, aby se to potom dalo pomocí cs2cs převést do
> druhého (to funguje 100%) a pak na tom vyzkoušet, jestli gdalwarp dělá
> co má
>
> no, GPSka mi dorazí tenhle pátek, já dorazím tuhle sobotu, pak se
> nejspíš rozkrájím, a pak tak za měsíc bych moh mít trochu času na
> výpravičku na brdy, naměření pár bodů tou gpskou (ať se trochu
> zahřeje) a oveření..
>
> jinak bych vám doporučil, pokud trochu můžete, zůstaňte ve wgs84,
> ušetříte si tím spoustu problémů
>
> jachym@mellon:~/data/cr$ gdalinfo QBsubset.TIF
> Driver: GTiff/GeoTIFF
> Size is 317, 317
> Coordinate System is `'
> Origin = (746558.400000000023283,5496931.200000000186265)
> Pixel Size = (2.400000000000000,-2.400000000000000)
> Corner Coordinates:
> Upper Left  (  746558.400, 5496931.200)
> Lower Left  (  746558.400, 5496170.400)
> Upper Right (  747319.200, 5496931.200)
> Lower Right (  747319.200, 5496170.400)
> Center      (  746938.800, 5496550.800)
> Band 1 Block=317x12 Type=UInt16, ColorInterp=Gray
> Band 2 Block=317x12 Type=UInt16, ColorInterp=Undefined
> Band 3 Block=317x12 Type=UInt16, ColorInterp=Undefined
> Band 4 Block=317x12 Type=UInt16, ColorInterp=Undefined
> jachym@mellon:~/data/cr$
>
> hmm, souřadný systém nedefinován můžeš mi prosím napsat, v jakym souř.
> systému to je? tohle je na mě trochu moc, hádat souř. systém ze
> souřadnic..
>
> dík
>
> j
>
>
>
> 11.4.07, HÁJEK Filip Ing. <Hajek.Filip AT uhul.cz>:
> > Nazdar,
> > dekuju za odpoved. posilam subset QucikBird, z nasi tesovaci oblasti v 
> Beskydech. Zkus!
> >
> > F.
> >
> > ::::::::::::::::::::::::::::::::::::::::::::::::::
> > Ing. Filip Hájek
> > zpracování dat DPZ
> > ÚHÚL Brandýs nad Labem, pobočka Olomouc
> > +420724935916
> >
> > -----Original Message-----
> > From: Jachym Cepicky [mailto:jachym.cepicky AT gmail.com]
> > Sent: Wednesday, April 11, 2007 4:47 PM
> > To: HÁJEK Filip Ing.
> > Subject: Re: reprojekce do krovaka
> >
> > ahoj,
> >
> > jsem rád, že jste to ještě nevzdali :-)
> >
> >
> > co se gdalwarp týče, je vždycky dobrý mít nějaký referenční data.
> > cs2cs |v.in.ascii je pro začátek dobrej krok.
> >
> > jak psal martin. od proj 4.5.0 by to mělo chodit, jak potřeba..
> >
> > co se gdalwarp týče, mělo byt to chodit asi takhle:
> >
> > gdalwarp  -s_srs +init="epsg:4326" -t_srs +init="epsg:2065"
> > +towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 soubor.tif
> > soubor-jtsk.tif
> >
> > no, ale jsem trochu v rozpacich, protoze ten parametr "towgs84"
> > nevypada, ze by delal to, co bych cekal.
> >
> > mas prosim nejaka data, na kterych bys to mohl overit?
> >
> > j
> > 11.4.07, Martin Landa <landa.martin AT gmail.com>:
> > > Ahoj,
> > >
> > > mrkni na
> > >
> > > http://grass.fsv.cvut.cz/wiki/index.php/S-JTSK
> > >
> > > od verze PROJ.4 4.5.0 by nemel byt problem.
> > >
> > > [CC Jachymovi]
> > >
> > > Martin
> > >
> > > PS: Lepsi napsat do konference, vic hlav vic vi:-)
> > >
> > > http://grass.fsv.cvut.cz/wiki/index.php/E-konference_FreeGeoCZ
> > >
> > > 11.4.07, HÁJEK Filip Ing. <Hajek.Filip AT uhul.cz>:
> > > >
> > > >
> > > >
> > > >
> > > > Ahoj Martine,
> > > >
> > > > byl jsem onehdy na kurzu GRASS co mel Jachym pro lidi z UHULu, pak jsem 
> te
> > > > videl prednaset na Gis Ostrava o QGisu. Pracuju na UHUL teprve cca 2 
> mesice
> > > > postupne se pokousim rozjet nejakou praci s druzicovymi daty. Standartni
> > > > nastroj je tady Topol, ja se však postupne ucim s ruznyma OSGIG 
> aplikacema
> > > > na praci s rastry - hlavne OpenEV a OSSIM ImageLinker (zatím pod 
> Windows).
> > > > OpenEV je skvely, poziva tez vyhod GDAL a Proj, ad...
> > > >
> > > >
> > > >
> > > > Ted resim, jestli uz je aktualne opravdu mozne prevest rastry ve WGS84 
> na
> > > > S-JTSK? Nova verze Proj by to mela podporovat, ze?
> > > >
> > > > Muzu to udelat prikazem gdalwarp, nebo doporucis něco jineho?
> > > >
> > > >
> > > >
> > > > Dekuju moc za odpoved.
> > > >
> > > > Filip
> > > >
> > > >
> > > >
> > > > ::::::::::::::::::::::::::::::::::::::::::::::::::
> > > >
> > > > Ing. Filip Hájek
> > > >
> > > > zpracování dat DPZ
> > > >
> > > > ÚHÚL Brandýs nad Labem, pobočka Olomouc
> > > >
> > > > +420724935916
> > > >
> > > >
> > >
> > >
> > > --
> > > Martin Landa <landa.martin AT gmail.com> * http://gama.fsv.cvut.cz/~landa *
> > >
> >
> >
> > --
> > Jachym Cepicky
> > e-mail: jachym.cepicky gmail com
> > URL: http://les-ejk.cz
> > GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub
> >
> >
>
>
> --
> Jachym Cepicky
> e-mail: jachym.cepicky gmail com
> URL: http://les-ejk.cz
> GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub


--
Jachym Cepicky
e-mail: jachym.cepicky gmail com
URL: http://les-ejk.cz
GPG: http://www.les-ejk.cz/pgp/jachym_cepicky-gpg.pub