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

Upper limit to number of GeoTiff tiles?

I'm using the GDAL SWIG python bindings to create a new GeoTiff raster for use in a blockwise exponential-decay function. kernel_size = 177764054 driver = gdal.GetDriverByName('GTiff') kernel_dataset ...
0
votes
2answers
32 views

Skip a specified input layer in ogr2ogr

I am using ogr2ogr to load data from multiple GML files into PostGIS. As shown in the ogrinfo output below, each GML file has multiple layers. Had to open data source read-only. INFO: Open of ...
3
votes
2answers
70 views

Can GDAL extract a raster window from disk without loading the whole file to memory?

I have a 500MB GeoTIFF with 3 bands on disk. I can easily read a very small, 25x25-pixel subwindow using GDAL like so: gdal_translate -srcwin 600 550 25 25 hugeimage.tif smallwindow.tif Does this ...
1
vote
2answers
122 views

Determining Min and Max elevation in Python using GDAL? [closed]

I have a DEM file and shape file of specific area ,and i want to Create a module in Python that writes name of each area which is shows in the picture, with the greatest differences in altitude on ...
3
votes
1answer
76 views

Optimizing Python GDAL ReadAsArray

I am using the GDAL ReadAsArray method to work with raster data using numpy (specifically reclassification). As my rasters are large, I process the arrays in blocks, iterating though each block and ...
1
vote
1answer
100 views

Python Gdal ERROR 6: GEOS support not enabled

I installed GDAL from here (https://pypi.python.org/pypi/GDAL/1.9.1) Before attempting the python GDAL I installed the prerequisite gdal libraries using Kyngchaos's description some time ago ...
1
vote
2answers
46 views

gdal_translate crashes on huge jp2

I have 44GB jp2 file, I want to translate it to tif file, but once it reaches 10% it crashes. No messages, just stopped. This is my command: gdal_translate.exe -of GTiff --config GDAL_CACHMAX 8000 ...
2
votes
1answer
26 views

Libraries for creating surfaces output in geotiff

I am looking for some libraries that can assist me in creating geotiffs from 3d polygons that becomes surfaces. Given a set of polygons where each node is [x,y,z] coordinates, i would like to create ...
6
votes
1answer
137 views

Using a geo-referenced image as an overlay to display markers on floorplan

I've created a tool based on the Google overlay tiler. This allows our users to upload a floor plan, adjust the position etc. using Google maps. Nothing fancy. We have a server-side job which creates ...
1
vote
0answers
43 views

gdalwarp UTM Zone 1 to WGS84 - fails to transform

I have a global dataset consisting of tiles in UTM projection that I want to mosaic and therefore need to reproject to a common projection. Can I solve my reprojection error with options to gdalwarp ...
1
vote
1answer
63 views

GDAL Python OGRFeature .GetGeometryRef vs .Clone().GetGeometryRef

Just recently I was trying to get the geometry of a polygon stored in an ESRI Shapefile. In the script I was trying to read the polygon as follows: import ogr path = "shape.shp" ds = ogr.Open(path, ...
1
vote
0answers
35 views

Handling warning about DX and DY in gdal_translate from .Gtiff to .ASC?

When I try to convert Gtiff to .asc file, QGIS shows the warning: Producing Golden Surfer style file with DX and DY instead of CELL SIZE since the input pixels are non-square. Use the FORCE CELL SIZE ...
2
votes
0answers
80 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
63 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
30 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
50 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
50 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
83 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
54 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
37 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
78 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
101 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
35 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
26 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
183 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
103 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
84 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
204 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
56 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
27 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
86 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 ...