gdalwarp WebGIS EPSG:3857 ke EPSG:32647 (UTM Zona 47 North)

Pekanbaru berada di zona 47 utara, sehingga kode EPSG (Masyarakat Geo Minyak Eropa) adalah 32647. Namun, warping gambar dari citra satelit “Google” ataupun webgis lain yang ber-EPSG 3857 terkadang menemui banyak kendala.

Sebagai contoh:

 

gdalwarp -overwrite -s_srs EPSG:3857 -t_srs EPSG:32647 -multi -of GTiff "D:/GIS/Pekanbaru/Pekanbaru Lv 17 World Imagery/17-102400_102495-65280_65343-24576x16384.jpg" D:/GIS/Pekanbaru/pekanbaru_32647/17-102400_102495-65280_65343.tif

Kode di atas untuk mengubah gambar dari JPG (yang sudah ada World File nya, wld) yang source_SRS nya adalah 3857 ke target kita, UTM Zona 47N, EPSG:32647

Namun, yang menjadi kendala adalah, selain prosesnya lama, karena memang ukuran citra-nya yang 24576×16384 pixels, hasilnya berkas GeoTIFF nya pun besar!

Berkas JPG awal hanya 123MB, GeoTIFF hasilnya sampai 1,17 Giga lebih!

 

Nyari-nyari solusi di Internet, ya … seharusnya GeoTIFF nya dikompress, NAMUN …

Microsoft Windows [Version 10.0.14393]
(c) 2016 Microsoft Corporation. All rights reserved.

C:\Users\USER>gdalwarp --help
'gdalwarp' is not recognized as an internal or external command,
operable program or batch file.

C:\Users\USER>C:\OSGeo4W\OSGeo4W.bat

OSGEO4W home is C:\OSGeo4W\

run o-help for a list of available commands
C:\Users\USER>gdalwarp --help
Usage: gdalwarp [--help-general] [--formats]
 [-s_srs srs_def] [-t_srs srs_def] [-to "NAME=VALUE"]
 [-order n | -tps | -rpc | -geoloc] [-et err_threshold]
 [-refine_gcps tolerance [minimum_gcps]]
 [-te xmin ymin xmax ymax] [-tr xres yres] [-tap] [-ts width height]
 [-ovr level|AUTO|AUTO-n|NONE] [-wo "NAME=VALUE"] [-ot Byte/Int16/...] [-wt Byte/Int16]
 [-srcnodata "value [value...]"] [-dstnodata "value [value...]"] -dstalpha
 [-r resampling_method] [-wm memory_in_mb] [-multi] [-q]
 [-cutline datasource] [-cl layer] [-cwhere expression]
 [-csql statement] [-cblend dist_in_pixels] [-crop_to_cutline]
 [-of format] [-co "NAME=VALUE"]* [-overwrite]
 [-nomd] [-cvmd meta_conflict_value] [-setci] [-oo NAME=VALUE]*
 [-doo NAME=VALUE]*
 srcfile* dstfile

Available resampling methods:
 near (default), bilinear, cubic, cubicspline, lanczos, average, mode, max, min, med, Q1, Q3.

C:\Users\USER>

 

WOW! ternyata gdalwarp dukung kompres!

tinggal pake -co “NAME=VALUE”

 

NAMUN

https://trac.osgeo.org/gdal/wiki/UserDocs/GdalWarp

[GeoTIFF output] -co COMPRESS= is broken!

waduh!

 

Katanya, harus di warp dulu uncompress, baru di translate dengan kompresi!

hm …

 

Iklan