[FreeGeoCZ] gdalwarp stavkuje asi na dateline

--===============1786772271==
Content-Type: multipart/alternative; boundary="001a113ff91e059cd005662f561c"

--001a113ff91e059cd005662f561c
Content-Type: text/plain; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

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=3DY=
ES -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.t=
if
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 sou=
boru 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=A1eno=
sti? Dateline je jedin=C3=A9,
co m=C4=9B napad=C3=A1.

D=C3=ADk

J

--001a113ff91e059cd005662f561c
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr">Ahoj,<div><br></div><div>zkou=C5=BE=C3=ADm warpovat data z=
 Arcticdem zdroje, kolem dateline</div><div><br></div><div><div>$ wget <a h=
ref=3D"http://data.pgc.umn.edu/elev/dem/setsm/ArcticDEM/mosaic/v2.0/56_26/5=
6_26_1_1_5m_v2.0.tar.gz">http://data.pgc.umn.edu/elev/dem/setsm/ArcticDEM/m=
osaic/v2.0/56_26/56_26_1_1_5m_v2.0.tar.gz</a></div><div><br></div><div>$ ta=
r xzf 56_26_1_1_5m_v2.0.tar.gz</div><div><br></div><div>$ ls=C2=A0</div><di=
v>...</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&#39;38.91&quo=
t;E, 70d16&#39;39.53&quot;N)</div><div>Lower Left=C2=A0 (-1500000.000, 1500=
000.000) (180d 0&#39; 0.00&quot;W, 70d35&#39;51.24&quot;N)</div><div>Upper =
Right (-1450000.000, 1550000.000) (178d 5&#39;27.05&quot;E, 70d35&#39;13.15=
&quot;N)</div><div>Lower Right (-1450000.000, 1500000.000) (179d 1&#39;44.3=
2&quot;E, 70d54&#39;44.94&quot;N)</div><div>....</div><div><br></div><div>C=
hci to p=C5=99etransformovat do wgs84, pou=C5=BE=C3=ADv=C3=A1m snad v=C5=A1=
echny triky kter=C3=A9 jsem na=C5=A1el, ale gdalwarp prost=C4=9B z=C5=AFsta=
ne viset</div><div><br></div><div>gdalwarp -wo SOURCE_EXTRA=3D1000 -wo SAMP=
LE_STEPS=3D100 -wo SAMPLE_GRID=3DYES -co COMPRESS=3DDEFLATE -co BIGTIFF=3DY=
ES -s_srs +init=3Depsg:3411 -t_srs &quot;+proj=3Dlonglat +datum=3DWGS84 +ov=
er +no_defs&quot; 56_26_1_1_5m_v2.0_reg_dem.tif 56_26_1_1_5m_v2.0_reg_dem-w=
arped.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_re=
g_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><b=
r></div><div>D=C3=ADk</div><div><br></div><div>J</div><div><br></div><div><=
br></div></div>

--001a113ff91e059cd005662f561c--

--===============1786772271==
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

--===============1786772271==--