freegeocz AT fsv.cvut.cz
Předmět: Svobodná geoinformační infrastruktura
List archive
- From: Antonin Orlik <antonin.orlik AT gmail.com>
- To: Svobodná geoinformační infrastruktura <freegeocz AT fsv.cvut.cz>, Jachym Cepicky <jachym.cepicky AT gmail.com>
- Cc:
- Subject: Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline
- Date: Thu, 1 Mar 2018 08:59:52 +0100
- List-archive: <http://mailman.fsv.cvut.cz/pipermail/freegeocz>
- List-id: Svobodná geoinformační infrastruktura <freegeocz.fsv.cvut.cz>
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 <jachym.cepicky AT gmail.com> Subject: [FreeGeoCZ] gdalwarp stavkuje asi na dateline Date: Tuesday, Feb 27, 2018 10:46 GMT To: FreeGeoCZ <freegeocz AT fsv.cvut.cz> 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 FreeGeoCZ AT fsv.cvut.cz http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz |
- Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline, Antonin Orlik, 03/01/2018
- Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline, Jachym Cepicky, 03/01/2018
- Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline, Antonin Orlik, 03/01/2018
- Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline, Jachym Cepicky, 03/01/2018
- Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline, Antonin Orlik, 03/01/2018
- Re: [FreeGeoCZ] gdalwarp stavkuje asi na dateline, Jachym Cepicky, 03/01/2018
Archivace běží na MHonArc 2.6.19+.