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
1answer
19 views

gdal_polygonize returns with a strange error

i'm trying to execute gdal_polygoinize.py command. I have a main python method that takes in input a raster file and a output shapefile, and then calls a method named polygonize.py that builds and try ...
3
votes
1answer
741 views

Extract data from netcdf using a shapefile polygon in python

I need to subset data from a NetCDF using a specific shapefile. The data are sea surface temperature and ocean color at 1/4 degree resolution. I have 4 polygons describing the US. Northeast ...
1
vote
1answer
27 views

Unable to add MrSID extension to GeoServer

I'm trying to serve some MrSID files with GeoServer following the geoserver.org guide, however, I can't get the MrSID Raster Data Source to show up when attempting to add a new GeoServer Store. This ...
2
votes
1answer
114 views

Implement r.cross / ArcGIS Combine with GDAL

I am using GDAL to process a lot of rasters. I then need to combine all the arrays to create a single output raster that represents all the unique combinations of values. I.E the same thing as ...
2
votes
1answer
21 views

Convert MODIS HDF file in sinusoidal projection into GeoTiff with

I'm looking for a way to convert MODIS LST data from HDF to Tiff format, and reproject from sinusoidal to WGS84 - all using python GDAL. An example file is here (I only care about band 0) : ...
2
votes
1answer
212 views

gdal_calc works but I get a Python error at the end of each process

I have a problem while using gdal_calc by a .bat file. Here is the code : call "C:\OSGeo4w64\bin\o4w_env.bat" gdal_calc ...
0
votes
0answers
9 views

Land-Water mask for MODIS 0.05deg grid

I need to create a worldwide land-water mask for a series of MODIS LST/VI images (0.05 deg). I know that each MODIS hdf file contains a land-water mask. However, this mask changes between different ...
2
votes
1answer
32 views

Changing the origin of a raster?

I have a raster file that is projected not detailed enough. Every pixel in the raster should be replaced, just a bit. With ArcGIS I calculated the offset per pixel. I couldn't find a way to do this ...
1
vote
0answers
43 views

How to Convert from GMT Line files to Shapfile for Qgis?

I have several line text-files of (used in GMT) saved in the format (this is a small sample): -111.29 42.28 -111.29 42.31 -111.29 42.33 > -111.37 41.99 -111.38 42.02 I am trying to convert this ...
2
votes
1answer
204 views

NASA SRTM: precipitous contour lines in TileMill

I was trying to make a nice looking print map of Bali (Indonesia) in TileMill. After creating a terrain layer (using this walkthrough) I noticed that the island's coastline was precipitous. Is ...
0
votes
0answers
8 views

NED (NAVD88) to WGS84 heights

I need to convert the 1/3 arc second DEM files from the NED into WGS84 ellipsoid heights. The documentation says that the NED is referenced to NAVD88. I've tried using gdalwarp as follows: gdalwarp ...
1
vote
1answer
181 views

Accessing bands in .pix (PCIDSK) file in python/arcpy

I am trying to perform zonal statistics on a single band from a multi-band raster in .pix format. I have things working great with a tiff file. I am able to access the bands using this code: desc = ...
1
vote
2answers
501 views

Import error: no module named _gdal_array

I am using python 2.7 and have gdal installed. I want to read an entire .jpg image at once and so I used statement import gdal from gdalconst import * driver = gdal.GetDriverByName('JPEG') ...
1
vote
2answers
427 views

gdalwarp changing value in nodata pixels

I'm using gdalwarp to convert MODIS NDVI mosaics from Sinusoidal to epsg:4326. To save space I'm using VRT files until I reach the final product. Here is my workflow. Mosaic NDVI Modis tiles using ...
0
votes
1answer
47 views

How to read MODIS hdf file into R [on hold]

I am having some problems with the processing MODIS HDF in R. I want to import MOD13A3 NDVI hdf file. Purpose is to read image cells value as a matrix. I have install MRT in my PC, I had planned to ...
1
vote
1answer
1k views

How do I add elevation data to a GeoJSON line from a GPX file?

I understand that elevation can be added as an optional third value in Coordinates (lat, long, elev) in a geojson file, but I can't find a way to actually do that from a gpx file. I've used ogr2ogr to ...
1
vote
1answer
1k views

Importing Census TIGER Shape Files to SQL Server 2012?

I've been doing my research on how to import shape files into Microsoft SQL Server 2012. and am still unable to figure out how to do it after hours of playing around with it. I'll document what I've ...
0
votes
1answer
526 views

Can't get gdal_merge Python script to run

I downloaded and installed GDAL. I am using Python 27. I had to use the following commands at the IDLE Python GUI: from osgeo import gdal,ogr sys.path.append('C:\Program Files\GDAL') import ...
1
vote
2answers
1k views

GDAL/proj4 equirectangular projection coordinate system units assumed by OGRSpatialReference and GDALDataset::GetGeoTransform()?

How do I configure my OGRSpatialReference object to specify an equirectangular projection system http://en.wikipedia.org/wiki/Equirectangular_projection where the affine geo-transform matrix (provided ...
-1
votes
1answer
29 views

Gdalbuildvrt dosen't output the vrt file

I am trying to use gdalbuildvrt to create a time series of landsat images. I used the following code but it is not outputting the vrt file. In the code inlist is the txt file with path to input ...
3
votes
1answer
58 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 ...
0
votes
0answers
9 views

Resampling using gdalwarp gives me no data

I have an image with the following information. Driver: GTiff/GeoTIFF Files:test_DEM.tiff Size is 91, 73 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", ...
3
votes
1answer
245 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 ...
0
votes
0answers
20 views

Gdal ogr simplify/intersects

Can anyone give an example: How to use the simplification? For example // read the coordinates of the polygon char * wgsRef = NULL; OGRPolygon *poPol = (OGRPolygon ...
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 ...
4
votes
1answer
62 views

Georeferencing Himawari-8 in GDAL (or other)

Does anyone know the appropriate project information to georeference (and hence reproject) Himawari-8 AHI data in the IR? The images are 5500x5500 pixels. Currently I try these: ulx=-5570248.832537 ...
4
votes
1answer
202 views

Tabulate Area using GDAL and Numpy

How could I perform a Tabulate Area using GDAL and Numpy? I am trying to summarize how many pixels of each class in a Landcover raster are within each polygon of a Census Blocks shapefile. I've ...
1
vote
1answer
67 views

Installing GDAL2.1

Here will be my sequence of command lines while trying to install gdal2.1 in a UBUNTU virtual machine. My virtual machine is a UBUNTU 16.04 LTS(64bit) I would need gdal2.1 and especially the Python ...
0
votes
1answer
372 views

Combine overlapping rasters in gdal_calc.py, taking the maximum value from each pixel?

I want to combine 3 overlapping rasters to create 1 raster. Each have the same dimensions, resolution and CRS. I want to populate the new raster with the maximum value from each of the layers. ...
1
vote
0answers
31 views

Gdalwarp UTM reprojection causes image distortions

I'm trying to re-project Sentinel-2 images to a common UTM zone in order to create a mosaic. For this I'm using gdalwarp with the following syntax: C:/OSGeo4W64/OSGeo4W.bat; gdalwarp -t_srs ...
1
vote
2answers
78 views

Reproject raster without projection info

I'm pretty new to gdal and got a presumably easy question for experienced gdal users. My problem: I have a raster ascii file of which I don't know the projection or epsg code. However, I want this ...
6
votes
2answers
5k views

What's the difference between gdalwarp and gdal_merge for mosaicing

I wanted to merge around 20 dem tiffs each around 100MB, and I used gdalwarp (I know about virtual file option): gdalwarp $(list_of_tiffs) merged.tiff However this command took very long time, and ...
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, ...
0
votes
1answer
191 views

How to get projection data (.prj) of raster image using GDAL/OGR?

I start with a boundary shape-file (think U.S. states) in lat/long that I convert to some projection, and then rasterize the result as a .tif. What I need now are the raster coordinates of the ...
1
vote
2answers
248 views

How to compile QGIS from source with Grass?

I recently decided to reinstall QGIS and build against the latest version of GDAL - something I have been putting off for ages because I suspected it would be a nightmare. Lo and behold it has been, ...
7
votes
1answer
88 views

How to index large VRT files in QGIS?

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 ...
6
votes
2answers
170 views

Transforming geostationary satellite image to lon/lat

The ICARE Data and Services Center distributes the SEV_AERUS-AEROSOL-D3 product in HDF5 format (http://www.icare.univ-lille1.fr/archive/?dir=GEO/MSG+0000/SEV_AERUS-AEROSOL-D3.v1.03/). The main ...
1
vote
0answers
34 views

Converting TIFF (DIMAP) to .xyz file format with GDAL

We are attempting to convert a .tif file into .xyz file. We have been using gdal_translate on a linux box and this is the line we've used: gdal_translate -strict -of XYZ ABC.tif ABC.xyz -mo ABC.dim ...
1
vote
1answer
41 views

tiff to shapefile using GDAL - slow speed

I am converting a 500MB .tif file to .shp file using GDAL. I am using the following: gdal_polygonize.py .tiff file -f “ESRI Shapefile” out.shp However, I abandoned the first attempt as in 12hrs it ...
3
votes
2answers
76 views

How to add support for FileGDB (Esri file gdb API) driver in fiona?

There has been a tons of questions on that but none of them solved my problem. I have a machine with: Windows 7 x64 Python 3.4.3 x64 gdal installed: gdalinfo --version --> GDAL 1.11.4, released ...
3
votes
2answers
439 views

GDAL with GEOS support, pandas and VTK using python on Windows?

So far I have been using OSGeo4W shell (or OSGEO install) as really convenient way to use python together with additional packages, namely gdal, shapely, numpy, scipy and matplotlib. I was happy till ...
0
votes
1answer
231 views

“Too many points” when using gdal2tiles.py in EPSG:4326

I want to put a custom map onto cartodb, which, I believes, requires EPSG 4326. When I try making tiles with gdal2tiles, using a 4326 georeferenced image, I get the following error: ERROR 1: Too many ...
3
votes
0answers
21 views

How to blend edges of aerial imagery with gdalwarp

I need to fade/blend the edges of a series of aerial images with gdal. The closest thing that I've been able to accomplish is using the -cutline and -cblend switches of gdalwarp to get the edges to ...
3
votes
0answers
47 views

GDAL RasterIO performance in C++ vs ReadAsArray in python/numpy

I'm currently writing optimisations for a raster viewer program which uses gdal as it's base. It's currently written purely in python, and has some major speed issues which cause problems when we are ...
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
1answer
16 views

Set PATH variables for ogr2ogr without admin privileges

In the office I use ogr2ogr to convert several file formats such as KML or ESRI shapefiles. Since the commands are the same almost every time I'd like to run ogr2ogr with a batch. I'm using the GDAL ...
1
vote
0answers
34 views

Is there any way to connect features (entities) of dxf to database in python?

I'm using gdal, ezdxf, sqlite python packages and sqlite ODBC driver. They work fine, and in AutoCAD I can easily manually connect to my dbase.sqlite (read only) using _DBCONFIGURE and _DBCONNECT ...
3
votes
0answers
24 views

GDAL gdal_polygonize.py TypeError

I am using GDAL to convert a .tif to .shp file. gdal_polygonize.py input.tif -f “ESRI Shapefile” out.shp However this does not work and results in: Traceback (most recent call last): File ...
4
votes
2answers
611 views

LSA-SAF data reprojection with GDAL and R

The Land Surface Analysis Satellite Applications Facility (LSA SAF) distributes several products derived from the EUMETSAT satellites data. I am working with R to read and process the data (the file ...
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 ...