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

0
votes
0answers
43 views

Relation between contour_interval and DN values of a georeferenced image in gdal.ContourGenerate() function

I am generating contours from a georeferenced image by using gdal.ContourGenerate() function in python. One of the parameters of the function is contour_interval. The value of contour_interval ...
1
vote
1answer
21 views

Georeferencing files that are orthogonal to a point on a sphere

I have about 65000 photos, each of them a square where an inscribed circle is radius 1 degree. Each of these photos has an associated latitude and longitude of the center. I'm trying to figure out ...
1
vote
1answer
26 views

QGIS GdalUtils runGdal function throws error for second argument

I am trying to use a script provided in the QGIS Python Cookbook by Joel Lawhead - its a simple script that georeferences images based on control points via the gdalwarp function. (p.124 in the book) ...
1
vote
0answers
44 views

Efficiently iterating through raster calculations with Rasterio or GDAL

I've got a very slow and round-about way to currently do this, but I struggle with the general concept of using windows in raster calculations and how to apply an affine or geotransformation between ...
6
votes
1answer
132 views

Adding GDAL as dependency to Python package?

I am trying to package a Python script that uses GDAL for PyPI. I started by including a direct reference in my setup.py: install_requires=['GDAL==1.11.2'], This way the package failed to install in ...
2
votes
1answer
47 views

MODIS data with .hdf extension not compiling in Ubuntu

I am trying to use GDAL to read .hdf file, but it does not seem to recognize the extension. It works fine in Windows system. The GDAL open command gives the following error message: ERROR 4: `...
5
votes
2answers
110 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 ...
4
votes
1answer
107 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
79 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 file....
1
vote
2answers
96 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
15 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
137 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
44 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 \"+proj=...
1
vote
0answers
39 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
69 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
1answer
157 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
27 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
55 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
78 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
52 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
26 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
51 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 ...
1
vote
1answer
38 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
32 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
37 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
82 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
30 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
234 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
64 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
40 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
43 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
46 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
61 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
50 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
73 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
69 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
58 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
72 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
140 views

Removing parts of raster layer with QGIS (equivalent of vector difference tool)?

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
80 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: http://e4ftl01.cr.usgs.gov/MOTA/MCD12Q1.051/2001....
0
votes
0answers
39 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
1answer
36 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
56 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 --calc="10*...
0
votes
0answers
27 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
79 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
197 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
35 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
37 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 'DN=...
0
votes
0answers
47 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 55....
1
vote
1answer
63 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, driver, ...