GDAL (Geospatial Data Abstraction Library) is an open source translator library & set of command-line utilities for processing raster spatial data formats. As a library, it presents a single abstract data model to the calling application for all supported formats. It also comes with a variety of ...

learn more… | top users | synonyms

5
votes
2answers
87 views

Why am I getting two different outputs from Python and C++ using GDAL?

I am trying to gather information from a *.tif file using the GDAL libraries. When I open the *.tif file using both the Python API, and C++ API, I encouter two completely different answers. Here are ...
1
vote
0answers
26 views

Georeferencing Himawari-8 in GDAL (or other)

Does anyone know the appropriate project information to georeference (and hence reproject) Himawari-8 AHI data in the IR? The images are 5500x5500 pixels. Currently I try these: ulx=-5570248.832537 ...
1
vote
0answers
64 views

How to make NoData transparent?

I have been using python to turn a netcdf into a geotiff. Once I have the geotiff, I am using gdalwarp to transfer it into the right projection. Then use gdal_translate to get the tif into a vrt ...
1
vote
2answers
82 views

Translating ArcPy Commands to GDAL Equivalents? [closed]

I'm trying to figure out GDAL equivalencies to the following arcpy functions: arcpy.env.overwriteOutput() arcpy.MakeFeatureLayer_management() arcpy.SpatialJoin_analysis() arcpy.da.SearchCursor() ...
1
vote
0answers
14 views

Is there a Guide to Translate ArcPy code to GDAL equivalent? [duplicate]

Due to cost-cutting measures I've been asked if I can migrate ArcPy code to a GDAL equivalent. I'm fairly new to GIS, ArcPy, and GDAL. Is there a website or guide which caters to providing ...
1
vote
2answers
59 views

Reproject raster without projection info

I'm pretty new to gdal and got a presumably easy question for experienced gdal users. My problem: I have a raster ascii file of which I don't know the projection or epsg code. However, I want this ...
1
vote
0answers
25 views

Gdalwarp UTM reprojection causes image distortions

I'm trying to re-project Sentinel-2 images to a common UTM zone in order to create a mosaic. For this I'm using gdalwarp with the following syntax: C:/OSGeo4W64/OSGeo4W.bat; gdalwarp -t_srs ...
1
vote
0answers
32 views

gdalwarp won't run with “cubic” resampling method

I am having an issue with running the gdalwarp command from command prompt. Basically I have a raster WGS84 file, which I would like to reproject to spatial reference which is using the Azimuthal ...
2
votes
1answer
31 views

Tiff: How to calculate Tiff Tag ModelPixelScaleTag (field 33550) and Field 34735

I have to manage a service which can return images in tiff format. This service is requesting another WCS Service from which it receives the Image File without geo-referenced informations. So i need ...
1
vote
0answers
56 views

mosaic rasters with gdal_merge.py

I want to mosaic 2 rasters "rasa" and "rasb" using gdal_merge.py . While running the code, it produces an error " No input files selected". Can somebody help me to find the error in the code? ...
0
votes
0answers
20 views

Installation and execution of GDAL

I am trying to install GDAL on the remote server. I do: pip install GDAL==1.9.1 It gets installed and then I try to do: gdalinfo --version But I get 'command not found'. I understand that I ...
1
vote
2answers
44 views

Create GeoJSON after processing.runalg('qgis:clip', …)

I have the following problem: I have two vector layers, foobar and clipLayer. About the python console I call the following command: processing.runalg ('qgis:clip', 'foobar', 'clipLayer', None) As a ...
1
vote
0answers
49 views

Troubling using MODIS package in R (mac). Error with gdal

I'm new to using MODIS data, and have been trying to set up a a way to use the HDF files produced by MODIS on R. I had trouble with gdal (could only download the package through system dialog), and ...
1
vote
0answers
41 views

How do I calculate the percentage of progress using raster data using Python?

How do I calculate the progress amount of two or more raster data being processed? I have seen examples of vector data but it becomes difficult to calculate the progress when no vector data are ...
4
votes
0answers
17 views

Is there a way to use the GDAL C++ API to find areas of overlap between two SAR images?

I have a bunch of SAR images (about 1 GB each) on which i need to do some image processing. Each image is georeferenced and there is significant overlap between each pair of images in the set ie. one ...
0
votes
1answer
27 views

How can I blend DEM rasters with GDAL?

I have several DEM rasters in different projections and at different resolutions. The rasters are not rectangular, but they have appropriate NoData values. They overlap enough that I should be able ...
0
votes
1answer
26 views

Install pygeoprocessing on ubuntu 14.04 and python3

I'm trying to install the pygeoprocessing package for Python 3.4 on an Ubuntu 14.04 VM. I have all the necessary versions of the required dependencies (GDAL, numpy, scipy, shapely, Cython). ...
0
votes
0answers
15 views

Keep NoData values when using gdal_merge

I'm using gdal_merge to join several BioClim rasters, which I will then feed into a Java program that "snaps" points to the nearest BioClim raster tile of non-NoData value. However, when I run ...
0
votes
1answer
26 views

use of gdaltransform in C# code

I have a problem to use correctly gdaltransform in a C# project to transform a set of coordinates from a crs to another. If first try in console command and this line is working good my inpufile ...
0
votes
1answer
38 views

How to use gdal_rasterize with postgis vector?

I'm trying to convert a multipolygon vector layer to a raster using the GDAL tool gdal_rasterize. The command I am using is: gdal_rasterize -b 1 -burn -9999 PG:"dbname=db password=password" \ ...
0
votes
0answers
12 views

ERROR: permission denied to set parameter “postgis.gdal_enabled_drivers”

I am trying to enable raster drivers and out-of-db raster storage for PostGIS 2.2+ using GUC (Grand Unified Custom) Variables. In particular, I was trying to set postgis.gdal_enabled_drivers following ...
1
vote
1answer
75 views

ImportError installing GDAL on Ubuntu 14.04 with Anaconda Python 3.5

I installed the GDAL library on Ubuntu 14.04 through the Anaconda (Python 3.5) distribution. After installation, which seems to work fine, I try importing it and I get this error: >>>from ...
0
votes
0answers
28 views

How to clip MrSID Files in R and export to GeoTiff?

I would like to use R to: Read in an MrSID raster (.sid format) Clip the extent using a shapefile Save the new clipped raster as a geotiff I am fine with steps #2 and 3, but I can't find any ...
1
vote
1answer
35 views

gdalwarp output shifted 15 KM

This is similar to this question regarding shifts in raster output, but the difference I'm experiencing is fairly large-- the output raster is shifted 15 KM to the south. Here's an image: When ...
0
votes
0answers
30 views

ogr2ogr - PostGIS to GML with nested attributes

I am attempting to create GML from PostGIS tables and having issues with nested tags. Before I waste too much time on this, I would like to know if what I am trying to do is possible with "ogr2ogr". I ...
1
vote
1answer
26 views

gdalwarp causing shift in pixels

Can somebody help me with gdalwarp issue? I am having an input raster and would like to clip it with polygons. After running the below code, while overlaying output raster on input raster, i see there ...
1
vote
1answer
30 views

How to get the pixels from a geotiff file in gdal python for a given polygon?

I use python and GDAL under apache to make a webservice. I have a GeoTIFF file and i want only a subset of the data, everything inside one polygon defined by a buffer of a line (coord is the input ...
0
votes
1answer
36 views

Converting files with ogr2ogr and dropping ALL attribute fields

I want to add a gpx file to a PostgreSQL/PostGIS table using ogr2ogr, but my table does not have any of the fields that are attributes of the GPX file. There are previous questions here dealing with ...
0
votes
1answer
53 views

Is there a tool to shrink/grow raster data edge by n pixels

Looking for an open source tool to grow/shrink raster edge by n pixels. I noticed Grass has r.grow but hoping there might be a way to do it with GDAL, PostGIS or python library. More information: I ...
-1
votes
1answer
60 views

The reason why dataset does not create layer

I have a plugin written in Python, and two versions of the QGIS, assembly there are differences, but the main suspect, that the version of GDAL==2.0.0dev (1) another 1.11.3 (2) for first ...
1
vote
0answers
31 views

ogr2ogr (or other) simplification of polygons with shared boundaries

I am trying to simplify a collection of 100,000 or so polygons containing nearly 4 million vertices. The data is stored in a 2GB shapefile so there is no explicit representation of the topology ...
3
votes
1answer
57 views

Convert Japan Administrative Boundary Data to Shapefile

I've downloaded Administrative Boundary data from the Japan MLIT (Ministry of Land, Infrastructure and Transport), but it is in an XML based format referred to as JPGIS or KSJ. Are there any scripts, ...
3
votes
2answers
87 views

How to remove parts of raster layer with QGIS (equivalent of difference tool but for raster)?

I'm trying to get the difference between two raster layers the same way we can do it with the difference tool with vector layers. I'm using QGIS 2.14. Anyone has already done that and can tell me how ...
0
votes
1answer
51 views

GDAL 2.0.1/QGIS 2.14.0 on OSX won't open MODIS HDF?

I'm on OSX 10.9.4, using GDAL 2.0.1 from the terminal. I want to mosaic MCD12Q1 HDFs together with GDAL but my build can't even read them. The file: ...
0
votes
0answers
27 views

How can I connect QGIS with compiled GDAL on linux?

I've compiled GDAL on Debian with FileGDB support. The output of ogr2ogr --formats | grep GDB is OpenFileGDB -vector- (rov): ESRI FileGDB FileGDB -vector- (rw+): ESRI FileGDB and gdalinfo ...
1
vote
0answers
12 views

How can I read the TIFFTAG_DATETIME using GDAL in Python

I am currently using GDAL from Python to read and transform coordinates and that works well. I would now like to be able to read some of the none geographic tags, for example, TIFFTAG_DATETIME using ...
0
votes
0answers
48 views

Create raster with decibel value from Geotiff Intensity using QGIS GDAL python

I have raster with intensity values and I want to create raster with decibel values, I work with QGIS GDAL, I use gdal_calc I do this from the terminal: gdal_calc.py -A input_image.tif ...
0
votes
0answers
26 views

Converting PostGIS data to image [duplicate]

I have an image showing the height information of a map. I have created a table using the 'raster2pgsql' command. This table stores the height values of the map. Now I want to create an image out of ...
3
votes
1answer
67 views

Understanding intersection of GDAL and QGIS?

Background I'm starting a project that will be done entirely through Python, except perhaps at the end for visualisation. So far I've been doing it with ArcGIS and arcpy, but we'd like to move it to ...
3
votes
1answer
149 views

transpose large netCDF with GDAL, NCL, NCO or CDO

I have a large netCDF file (444GB) that I would like to transpose. The file is flipped with latitude on the x-axis and lon on the y-axis. See pictures of one band below. Not correct Correct: I ...
0
votes
0answers
23 views

Best file format to manipulate raster data with big length and small width

How is the best way to manipulate data from specific distance from an axis? Problem: I have a road axis; I want to manipulate 1Km of data from each side of this axis; Pixel size is 1m X 1m Axis Length ...
2
votes
0answers
35 views

How to use GDAL RasterizeLayer (python binding) -where flag

I am trying to use GDAL RasterizeLayer python binding, but I can not find a way to use the -where flag. The 'pure' GDAL command I am trying to run is: gdal_rasterize -b 1 -burn 0.00001 -where ...
0
votes
0answers
26 views

How to improve poor image quality of “gdalwarp” operation result?

I'm trying to improve the result of gdalwarp operation. The image width and height are 6732x4721 GCP's are: -gcp 0 0 38.005301407753 55.048725423557 -gcp 6732 0 38.184943751144 ...
1
vote
1answer
45 views

Select by attribute using fiona

with fiona.open(constants.file_combined, 'r') as source: filtered = filter(lambda f: f['properties']['NAME'] == name_dm1, source) # **source.meta is a shortcut to get the crs, ...
3
votes
0answers
40 views

Gdal and module import error when creating an .exe application with pyinstaller

I’ve been struggling since three days in a row to create a windows executable application from my python application. I have managed to get results but I’m not at all satisfied with the workarounds I ...
2
votes
2answers
35 views

Assigning column names dynamically

I have an interesting problem: I need to batchprocess a few files. I do that in a for-loop, where some operations are performed. Later on I want to write the outcome into an existing shapefile which ...
3
votes
1answer
52 views

Table of statistics of rasters input in QGIS

I work with QGIS GDAL,I have many raster image 'geotiff' with one band for many dates, I want to create a table that contains statistics for all rasters that I have. I tried the algorithm 'raster ...
1
vote
1answer
38 views

gdal/python WriteArray working for 8-bit but not 16-bit

I'm trying to use GDAL/python API to open a raster, do some math, then write the two resulting arrays as rastes. One is to be an 8-bit raster and the other 16-bit signed. However, only the 8-bit ...
0
votes
0answers
28 views

SharpMap.Layers.GdalRasterLayer threw an exception

I am trying to load a raster layer from PostGIS database. I used this code var rlay = new SharpMap.Layers.GdalRasterLayer("demelevation", "PG\":dbname=Bahram host=localhost user=postgres ...
1
vote
0answers
39 views

gdal_merge.py is not working with NODATA values

I am trying to merge a couple of GeoTiff files with gdal_merge.py. Example input: In the input data the nodata value is -32678. Hence I tried: gdal_merge.py -n -32678 -o DEM.tif -of GTiff ...