Hi,

I would expect these 2 ways of running GDAL warp should produce the same
results.

gdalwarp -t_srs EPSG:3857 \
         -te 558728 5789851 1242296 7544030 \
         -ts 500 500 \
         "WMS:
http://127.0.0.1:8080/api/resource/5840/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=test_rounds&SRS=EPSG:4326&BBOX=-180.000000,-85.051129,180.000000,85.051129";
\
         /tmp/result1.tif

and

from osgeo import gdal
from osgeo import gdalconst

width = height = 500
url = "WMS:
http://127.0.0.1:8080/api/resource/5840/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=test_rounds&SRS=EPSG:4326&BBOX=-180.000000,-85.051129,180.000000,85.051129
"
ds = gdal.Open(url, gdalconst.GA_ReadOnly)
ds_img = gdal.Warp(
    "/tmp/result2.tif",
    ds,
    options=gdal.WarpOptions(
        width=width,
        height=height,
        outputBounds=(558728, 5789851, 1242296, 7544030),
        dstSRS="EPSG:3857",
    ),
)

But in the first case I get what I expect but in the second one I get a
raster with right spatial extent but 0 in all bands.

What am I missing?
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to