Questions tagged [gdal]

GDAL (Geospatial Data Abstraction Library) is an open source translator library for raster and vector geospatial data formats.

Filter by
Sorted by
Tagged with
0 votes
0 answers
13 views

NED1 data and GDAL

I have a full dataset of the older USGS data that was formerly referred to as "NED1". (One arcsecond data, 1 degree tiles, raster of 3612x3612 floats in a .flt file, .hdr and .prj files ...
Dan H's user avatar
  • 101
0 votes
0 answers
5 views

Parsing complex GML with GDAL

We have some complex GML that looks like <?xml version="1.0" encoding="UTF-8"?> <abpr:AddressBaseSupplySet xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance&...
tjmgis's user avatar
  • 3,320
0 votes
0 answers
13 views

Resize NITF while maintaining metadata

What I want to do: I'm working on enhancing a legacy NITF reader in C++. Before I touch it, I want to create a battery of unit tests to make sure I don't break things while I'm trying to improve them. ...
user55937's user avatar
  • 1,281
0 votes
0 answers
25 views

GDAL - Unable to open PDF or SVG

I'm trying to import a PDF of a floor plan of a school into QGIS. I was able to open this PDF in QGIS georeferencer and it recognized points, lines, and polygons and I was able to georeference these. ...
Evan S's user avatar
  • 1
0 votes
0 answers
22 views

gdal_proximity - setting processing extent using -te parameter in QGIS not working

the default processing extent of many tools in QGIS is the input layer's extent. GDAL Proximity (Raster Distance) is no different. According to the documentation (https://gdal.org/programs/...
Guziko's user avatar
  • 1
0 votes
1 answer
25 views

gdal_retile.py - Set Resolution of finest pyramid

I have a raster that i want to publish using geoserver. The rasters resolution is 10x10m. However even the smallest tiles in the targetDir have 20x20m resolution and thus my wms has a resolution fo ...
Felix's user avatar
  • 165
0 votes
0 answers
33 views

Opening .png.aux.xml file in QGIS

I received a .png.aux.xml file. I have tried to read it with the GDAL translator and by making a URL from the code. I haven't been able to open the file in QGIS. Any ideas?
Pedro Castro Gomez's user avatar
1 vote
0 answers
32 views

Rasterio convert vertical datum of NAVD88 to WSG84

I've been asked to provide an elevation data for projected to GPS. To achieve that i've tried to convert the NAD83 geotiff with NAVD88 vertical datum (downloaded from USGS) into WSG84. To do the ...
Szymson's user avatar
  • 11
0 votes
0 answers
29 views

Pred and zlevel arguments from deprecated gdalUtils R library mosaic_rasters function [closed]

I'm trying to run some older R code that relies on the now deprecated gdalUtils library. I'm trying to find a way to run the following function in an equivalent way using a library like raster or sf. ...
Samuel T's user avatar
1 vote
1 answer
36 views

How to use the 'sqlite' dialect in ogr2ogr queries when one of the GeoJSON property fields is "id"

I have an example GeoJSON, and I would like to use it as input to an ogr2ogr command. I aim to generate a new file with positive integer feature ids, as described here. The example GeoJSON, src....
zabop's user avatar
  • 1,810
0 votes
0 answers
32 views

Opening Esri .gdb database with values coded or stored in domains

I have a .gdb file with multiple layers of different geometries and a table called "Toolbox". When i import it into QGIS it shows me all the layers and the "Toolbox" as a table. ...
m4D_guY's user avatar
  • 21
1 vote
0 answers
44 views

How to combine ODA File Converter with OGR (GDAL) so that dwg is a format option

I have downloaded ogr/gdal with QGIS 3.34.1 and I have downloaded the ODA file converter from https://www.opendesign.com/guestfiles/oda_file_converter However, I am not sure how to combine the two so ...
guest_231213's user avatar
0 votes
0 answers
32 views

Strange result from GDAL in python from ArcGIS Pro

I am in the process of converting my Desktop scripts to be Pro compatable, one script uses GDAL to calculate a sum of overlapping rasters generated previously in the script (because Cell Statistics ...
Michael Stimson's user avatar
1 vote
0 answers
32 views

Writing COG using GDAL C++ API fails with block order warning

I'm using the GDAL C++ API to write a COG which is constructed by the rest of the code. If I use the Geotiff driver, there are no problems. But when I switch over to the COG driver, I'm not able to ...
goalie39's user avatar
0 votes
1 answer
31 views

GDAL /vsicurl/ does not use range requests

I'm trying to extract a small part of a largish COG GeoTIFF using gdal_translate. The source file is stored on S3 and I have a presigned URL to it, and use /vsicurl/ to access it. For whatever reason, ...
Liedman's user avatar
  • 1,073
2 votes
1 answer
47 views

Why doesn't `gdalinfo` with `json` flag report searched but not found attributes?

Motivation I have satellite .tiff data files. I want to understand what information exists in these files. When I run gdalinfo on them, with the -json flag, I don't get any null values from the output,...
Louie_the_unsolver's user avatar
0 votes
0 answers
26 views

gdal_fillnodata strange artifacts

I'm seeing an odd artifact when using gdal_fillnodata.py. I have a raster (purple is no data, white and black are two different numeric values). Using a mask exactly corresponding to the no data ...
Stev_k's user avatar
  • 6,679
1 vote
1 answer
132 views

Create copies of rasters in R via GDAL

I am looking for an efficient way to create copies of rasters in R via GDAL. Since the advantages of such an approach might not be directly obvious, here is some more context. Via GDAL we can read ...
darius.goergen's user avatar
1 vote
1 answer
26 views

Zoom level with gdaladdo

I am trying to convert a .las file into mbtiles. Here are the steps I followed: Convert .las file into .tif in cloud compare. Reproject the .tif file with gdalwarp. Creating the tif file into ...
Tahmid's user avatar
  • 21
1 vote
0 answers
24 views

Loading large Rasters into PostGIS

I've been trying to load large Rasters into PostGIS using raster2pgsql command, however since it executes under one process, it always seems to be bound by the CPU of the vCore it's scheduled under. ...
Nishant Sondhi's user avatar
3 votes
1 answer
91 views

QGIS 3.28/GDAL Unable to export ESRI file geodatabase layer to PostgreSQL

In the past, using the OSGEO4W Shell installed as part of QGIS, I was able to run an ogr2ogr command to transfer a feature layer from an ESRI File Geodatabase to a table in PostgreSQL. However, after ...
Mike Smith's user avatar
0 votes
3 answers
266 views

Installing GDAL for Python

I'm looking to install GDAL for Python. The first thing I did is installing Python, then I install OSGeo including its library GDAL, then in CMD when I type pip install gdal I get this error ERROR: ...
hally-mio's user avatar
1 vote
0 answers
21 views

Clipping WMS with gdal_translate [closed]

I want to clip this WMS https://www.ign.es/wms-inspire/pnoa-ma?SERVICE=WMS& to a specific extent. I have tried what was suggested in these two posts, but they did not work: Trying to use ...
Isa's user avatar
  • 175
0 votes
1 answer
38 views

How do I set the result of a SELECT query equal to a new column in a GPKG file using ogrinfo?

This thread: How do I add a new column to a GPKG file using ogrinfo? explains how to add a new column to a GPKG file and set all rows in it to a specific constant value. Here, I would like to add a ...
zabop's user avatar
  • 1,810
0 votes
0 answers
39 views

Using GDAL and OSR to convert from latitude/longitude to DEM coordinates

I'm trying to convert back and for between DEM and "world" (aka latitude/longitude) coordinates using GDAL and I think I might be misunderstanding something which is causing me to get "...
Andrew's user avatar
  • 101
0 votes
0 answers
22 views

Translating (cropping) a GeoTIFF file to a specific extent but ending up with a different final extent

I'm trying to crop a larger GeoTIFF file to a smaller extent. The darker GeoTIFF file in the middle has the extent I want to reproject to. I've used the following command to do the translation: ...
Jamess11's user avatar
2 votes
1 answer
93 views

Why can ogr2ogr upload data to a not previously existing table faster than to an existing empty one?

Recently I observed that it takes much less time for ogr2ogr to upload data to a PostGIS database if I do not create the destination table first. I would like to understand why this happens. Below is ...
zabop's user avatar
  • 1,810
0 votes
0 answers
14 views

Produce ISEA "Icosahedral Snyder Equal Area" GEOTIFF/raster with gdalwarp

Is an icosahedral projection impossible to produce with PROJ/GDAL ? Produce ISEA "Icosahedral Snyder Equal Area" GEOTIFF/raster with gdalwarp. I am trying to use gdal tools to convert Plate-...
Molly's user avatar
  • 35
4 votes
1 answer
48 views

ST_CLIP with one raster and one polygon

I don't understand how I can select a raster and polygon by ID and then clipping them. Something like: SELECT ST_CLIP(selected_raser.rast, 1, selected_polygon.geometry) From raster WHERE ID = 99 ...
Squeezie's user avatar
  • 240
0 votes
0 answers
43 views

Reproject pseudo Plate-Carree to "Icosahedral Snyder Equal Area" using GDALWARP

Reproject pseudo Plate-Carree to "Icosahedral Snyder Equal Area" with gdalwarp. I am trying to reproject a plate-carree to ISEA "Icosahedral Snyder Equal Area". I suspect ISEA ...
Molly's user avatar
  • 35
3 votes
0 answers
49 views

How to close a GDAL PostGIS connection?

This issue has been bugging me for years, desperation got the better of me so now I'm here to seek the help of the public. I'm storing raster data in a PostGIS database, and I'm using GDAL from Python ...
Jeroen's user avatar
  • 31
0 votes
0 answers
107 views

Correctly merge Landsat images in a median mosaic using different tiles

I have several Landsat 8 and 9 images, from two different tiles that I want to merge together in a median mosaic. The tiles are tiles that overlap, with two images for each tile. What I'm doing to ...
rrwork___'s user avatar
1 vote
0 answers
31 views

ogr2ogr MSSQL to gpkg encoding issue

I am trying to write data from SQL server to geopackage but I have met some encoding issues. Warning 1: Value of field 'IDENTIFICATION' is not a valid UTF-8 string. I dont understand where in the ...
geogrow's user avatar
  • 1,663
0 votes
0 answers
18 views

Create GDAL dataset from scattered points in python

I have a set of points in XYZ format. As they are not gridded data (no guarantee), I want to use GDAL's algorithms to create a regular grid out of the scattered data points. I created a minimal ...
Aristotelis's user avatar
1 vote
1 answer
45 views

Using GDAL (MaxRev.GDAL) to write TAB files with mixed geometry

I'm migrating from calling native MiTAB operations for exporting features to a MapInfo TAB file to using GDAL and the MapInfo driver. I'm using C# and MaxRev.Gdal.Core package (and MaxRev.Gdal....
Mpalle's user avatar
  • 11
1 vote
0 answers
27 views

Double geometries in GML with different EPSG using GDAL

In our GML we have two geometries (a British National Grid and a ETRS89) so to try and extract these we use a GFS file with the following additional XML to extract them. <GeomPropertyDefn> <...
tjmgis's user avatar
  • 3,320
0 votes
1 answer
28 views

GDAL Geo-Transform incorrect after I using Register Raster of ArcMap

I have a .tiff raster and a .shp vector. And I found they are not aligned after loading them into ArcMap. So I used the registration function of ArcMap to add control points to align the raster with ...
hoca's user avatar
  • 3
1 vote
2 answers
57 views

Use gdal_translate on an orthographic PNG image of the Earth centered at Lat=+50N Lon=+5E to proper orthographic GeoTIFF format

How do I gdal_translate an orthographic PNG image of the Earth centered at Lat=+50N Lon=+5E to proper GeoTIFFf format ? I cannot work out how to specify the limits of a full disc. (diameter of Earth ...
Molly's user avatar
  • 35
0 votes
2 answers
44 views

Enlarge numpy array to fill/pad bounds in rasterio whilst keeping existing data intact

I'm looking to do what this question and answer is doing using gdal_translate / gdal_warp, but in Rasterio. Does anyone know of a way? I can't find anything in the rasterio.warp.reproject module that ...
Stev_k's user avatar
  • 6,679
1 vote
1 answer
53 views

Find coordinates corresponding to max and min values using Python

I'm using GDAL and numpy so far in my Python program. I am able to get the max and min values of a band within a .grib2 file with: def getDataStats(grib2Filepath, bandNr): local_minimum = None ...
interwebjill's user avatar
1 vote
1 answer
23 views

How to assign a layer to a group in GDB datasets using gdal python

I'm currently working with ESRI GDB datasets and I can easily open hierarchical geodatabases and walk the tree using the API that gdal provides. root_group = ds.GetRootGroup() for g in root_group....
Alireza S.N's user avatar
1 vote
0 answers
44 views

Is there any way to merge multiple .cog files into one using python?

I need to merge several .cog files into one. I have tried using rasterio, rioxarray, gdal, etc but I have not been able to get it to work. The closest I have been able to get was through this function ...
Guzmán R's user avatar
0 votes
0 answers
36 views

How to load Parquet driver in GDAL with c# wrappers?

I am using .Net nugets for GDAL and GDAL.Plugins. In GDAL documentation, in order for Parquet drivers to work, GDAL needs the Apache Arrow C++ library (https://gdal.org/drivers/vector/parquet.html). ...
gktval's user avatar
  • 1
0 votes
1 answer
40 views

Converting GML to GeoJSON - offset is out of bounds

I need to convert a GML file to GeoJSON from this website (https://use-land-property-data.service.gov.uk/datasets/inspire/download) and store it in my MongoDB database. I have read a few posts on here ...
GentleCynic's user avatar
0 votes
1 answer
35 views

Using gdal_translate to remove 4th band from command line

I have quite a few geotiff files that I need to process with gdal_translate to remove the 4th (alpha) band. They are all in the same directory (for example, /imagery). The original names are in the ...
Stevenj279's user avatar
1 vote
0 answers
24 views

Cropping VRT with Shapefile to create new VRT with no dependecy from input

I'm new to work with gdal and I'm trying to crop a VRT-File with a shapefile to create a new VRT-File that is independent from the input. I use gdalbuildvrt [...] to create the base VRT-File from ...
Martin Reger's user avatar
2 votes
1 answer
94 views

How to use ogr2ogr's spat_extent keyword to clip a layer of a GPKG file?

I am trying to clip geometries within a GPKG file using a bounding box. I write an L-shaped line to a file: import geopandas as gpd import shapely.geometry l = shapely.geometry.LineString([[0,2],[0,-...
zabop's user avatar
  • 1,810
1 vote
1 answer
50 views

What 'options' are available in osgeo.gdal.ContourGenerateEx?

I am working on generating contour polygons with gdal in python. I have a working script, which utilizes the osgeo.gdal.ContourGenerateEx function from the gdal library. However, I am having trouble ...
MarcusR's user avatar
  • 113
0 votes
0 answers
46 views

Reprojected image is rotated, can't match it with DEM

Why would reprojected sat image be rotated? I can't match it with source DEM map, if there is rotation. gdalwarp -t_srs EPSG:31287 input_satmap.tif output_reprojected.tif I use that line in gdal. ...
Alex J.'s user avatar
  • 119
1 vote
1 answer
62 views

Change projection from Unknown to given EPSG with GDAL

I've got a TIFF in projection (proj4: +proj=stere +lat_0=90 +lat_ts=60 +lon_0=-80 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs). Would like to convert it to EPSG:3857. ArcMap does it correctly, but ...
Flash Thunder's user avatar

1
2 3 4 5
124