Re: [gdal-dev] [Question]Output from gdal_translate projwin is different from geoserver WCS GetCoverage request?

It completes with warning: > Warning 1: Computed -srcwin -1 -1 17 4 falls partially outside raster extent. > Going on however. > Try again with your -projwin values within your geoserver.tif extent and see > what happens. On Mon, 3 Apr 2023 at 14:29, Pham Huu Bang wrote: >

[gdal-dev] [Question]Output from gdal_translate projwin is different from geoserver WCS GetCoverage request?

Hello, I've a hard time understanding which one should be the correct result by gdal_translate or GeoServer. My system (Ubuntu 22.04) has: - GDAL 3.4.1 - Geoserver 2.22 running on Tomcat on port 9090 I have a small tiff file (link to download at the bottom of this thread) in EPSG:4326, then I

Re: [gdal-dev] Process Sentinel-5p netCDF to tiff files but overlapping appears clearly between 2 files?

=sharing Best, From: Pham Huu, Bang Sent: Monday, November 8, 2021 07:09 To: gdal-dev@lists.osgeo.org Subject: Process Sentinel-5p netCDF to tiff files but overlapping appears clearly between 2 files? Hello, I've used 2 Sentinel-5p scenes in hdf5 / netCDF to geoTiff

Re: [gdal-dev] gdalwarp from EPSG:25832 to EPSG:4326 with error latitude or longitude exceeded limits?

Thank to Jukka, so the problem was with my gdal 1.11.4 on Centos 7. When I translated it to EPSG:4647, it doesn't contain proper WKT for EPSG, e.g: Coordinate System is: LOCAL_CS["ETRS89 / UTM zone 32N (zE-N)", GEOGCS["ETRS89", DATUM["unknown", SPHEROID["unretrievable -

Re: [gdal-dev] gdalwarp from EPSG:25832 to EPSG:4326 with error latitude or longitude exceeded limits?

The tiff file you can download from https://files.fm/f/5ttjhyh9p The xyz file I already shared in the beginning: https://pastebin.com/40xSNqm4 Thanks, On Wed, 25 Nov 2020 at 11:10, Tobias Wendorff < tobias.wendo...@tu-dortmund.de> wrote: > Am 25.11.2020 um 11:00 schrieb Pham

Re: [gdal-dev] gdalwarp from EPSG:25832 to EPSG:4326 with error latitude or longitude exceeded limits?

SG:4326 test2.tif ERROR 1: No PROJ.4 translation for source SRS, coordinate transformation initialization has failed. On Wed, 25 Nov 2020 at 10:55, Tobias Wendorff < tobias.wendo...@tu-dortmund.de> wrote: > Am 25.11.2020 um 10:17 schrieb Pham Huu Bang: > > Corner Coordinates: > > Up

[gdal-dev] gdalwarp from EPSG:25832 to EPSG:4326 with error latitude or longitude exceeded limits?

Hello, I have a DGM 10m tiff file (translated from xyz file: e.g:dgm10_32572_6020_1_sh.xyz - example: https://pastebin.com/40xSNqm4) in EPSG:25832 CRS with this gdalinfo as below. I tried to warp it to EPSG:4326 but I got error. Do I miss something? Thanks, gdalwarp test.tif -t_srs EPSG:4326

Re: [gdal-dev] Python/Java GDAL - How to estimate the geo coordinates and width / height without warping dataset?

@Even Rouault Unfortunately, this method is not enough to solve my problem as I *also have a given geo X, Y resolutions *in target CRS (EPSG:4326) as input. Below is what I've with Python GDAL: src_ds = gdal.GetDriverByName("VRT").Create("", 1830, 1830)

[gdal-dev] Python/Java GDAL - How to estimate the geo coordinates and width / height without warping dataset?

Hello, I have a Sentinel 2 in UTM32 with these metadata from gdalinfo: Size is 1830, 1830 Pixel Size = (60.000,-60.000) Lower Left ( 60.000, 5790240.000) Upper Left ( 60.000, 5900040.000) How can I have these same metadata* by GDAL Python/Java binding* but in

Re: [gdal-dev] Results of gdalwarp from UTM32 to EPSG:4326 are not aligned?

T32UQD_20200421T102021_TCI_60m.jp2 full_warp.tif -t_srs EPSG:4326 -dstnodata 0 -tap gdalwarp -tr 0.0007161230037719291 -0.0007161230037719291 crop.tif small_warp.tif -t_srs EPSG:4326 -dstnodata 0 -tap On Tue, 12 May 2020 at 18:12, Even Rouault wrote: > On mardi 12 mai 2020 16:19:43 CEST Pham Huu Bang wr

[gdal-dev] Results of gdalwarp from UTM32 to EPSG:4326 are not aligned?

Hello, I've a Sentinel 2 UTM32 scene (JPEG2000 format). What I want to do is: + Warp this scene from UTM32 to EPSG:4326 by gdalwarp gdalwarp T32UQD_20200421T102021_TCI_60m.jp2 -t_srs EPSG:4326 full_warp.tif + Clip a small subset of this file by gdal_translate gdal_translate -projwin

[gdal-dev] Difference results from QGIS project on fly and gdalwarp for EPSG:4326 of Sentinel 2 UTM zone 32?

Hello, I always thought QGIS and gdalwarp should return same result for projecting a 2D Sentinel 2 image (e.g: T32UQD_20200421T102021_TCI_60m.jp2) from UTM 32 to EPSG:4326. What I did was: - Open the .jp2 in QGis 2.14.9 with 'Enable on the fly CRS transformation'. - Warp the .jp2

[gdal-dev] How to calculate affine transformation by Python gdal as same as gdal_translate?

Hello, I've tried to reproduce a problem to subset (clip) an UTM 32 Sentinel 2 GeoTiff by a given Bounding box (BBOX) in EPSG:4326. I've created a python gdal script for it, but the result is always contain few more pixels than what I've tried with gdal_translate binary. All information of the

[gdal-dev] How to find correct 2D bounding box to subset in different CRS with gdal (different outputs for gdalwarp)?

Hello, I'm using GDAL version 1.11.4 and I have a problem how to use gdal to translate a bounding box correctly [Lat(62:63), Long(10:11)] from source CRS EPSG:4326 to target CRS EPSG:32632. I've tried so far with gdaltransform to translate coordinates: Long min, Lat min and Long max, Lat max