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
0answers
75 views

NetCDF to Shapefile with Python and GDAL?

I need to save my outputs, which are netcdf files or panda python tables, to shapefiles. They cannot be rasters because the software I am loading them into accepts shapefiles only. I have googled ...
0
votes
1answer
58 views

How to open SMAP's Surface Moisture data files

I downloaded surface moisture data from Nasa's Soil Moisture Active Passive (SMAP) sensor website (https://nsidc.org/data/smap/smap-data.html). For each processing level (1-4), three types of files ...
0
votes
1answer
28 views

gdaldem hillshade using command to calculate hillshade of whole folder

I need to calculate hillshade of folder using gdaldem hillshade which has 60 files in it. How can i do it using command prompt window.
0
votes
0answers
19 views

USGS Raster Conversion with GDAL Returns Blank .TIF File

I am new to GIS systems. I installed GDAL to convert .ADF files to .TIF files to obtain image heightmaps. It outputs a GeoTIFF .TFW file that I can't open, and a .TIF file that is blank (all white). ...
1
vote
1answer
38 views

raster2pgsql size problem

I am about to import a raster (.tif with 50mb, size 3301x21676, pixel size 10x10) into postgres with raster2pgsql. With small extracts of the image it works but ... raster2pgsql [myfilename].tif -F ...
0
votes
1answer
47 views

Creating point shapefile from GeoTIFF using GDAL/OGR?

I have a 0.05 deg global maps (~8000 files) in GeoTiff format. It's a kind of a mask containig "1" and "0". I need to create a point shapefile for each raster with points located over "1" values of ...
1
vote
1answer
49 views

Dissolve in ogr / python

I have a set of ~100,000 polygons (.shp format). Using python and ogr I would like to dissolve them to one feature class to speed up other downstream operations. The command I am using is: ogr2ogr -f ...
2
votes
1answer
61 views

Create new shapefile based on values ranges in existing shapefile using gdal/ogr/python

I have a shapefile containing many thousands of polygons with values stored in an attribute field and ranging 0-105. I would like to create a much simplified file that contains classes for value ...
1
vote
2answers
78 views

ESRI Shapefile not loading into SQL Server

I have been loading a number of shapefiles into SQL server 2012 using ogr2ogr.exe and the MSSQLSpatial driver. Up until this point, I have been very successful with doing so. However, I have now run ...
2
votes
0answers
13 views

It is possible to use a GDAL in plugin qgis? [duplicate]

I work on qgis plugin, and I say why not not work with GDAL, it'll be too useful. it is possible to integrate the GDAL library to develop QGIS plugin?
1
vote
1answer
50 views

How can a make a netcdf with subgroups smaller?

I'm trying to write a test for code that works with hycom netcdf global ocean circulation data. The source file 2 subdatasets: gdalinfo hycom_glb_911_2015110200_t000_uv3z.nc | grep _NAME ...
0
votes
1answer
31 views

GDAL “dll not found” error

After following step by step the instructions here to install GDAL: http://sandbox.idre.ucla.edu/sandbox/tutorials/installing-gdal-for-windows When I run from osgeo import gdal, I get the following: ...
0
votes
0answers
19 views

OGR_LINK_FIELD to filename in a Geospatial PDF

Im using gdal_translate to create a Geospatial PDF. im trying to add file path to my feature when clicking on them by using: -co OGR_LINK_FIELD=<column whit path> i creating my pdf in the ...
0
votes
0answers
33 views

GDAL error in QGIS application displaying VRT files

I am working with a standalone app build with QGIS 2.8.1 on a CentOs7 machine. I got a new DTED data set which raises errors "somewhere" in the app. The app doesn't read the .dt? files directly, but ...
0
votes
1answer
73 views

Best resampling method for DEM?

I have a task to slice 10m elevation model into the 0.1 x 0.1 degree tiles. I also have to change the projection method used in the TIFF - the WGS84 (Mercator) is a must. I have implemented the ...
1
vote
0answers
25 views

Raster statistics change after re-projection — which one is correct

I'm looking into some global DEMs (like ETOPO1 - http://www.ngdc.noaa.gov/mgg/global/global.html), and noticed that the statistics for the elevation values can change quite a bit after reprojecting ...
0
votes
0answers
17 views

using import ogr in Arcgis python [duplicate]

I need to use import ogr in my python window provided by ArcGIS 10.2, but i can't configure ogr and gdal on my system. Please tell me how i can fix this issue and import ogr in my python window ...
0
votes
0answers
12 views

Scripting procedure to clip vector lat-long file to those that fall within a raster? [duplicate]

I am looking for something similar to this thread, but in a scripting procedure which can be run through a Bash shell. In essence, I have a vector of latitudes and longitudes that span most of ...
5
votes
1answer
92 views

Reprojecting MODIS-Sphere ESRI GRID file to WGS84 using QGIS/GDAL

I have an ESRI GRID file projected in MODIS-Sphere (here's the dataset http://urban.yale.edu/sites/default/files/data/setoguneralphutyra_uf.zip) that I need to convert to a WGS84 GeoTIFF. I tried ...
0
votes
0answers
23 views

Python pyModis installation give error on windows 10 using pip install pyModis [duplicate]

I am using windows 10 operating system on my HP i7 machine with Python 2.7 installed. For windows, it is possible to install pyModis using pip by writing the following syntax in Windows PowerShell: ...
0
votes
0answers
27 views

gdal2tiles_parallel seems to be signiciantly slowing down for super high-res images

I have 3 images that I've been trying to 'tileize' 1 86,000*41,000 2 86,000*21,000 both are simple .pngs I've been using the parallel version of gdal2tiles on a DigitalOcean box with 20cores and ...
0
votes
1answer
51 views

Improving performance of Clip using OGR?

I have a set of points which I am trying to clip to a set of polygons. The files have many features (150,000 polys, 8,000 points). I'd like to automate this process carry it out once a day, hence I am ...
0
votes
0answers
38 views

Using gdal in ironpython

I would like to replicate the following two statements GDAL python has: # project the decimal degrees to Mercator: gdalwarp -s_srs EPSG:4326 -t_srs EPSG:32631 -r cubic original.asc ...
1
vote
0answers
34 views

Extract raster values to multipoint shapefile and write to new field using gdal/python [duplicate]

I have a raster with one band and I would like to query a set of overlying points (.shp) to determine how many points as well as which points points fall under a range of certain pixel values ...
0
votes
0answers
18 views

Why does gdal2tiles drop lowest values of source image?

I am trying to use gdal2tiles to batch-create a bunch of KML superoverlays from a set big geotiff files; problem is that the gdal2tiles output seems to be stripped of its lowest data. Here is an ...
1
vote
0answers
26 views

Theory notions on exporting bufr data in geotiff format

I'm trying to export some bufr data from polar satellites/sensors like DMSP, ASCAT, ATMS etc in geotiff image. When my program rasterizes satellite data (brigthness temperature for example) by ...
1
vote
1answer
25 views

How to make transparent the Index 0 of the color table of a Tif/GTiff/GeoTiff using gdal?

With a Tif file, I'm trying to generate tiles using gdal tools. gdalinfo gdalinfo I_III_IVA_MRA_WGS.tif | more Output Driver: GTiff/GeoTIFF Files: I_III_IVA_MRA_WGS.tif Size is 17119, 23309 ...
1
vote
3answers
176 views

Open LiDAR DEM of the UK with *.asc extension

I've downloaded a tile of a freely accessible LiDAR DEM of the UK from http://environment.data.gov.uk/ds/survey#/download. But I don't know how to open it. The files inside have the *.asc extension. ...
0
votes
0answers
34 views

How to create a geospatial PDF using GDAL (in python)? [duplicate]

Can anyone suggest a way to create a geospatial PDF instead of geoTiff using GDAL in python? Aside, I also tried creating first a geoTiff and then used the 'gdal_transform' command to try and create ...
1
vote
1answer
101 views

Converting large jp2 files using gdal

I have a bunch of jp2 files with size going from 40MB to 1GB and I need to convert them all to a .tif format. When using gdal_translate for the small files (having the jasper driver to support jp2 ...
0
votes
0answers
10 views

How to change the pixel aspect ratio of a raster? [duplicate]

I hva a raster where the pixel size is (0.000000610870000,-0.000000273755000), but I'd like them square (e.g. (0.000000234155000,-0.000000234155000)). Is there a way to accomplish this using GDAL? ...
2
votes
2answers
80 views

Getting origin as long, lat in GDAL without knowing EPSG

My goal is to get the long, lat coordinates for the origin of a particular geotif. gdalinfo outputs: Upper Left (-2764486.928, 3232110.510) (130d 7'22.42"W, 49d56'48.81"N) But I don't understand the ...
2
votes
1answer
190 views

Calculate area of polygons using OGR in python script

I have a shapefile of many thousands of polygons. I'm trying to append an area field to the attribute table and calculate the area of each polygon (in sq. meters) driver = ogr.GetDriverByName("ESRI ...
6
votes
0answers
58 views

QGIS: large .vrt indexing

I have a large number of small tiles merged into a single VRT. The rasters are only viewed at large scale, so only a few of them are visible at the same time. The map refresh in qgis is slowing down ...
3
votes
2answers
55 views

How to open and save QGIS datasource with/to lower case fields?

I want to be able to change all field names to lower case from a data source when loading them in QGIS (SHP in particular). Is it possible? Via the GUI, via Python. Is it possible to convert field ...
1
vote
0answers
26 views

WMS layer throws “IReadBlock failed at X offset 0, Y offset 1” error. Possible XML issues

I am working with a few WMS layers to try and see which one would work best with an application I'm developing out of QGIS 2.8. I need the CRS to be in EPSG 4326 but, as many of you know, most webmaps ...
1
vote
0answers
20 views

Should I set the projection if the data has already been projected?

I need to plot some data which was produced on a native grid with a lambert conformal conic projection. I need to create a geotiff (using python and gdal) and was wondering how to go about and do ...
2
votes
2answers
79 views

gdal/python: extracting projection info from hdf file

I'm trying to do something fairly simple. Take two subdatasets from an hdf file, read them as an array, do some math and write the new array back to a geotiff, using the same projection as the input ...
2
votes
1answer
48 views

Reading GRIB band metadata with GDAL

I'm trying to use GDAL in python to get only the temperature band from this GRIB file that has 2 bands. I can't just do temp_band = ds.GetRasterBand(1) because sometimes, but not always, the grib ...
1
vote
0answers
76 views

how to set the geotransform using GDAL in python

I am trying to create an image from a numpy array and georefenrece it so that I can import it into qgis. I have the latitude and longitude in degrees for the top left corner of the array as well as ...
1
vote
2answers
119 views

ogr2ogr without spatialite support - How to fix it?

I am trying to use ogr2ogr to convert OSM data from pbf to spatialite. I used to do so with the following command: ogr2ogr -f "SQLite" -dsco SPATIALITE=YES png-full-last.db ...
0
votes
0answers
34 views

How to create raster with origin in degrees and dx dy in metres

I need to create a raster from array. In gdal if i understand correctly i need to specify the top left lat lon for the SetGeoTransform((originX, pixelWidth, 0, originY, 0, pixelHeight)) however the ...
1
vote
0answers
23 views

Projected Coordinate System issue using GDAL and qgis when defining x and y

I have a numpy array which I would like to save as GeoTiff using GDAL in python. The array was created using the following PCS in mind: +proj=lcc +lat_1=48 +lat_2=48 +lat_0=48 +lon_0=0 +x_0=0 ...
0
votes
0answers
37 views

How to transform any GTiff to Lambert Azimuthal Equal Area with GDAL?

I am trying to transform a GTiff EPSG:4326 +proj=longlat +datum=WGS84 +no_defs and a GTiff EPSG:31982 +proj=utm +zone=22 +south +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs both to ...
1
vote
1answer
59 views

No GDAL/OGR package for Ubuntu 15.04?

Using Ubuntu 15.04, I'm trying to install GDAL/OGR. I did: [] sudo add-apt-repository ppa:ubuntugis/ubuntugis-unstable [] sudo apt-get update However, I found this error message W: Failed to ...
0
votes
1answer
178 views

gdal_merge.py can't find libgdal.so.20

I am having some trouble with GDAL (v2.1-dev). I build it from source today using Poppler and OpenJPEG. All the command line utilities work fine, but none of the Python plugins work at all, ...
0
votes
1answer
15 views

Rectifying a Tiff Given One Corner Coordinate

I'm trying to rectify an image that I don't know much about other than the coordinate of the southwest corner, which I was able to pull using ExifTool. Is there any programming solution to ...
0
votes
1answer
65 views

Calculating the difference between two grids in percent with QGIS

I'm using the QGIS GDAL Raster Calculator to get the difference in percent between two grids. Grid A contains an integer, and grid B contains an integer which may be the same as or larger than the ...
1
vote
1answer
55 views

GDAL install on Ubuntu 12.04 works for python 2.x but not python3

I need to install gdal 1.10 or greater to be compatible with python3 on ubuntu 12.04. The machine also has python2 installed. I followed instructions here Error While installing GDAL in Ubuntu ...
1
vote
0answers
92 views

How to handle GDAL dataset error when accessing Arcgis Online world imagery data?

I would like to use the world imagery data set for loading base maps to my QGIS projects. Up until now, I have had no problem by typing the following command in the Python console: ...