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

0
votes
0answers
34 views

gdal_grid performance issue — alternatives?

I use gdal_grid within a Python script, but for my large files the calculation time is enormous. Apparently this is a gdal_grid issue http://trac.osgeo.org/gdal/ticket/2411 Does anyone have an advice ...
0
votes
1answer
48 views

going from hillshade to percent shaded (gdaldem)

I assumed hillshades output percent shaded (scaled to 255) at the given azimuth and altitude but this is incorrect. A flat raster at 80 degrees altitude has a maximum value of 251 instead of 255. How ...
0
votes
1answer
44 views

Gdal_sieve not working?

I'm using the Sieve command within QGIS (Raster > Analysis > Sieve), but when I run the command with any variation of parameters, it completely changes the histogram of my raster data from between ...
1
vote
0answers
29 views

Sending commands into FGDB

I have a large .gdb of all blocks in the United States, available from Census. I would like to create a point shape file of all census tract centroids in the US. I have QGIS with FileGDB enabled, as ...
8
votes
2answers
125 views

How can I use data from PDF map?

I have a map over Sweden with plenty of data represented as circles of various sizes. The map is saved as PDF from Illustrator, fully editable. I wonder if there is any way I could get these circles ...
3
votes
3answers
162 views

Fast and robust open source portable script or tool for interpolation - triangulation or similar

I am looking for some fast and robust open source interpolation library/script/tool. The goal is to re-interpolate relatively big (1M+ points) unstructured (xyz) grids into high resolution cartesian ...
1
vote
2answers
74 views

Using very large image formats with GDAL

I would like to use gdalwarp and gdal_translate with some very large raster images (313,040 x 313,040 px) and I'm wondering what the best format or method is to do this, if it's even possible. I have ...
1
vote
0answers
62 views

GDAL WMS in QGIS 2.4

I'm trying to add in QGIS a layer shown here : propluvia When I watch in my console I see this request corresponding to the layer I need. The server is not able to respond to a getCapabilities ...
1
vote
0answers
43 views

Converting MDB File With Unicode Characters

I'm using C# wrapper of GDAL and I'm having a problem converting MDB file that contains unicode characters. Lets take for example conversion to KML file, non-English characters are replaced by '?' ...
2
votes
0answers
45 views

Error in polygonising raster data with QGIS on Mac

I'm using QGis Chugiak on a Mac and I'm trying to converting a raster layer in a vector data by: Raster/Conversion/Polygonize. So, I select the layer needed from the drop-down menu, I give a name to ...
3
votes
2answers
117 views

GDAL (Python) Aggregate Raster Into Lower Resolution

I have a (continent-scale) raster with 30m resolution in an Albers equal area projection. Each 30m cell has a value 1-34. I would like to aggregate information from this raster into a 1km grid such ...
4
votes
3answers
233 views

Install gdal in CentOS without Root

I have been struggling to get python setup to do some image processing, so I really need the http://www.gdal.org/ GDAL libraries installed on CentOS 6.5. I want to install the package in a local ...
2
votes
1answer
266 views

Problem in displaying tile images using openlayers3

I have problem with displaying my TMS layer using OL3. This is my first time working with tile images this way so maybe it doesn't have to do anything with ol3, maybe I just don't understand it right. ...
2
votes
1answer
47 views

Merging h5 files from VIIRS to gtiff

I am trying to mosaic all the data from VIIRS of a day into a single gtiff. I've already done it with MODIS data using GDAL. The problem it seems is that the coordinates of each h5 files doesn't seem ...
1
vote
2answers
66 views

Is there any efficient way to create a raster with pixel value representing a unique id calculated based on pixel lat/lon?

I have 10m raster for each state in US. Current pixel value represents something. I want to have 'same' raster but pixel size representing a unique id calculated based on pixel lat and lon. I am ...
1
vote
1answer
76 views

raster2pgsql does nothing

I try to load rasters through raster2pgsql. I try to call it through pgAdmin command line, but when I call any raster2pgsql command (even raster2pgsql -G), it does nothing, I just get another line to ...
0
votes
0answers
11 views

How to read *.map file by using geotools api in Java? [duplicate]

How to read *.map file by using geotools api? String filePath ="C:\\CTO-New\\SOCLOMO\\Location\\Map Info FIles\\New folder\\18-july-2014\\1301\\filename"; File file = new File(filePath); ...
1
vote
2answers
66 views

How to restore gdalinfo metadata after ImageMagick processing?

I first print out the file info : $gdalinfo crop.origin.tif result: + gdalinfo crop.origin.tif Driver: GTiff/GeoTIFF Files: crop.origin.tif Size is 1920, 1950 Coordinate System is `' Origin = ...
2
votes
0answers
44 views

(How) can I write string data to a raster with Python/GDAL?

In Python, I have a NumPy array of string values (categories of some variable) that I'd like to write to a raster file (.tif or .img) using GDAL. The impression I get is that it's not possible to ...
1
vote
1answer
105 views

Convert GML to DXF using ogr2ogr

I cannot make ogr2ogr.exe work. I'm using GDAL 1.11. I keep getting the following error message: ERROR 1: Did not get at least 2 values or invalid number of set of coordinates 709993.98 ...
1
vote
1answer
54 views

Convert list of coordinates to OGRGeometry or WKT

I am using a Visvalingam-Whyatt simplification script and the output is a list of coordinates such: [[[[-10885205.10690245 3525662.26531131] [-10885369.01690829 3525374.0741439 ] ...
2
votes
1answer
45 views

Error reading large virtual raster with rasterio: “Cannot invert degenerate transform”

I've created a virtual raster using gdalbuildvrt, which covers a large extent (England) at a high resolution. The data is fairly sparse, focused around the coastlines with no files covering inland ...
1
vote
0answers
33 views

Rendering a geotiff block with gdal and gdk-pixbuf - Some flaws in the result

I currently try to render a geotiff file with GTKmm\GDKmm (gdk-pixbuf). The file is a tiled geotiff with a block size of 256*256 and three rasterbands. The shortened gdalinfo output is: Driver: ...
2
votes
1answer
76 views

Are there examples of high-level interfaces to NetCDF

The R raster package and GDAL both have independent wrappers around NetCDF-4* that treat an array variable as a "complete" dataset, i.e. they provide objects or interfaces that seamlessly give access ...
-1
votes
2answers
74 views

Converting *.kml file to *.lyr file so Google Earth can create spatial data to use in ArcGIS for Desktop? [duplicate]

Want to convert the .kml file into .shp file so that I can create new the point, line and polygons data to work in ArcGIS. Actually I want to use Google Earth to make data for ArcGIS.
2
votes
0answers
50 views

Memory leak with GDALOpen?

I'm very new in using GDAL, and trying the tutorial given here. But with a piece of code with only GDALRegister_GTiff(); GDALDataset *poDataset = (GDALDataset *) GDALOpen(argv[1], GA_ReadOnly ); ...
2
votes
1answer
224 views

Is there any way to fast read pixel value using gdal in C#

I create with C# a series of lat/lon value pair in the code(using i, j loop) and I want to extract the pixel value based on either lat/lon pair or pixel location (i, j). I know there is one gdal ...
1
vote
0answers
97 views

Correct way to extract pixel data using GDAL. Current method gives repeated but changing value.

I am trying to open a spectral file and move bands of data from GDAL to OpenCV's Mat data structures. As of now I am just trying to make sure that I can correctly access each band using GDAL and pull ...
1
vote
0answers
64 views

ogr2ogr “Error: pointer being freed was not allocated” on mac mavericks

I've successfully installed GDAL via homebrew on mac But while I'm using command ogr2ogr, keep getting error below $ ogr2ogr -f geoJSON tl_2010_20_tract10.json tl_2010_20_tract10.shp ...
0
votes
1answer
31 views

Gdal_grid gives blank .tiff file after conversion

After using gdal_grid I get blank .tiff file. I use the following .vrt file : <OGRVRTDataSource> <OGRVRTLayer name="p"> <SrcDataSource>p.csv</SrcDataSource> ...
1
vote
1answer
186 views

how to use gdal_rasterize? how to know which attribute to burn with ogrinfo?

I have a vector layer I need to rasterize. I guess this is a basic operation, but I haven't found the information I need yet. First, I have reprojected the vector shapefile: ogr2ogr -s_srs ...
1
vote
1answer
72 views

geotiff conversion to coordinate system with gdal

I would like to get the coordinates X Y Z from a geotiff file. I am using the gdal translate for the conversion. Which output format is the more adequate for containing a 3 numpay array X Y Z ?
2
votes
1answer
72 views

How to convert ASCII in x,y,z format using gdal_grid?

I have a map saved in (AFAIK) Arc/Info ASCII format and I would like to convert it with gdal grid but I am getting errors. My data (p.csv) is in the following format: 319053.25,319301.12,544.28 ...
0
votes
0answers
70 views

Merge/mosaic categorical rasters in QGIS/GDAL

I am attempting to make yearly bins of some monthly MODIS data for all of the US. The cells of the raster are defined with Julian dates indicating the date at which a burn event took place. After ...
0
votes
0answers
43 views

Read esri shape point data and spline interpolation using python shell scripts?

I am quite new in python scripting. I am a dedicated user of arcgis but this time I have to work without arcgis. I have a point shapefile with 10 different attribute fields of different peat layer ...
1
vote
0answers
58 views

How to read .adf (TIN) map files?

I have a set of maps in few formats : TIN ├── [4.0K] info │   └── [2.6K] arc.dir ├── [ 708] log ├── [374K] M-33-45-D-d-4.asc └── [4.0K] m-33-45-d-d-4tin ├── [ 96] tdenv.adf ├── ...
0
votes
1answer
77 views

Does gdal_calc.py not calculate arithmetic with in the same file?

I have one GeoTIFF file with 2 bands. The data type is Byte. I want to calculate (Band2-Band1)/(Band2+Band1). But the following command provided a dark image: gdal_calc.py -A infile --A_band 1 -B ...
0
votes
0answers
57 views

Read large raster into lat/lon value pair to csv file using python and/or gdal?

I have rasters for all of the states in US and they are in 10m resolution. So they can be very large, containing a large number of pixels if you think of TX or AK. I want to create a csv file with ...
2
votes
3answers
114 views

Use GDAL command line to copy projections

I have two images, a georeferenced Tiff and a second Tiff without a reference. The second image is a copy of the first, but it lost its reference. Now, I want to copy the reference from the first Tiff ...
1
vote
0answers
39 views

Clip Aerial Image Pair by Overlapping Area?

I have a set of aerial images which I am trying to clip by their overlapping areas, for use in a more involved program I am writing in Python with Python Imaging Library (PIL). Is there a command ...
0
votes
1answer
41 views

Trouble with r.out.gdal in GRASS GIS

I'm processing some landsat 8 files. I'm using this scene LC82260982013325LGN00 and I have to apply i.landsat.toar and export the result as a GTiff. However when I use r.out.gdal it takes a very long ...
3
votes
1answer
109 views

Join CSV using gdal/ogr

I'm trying to use gdal/ogr to join a shapefile called myshape.shp with a csv called mycsv.csv. The join field in the shape in joinshp and the join field in the csv is joincsv. I've found this command ...
0
votes
1answer
31 views

GDAL: Update spreadsheet from polygons in a shapefile

I am a beginner with GIS data, and have been looking for an answer but have not yet found it. I have a spreadsheet with lat/long points. I also have a shapefile with several numbered polygons over a ...
1
vote
0answers
85 views

Rasterize shapefile using another raster's size and resultion

Rasterizing a vector (shp) I face the following problem: I want my output file (the rasterized vector) to have the same resolution and dimensions as an existing raster. If I understand the interface ...
0
votes
1answer
39 views

Adding a Vertex To Exisiting Feature(linear) Using OGR Python

I am working on a project that involves finding the intersection of linears in a shapefile, if that vertex(the found intersection) doesn't already exist for each(or one) of the features, then that ...
0
votes
2answers
106 views

gdal2Tiles output differs from MapTiler

I'm in need to process a large batch of Landsat images with MapTiler and since the version I'm using (1.0 Beta 2) does not support multiple input files, I created a simple python script that calls ...
2
votes
1answer
38 views

Invert gdalwrap cutline

I have a big raster image with medium details, but for some regions I do have far more details as a seperate raster image. Now I'd like to display both of them in my map that I render with Tilemill. ...
1
vote
1answer
63 views

Holes between imagery tiles projected using gdalwarp

I am having the following issue when reprojecting imagery tiles using gdalwarp. These are 4096x4096 imagery tiles, which I'm gdalwarp'ing from WGS84 to a State Plane projection. When reprojected, ...
0
votes
0answers
61 views

Polygonize only pixel values of 1 from a raster

I have a raster mask that I have created in QGIS (2.2.0) that consist of 0, 1 and nodata values. I would like to end up with a polygon dataset that consist of only the value of one pixels from this ...
3
votes
1answer
100 views

Projection mismatch when plotting precipitation with ggplot

I'm trying to plot accumulated precipitation in Brazil for a given period, see code below: Load packages library(raincpc) library(SDMTools) library(raster) library(ggplot2) library(rgdal) ...