Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline
Antonin Orlik <[email protected]> Thu, 1 Mar 2018 10:27:27 +0100
| Newsgroups | gmane.comp.gis.grass.cz |
|---|---|
| Message-ID | <[email protected]> |
This is a multi-part message in MIME format. --===============0531825127== Content-Type: multipart/alternative; boundary="------------965FDA83C9D8144EC48D8F9D" Content-Language: cs-CZ This is a multi-part message in MIME format. --------------965FDA83C9D8144EC48D8F9D Content-Type: text/plain; charset=utf-8; format=flowed Content-Transfer-Encoding: 8bit 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 <[email protected]> *Subject:* [FreeGeoCZ] gdalwarp stavkuje asi na dateline *Date:* Thursday, Mar 1, 2018 8:35 GMT *To:* Antonin Orlik <[email protected]> *Cc:* Svobodná geoinformační infrastruktura <[email protected]> > 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 <[email protected] > <mailto:[email protected]>> 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 <[email protected]> > <mailto:[email protected]> > *Subject:* [FreeGeoCZ] gdalwarp stavkuje asi na dateline > *Date:* Tuesday, Feb 27, 2018 10:46 GMT > *To:* FreeGeoCZ <[email protected]> <mailto:[email protected]> > >> 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 list >> [email protected] <mailto:[email protected]> >> http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz > --------------965FDA83C9D8144EC48D8F9D Content-Type: text/html; charset=utf-8 Content-Transfer-Encoding: 8bit <html> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> </head> <body text="#000000" bgcolor="#FFFFFF"> 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.<br> Nemyslím si proto, že by to byla chyba GDALu - on nemůže vedět, ze které strany timeline chceš hranici extentu...<br> <br> <div class="moz-cite-prefix"> <div id="rwhMsgHeader">Tonda<br> <br> <br> <hr id="rwhMsgHdrDivider" style="border:0;border-top:1px solid #B5C4DF;padding:0;margin:10px 0 5px 0;width:100%;"><span style="margin: -1.3px 0 0 0 !important;"><font style="font: 13px Tahoma !important; color: #000000 !important;" face="Tahoma" color="#000000"><b>From:</b> Jachym Cepicky <a class="moz-txt-link-rfc2396E" href="mailto:[email protected]"><[email protected]></a></font></span><br> <span style="margin: -1.3px 0 0 0 !important;"><font style="font: 13px Tahoma !important; color: #000000 !important;" face="Tahoma" color="#000000"><b>Subject:</b> [FreeGeoCZ] gdalwarp stavkuje asi na dateline</font></span><br> <span style="margin: -1.3px 0 0 0 !important;"><font style="font: 13px Tahoma !important; color: #000000 !important;" face="Tahoma" color="#000000"><b>Date:</b> Thursday, Mar 1, 2018 8:35 GMT</font></span><br> <span style="margin: -1.3px 0 0 0 !important;"><font style="font: 13px Tahoma !important; color: #000000 !important;" face="Tahoma" color="#000000"><b>To:</b> Antonin Orlik <a class="moz-txt-link-rfc2396E" href="mailto:[email protected]"><[email protected]></a></font></span><br> <span style="margin: -1.3px 0 0 0 !important;"><font style="font: 13px Tahoma !important; color: #000000 !important;" face="Tahoma" color="#000000"><b>Cc:</b> Svobodná geoinformační infrastruktura <a class="moz-txt-link-rfc2396E" href="mailto:[email protected]"><[email protected]></a></font></span><br> <br> </div> </div> <blockquote type="cite" cite="mid:CAAZUH4F6Ppu0nKJmqUCPqRWZtxN76f7rkZgF2w6DYBY64yPSQg@mail.gmail.com" style="border:none !important; margin-left:0px !important; margin-right:0px !important; margin-top:0px !important; padding-left:0px !important; padding-right:0px !important"> <div dir="ltr">Vidíš, to je zajímavej typ, <br> <br> mezitím tam něco patchnul i Even <a href="https://trac.osgeo.org/osgeo/ticket/2096" moz-do-not-send="true">https://trac.osgeo.org/osgeo/ticket/2096</a> ale zatím se mi to nepodařilo rozběhnout <div><br> </div> <div>J</div> </div> <br> <div class="gmail_quote"> <div dir="ltr">čt 1. 3. 2018 v 8:59 odesílatel Antonin Orlik <<a href="mailto:[email protected]" moz-do-not-send="true">[email protected]</a>> napsal:<br> </div> <blockquote class="gmail_quote" style="border:none !important; margin-left:0px !important; margin-right:0px !important; margin-top:0px !important; padding-left:0px !important; padding-right:0px !important"> <div text="#000000" bgcolor="#FFFFFF"> Ahoj Jáchyme,<br> <br> 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:<br> <br> 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<br> <br> Creating output file that is 42477P x 14125L.</div> <div text="#000000" bgcolor="#FFFFFF"><br> Processing input file 56_26_1_1_5m_v2.0_reg_dem.tif.<br> Using internal nodata values (e.g. -9999) for image 56_26_1_1_5m_v2.0_reg_dem.tif.<br> 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.<br> </div> <div text="#000000" bgcolor="#FFFFFF"> 0...10...20...30...40...50...60...70...80...90...100 - done.<br> <br> Konverze trvá pár vteřin a výsledek sedí dobře...<br> <br> Měj se!<br> Tonda <div class="m_-7111605648087380987moz-cite-prefix"> <div id="m_-7111605648087380987rwhMsgHeader"><br> <hr id="m_-7111605648087380987rwhMsgHdrDivider" style="border:0;border-top:1px solid #b5c4df;padding:0;margin:10px 0 5px 0;width:100%"><span><font style="font:13px Tahoma!important;color:#000000!important" face="Tahoma" color="#000000"><b>From:</b> Jachym Cepicky <a class="m_-7111605648087380987moz-txt-link-rfc2396E" href="mailto:[email protected]" target="_blank" moz-do-not-send="true"><[email protected]></a></font></span><br> <span><font style="font:13px Tahoma!important;color:#000000!important" face="Tahoma" color="#000000"><b>Subject:</b> [FreeGeoCZ] gdalwarp stavkuje asi na dateline</font></span><br> <span><font style="font:13px Tahoma!important;color:#000000!important" face="Tahoma" color="#000000"><b>Date:</b> Tuesday, Feb 27, 2018 10:46 GMT</font></span><br> <span><font style="font:13px Tahoma!important;color:#000000!important" face="Tahoma" color="#000000"><b>To:</b> FreeGeoCZ <a class="m_-7111605648087380987moz-txt-link-rfc2396E" href="mailto:[email protected]" target="_blank" moz-do-not-send="true"><[email protected]></a></font></span><br> <br> </div> </div> </div> <div text="#000000" bgcolor="#FFFFFF"> <blockquote type="cite" style="border:none !important; margin-left:0px !important; margin-right:0px !important; margin-top:0px !important; padding-left:0px !important; padding-right:0px !important"> <div dir="ltr">Ahoj, <div><br> </div> <div>zkoužím warpovat data z Arcticdem zdroje, kolem dateline</div> <div><br> </div> <div> <div>$ wget <a href="http://data.pgc.umn.edu/elev/dem/setsm/ArcticDEM/mosaic/v2.0/56_26/56_26_1_1_5m_v2.0.tar.gz" target="_blank" moz-do-not-send="true">http://data.pgc.umn.edu/elev/dem/setsm/ArcticDEM/mosaic/v2.0/56_26/56_26_1_1_5m_v2.0.tar.gz</a></div> <div><br> </div> <div>$ tar xzf 56_26_1_1_5m_v2.0.tar.gz</div> <div><br> </div> <div>$ ls </div> <div>...</div> <div>56_26_1_1_5m_v2.0_reg_dem.tif<br> </div> <div>...</div> <div><br> </div> <div>$ gdalinfo 56_26_1_1_5m_v2.0_reg_dem.tif</div> <div>...</div> <div>Upper Left (-1500000.000, 1550000.000) (179d 3'38.91"E, 70d16'39.53"N)</div> <div>Lower Left (-1500000.000, 1500000.000) (180d 0' 0.00"W, 70d35'51.24"N)</div> <div>Upper Right (-1450000.000, 1550000.000) (178d 5'27.05"E, 70d35'13.15"N)</div> <div>Lower Right (-1450000.000, 1500000.000) (179d 1'44.32"E, 70d54'44.94"N)</div> <div>....</div> <div><br> </div> <div>Chci to přetransformovat do wgs84, používám snad všechny triky které jsem našel, ale gdalwarp prostě zůstane viset</div> <div><br> </div> <div>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<br> </div> </div> <div><br> </div> <div> <div>Creating output file that is 8009770P x 14124L.</div> <div>Processing input file 56_26_1_1_5m_v2.0_reg_dem.tif.</div> <div>Using internal nodata values (e.g. -9999) for image 56_26_1_1_5m_v2.0_reg_dem.tif.</div> <div>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.</div> </div> <div>0.</div> <div><br> </div> <div>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á.</div> <div><br> </div> <div>Dík</div> <div><br> </div> <div>J</div> <div><br> </div> <div><br> </div> </div> <br> <fieldset class="m_-7111605648087380987mimeAttachmentHeader"></fieldset> <br> </blockquote> </div> <div text="#000000" bgcolor="#FFFFFF"> <blockquote type="cite" style="border:none !important; margin-left:0px !important; margin-right:0px !important; margin-top:0px !important; padding-left:0px !important; padding-right:0px !important"> <pre>_______________________________________________ FreeGeoCZ mailing list <a class="m_-7111605648087380987moz-txt-link-abbreviated" href="mailto:[email protected]" target="_blank" moz-do-not-send="true">[email protected]</a> <a class="m_-7111605648087380987moz-txt-link-freetext" href="http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz" target="_blank" moz-do-not-send="true">http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz</a> </pre> </blockquote> <br> </div> </blockquote> </div> </blockquote> <br> </body> </html> --------------965FDA83C9D8144EC48D8F9D-- --===============0531825127== Content-Type: text/plain; charset="iso-8859-2" MIME-Version: 1.0 Content-Transfer-Encoding: quoted-printable Content-Disposition: inline _______________________________________________ FreeGeoCZ mailing list [email protected] http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz --===============0531825127==--