Re: [gdal-dev] HDF5 and geolocation arrays

2024-04-29 Thread Michael Sumner via gdal-dev
you need gdalinfo --config GDAL_HTTP_HEADERS="Authorization: Bearer " "/vsicurl/https://n5eil01u.ecs.nsidc.org/… Where "" is an Earthdata access token from your account. https://urs.earthdata.nasa.gov/documentation/for_users/user_token HTH (there are other methods like setting GDAL_HTTP_HEAD

Re: [gdal-dev] Re-post, user-defined

2024-04-29 Thread Liyuneh Gebre via gdal-dev
Well the arcpy made transformation before re-projection. Appended the sample script in the shared folder as txt. Anyway how can I achieve this ArcPy output re-projection from gdal with the given gdalinfo on the input file and by editing the PROJ.4 string (target projection )? PARAMETER["Latitude o

Re: [gdal-dev] Re-post, user-defined

2024-04-29 Thread Even Rouault via gdal-dev
At first sight, I'd say that the ArcPy output is incorrect. The SRS definition of its output doesn't match what you want since it has PARAMETER["Latitude of natural origin",20, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8801]], PARAMETER["Longitude of natural origin",6378137, ANGLEUNIT["d

Re: [gdal-dev] Re-post, user-defined

2024-04-29 Thread Liyuneh Gebre via gdal-dev
Hi Even, thank you for your response, gdal version is: GDAL 3.6.4, released 2023/04/17 I got the defined target projection included in the link. yes, it's geocentric and for further exploration the custom project transformation file .gtf extension is also included. https://u.pcloud.link/publink/sho

Re: [gdal-dev] Test problem After Merge

2024-04-29 Thread Even Rouault via gdal-dev
Andrew, Perhaps you're running updated GDAL Python bindings against a libgdal that hasn't been rebuilt recently? OSRIsDerivedProjected() has been added recently both in libgdal and the bindings Even Le 29/04/2024 à 23:27, Andrew Bell via gdal-dev a écrit : Hi, I just merged master into my

Re: [gdal-dev] HDF5 and geolocation arrays

2024-04-29 Thread Joaquim Manuel Freire Luís via gdal-dev
> This HDF5 (requires earthdata credentials your "Authorization: Bearer > " in GDAL_HTTP_HEADERS, or equiv) presents without geolocation arrays. gdalinfo "/vsicurl/https://n5eil01u.ecs.nsidc.org/AMSA/AU_SI12.001/2012.07.02/AMSR_U2_L3_SeaIce12km_B04_20120702.he5"; -sd 26 Excuse for this little

[gdal-dev] Test problem After Merge

2024-04-29 Thread Andrew Bell via gdal-dev
Hi, I just merged master into my branch and I'm now getting an error when trying to run a test. If someone might know what happened, I'd appreciate a tip. (gdal) [master] $ pytest -v autotest/utilities/test_gdal_viewshed.py ImportError while loading conftest '/Users/abell/gdal.2/build/autotest/co

Re: [gdal-dev] gdal2tiles for floating point images

2024-04-29 Thread Michael Sumner via gdal-dev
On Fri, Apr 26, 2024 at 5:37 AM lefsky--- via gdal-dev < gdal-dev@lists.osgeo.org> wrote: > I'd like to have a version of gdal2tiles that handles image types other > than uint8. Is there a reason why it doesn't handle those types of images? > It appears I'd have to modify the output format from pn

Re: [gdal-dev] HDF5 and geolocation arrays

2024-04-29 Thread Michael Sumner via gdal-dev
Oh wow, that's fantastic - thanks so much. Lots for me to learn from there too, with the HDF5EOS GRID details. Cheers, Mike On Tue, Apr 30, 2024 at 6:08 AM Even Rouault wrote: > Michael, > > actually support for those products was already there at 95% since they > use the HDF5EOS GRID formalis

Re: [gdal-dev] HDF5 and geolocation arrays

2024-04-29 Thread Even Rouault via gdal-dev
Michael, actually support for those products was already there at 95% since they use the HDF5EOS GRID formalism which we support since 3.7. But that particular type of products has an oddity. They have an empty GROUP=Dimension within the GROUP=GridStructure in the HDF5EOS metadata, which our

[gdal-dev] HDF5 and geolocation arrays

2024-04-29 Thread Michael Sumner via gdal-dev
This HDF5 (requires earthdata credentials your "Authorization: Bearer " in GDAL_HTTP_HEADERS, or equiv) presents without geolocation arrays. gdalinfo "/vsicurl/ https://n5eil01u.ecs.nsidc.org/AMSA/AU_SI12.001/2012.07.02/AMSR_U2_L3_SeaIce12km_B04_20120702.he5"; -sd 26 Driver: HDF5Image/HDF5 Dataset

Re: [gdal-dev] Reading interpolated values on DSM

2024-04-29 Thread Javier Jimenez Shaw via gdal-dev
Hi all Thanks for your answers Just to get some values, the VRT was enough to get an idea. Thanks Jukka. The questions was more like "this should be somewhere, and I do not find it", but apparently "it was not there". I see that GDALRasterInterpolateAtPoint could be useful (I am surprised it was

Re: [gdal-dev] Inquire mapped column names for conversion to shapefile

2024-04-29 Thread Jan Heckman via gdal-dev
That's pretty cool, you just saved me an awful lot of time! Thanks, Jan On Mon, Apr 29, 2024 at 4:31 PM Even Rouault wrote: > Hi, > > Something like: > > > from osgeo import ogr > original_field_names = [ "toolongforshapefile1", "toolongforshapefile2"] > map_to_shp = {} > tmpfilename = "/vsimem/

Re: [gdal-dev] Inquire mapped column names for conversion to shapefile

2024-04-29 Thread Even Rouault via gdal-dev
Hi, Something like: from osgeo import ogr original_field_names = [ "toolongforshapefile1", "toolongforshapefile2"] map_to_shp = {} tmpfilename = "/vsimem/temp.shp" ds = ogr.GetDriverByName("ESRI Shapefile").CreateDataSource(tmpfilename) lyr = ds.CreateLayer("temp") lyr_defn = lyr.GetLayerDefn()

[gdal-dev] Inquire mapped column names for conversion to shapefile

2024-04-29 Thread Jan Heckman via gdal-dev
Hi everyone, I want to make a map of original (let's say postgresql/postgis) column names to those used in a shapefile after conversion; assuming that at least some column names in postgis will be over 10 characters in length, and shortening may produce conflicts. Basically I know how this is done,

Re: [gdal-dev] Re-post, user-defined

2024-04-29 Thread Even Rouault via gdal-dev
Hi, it is difficult to help you with just the elements you've provided. What could be helpful is: - the output of gdalinfo on the input file - the output of gdalinfo on the output file generate by ArcPy - how you define the target projection in ArcPy - the output of gdalinfo on the output fi

[gdal-dev] Re-post, user-defined

2024-04-29 Thread Liyuneh Gebre via gdal-dev
Dear all, I am requesting again if you kind have some responds or direction, thank you I have difficulties in accurately reprojecting a raster file to a user-defined projection. I used the following commands in gdal while the output is not similar to arcpy results. kindly offer your guidance to ac

[gdal-dev] Reading DXF block attributes

2024-04-29 Thread Simon Gröchenig via gdal-dev
Hi all, I attempt to copy DXF data to GPKG using the following command: ogr2ogr --config DXF_MERGE_BLOCK_GEOMETRIES FALSE --config DXF_INLINE_BLOCKS FALSE -f "GPKG" D:/tmp/632_20-2N.gpkg D:/tmp/632_20-2N.dxf --debug ON My DXF contains some blocks with several block attributes. According to htt

[gdal-dev] GDAL 3.9.0beta2 available for testing

2024-04-29 Thread Even Rouault via gdal-dev
Hi, I've prepared a beta2 of GDAL 3.9.0 as there have been a few non-trivial changes since beta1 The major ones are: * pyproject.toml: use numpy>=2.0.0rc1 for python >=3.9 (#9751) * Docker: update ubuntu-small and ubuntu-full to 24.04 * TileDB: migrate away from deprecated APIs, and bump mini

Re: [gdal-dev] Total Viewshed and GDAL

2024-04-29 Thread Rahkonen Jukka via gdal-dev
Hi, Can this tool utilize digital surface model and deal with obstacles like buildings and trees? -Jukka Rahkonen- -Alkuperäinen viesti- Lähettäjä: gdal-dev Puolesta Luis Felipe Romero via gdal-dev Lähetetty: maanantai 29. huhtikuuta 2024 12.16 Vastaanottaja: gdal-dev@lists.osgeo.org

Re: [gdal-dev] Total Viewshed and GDAL

2024-04-29 Thread Luis Felipe Romero via gdal-dev
Hi, After several months working on other things, I have returned to an application of the total viewshed. What I have done is use Tamas Szekeres' C++ code as a template and, at least on Windows, I have managed to compile and run. With just a GPU it is quite fast: on a 6000x6000 raster it does

[gdal-dev] Transparent background in hillshade map

2024-04-29 Thread Elena Ruiz via gdal-dev
Hello, I am using GDAL version 3.6.2 to make a Hillshade map from an elevation map (tif) using the gdaldem command. My question is because I need the background of the image generated with the previous command to be transparent, I have tried different ways but it usually gives me this error: ERR