Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline
Jachym Cepicky <[email protected]> Thu, 01 Mar 2018 08:35:30 +0000
| Newsgroups | gmane.comp.gis.grass.cz |
|---|---|
| Message-ID | <CAAZUH4F6Ppu0nKJmqUCPqRWZtxN76f7rkZgF2w6DYBY64yPSQg@mail.gmail.com> |
--===============1551614901== Content-Type: multipart/alternative; boundary="001a1141aeecbd4a70056655bed3" --001a1141aeecbd4a70056655bed3 Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: quoted-printable Vid=C3=AD=C5=A1, to je zaj=C3=ADmavej typ, mezit=C3=ADm tam n=C4=9Bco patchnul i Even https://trac.osgeo.org/osgeo/tic= ket/2096 ale zat=C3=ADm se mi to nepoda=C5=99ilo rozb=C4=9Bhnout J =C4=8Dt 1. 3. 2018 v 8:59 odes=C3=ADlatel Antonin Orlik <antonin.orlik@gmai= l.com> napsal: > Ahoj J=C3=A1chyme, > > probl=C3=A9m je skute=C4=8Dn=C4=9B v dateline. Kdy=C5=BE se pod=C3=ADv=C3= =A1=C5=A1 na sou=C5=99adnice extentu, > kter=C3=A9 jsi poslal, tak minx v lat/lon (resp. ve WGS 84 do kter=C3=A9h= o > transformuje=C5=A1) je 180d 0' 0.00"W a maxx je 179d 3'38.91"E. Rozsah > v=C3=BDstupn=C3=ADho extentu tedy za=C4=8D=C3=ADn=C3=A1 na -180=C2=B0 a k= on=C4=8D=C3=AD na n=C4=9Bjak=C3=BDch +179. Extent > pak tedy obep=C3=ADn=C3=A1 celou zem=C4=9Bkouli a velikost v=C3=BDstupn= =C3=ADho rastru je pak tedy > neuv=C4=9B=C5=99iteln=C3=BDch 8009770 x 14124 pixel=C5=AF a transformace = je tedy prakticky > nemo=C5=BEn=C3=A1. Proto je nutn=C3=A9 omezit v=C3=BDstupn=C3=AD extent p= omoc=C3=AD parametru -te - ve > tv=C3=A9m p=C5=99=C3=ADpad=C4=9B konkr=C3=A9tn=C4=9B takto: > > gdalwarp -wo SOURCE_EXTRA=3D1000 -wo SAMPLE_STEPS=3D100 -wo SAMPLE_GRID= =3DYES > -te 178.090847222 70.2776472222 179.99999 70.9124833333 -co > COMPRESS=3DDEFLATE -co BIGTIFF=3DYES -s_srs +init=3Depsg:3411 -t_srs > "+proj=3Dlonglat +datum=3DWGS84 +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=C3=A1 p=C3=A1r vte=C5=99in a v=C3=BDsledek sed=C3=AD dob=C5= =99e... > > M=C4=9Bj se! > Tonda > > ------------------------------ > *From:* Jachym Cepicky <[email protected]> > <[email protected]> > *Subject:* [FreeGeoCZ] gdalwarp stavkuje asi na dateline > *Date:* Tuesday, Feb 27, 2018 10:46 GMT > *To:* FreeGeoCZ <[email protected]> <[email protected]> > > Ahoj, > > zkou=C5=BE=C3=ADm 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=C5=99etransformovat do wgs84, pou=C5=BE=C3=ADv=C3=A1m snad v=C5= =A1echny triky kter=C3=A9 jsem > na=C5=A1el, ale gdalwarp prost=C4=9B z=C5=AFstane viset > > gdalwarp -wo SOURCE_EXTRA=3D1000 -wo SAMPLE_STEPS=3D100 -wo SAMPLE_GRID= =3DYES > -co COMPRESS=3DDEFLATE -co BIGTIFF=3DYES -s_srs +init=3Depsg:3411 -t_srs > "+proj=3Dlonglat +datum=3DWGS84 +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=C5=AFstane viset na 0% a d=C3=A1l se ani nehne. Do c=C3=ADlov=C3=A9ho s= ouboru i n=C4=9Bco zap=C3=AD=C5=A1e, > ale je to jenom p=C3=A1r mega. Nem=C3=A1 s t=C3=ADm n=C4=9Bkdo zku=C5=A1e= nosti? Dateline je jedin=C3=A9, > co m=C4=9B napad=C3=A1. > > D=C3=ADk > > J > > > > > _______________________________________________ > FreeGeoCZ mailing [email protected]://mailman.fsv.cvut.cz/mai= lman/listinfo/freegeocz > > > --001a1141aeecbd4a70056655bed3 Content-Type: text/html; charset="UTF-8" Content-Transfer-Encoding: quoted-printable <div dir=3D"ltr">Vid=C3=AD=C5=A1, to je zaj=C3=ADmavej typ,=C2=A0<br><br>me= zit=C3=ADm tam n=C4=9Bco patchnul i Even=C2=A0<a href=3D"https://trac.osgeo= .org/osgeo/ticket/2096">https://trac.osgeo.org/osgeo/ticket/2096</a>=C2=A0a= le zat=C3=ADm se mi to nepoda=C5=99ilo rozb=C4=9Bhnout<div><br></div><div>J= </div></div><br><div class=3D"gmail_quote"><div dir=3D"ltr">=C4=8Dt 1. 3. 2= 018 v 8:59 odes=C3=ADlatel Antonin Orlik <<a href=3D"mailto:antonin.orli= [email protected]">[email protected]</a>> napsal:<br></div><blockquote c= lass=3D"gmail_quote" style=3D"margin:0 0 0 .8ex;border-left:1px #ccc solid;= padding-left:1ex"> =20 =20 =20 <div text=3D"#000000" bgcolor=3D"#FFFFFF"> Ahoj J=C3=A1chyme,<br> <br> probl=C3=A9m je skute=C4=8Dn=C4=9B v dateline. Kdy=C5=BE se pod=C3=ADv= =C3=A1=C5=A1 na sou=C5=99adnice extentu, kter=C3=A9 jsi poslal, tak minx v lat/lon (resp. ve WGS 84 do kter=C3=A9ho transformuje=C5=A1) je 180d 0' 0.00"W a maxx je 1= 79d 3'38.91"E. Rozsah v=C3=BDstupn=C3=ADho extentu tedy za=C4=8D=C3=ADn=C3=A1 na -180= =C2=B0 a kon=C4=8D=C3=AD na n=C4=9Bjak=C3=BDch +179. Extent pak tedy obep=C3=ADn=C3=A1 celou zem=C4=9Bkouli a velikost= v=C3=BDstupn=C3=ADho rastru je pak tedy neuv=C4=9B=C5=99iteln=C3=BDch 8009770 x 14124 pixel= =C5=AF a transformace je tedy prakticky nemo=C5=BEn=C3=A1. Proto je nutn=C3=A9 o= mezit v=C3=BDstupn=C3=AD extent pomoc=C3=AD parametru -te - ve tv=C3=A9m p=C5= =99=C3=ADpad=C4=9B konkr=C3=A9tn=C4=9B takto:<br> <br> gdalwarp -wo SOURCE_EXTRA=3D1000 -wo SAMPLE_STEPS=3D100 -wo SAMPLE_GRID=3DYES -te 178.090847222 70.2776472222 179.99999 70.9124833333 -co COMPRESS=3DDEFLATE -co BIGTIFF=3DYES -s_srs +init=3Depsg:3411 -t_srs "+proj=3Dlonglat +datum=3DWGS84 +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=3D"#000000= " bgcolor=3D"#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=3D= "#000000" bgcolor=3D"#FFFFFF"> 0...10...20...30...40...50...60...70...80...90...100 - done.<br> <br> Konverze trv=C3=A1 p=C3=A1r vte=C5=99in a v=C3=BDsledek sed=C3=AD dob= =C5=99e...<br> <br> M=C4=9Bj se!<br> Tonda <div class=3D"m_-7111605648087380987moz-cite-prefix"> <div id=3D"m_-7111605648087380987rwhMsgHeader"><br> <hr id=3D"m_-7111605648087380987rwhMsgHdrDivider" style=3D"border:0= ;border-top:1px solid #b5c4df;padding:0;margin:10px 0 5px 0;width:100%"><sp= an><font style=3D"font:13px Tahoma!important;color:#000000!important" face= =3D"Tahoma" color=3D"#000000"><b>From:</b> Jachym Cepicky <a class=3D"m_-7111605648087380987moz-txt-link-rfc2396E" href= =3D"mailto:[email protected]" target=3D"_blank"><jachym.cepicky@g= mail.com></a></font></span><br> <span><font style=3D"font:13px Tahoma!important;color:#000000!impor= tant" face=3D"Tahoma" color=3D"#000000"><b>Subject:</b> [FreeGeoCZ] gdalwarp stavkuje asi na dateline</font></span><br> <span><font style=3D"font:13px Tahoma!important;color:#000000!impor= tant" face=3D"Tahoma" color=3D"#000000"><b>Date:</b> Tuesday, Feb 27, 2018 10:46 GMT</font></span><br> <span><font style=3D"font:13px Tahoma!important;color:#000000!impor= tant" face=3D"Tahoma" color=3D"#000000"><b>To:</b> FreeGeoCZ <a class=3D"m_-7111605648087380987moz-txt-link-rfc239= 6E" href=3D"mailto:[email protected]" target=3D"_blank"><freegeocz@f= sv.cvut.cz></a></font></span><br> <br> </div> </div> <blockquote type=3D"cite" style=3D"border:none!important;margin-left:0p= x!important;margin-right:0px!important;margin-top:0px!important;padding-lef= t:0px!important;padding-right:0px!important"></blockquote></div><div text= =3D"#000000" bgcolor=3D"#FFFFFF"><blockquote type=3D"cite" style=3D"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=3D"ltr">Ahoj, <div><br> </div> <div>zkou=C5=BE=C3=ADm warpovat data z Arcticdem zdroje, kolem date= line</div> <div><br> </div> <div> <div>$ wget <a href=3D"http://data.pgc.umn.edu/elev/dem/setsm/Arc= ticDEM/mosaic/v2.0/56_26/56_26_1_1_5m_v2.0.tar.gz" target=3D"_blank">http:/= /data.pgc.umn.edu/elev/dem/setsm/ArcticDEM/mosaic/v2.0/56_26/56_26_1_1_5m_v= 2.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=C2=A0</div> <div>...</div> <div>56_26_1_1_5m_v2.0_reg_dem.tif<br> </div> <div>...</div> <div><br> </div> <div>$ gdalinfo=C2=A056_26_1_1_5m_v2.0_reg_dem.tif</div> <div>...</div> <div>Upper Left=C2=A0 (-1500000.000, 1550000.000) (179d 3'38.= 91"E, 70d16'39.53"N)</div> <div>Lower Left=C2=A0 (-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&qu= ot;E, 70d35'13.15"N)</div> <div>Lower Right (-1450000.000, 1500000.000) (179d 1'44.32&qu= ot;E, 70d54'44.94"N)</div> <div>....</div> <div><br> </div> <div>Chci to p=C5=99etransformovat do wgs84, pou=C5=BE=C3=ADv=C3= =A1m snad v=C5=A1echny triky kter=C3=A9 jsem na=C5=A1el, ale gdalwarp prost=C4=9B z=C5= =AFstane viset</div> <div><br> </div> <div>gdalwarp -wo SOURCE_EXTRA=3D1000 -wo SAMPLE_STEPS=3D100 -wo SAMPLE_GRID=3DYES -co COMPRESS=3DDEFLATE -co BIGTIFF=3DYES -s_s= rs +init=3Depsg:3411 -t_srs "+proj=3Dlonglat +datum=3DWGS84 += 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=C5=AFstane viset na 0% a d=C3=A1l se ani nehne. Do c=C3=ADlo= v=C3=A9ho souboru i n=C4=9Bco zap=C3=AD=C5=A1e, ale je to jenom p=C3=A1r mega. Nem= =C3=A1 s t=C3=ADm n=C4=9Bkdo zku=C5=A1enosti? Dateline je jedin=C3=A9, co m=C4=9B napad=C3=A1.= </div> <div><br> </div> <div>D=C3=ADk</div> <div><br> </div> <div>J</div> <div><br> </div> <div><br> </div> </div> <br> <fieldset class=3D"m_-7111605648087380987mimeAttachmentHeader"></fiel= dset> <br> </blockquote></div><div text=3D"#000000" bgcolor=3D"#FFFFFF"><blockqu= ote type=3D"cite" style=3D"border:none!important;margin-left:0px!important;= margin-right:0px!important;margin-top:0px!important;padding-left:0px!import= ant;padding-right:0px!important"><pre>_____________________________________= __________ FreeGeoCZ mailing list <a class=3D"m_-7111605648087380987moz-txt-link-abbreviated" href=3D"mailto:= [email protected]" target=3D"_blank">[email protected]</a> <a class=3D"m_-7111605648087380987moz-txt-link-freetext" href=3D"http://mai= lman.fsv.cvut.cz/mailman/listinfo/freegeocz" target=3D"_blank">http://mailm= an.fsv.cvut.cz/mailman/listinfo/freegeocz</a> </pre> </blockquote> <br> </div> </blockquote></div> --001a1141aeecbd4a70056655bed3-- --===============1551614901== 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 --===============1551614901==--