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
141 views

GeoTiff image not aligned after using gdal_translate and gdalwarp

I don't see what I am doing wrong. The coordinate system was assigned to the GeoTiff using the following gdal_translate command: (Usage: gdal_translate -of GTiff -a_srs EPSG#### -a_ullr upper left ...
0
votes
1answer
139 views

Building GDAL with libcurl support on Debian Linux

Does anyone know how to build GDAL with libcurl support on Debian Linux (specifically LinuxMint and/or Ubuntu)? After following build instructions here and tip from @EvilGenius below, I'm specifying ...
1
vote
0answers
27 views

How Do I Split a Raster into parts using GDAL? [duplicate]

I have a large Raster with High resolution DEM of a large area, and I need to split it into smaller squares of 500mx500m, so that I can run other raster functions on it. Basically I'm looking for the ...
4
votes
0answers
249 views

How to remove nodata values from XYZ output with GDAL?

I have a buffered line route that is overlayed on a DEM: I've clipped the DEM to the buffered line, and exported it as an ADF ESRI grid file. I then used gdal_translate to convert it from the Arc/...
4
votes
2answers
1k views

Convert multi-layer DWG into several KML or geoJSON files (one per layer)

I have some big DWG files, with 5 to 50 layers, which I draw and geolocated using AutoCAD (with some external modules). I need to convert them to KML or geoJSON so that the JS library can load them ...
0
votes
2answers
249 views

Gdal driver.Create TypeError

I have been writing a python script to convert an ascii grid into a geotiff file and when running the Create method in gdal i am receiving the following error message: william@will-HP-EliteBook-850-...
1
vote
0answers
117 views

WKID issue when adding raster layer from ArcGIS “mapserver” to QGIS via REST

QGIS docs here describe how to do this. I am trying to add a layer from this MapServer: gps.digimap.gg/arcgis/rest/services/JerseyBaseMap/MapServer When I try the following command from the Python ...
2
votes
1answer
164 views

HDF5 and georeferencing?

I'm doing quite heavy processing of Earth Observation data, and it appears that a lot of my workflow can be speeded up by using HDF5 for local storage. I mostly work with large stacks of raster files, ...
3
votes
1answer
65 views

Display or export gdal_translate exponential scaling equation

I have an GeoTiff raster that contains elevation information at each pixel. As OpenLayers 2 can't display GeoTiff, I convert it to a PNG image and have to apply an exponential scale so the resulting ...
1
vote
1answer
998 views

SVG To Shapefile [duplicate]

Is it possible to draw out some SVG shapes in Illustrator or Inkscape and convert them to Shapefile using GDAL or QGIS (or something else perhaps?). I have an image that I can easily georeference, ...
1
vote
0answers
126 views

32Bit QGIS 2.10 Problems with Gdal WMS TMS google maps layer

I have been using gdal wms layers to view various basemaps in QGIS. I am using this technique rather than open layers because this I can reproject gdal WMS layers as needed. When I get close to the ...
2
votes
1answer
62 views

Change number of cells in raster using GDAL?

i have 5 raster files and i want to merge them to one file with exact 50x50 cells because i have an R script to process it later. I tried the ArcGIS function "mosaic to raster" but it delivers a ...
2
votes
0answers
344 views

Comparing the use of GDAL to build raster pyramids or overviews versus Arcmap 10.1

I have a large set of imagery that I created a mosaic from in arcmap 10.1. I want to use GDAL to create the overviews because arcmap has issues building the pyramids and never seems to complete. I ...
1
vote
2answers
240 views

Extracting all vertices from line shapefile in a simple way? (no python scripting)

I found that shpdump is a very good tool to extract all vertices from an ESRI polyline shapefile, but I need a single triplet (X,Y,Z) file (ie. shpdump provides all coordinates by parts). Is there a ...
2
votes
1answer
217 views

Mosaic Imagery using GDAL Edited

I have a couple of large sets of imagery I need to mosaic using GDAL. The sets are currently in tiff format with 18 tiles totaling about 20gb in size. I have attempted to mosaic using both the ...
3
votes
1answer
312 views

raster tiles from vector data with GDAL - how to avoid aliasing artefacts?

Prequesites I have vector data in a Shapefile N47E007_CONTOURS.shp produced with wget -O N47E007.zip http://e4ftl01.cr.usgs.gov/SRTM/SRTMGL1.003/2000.02.11/N47E007.SRTMGL1.hgt.zip unzip N47E007.zip ...
2
votes
0answers
104 views

Output formats not showing in heatmap plugin (QGis 2.10.1 in ubuntu 14.04)

I just upgraded to QGIS 2.10 (using the debian repos) and the heatmap plugin doesn't show any formats in the dropdown menu. I'm able to open every raster format I have at hand and, of course, the ...
2
votes
1answer
145 views

How to define a projection for a file with geo points?

I have a HDF file that I want to convert to tif or ENVI file. I know that the +Proj of the original file is: "+proj=cea +lon_0=0 +lat_ts=30 +ellps=WGS84 +units=m" so: gdal_translate -of "ENVI" -...
0
votes
1answer
43 views

Dtm format save only heights

I need dtm format that store only heights and scaling and left up position for saving storage I glad to know if any stadard save heights that way
0
votes
1answer
91 views

GDAL -SQL or -where options

ogr2ogr -overwrite -t_srs EPSG:3978 -f "ESRI Shapefile" -sql "SELECT * FROM 'Layer Name' WHERE FIELD_NAME = 'A CLASS'" dst.shp src.shp help me out here - looking to 1) reproject a shape and 2) ...
1
vote
1answer
69 views

Why does Spatialite layer prompt for CRS if SRID was not explicitly defined in Layer Options?

In QGIS, if I Save As... a SpatiaLite file and don't enter an SRID in Layer Options, I will be prompted for a CRS upon opening the file. Yet GDAL documentation says that this option is only for ...
3
votes
2answers
274 views

Finding the pixel coordinates using GDAL/numpy

Is there a way to get the XY cordinates of a pixel based on the array index using numpy/gdal? I have a new set of extents for a raster and, based on these newer extents, want to determine the ...
0
votes
0answers
137 views

find pairs of near features with geopandas, fiona, gdal etc

I have a feature class of points, a few thousands I have now, but I expect to deal with hundreds of thousands of points. I would like find pair of points that is within a specific distance. The ...
2
votes
1answer
128 views

Image footprint in Landsat mosaic using GDAL?

I need to produce a Landsat mosaic over part of Indonesia and an extra shape file containing the dates of each part of the mosaic. The steps I used are: Stack 3 bands to produce RGB image Project ...
0
votes
1answer
32 views

C# ogrinfo with virtual vector source doesn't work GDAL 1.10

I am trying to use GDAL 1.10 with a vector virtual source which is nothing more than the example union of two shapefile layers. from the commandline, ogrinfo works fine on the source. But the C# ...
2
votes
1answer
321 views

How to create GTiff from Tiff using GDAL

I have a tiff image and I would like to georeference it, and so create a Geotiff. To do that, I have to use GDAL. I use in particular, the executable store in the QGIS file (gdalwarp.exe, ...
1
vote
1answer
295 views

Geoserver ecw plugin problem

I am trying to publish some ecw files with geoserver 2.7.2 on ubuntu 15.04 x64. I downloaded gdal plugin, gdal_data and native libraries for image-io 1.1.11 and set environment variables; GDAL_DATA, ...
2
votes
1answer
218 views

gdal2tiles georeferenced image as overlay not working

Ive been looking for a solution to my problem for some time now. I'm using gdal2tiles-leaflet to create tiles for my Leaflet application. The images i use are GeoTiff's with the embedded projection ...
2
votes
1answer
491 views

Converting Lambert Conformal Conic (LCC) to UTM

I have two numpy meshgrids (X, Y both 2D) in Lambert Conformal conic with the following information: I'd like to directly (or indirectly if not possible) project them to UTM. I have both pyproj and ...
1
vote
0answers
35 views

How to Divide geo spatial Raster images into tiles and combine back into a complete image

I am reading a 2048 x 2048 image in the form of (204 x204) tiles and bilinear interpolating them to get (819x819) tiles. Then combining several 819x819 tiles to form a 8190x8190 image. Now as you can ...
1
vote
2answers
103 views

Use gdalwarp to crop area with simple coordinate system

I'm new to GIS and trying to cut image with gdalwarp I need image projected like google map and I need to input coordinates of area of interest in latidute longitute system, like: from 45.3 10.2 to ...
1
vote
1answer
210 views

ogr2ogr problem mapping shapefile to geojson ETRS89 / Poland CS92

I'm trying to convert official boundaries of administrative regions of Poland from official governmental site (codgik.gov.pl). I want to use those data because they are official ones, and set contains ...
0
votes
0answers
725 views

GeoPandas save to shapefile issue: Unable to open EPSG support file

When I try to use to_file() to write a GeoPandas GeoDataFrame to a shapefile I get the following error. df3.to_file('~/Desktop/testout.shp', driver = 'ESRI Shapefile'); ERROR:Fiona:OGR Error 4: ...
1
vote
2answers
142 views

Suitable raster format for elevation data with sub-meter accuracy?

I've got some elevation data represented in ascii gridded xyz-files which I want to convert to a conventional raster format. What's a bit special is that the elevation values have sub-meter accuracy, ...
3
votes
1answer
391 views

Is it possible to read raster files by block with rasterio?

I find rasterio much easier to use that plain ogr/gdal, so I'd prefer using it. However, I need to process large rasters, so reading rasters by block as explained in the Python GDAL course is ...
2
votes
3answers
381 views

-9999 (no data value) becomes 0 when writing array to GDAL memory file

I am trying to reproject and resample (1000m ->> 30m) a raster image (shape 238X406) using gdal.ReprojectImage() function in Python. The input file has -9999 cells as no data value. The result is an ...
-1
votes
2answers
92 views

Perform band calc in numpy array

I'd like to perform a band calc in some raster bands. What I am trying to do is import the raster (ENVI format) like array with gdal, then I do the calc but when I try to save the raster I get an ...
0
votes
0answers
27 views

GDAL data - change top left coordinate

I have GDAL raster dataset. I want to extract raw values from it. All is working correctly, but problem is with top left coorinates. If I extract data from GDALRasterBand, point[0,0] is at topLeft. ...
6
votes
1answer
202 views

Extracting the pixelwise proportion of classes of a finer raster in a coarser one

I have two rasters, one has a 1 km resolution. The other has a spatial resolution of 30 m. This latter pixelwise is a classification. I want to overlay the 1 km raster on the classification raster and ...
3
votes
2answers
210 views

How can I convert points in EPSG:2956 to EPSG:4326 with pyproj?

I have some points in EPSG:4326 -110.28879750 64.46310187 405.49938863 -110.28879359 64.46310086 405.45967110 -110.28880175 64.46310446 405.93823604 -110.28880245 64.46310258 405.48060850 -110....
2
votes
1answer
292 views

Python-gdal create geotiff from array with colormapping

I am creating grayscale, single channel geotiff files from numpy arrays with the below module, but I would like to apply a fairly standard colormap to the data. I've tried three channels of R,B and G ...
2
votes
1answer
108 views

How can I export a shapefile from PostGIS with a GeoJSON query?

I tried to export a shapefile from my PostGIS installation using pgsql2shp as follow pgsql2shp -f tiles.shp -h <host> -u <user> -P <pass> <database> "SELECT id, the_geom FROM &...
0
votes
0answers
58 views

How to transform ENVI Hdr to PNG using GDAL?

I am trying to transform a ENVI Hdr image into a PNG (or other RGB image format) in order to visualize this image. I am selecting the three default bands (55, 41, and 12). I am using GDAL by using ...
0
votes
0answers
61 views

Aggregate/merge data upon insert into PostGIS (via gdal)

I am working on a project where I take a large raster file (1.5GB+), put it through a number of image-processing operations, and finally convert it to vector data. Being that the raster is large, I am ...
2
votes
1answer
118 views

python crashes with returned gdal band object from a function

I wrote these two functions in a module when I want to read the band in NetCDF or HDF files and return the band outside the function. Each HDF and NetCDF has many subdatasets. Each subdataset in a ...
2
votes
0answers
124 views

gdal_merge shifts rasters after processing, layer stack is slightly off

I have experienced that gdal shifts rasters. The outputs of gdal_warp are often skewed, have different pixel sizes and so forth. A known error which has also been discussed here. But now, I had to ...
1
vote
0answers
37 views

Python crashes when I try to import GDAL after re-installing it

A while back I went through the installation process for GDAL by installing the core components and then the Python/GDAL bindings. Everything was fine when I would type in from osgeo import gdal. ...
3
votes
1answer
401 views

How to re-compile ogrinfo against gdal used by QGIS?

I just solved a problem with missing libraries blocking use of ogrinfo and ogr2ogr, only to find that these utilities are garbled in one more way: whenever I run any of these, I get following error (...
0
votes
2answers
45 views

How do I run gdalwarp on a number of csv of coordinates?

I have a CSV of Ground Control Points in NAD83 UTM 12N and I need to convert them to WSG84. Is it possible to do this with gdalwarp or something similar?
1
vote
1answer
95 views

Erase shapfile (polygon) area from GeoTIFF

I've been reading several posts to do so and one option would be to use gdalwarp -cutline, but I would need to keep what it is NOT inside the polygon area. GDAL or any command-line tool does that? ...