no zejo, to by melo vyplyvat z tech dat, ale dik za tip, to "-te" by melo zabrat
čt 1. 3. 2018 v 10:27 odesílatel Antonin Orlik <antonin.or...@gmail.com> napsal: > Abych to doplnil - problém je v tom, že hranice timeline může být uvedena > jako -180 (180W) i jako 180 (180E) a je možno mezi ně dát rovnítko. > gdalwarp to bohužel vzal z té opačné strany než bylo potřeba a tomu > uzpůsobil výsledný extent - a tak mu bylo potřeb pomoci a v -te uvést místo > -180 hodnotu 179.9999999. > Nemyslím si proto, že by to byla chyba GDALu - on nemůže vedět, ze které > strany timeline chceš hranici extentu... > > Tonda > > > ------------------------------ > *From:* Jachym Cepicky <jachym.cepi...@gmail.com> > <jachym.cepi...@gmail.com> > *Subject:* [FreeGeoCZ] gdalwarp stavkuje asi na dateline > *Date:* Thursday, Mar 1, 2018 8:35 GMT > *To:* Antonin Orlik <antonin.or...@gmail.com> <antonin.or...@gmail.com> > *Cc:* Svobodná geoinformační infrastruktura <freegeocz@fsv.cvut.cz> > <freegeocz@fsv.cvut.cz> > > Vidíš, to je zajímavej typ, > > mezitím tam něco patchnul i Even https://trac.osgeo.org/osgeo/ticket/2096 ale > zatím se mi to nepodařilo rozběhnout > > J > > čt 1. 3. 2018 v 8:59 odesílatel Antonin Orlik <antonin.or...@gmail.com> > napsal: > >> Ahoj Jáchyme, >> >> problém je skutečně v dateline. Když se podíváš na souřadnice extentu, >> které jsi poslal, tak minx v lat/lon (resp. ve WGS 84 do kterého >> transformuješ) je 180d 0' 0.00"W a maxx je 179d 3'38.91"E. Rozsah >> výstupního extentu tedy začíná na -180° a končí na nějakých +179. Extent >> pak tedy obepíná celou zeměkouli a velikost výstupního rastru je pak tedy >> neuvěřitelných 8009770 x 14124 pixelů a transformace je tedy prakticky >> nemožná. Proto je nutné omezit výstupní extent pomocí parametru -te - ve >> tvém případě konkrétně takto: >> >> gdalwarp -wo SOURCE_EXTRA=1000 -wo SAMPLE_STEPS=100 -wo SAMPLE_GRID=YES >> -te 178.090847222 70.2776472222 179.99999 70.9124833333 -co >> COMPRESS=DEFLATE -co BIGTIFF=YES -s_srs +init=epsg:3411 -t_srs >> "+proj=longlat +datum=WGS84 +over +no_defs" 56_26_1_1_5m_v2.0_reg_dem.tif >> 56_26_1_1_5m_v2.0_reg_dem-warped.tif >> >> Creating output file that is 42477P x 14125L. >> >> Processing input file 56_26_1_1_5m_v2.0_reg_dem.tif. >> Using internal nodata values (e.g. -9999) for image >> 56_26_1_1_5m_v2.0_reg_dem.tif. >> Copying nodata values from source 56_26_1_1_5m_v2.0_reg_dem.tif to >> destination 56_26_1_1_5m_v2.0_reg_dem-warped.tif. >> 0...10...20...30...40...50...60...70...80...90...100 - done. >> >> Konverze trvá pár vteřin a výsledek sedí dobře... >> >> Měj se! >> Tonda >> >> ------------------------------ >> *From:* Jachym Cepicky <jachym.cepi...@gmail.com> >> <jachym.cepi...@gmail.com> >> *Subject:* [FreeGeoCZ] gdalwarp stavkuje asi na dateline >> *Date:* Tuesday, Feb 27, 2018 10:46 GMT >> *To:* FreeGeoCZ <freegeocz@fsv.cvut.cz> <freegeocz@fsv.cvut.cz> >> >> Ahoj, >> >> zkoužím warpovat data z Arcticdem zdroje, kolem dateline >> >> $ wget >> http://data.pgc.umn.edu/elev/dem/setsm/ArcticDEM/mosaic/v2.0/56_26/56_26_1_1_5m_v2.0.tar.gz >> >> $ tar xzf 56_26_1_1_5m_v2.0.tar.gz >> >> $ ls >> ... >> 56_26_1_1_5m_v2.0_reg_dem.tif >> ... >> >> $ gdalinfo 56_26_1_1_5m_v2.0_reg_dem.tif >> ... >> Upper Left (-1500000.000, 1550000.000) (179d 3'38.91"E, 70d16'39.53"N) >> Lower Left (-1500000.000, 1500000.000) (180d 0' 0.00"W, 70d35'51.24"N) >> Upper Right (-1450000.000, 1550000.000) (178d 5'27.05"E, 70d35'13.15"N) >> Lower Right (-1450000.000, 1500000.000) (179d 1'44.32"E, 70d54'44.94"N) >> .... >> >> Chci to přetransformovat do wgs84, používám snad všechny triky které jsem >> našel, ale gdalwarp prostě zůstane viset >> >> gdalwarp -wo SOURCE_EXTRA=1000 -wo SAMPLE_STEPS=100 -wo SAMPLE_GRID=YES >> -co COMPRESS=DEFLATE -co BIGTIFF=YES -s_srs +init=epsg:3411 -t_srs >> "+proj=longlat +datum=WGS84 +over +no_defs" 56_26_1_1_5m_v2.0_reg_dem.tif >> 56_26_1_1_5m_v2.0_reg_dem-warped.tif >> >> Creating output file that is 8009770P x 14124L. >> Processing input file 56_26_1_1_5m_v2.0_reg_dem.tif. >> Using internal nodata values (e.g. -9999) for image >> 56_26_1_1_5m_v2.0_reg_dem.tif. >> Copying nodata values from source 56_26_1_1_5m_v2.0_reg_dem.tif to >> destination56_26_1_1_5m_v2.0_reg_dem-warped.tif. >> 0. >> >> Zůstane viset na 0% a dál se ani nehne. Do cílového souboru i něco >> zapíše, ale je to jenom pár mega. Nemá s tím někdo zkušenosti? Dateline je >> jediné, co mě napadá. >> >> Dík >> >> J >> >> >> >> >> _______________________________________________ >> FreeGeoCZ mailing >> listfreege...@fsv.cvut.czhttp://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz >> >> >> >
_______________________________________________ FreeGeoCZ mailing list FreeGeoCZ@fsv.cvut.cz http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz
