Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline
Jachym Cepicky <[email protected]> Thu, 01 Mar 2018 20:03:44 +0000
| Newsgroups | gmane.comp.gis.grass.cz |
|---|---|
| Message-ID | <CAAZUH4GP61cQzygdGTJB0tYiEz61JvR_FKGtHXjmZkQyqc2ZSQ@mail.gmail.com> |
--===============0795777287== Content-Type: multipart/alternative; boundary="001a114198eefcc26f05665f5b96" --001a114198eefcc26f05665f5b96 Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: quoted-printable no zejo, to by melo vyplyvat z tech dat, ale dik za tip, to "-te" by melo zabrat =C4=8Dt 1. 3. 2018 v 10:27 odes=C3=ADlatel Antonin Orlik <antonin.orlik@gma= il.com> napsal: > Abych to doplnil - probl=C3=A9m je v tom, =C5=BEe hranice timeline m=C5= =AF=C5=BEe b=C3=BDt uvedena > jako -180 (180W) i jako 180 (180E) a je mo=C5=BEno mezi n=C4=9B d=C3=A1t = rovn=C3=ADtko. > gdalwarp to bohu=C5=BEel vzal z t=C3=A9 opa=C4=8Dn=C3=A9 strany ne=C5=BE = bylo pot=C5=99eba a tomu > uzp=C5=AFsobil v=C3=BDsledn=C3=BD extent - a tak mu bylo pot=C5=99eb pomo= ci a v -te uv=C3=A9st m=C3=ADsto > -180 hodnotu 179.9999999. > Nemysl=C3=ADm si proto, =C5=BEe by to byla chyba GDALu - on nem=C5=AF=C5= =BEe ved=C4=9Bt, ze kter=C3=A9 > strany timeline chce=C5=A1 hranici extentu... > > Tonda > > > ------------------------------ > *From:* Jachym Cepicky <[email protected]> > <[email protected]> > *Subject:* [FreeGeoCZ] gdalwarp stavkuje asi na dateline > *Date:* Thursday, Mar 1, 2018 8:35 GMT > *To:* Antonin Orlik <[email protected]> <[email protected]> > *Cc:* Svobodn=C3=A1 geoinforma=C4=8Dn=C3=AD infrastruktura <freegeocz@fsv= .cvut.cz> > <[email protected]> > > 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/t= icket/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@gm= ail.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=A9= ho >> 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 = 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 omezit 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: >> >> 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_de= m.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_de= m.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 = 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. >> >> D=C3=ADk >> >> J >> >> >> >> >> _______________________________________________ >> FreeGeoCZ mailing [email protected]://mailman.fsv.cvut.cz/ma= ilman/listinfo/freegeocz >> >> >> > --001a114198eefcc26f05665f5b96 Content-Type: text/html; charset="UTF-8" Content-Transfer-Encoding: quoted-printable <div dir=3D"ltr">no zejo, to by melo vyplyvat z tech dat, ale dik za tip, t= o "-te" by melo zabrat</div><br><div class=3D"gmail_quote"><div d= ir=3D"ltr">=C4=8Dt 1. 3. 2018 v 10:27 odes=C3=ADlatel Antonin Orlik <<a = href=3D"mailto:[email protected]">[email protected]</a>> nap= sal:<br></div><blockquote class=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"> Abych to doplnil - probl=C3=A9m je v tom, =C5=BEe hranice timeline m=C5= =AF=C5=BEe b=C3=BDt uvedena jako -180 (180W) i jako 180 (180E) a je mo=C5=BEno mezi n=C4=9B= d=C3=A1t rovn=C3=ADtko. gdalwarp to bohu=C5=BEel vzal z t=C3=A9 opa=C4=8Dn=C3=A9= strany ne=C5=BE bylo pot=C5=99eba a tomu uzp=C5=AFsobil v=C3=BDsledn=C3=BD extent - a tak mu= bylo pot=C5=99eb pomoci a v -te uv=C3=A9st m=C3=ADsto -180 hodnotu 179.9999999.<br> Nemysl=C3=ADm si proto, =C5=BEe by to byla chyba GDALu - on nem=C5=AF= =C5=BEe ved=C4=9Bt, ze kter=C3=A9 strany timeline chce=C5=A1 hranici extentu...<br> <br> <div class=3D"m_-8336240161457324689moz-cite-prefix"> <div id=3D"m_-8336240161457324689rwhMsgHeader"></div></div></div><div= text=3D"#000000" bgcolor=3D"#FFFFFF"><div class=3D"m_-8336240161457324689m= oz-cite-prefix"><div id=3D"m_-8336240161457324689rwhMsgHeader">Tonda<br> <br> <br> <hr id=3D"m_-8336240161457324689rwhMsgHdrDivider" 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_-8336240161457324689moz-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> </div></div></div><div text=3D"#000000" bgcolor=3D"#FFFFFF"><div cl= ass=3D"m_-8336240161457324689moz-cite-prefix"><div id=3D"m_-833624016145732= 4689rwhMsgHeader"><span><font style=3D"font:13px Tahoma!important;color:#00= 0000!important" face=3D"Tahoma" color=3D"#000000"><b>Date:</b> Thursday, Mar 1, 2018 8:35 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> Antonin Orlik <a class=3D"m_-8336240161457324689moz-txt-link-rf= c2396E" href=3D"mailto:[email protected]" target=3D"_blank"><anton= [email protected]></a></font></span><br> <span><font style=3D"font:13px Tahoma!important;color:#000000!impor= tant" face=3D"Tahoma" color=3D"#000000"><b>Cc:</b> Svobodn=C3=A1 geoinforma=C4=8Dn=C3=AD infrastruktura <a class=3D"m_-8336240161457324689moz-txt-link-rfc2396E" href= =3D"mailto:[email protected]" target=3D"_blank"><[email protected].= cz></a></font></span><br> <br> </div> </div></div><div text=3D"#000000" bgcolor=3D"#FFFFFF"> <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"> <div dir=3D"ltr">Vid=C3=AD=C5=A1, to je zaj=C3=ADmavej typ,=C2=A0<br> <br> mezit=C3=ADm tam n=C4=9Bco patchnul i Even=C2=A0<a href=3D"https://= trac.osgeo.org/osgeo/ticket/2096" target=3D"_blank">https://trac.osgeo.org/= osgeo/ticket/2096</a>=C2=A0ale 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. 2018 v 8:59 odes=C3=ADlatel Antonin = Orlik <<a href=3D"mailto:[email protected]" target=3D"_blank">= [email protected]</a>> napsal:<br> </div> <blockquote class=3D"gmail_quote" style=3D"border:none!important;ma= rgin-left:0px!important;margin-right:0px!important;margin-top:0px!important= ;padding-left:0px!important;padding-right:0px!important"> <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&= #39; 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 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 omezit v= =C3=BDstupn=C3=AD extent pomoc=C3=AD parametru -te - ve tv=C3=A9m p=C5=99=C3=ADpa= d=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_-8336240161457324689m_-7111605648087380987moz-c= ite-prefix"> <div id=3D"m_-8336240161457324689m_-7111605648087380987rwhMsg= Header"><br> <hr id=3D"m_-8336240161457324689m_-7111605648087380987rwhMs= gHdrDivider" style=3D"border:0;border-top:1px solid #b5c4df;padding:0;margi= n:10px 0 5px 0;width:100%"><span><font style=3D"font:13px Tahoma!important;= color:#000000!important" face=3D"Tahoma" color=3D"#000000"><b>From:</b> Jac= hym Cepicky <a class=3D"m_-8336240161457324689m_-7111605648= 087380987moz-txt-link-rfc2396E" href=3D"mailto:[email protected]" ta= rget=3D"_blank"><[email protected]></a></font></span><br> <span><font style=3D"font:13px Tahoma!important;color:#0000= 00!important" face=3D"Tahoma" color=3D"#000000"><b>Subject:</b> [FreeGeoCZ] gdalwarp stavkuje asi na dateline</font></s= pan><br> <span><font style=3D"font:13px Tahoma!important;color:#0000= 00!important" 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:#0000= 00!important" face=3D"Tahoma" color=3D"#000000"><b>To:</b> FreeGeoCZ <a cla= ss=3D"m_-8336240161457324689m_-7111605648087380987moz-txt-link-rfc2396E" hr= ef=3D"mailto:[email protected]" target=3D"_blank"><[email protected]= t.cz></a></font></span><br> <br> </div> </div> </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;pad= ding-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, ko= lem dateline</div> <div><br> </div> <div> <div>$ wget <a href=3D"http://data.pgc.umn.edu/elev/dem/s= etsm/ArcticDEM/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_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=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"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=C5=99etransformovat do wgs84, pou=C5=BE=C3= =ADv=C3=A1m snad v=C5=A1echny triky kter=C3=A9 jsem na=C5=A1el, ale gdal= warp 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_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> </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=ADlov=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_-8336240161457324689m_-71116056480873809= 87mimeAttachmentHeader"></fieldset> <br> </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;pad= ding-left:0px!important;padding-right:0px!important"> <pre>_______________________________________________ FreeGeoCZ mailing list <a class=3D"m_-8336240161457324689m_-7111605648087380987moz-txt-link-abbrev= iated" href=3D"mailto:[email protected]" target=3D"_blank">FreeGeoCZ@fs= v.cvut.cz</a> <a class=3D"m_-8336240161457324689m_-7111605648087380987moz-txt-link-freete= xt" href=3D"http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz" target= =3D"_blank">http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz</a> </pre> </blockquote> <br> </div> </blockquote> </div> </blockquote> <br> </div></blockquote></div> --001a114198eefcc26f05665f5b96-- --===============0795777287== 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 --===============0795777287==--