Vincent Sarago vincentsarago

the test.py is equivalent to : gdalwarp -t_srs EPSG:3857 -te -8922952.933898335 5244191.636589372 -8883817.175416324 5283327.395071382 -te_srs EPSG:3857 -ts 128 128 -r bilinear /vsis3/landsat-pds/c1/L8/017/030/LC08_L1TP_017030_20180228_20180308_01_T1/LC08_L1TP_017030_20180228_20180308_01_T1_B1.TIF test.tif

VSI_CACHE=FALSE VSI_CACHE_SIZE=0 CPL_DEBUG=ON CPL_CURL_VERBOSE=YES gdalwarp -t_srs EPSG:3857 -te -8922952.933898335 5244191.636589372 -8883817.175416324 5283327.395071382  -te_srs EPSG:3857 -ts 128 128 -r bilinear /vsis3/landsat-pds/c1/L8/017/030/LC08_L1TP_017030_20180228_20180308_01_T1/LC08_L1TP_017030_20180228_20180308_01_T1_B1.TIF test.tif 2>&1 >/dev/null | grep "Content-Length:"
< Content-Length: 16384
< Content-Length: 8700
< Content-Length: 8700
< Content-Length: 16384

All usefull data must be downloaded from https://s3-us-west-2.amazonaws.com/remotepixel-pub/data/IMW/IMW2018_QGIS_data_processing.zip

Data

/bin/sh /usr/local/src/gdal-2.3.1/libtool --mode=link --silent g++ -std=c++11   -lcrypto -lproj -L/usr/lib -lgeos_c -L/usr/local/lib -lwebp -lexpat -ljpeg -lpng -L/usr/local/lib -lzstd -lz -L/usr -L/usr/lib -lpthread -lm -lrt -ldl  -L/usr/lib -lspatialite   -L/usr/lib/x86_64-linux-gnu -lcurl -L/usr/local/src/MrSID_DSDK-9.1.0.4045-linux.x86-64.gcc44/Raster_DSDK/3rd-party/lib -L/usr/local/src/MrSID_DSDK-9.1.0.4045-linux.x86-64.gcc44/Raster_DSDK/lib -lltidsdk -lpthread -L/usr/local/src/MrSID_DSDK-9.1.0.4045-linux.x86-64.gcc44/Lidar_DSDK/lib -llti_lidar_dsdk -L/usr/local/src/lib/x64/release -lNCSEcw              -o libgdal.la ./ogr/gml2ogrgeometry.lo ./ogr/ogr2gmlgeometry.lo ./ogr/ogr_api.lo ./ogr/ogr_expat.lo ./ogr/ogr_fromepsg.lo ./ogr/ogr_geo_utils.lo ./ogr/ogr_geocoding.lo ./ogr/ogr_opt.lo ./ogr/ogr_srs_dict.lo ./ogr/ogr_srs_erm.lo ./ogr/ogr_srs_esri.lo ./ogr/ogr_srs_ozi.lo ./ogr/ogr_srs_panorama.lo ./ogr/ogr_srs_pci.lo ./ogr/ogr_srs_proj4.lo ./ogr/ogr_srs_usgs.lo ./ogr/ogr_srs_validate.lo ./ogr/ogr_srs_x

	# download http://soliton.vm.bytemark.co.uk/pub/cpt-city/wkp/schwarzwald/wiki-schwarzwald-cont.pg

	with open('schwarzwald.txt') as cmap:
	lines = cmap.read().splitlines()
	colormap = [list(map(int, line.split())) for line in lines if not line.startswith('#')][1:]

	new_cmap = []
	nb = len(colormap)
	cmap = [[round((ix / nb) * 257), cp] for ix, cp in enumerate(colormap[::-1])]