Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline
Antonin Orlik <[email protected]> Thu, 1 Mar 2018 08:59:52 +0100
| Newsgroups | gmane.comp.gis.grass.cz |
|---|---|
| Message-ID | <[email protected]> |
This is a multi-part message in MIME format. --===============0281303012== Content-Type: multipart/alternative; boundary="------------9A0ACCE88B6B835855E3594F" Content-Language: cs-CZ This is a multi-part message in MIME format. --------------9A0ACCE88B6B835855E3594F Content-Type: text/plain; charset=utf-8; format=flowed Content-Transfer-Encoding: 8bit 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]> *Subject:* [FreeGeoCZ] gdalwarp stavkuje asi na dateline *Date:* Tuesday, Feb 27, 2018 10:46 GMT *To:* FreeGeoCZ <[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] > http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz --------------9A0ACCE88B6B835855E3594F 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"> 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.<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> 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="moz-cite-prefix"> <div id="rwhMsgHeader"><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> Tuesday, Feb 27, 2018 10:46 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> FreeGeoCZ <a class="moz-txt-link-rfc2396E" href="mailto:[email protected]"><[email protected]></a></font></span><br> <br> </div> </div> <blockquote type="cite" cite="mid:CAAZUH4FdUxUK2X8zJhw7iSUxsc7tu+v0R3UTzSH_tHKbqCw0fA@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">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" 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="mimeAttachmentHeader"></fieldset> <br> <pre wrap="">_______________________________________________ FreeGeoCZ mailing list <a class="moz-txt-link-abbreviated" href="mailto:[email protected]">[email protected]</a> <a class="moz-txt-link-freetext" href="http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz">http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz</a> </pre> </blockquote> <br> </body> </html> --------------9A0ACCE88B6B835855E3594F-- --===============0281303012== 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 --===============0281303012==--