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

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
73 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
162 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
53 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 ...
0
votes
0answers
35 views

Overwrite the mosaic raster generated using gdal_merge.py?

Is it possible to overwrite the raster generated using gdal_merge.py? When I ran the command: % gdal_merge.py -init 255 -o out.tif in1.tif in2.tif it works fine, but if I ran the same command and ...
3
votes
1answer
63 views

Importing ESRI shp / gdb files with M value into Postgis

Has anyone done this successfully? I have been using DB Manager in QGIS 2.12 and my MutilinestringM values are converted to Multilinestring in PostGIS. When I use the PostGIS shapefile importer, this ...
0
votes
0answers
25 views

gdal.SaveArray returns Received a NULL pointer

I am saving my 2D array into a raster data using GDAL but could not get it to work properly. I am using the gdal.SaveArray option and it returns ValueError: Received a NULL pointer. I think I used it ...
3
votes
1answer
66 views

Build a table with columns of raster values from multiple raster datasets, using Python, GDAL, or PyQGIS?

I have a series of raster datasets, clipped to the same extent, with matching cell resolutions, each of which represent land surface temperatures (LST) captured by the MODIS sensor. Each raster ...
2
votes
1answer
77 views

aggregating using 'gdalwarp -average' resulting in incorrect values

I'm trying to aggregate data up from 0.05 to a coarser 0.25 resolution and am having trouble getting gdalwarp to work correctly. From comparing the original input data with the aggregated data and ...
3
votes
0answers
29 views

gdal_translate multiple GDAL_PDF_LAYERS_OFF entries gives error

The documentation says that GDAL_PDF_LAYERS_OFF takes a list of comma separated layers. Examples look just like what I am trying to do. Executing the following gives an error FAILURE: Too many ...
3
votes
1answer
59 views

ogr2ogr batch load shapefiles into geodatabase

I'm trying to load a directory of shapefiles into a geodatabase with ogr2ogr. It works fine if I specify a single shapefile, but can't open the datasources if I specify a wildcard. Is there a wildcard ...
2
votes
1answer
34 views

Using gdalwarp to remove one row of pixels or add one row of pixels

I am trying to use GDAL_calc but my two rasters are off by one pixel in the y dimension. One is 256005 the other is 256004 so GDAL_calc fails with a multiplication. Is it possible in gdalwarp to ...
1
vote
1answer
52 views

gdal_merge not working in GDAL 201 windows install

I realise this looks like a duplicate of http://stackoverflow.com/questions/31174320/gdal-merge-py-not-working-after-gdal-and-it-python-bindings-have-been-installed, but I think the problem is a bit ...
2
votes
2answers
65 views

MODIS MOD11A1 Reproject Change Data Value

I'm import and reprojected MOD11A1 with GDAL 1.10....my shell script is easy.. for file in *.hdf; do echo "map: $file" gdal_translate -of Gtiff -ot Int16 ...
1
vote
0answers
28 views

Flipped (Y-axis) Raster from point cloud

I have been working on a project where I need to create rasters from irregularly spaced point clouds. I am resampling my point clouds to a regular 0.25 cm grid using a nearest neighbor method ...
1
vote
0answers
42 views

How to correctly rasterize a polygon which crosses 180E longitude?

I want to rasterize a polygon wich crosses 180E lon with gdal_rasterize. The polygon is in the WGS84 projection. I used the following command: gdal_rasterize -burn 1 -of GTiff -ot Byte -ts 1500 ...
2
votes
0answers
27 views

failed gdal script from QGIS Georeferencer using Thin Plate Spline (TPS)

I'm trying to georeference a photo using QGIS Georeferencer plugin (this is an oblique photograph. I know I should be using a proper georectification routine, but my current work now is just to simply ...
2
votes
0answers
128 views

gdal/python: tif not recognised a supported file format

I'm trying to open a tif using gdal, read it as a array, do some math, and write it back as a tif. However, when I try to gdal.Open(tif) in python, I get the error: RuntimeError: ...
1
vote
1answer
42 views

how to read a specific overview layer from GeoTIFF file using GDAL / Python

I know how to create the overviews to a specific geoTIFF file using GDAL ( i have tested it in QGIS), but i cannot find any good information on how to read these overviews individually in python. ...
-1
votes
1answer
35 views

Convert all geo TIFF in a directory in ASCII [duplicate]

I would like to convert all geoTIFF in a directory in ASCII with the use of a command like gdal_translate -of GTiff... in OSGeo4W Shell. Is this possible?
1
vote
0answers
25 views

GDAL with PostgreSql LIST_ALL_TABLES fails (no it doesn't) - SOLVED [closed]

I have compiled GDAL 2.0.0 and am accessing a PG database. I note that the online documentation has two options available with the word "TABLES" in it. I think the logic in the code isn't handling the ...
0
votes
0answers
38 views

Mapserver mapfile reading problem at maximum zoom level

I have created a mapfile for cartosat Raster data(B/W) which is working in browser fine but after zooming to maximum it is throwing an error as follow drawGDAL():unable to acess file ...
1
vote
1answer
75 views

GRIB to NetCDF with GDAL, specifying dimensions

I have a GRIB dataset representing temperature over North America, in Lambert Conformal Conic (with specified central meridian, latitude of origin, etc.). It has what should be understood as a ...
2
votes
1answer
45 views

R gdalUtils::ogrinfo error or very slow

I am pulling WFS data into R and using gdalUtils::ogrinfo to look at the available layer names on a WFS server. Using the same url in QGIS, this is very fast and gives complete information. With ...
2
votes
1answer
37 views

org2ogr PostgreSQL to GPX

I already have found how to import GPX files into PostgreSQL using gdal/org2org. Now I am wondering of dumping into GPX files the results of my modification on dataset stored into PostgreSQL. But I ...
4
votes
0answers
87 views

QGIS/OGR shp to tab export - crashing Mapinfo due to number format issues

I created a vector grid in QGIS and it exports as SHP -no issues. When I save as tab and try to open the file in MapInfo it crashes. According to ...
0
votes
0answers
29 views

Image pyramid with transparent pixels

I have a GeoTiff image which I wan't to use to build tile pyramid and serve it via geoserver. In most cases it works well but fails to publish a layer if tiled images have NO_DATA (transparent) ...
1
vote
0answers
26 views

How delete feature from layer using Gdal binding java

I am trying delete one feature from layer using GDAL in java https://trac.osgeo.org/gdal/wiki/GdalOgrInJava But nothing happens. ogr.RegisterAll(); Driver driver = ogr.GetDriverByName("ESRI ...
0
votes
1answer
42 views

GDAL Python: how to extract a cutline from VRT file in world coordinates?

I'm trying to extract the world coordinates of a cutline (clipping polygon) in a GDAL VRT file. Unfortunately, the cutline is stored in image coordinates (pixels), and i could not find the right ...
1
vote
1answer
17 views

How can I create an aggregate GDAL dataset with a specific ordering for overlapping regions?

I have a number of raster datasets (with common projections, pixel sizes, and pixel alignments) that I want to mosaic into a single dataset. The source datasets overlap in a number of areas; in those ...
2
votes
1answer
38 views

Encode an azimuthal equidistant projection in NetCDF conversion

I'm trying to do Azimuthal Equidistant in CF-compliant NetCDF, but I didn't get the correct projection. If you compare the first two information (gdalinfo) you see there is no AUTHORITY and also ...
1
vote
0answers
20 views

gdal_retile — force full tiles?

I have a set of large geotifs (DEM + orthoimagery) that I need to generate image pyramids from. So far gdal_retile from osgeo4w is working very well for my use case with one exception: I need all my ...
1
vote
1answer
92 views

Can't write decimal data into shape (shp) attribute table with Python

I am trying to create a shapefile with Python GDAL/OGR library. I add to it a real number field, but when I write into it doesn't appear as it should but rounded. I will explain better with the code: ...
1
vote
0answers
51 views

undefined reference to `TIFFDefaultStripSize@LIBTIFF_4.0'

I am trying for the first time to use gdal in Ubuntu 15.10 with gcc. I call gdal.h in my program. At first I got the error "undefined reference to `GDALOpen'", so I check online and I found that I ...
0
votes
1answer
60 views

Analyze height data using a digital elevation map

i ran the following code on python 2.7: import sys, struct from osgeo import gdal from osgeo import gdalconst minLat = -48 maxLat = -33 minLong = 165 maxLong = 179 dataset = ...
1
vote
1answer
33 views

Keep “no data” tiles as “no data” w/ proximity (raster distance) operation

Working in QGIS 2.8.1 Wien. I have a raster file with three classes: "No data" for tiles outside of the country of interest's bounds "0" for tiles that correspond to roads "1" for all other tiles. ...
1
vote
1answer
130 views

Convert to lat-lon coordinates GDAL

I have a netCDF file and would like to convert it to LAT-LON coordinates. 1) I do not know the projection in the file but from the website of the data provider, the data is on a polar-stereographic ...
0
votes
1answer
90 views

Project NetCDF file using GDAL or QGis

Can I project a NetCDF file with the help of GDAL or QGis? The file info is: gdalinfo c:\tomcat-6.0.36\webapps\ncWMS\layers\RAIN_MM_1hour_ac.nc Driver: netCDF/Network Common Data Format Files: ...
4
votes
3answers
33 views

gdalwarp to geographic coordinate system

I'm getting the folllowing Error when trying to re-project a raster from NAD_1983_Albers to GCS_North_American_1983. I used the proj4 defenition from going to Properties in QGIS of the layer whose ...
1
vote
1answer
51 views

GDAL: gdalwarp wrong transformation to “South_Pole_Lambert_Azimuthal_Equal_Area” EPSG:102020

I am trying to reproject a raster that is defined in EPSG:4326 to EPSG:102020. Command used gdalwarp -ot Float32 -s_srs EPSG:4326 -t_srs EPSG:102020 -r near -of GTiff -co COMPRESS=DEFLATE -co ...
3
votes
1answer
81 views

How to make a network /tower coverage map in in mapbox

I am new to mapping and currently I am trying to make a map showing coverage of a network provider here in my country. The map I am trying to achieve should look similiar to that at ...
2
votes
0answers
59 views

Why azimuth result is wrong for this pyqgis code?

In this code I want to obtain azimuth from a single point to each pixel center, but when run it, in raster created azimuth result for each pixel is wrong, what is correct code? from osgeo import ...
2
votes
2answers
95 views

GDAL: *.nc rasters, opened as “netCDF/Network Common Data Format” under Linux and as “ HDF5/Hierarchical Data Format Release 5” under Windows

I have a set of rasters that are in the "*.nc" format and I noticed that gdalinfo (and subsequently other GIS packages based on GDAL like QGIS) is opening them with two different drivers. Linux: ...
1
vote
1answer
75 views

How to extract a 5 bands geotiff file into RGB file?

I have a five bands geotiff file for Red, Green, Blue, Red-edge and Infra-red, and want to use gdal_translate to extract bands Red, Green and Blue gdal_translate -b 1 -b 2 -b 3 -mask "none" ...
1
vote
1answer
111 views

ogr2ogr Shapefiles to Postgres SRS problem

I am trying to import to my Postgresql (9.1+129ubuntu1) + PostGIS (1.5.3-2) database a Shapefile with polygons with the ogr2ogr tool from Gdal-Bin package (1.7.3-6ubuntu3). I know the origin SRS of ...
1
vote
1answer
67 views

Why is GDAL/QGIS creating Geospatial PDF's with erroneous spatial data?

I have a GeoTiff that I wish to convert into a geospatial PDF for use in Avenza. I attempted to use gdal_translate, and also used the QGIS's Raster->Conversion->Translate which apparently just builds ...
3
votes
1answer
76 views

Tearing of polygons using ggmap and readOGR

I'm trying to incorporate a shape file of US states with ggmap and keep getting "tearing" of my polygons. I tried changing group=id to group=group which helped, but didn't seem to resolve the ...
0
votes
0answers
43 views

Copy specfic layer to new datasource using OGR Swig Java bindings

Is it possible to efficiently copy a specific layer from one datasource to a different datasource using GDAL/OGR Java bindings? I know it is possible to create a new layer in the target data source ...