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
103 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 ...
1
vote
1answer
107 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
29 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
244 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
250 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
175 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
406 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
31 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
95 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
178 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
578 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
126 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
338 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
305 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
76 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
25 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
188 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
168 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 ...
2
votes
1answer
216 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
90 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
52 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
58 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
102 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
110 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
34 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
294 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
81 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? ...
0
votes
1answer
183 views

Problem doing a coordinate transform with GDAL python

I'm trying to do a concordant transform from the format of the shapefile (3309) to google maps(4326). This code works on one Centos 6.6 box with the SCL python 2.7 and a custom gdal 2.0 build. ...
-1
votes
1answer
40 views

ASC file to GeoTIFF

I would like to convert using Arcpy an ASC file to GeoTiff but I can't find nothing about it. I tried : arcpy.ProjectRaster_management() changing .tiff to .geotiff but it didn't work.
1
vote
0answers
76 views

Serving GeoTIFFs of different resolutions as single layer

I'm trying to serve a 30 arc second DEM with global coverage combined with high resolution 1/3 arc second DEM for a very small area. I believed at first I could combine them into a single file, ...
0
votes
0answers
248 views

WCS elevation with Cesium + Geoserver

I'm trying to develop an offline demonstrator showing hi-resolution bathymetry data using Cesium 1.11 and Geoserver. I've merged a number of 1/3 arc second areas around Puerto-Rico from NOAA ...
2
votes
2answers
68 views

Fusion table connexion issue

I'm trying to import some data from a SHP to GFT (Google Fusion Table). I get a 401 error. I'm using an Ubuntu 14.04 LTS Linux 64 bits with GDAL 1.11.0, released 2014/04/16. GFT is available (tested ...
1
vote
3answers
121 views

OSGeo4W GDAL Polygonize Batch Convert?

I'm trying to use GDAL polygonize to run a batch convert but I cannot for the life of my figure out a way to get it to read a text file as input that contains all of the commands. Do any of you know ...
0
votes
2answers
48 views

I cannot find osr after installing OSGeo4W package

I wanted to be able to import OSR from GDAL on my Python 2.7 (Anaconda) Windows 7 64-bit. I installed OSGeo4W after reading this post and eventually I'd like to import osr from osgeo import osr ...
1
vote
1answer
362 views

How do I convert WGS84 to Lat, Long to Google Maps in R?

I have this coordinate which I believe it's in WGS84 format x=1005740.867110, y=243957.356623 I want to plot that in Google Maps which I need to convert that to Lon, Lat. I've found this script ...
3
votes
2answers
313 views

How to have GDAL print layers of GeoPDF AND say which are raster vs vector

My Objective: I would like to use GDAL to convert a GeoPDF. I want the vector layers as shp files and the raster layers as tif files. I want to do this in a programmatic way. Edit: In reality, I ...
0
votes
0answers
19 views

Using GDAL Convert .SHP to GEO TIFF in Command Line Python [duplicate]

I have dozen Shape files (.SHP) with Georeferenced. I need to convert to GeoTiff or other format that GDAL primarly gdal2tiles.py to convert to MBTiles format. All I want is convert the .shp files ...
0
votes
0answers
29 views

Installing GDAL Python bindings on Windows [duplicate]

This has been addressed before (How to install GDAL with Python on windows?), but I wanted to bring the question up again, to see if anybody has come up with a better method in the last 5 years, and ...
0
votes
1answer
187 views

How to get elevation from dted only USGEO format using gdal apps

I need to get hieght by lat long position using gdal utilities
0
votes
0answers
74 views

What are the units in GMTED2010 elevation data?

I have downloaded the GMTED2010 elevation data from http://topotools.cr.usgs.gov/gmted_viewer/viewer.htm I am merging all the files/tiles (to create a global dataset) with: gdal_merge.py There is ...
1
vote
0answers
48 views

Non-Square Pixel on georeferencing of images using GDAL

I have a image in .jpg format. I have some GCP and also know the pixels values of these GCp in that image. Now i trying to convert this image into geotiff. for that I have coded the following lines. ...
2
votes
0answers
71 views

Multiple layer intersection

I am trying to find a way with qgis/python to get the intersection area of more than 2 layers. All the tools in Qgis take just the 2 input layers. ='( For example in the image below, I have 4 ...
3
votes
0answers
130 views

How do I specify an outsize for gdal translations to an ENVI file?

I'm trying to create an elevation map using the output from gdal_translate. When I run gdal I can convert to both PNG and ENVI and combine the results to make something like the following image. ...
2
votes
2answers
671 views

Install GDAL 2.0 on a MacOSX

Does anybody know what's the easiest way to install GDAL 2.0 on a MacOSX? The usual suspects kingchaos binary distribution and homebrew seem to be stuck on GDAL 1.11 I need support for TopoJSON ...
2
votes
1answer
251 views

Reprojecting MODIS Sea Ice Extent and IST data using GDAL

This is an additional exploration of the problem described at Reprojecting MODIS Ice Sea Temperature products with gdalwarp which was never answered in a satisfying way (I apologise for the length of ...
4
votes
1answer
71 views

Error trying to import data to CartoDB using ogr2ogr

I have a local install of CartoDB in my Ubuntu 14.04. I have created a CartoDB account called "development" and been playing with it for a while. Now I'd like to import some data from a shapefile ...
1
vote
1answer
134 views

Using GDAL library in QGIS plugin?

i've been working on my project to create QGIS plugin for raster analysis purpose. It will be easier if we use gdal library. Are we able to use gdal library to develop QGIS plugin? If so, how do ...
-1
votes
1answer
98 views

How to cut image into Map Sheets with GDAL?

I want to divide a whole raster map image into map sheets. The gdal tiling or clipping functionality doesn't seem to do it.
1
vote
1answer
316 views

How to Install GDAL with ECW Support on Windows

What do I download from http://www.gisinternals.com/stable.php for Win7-32 to cut 1:1M map sheets out of Australia.ecw file in GDA94/Geographic and save them as geotiff in WGS84/Geographic? I would ...