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
1answer
46 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
17 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
118 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
40 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
37 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
33 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
33 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
32 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
39 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
63 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
62 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
36 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
62 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
97 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
57 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
28 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
16 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
49 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
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
71 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
161 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
26 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
33 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
48 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
47 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
55 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
45 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
40 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
48 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 ...
6
votes
3answers
109 views

Variable window (“kernel”) size for raster processing

I'd like to use differently sized windows as I process (a.k.a. filter) a raster, with the size varying according to latitude. There are tools that allow a user to specify custom windows (e.g., ...
3
votes
1answer
36 views

gdal2tiles generating nosensical data

I'm trying to convert a global irradiance map - generated by r.sun - into tiles to be layered onto OpenLayers using gdal2tiles. However, every single one of my outputs is turning into white noise and ...
0
votes
0answers
26 views

changing shape file polygon values on the command line or in a library

I am using gdal_polygonize to create shapefiles. As gdal_polygonize only seems to work for integer values, I gave the different value ranges in my original raster data individual IDs. Now in the ...
1
vote
0answers
51 views

Storing non-orthorectified images in GeoTIFF

I'm reading GeoTIFF and GDAL documentations and I'm trying to understand what type of transformations are supported. GeoTIFF documentation says you can use ModelTransformationTag to define a ...
0
votes
0answers
16 views

Converting USGS (2010) GeoPDF to GeoTIFF

I'm having problems on coming up with Topo maps by converting USGS 2010 GeoPDF Maps to GeoTiffs (while removing the Orthoimage). I'm currently using GDAL in the conversion process. What happens is ...
0
votes
1answer
36 views

Converting Shapefiles to KML and extruding on height

I am working on a free/open source workflow to ingest OSM data and convert it to a KML so I can view it in 3D with Google Earth. I used ogr2ogr to convert to a shapefile and have extracted the height ...
3
votes
0answers
89 views

MODIS R Package: GDAL Not Available

I am trying to analyze the MODIS EVI 250m, 16day dataset. I found the R MODIS package and am following the instructions here: ...
2
votes
1answer
51 views

How to display this map in QGIS

I downloaded a map with a .tif, .tfw, and .xml file. I can't get it to display in QGIS. I used geotifcp -e to create a GeoTIFF. I tried Lambert Conformal Conic since it was mentioned in the .xml. Hard ...
1
vote
0answers
18 views

Local Maxima (or Pour Point) Pixel Value extract in Raster

I am doing watershed & pourpoint using python GDAL. I have calculated Fill, Flow-Direction & Flow-Accumulation using Python & GDAL. The next is to get Pour Point in certain window ...
0
votes
0answers
25 views

intersect geotiff with bbox

In a project using gdal I need to crop a geotiff file with bbox in pixel or coordinates epsg:4326. This step is ok. The next step is to manage cases where the image size is "lower" than the bbox size ...
0
votes
0answers
17 views

mosaic/reprojection/conversion order of operations

I'm writing a python/gdal program that ingests HDF4 files and: 1) Converts subdatasets from HDF4 to GTiff using gdal_translate 2) Mosaic's the new Gtiffs using gdal_merge 3) Defines the ...
2
votes
0answers
51 views

Is there a good example of using GDAL_WMS to add an authenticated TMS layer to QGIS?

I've been trying to add a TMS layer in QGIS using a GDAL_WMS xml configuration file. I added UserPwd to the xml and this seems to provide the authentication required by GeoServer but I'm struggling ...
0
votes
1answer
43 views

Applying Attribute Filter and Spatial Filter simultaneously (GDAL - Python)

quick question on something that has to be simple, but I'm getting stuck on. I want to use an Attribute Filter and a Spatial Filter simultaneously. I have two vector files (an overlaying grid and a ...
2
votes
1answer
32 views

Transformation from EPSG 2178 to EPSG 4326 gml.file?

How can I transform .gml file from EPSG 2178 to EPSG 4326? I used the following without success: gdaltransform -s_srs EPSG:2178 -t_srs EPSG:4326 "file1.gml"
0
votes
0answers
71 views

Error installing GDAL for Python on OSX

I'm trying to install GDAL for Python on my mac (OSX 10.9.2). I did: >>> brew install GDAL and got: >>> 🍺 /usr/local/Cellar/gdal/1.11.3_1: 230 files, 33.1M now: ...
1
vote
1answer
157 views

Installing QGIS using Homebrew

I've been unsuccessfully trying to install QGIS using Homebrew (used Homebrew to install GDAL and GSL) as outlined here My OS is Mac Yosemite, 10.10.5. The installation sort of finishes, but with a ...
0
votes
1answer
19 views

How to set coordinate dimensions in Django GDAL API?

The Django GDAL documentation for coord_dim says: coord_dim Returns or sets the coordinate dimension of this geometry. For example, the value would be 2 for two-dimensional geometries. But how do ...
0
votes
0answers
55 views

Converting lat/long to image pixel coordinate and vice versa using gdaltransform with image and its RPC file?

I am a new user of gdal. I am having a image (.tif) and its corresponding RPC file (.txt) How to convert lat/long to image pixel coordinate and vise versa using gdaltransform with the image and its ...
1
vote
0answers
52 views

error using gdal_translate with -projwin option

I have to crop geotiff image, I'm using gdal assembly and I'm trying to use the gdal_translate method If I crop my image using the -srcwin option all is fine gdal_translate -srcwin 0 0 400 400 ...