Hi,
I would guess that the appropriate zoom level is the one with pixel size
closest to the pixel size of the source file after it has been warped into
EPSG:3857.
Please add gdalinfo report about the source image.
-Jukka Rahkonen-
Evert Etienne (SITEMARK) wrote
> Dear,
>
> I have a question
Another option is the GMT program grdtrack
(https://docs.generic-mapping-tools.org/dev/grdtrack.html). It's pure C and
extremely fast.
-Original Message-
From: gdal-dev On Behalf Of Felix
Sent: Thursday, December 3, 2020 1:16 PM
To: gdal-dev@lists.osgeo.org
Subject: [gdal-dev] Sampling
Dear,
I have a question about the following extract of the MBTiles driver docs
(https://gdal.org/drivers/raster/mbtiles.html#raster-creation-issues)
> When using the CreateCopy() API (such as with gdal_translate), automatic
> reprojection of the input dataset to EPSG:3857 (WebMercator) will be
Hi Felix,
I'm really not sure I understand what exactly you are trying to do.
But if you have not tried pyresample yet it might do what you want. See
https://pyresample.readthedocs.io/en/latest/
Christoph
On Thu, 3 Dec 2020 at 14:16, Felix
wrote:
> Hello All,
>
> I'm here for some suggestions
Hello All,
I'm here for some suggestions on how to sample from a raster band at
_non-regular_ points. Non-regular in this context means that the points
have arbitrary coordinates and do not lie on a regular lattice/mesh like
all pairs of integers would do.
In my application, there is an irre
Dear Even et al.,
To answer my own question. Based on your comment on the units, I was able
to make it. By changing the Affine transformation as well. So the
gdaltransform example command becomes like this:
echo "0 -365" | gdaltransform -s_srs "+proj=stere +x_0=0 +y_0=0
+lat_0=90 +lon_0=0 +lat
Dear Even,
Thanks for your response! It's much appreciated. I control the creation of
the source file, so I can adjust the unit of the ellipsoid axis in the
source. Would that help with avoiding the need of a -ct pipeline?
Regards,
Mattijn
(forgot the Replay All, to send to the full list, inste
Mattijn,
> Question: how do I define my source raster instance, to avoid using the
> advance `-ct` operator in gdalwarp in GDAL3 during reprojection?
I don't think you can avoid it here. There seems to be an issue in either your
source file or in the HDF5 driver regarding the unit of the ellipso
