Tagged Questions

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
2answers
77 views

How to overwrite metadata in a tif file with GDAL?

I have a set of tif files generated from a process with multiple steps. Some metadata tags are added at different steps. I would like to update the metadata of those file before distributing the ...
0
votes
2answers
30 views

gdalcompare.py not found in GDAL 1.10.1

I cannot find the script gdalcompare.py within the osGeo distribution of GDAL by http://www.osgeo.org/. It is the 1.10.1 though. Is there any way of downloading that script separately from any other ...
2
votes
0answers
30 views

Improve gdaldem hillshade output (lighten up flat parts)

I have a hill-shading that I made using gdaldem hillshade which is already pretty good. I used the option -alt 60 to get nice contrast, but the problem is, that with lowering the light-source, also ...
0
votes
2answers
96 views

Conflicts while installing postgresql9.3-postgis2.1,QGIS and gdal-bin

I want to use postgresql9.3-postgis2.1 and gdal-bin parallel . I have refer instruction here for installing postgresql-9.3 on ubuntu 12.04. I have install qgis from here While installing gdal-bin it ...
2
votes
1answer
113 views

How to correctly resize raster GIS images to a given px width?

I work on ETOPO1, a 466MB file, and want to move to SRTM (~20GB Earth). My workflow involve several time-consuming processing steps so I strongly want to work on a minimum raster. It usually want to ...
1
vote
1answer
19 views

why gdaldem is not producing any result?

I have a geotiff file with class values (0, 1, 2, 3, and 4). I want to show each class as different rgb color. So, I created the text file (color_relief.txt) which contans the following: 0 55 78 94 ...
0
votes
0answers
32 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
41 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
31 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
121 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
148 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
65 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
55 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
40 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
44 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
92 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
138 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
203 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
39 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
58 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
64 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
54 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
41 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
82 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
51 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
43 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
28 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
72 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
61 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
45 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
154 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
82 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
37 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
26 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
116 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
60 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
61 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
59 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
38 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
48 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
59 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
53 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
92 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
35 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
40 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
98 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
30 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
62 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 ...