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

1
vote
0answers
11 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
18 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, ...
0
votes
0answers
19 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
40 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
20 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
94 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
39 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 ...
2
votes
1answer
42 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 ...
1
vote
0answers
22 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
66 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
19 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
38 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) ...
0
votes
1answer
18 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 ...
2
votes
2answers
44 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
15 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
21 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
10 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# ...
1
vote
1answer
65 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, ...
0
votes
1answer
54 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, ...
1
vote
1answer
38 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
67 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
14 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
25 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
36 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
34 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
86 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, ...
2
votes
1answer
143 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 ...
1
vote
1answer
64 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 ...
0
votes
2answers
29 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
18 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. ...
0
votes
0answers
43 views

gdal2tiles: how to make filenames of tiles to represent xyz coordinates?

I'm converting ECW file to tiles using gdal2tiles. Right now the tiles and folders generated by gdal2tiles are named as numbers starting from 0 successively. The folder structure looks like: How to ...
6
votes
1answer
138 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 ...
2
votes
2answers
49 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 ...
1
vote
0answers
25 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
39 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
29 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
27 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
50 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 ...
1
vote
0answers
37 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 ...
0
votes
0answers
24 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
48 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
33 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
56 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
64 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
24 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
32 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
39 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
12 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 ...
0
votes
1answer
30 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
26 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 ...