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

How can I maintain MapInfo GRC file attributes in QGIS/GDAL?

I am trying to manage a classified MapInfo grid file in QGIS - I'm wondering if anyone has any experience with maintaining the attributes of the classified grid? QGIS seems to assign them just a ...
1
vote
0answers
336 views

Topography: How to convert raster into vector layers?

(nb:This project is part of the #Wikipedia #wikimaps project) 1. Question: How to create vector topographic maps from raster image ? How to convert raster topographic image into closed vector ...
1
vote
0answers
263 views

gdal2tiles.py problem - creating mini tiles

I have a similar problem to the issue posted here (gdal2tiles problem). I have 80 GeoTIFFs (4000x4000) with EPSG:27700 projection embedded and TFW and PRJ files. These I have built into a virtual ...
1
vote
0answers
179 views

problem with WFS Data, GDAL, Python, ArcGIS, & Defining Projection

I'm pulling in data in gml2 format with a EPSG#4267 from a wfs. I create a datastream in python and then run the following from the python window in ArcMap: data = datastream.read() #write the gml ...
1
vote
0answers
185 views

QGIS 1.8.0 on Linux treats NetCDF 4 as HDF5

I've been generating some large (28800x13920) grids of BYTE data representing the NALCMS2005 North American Land Cover dataset (original 250m projected data resampled to a WGS84 grid at 15 arc sec). I ...
1
vote
0answers
146 views

How to create multi-image single band single file geotiff?

I have a geotiff file that is single band but contains multiple subdatasets. I use gdal with numpy to open a specific subdataset and mask some of its values. Example: $ gdalinfo sample.tif ... ...
1
vote
0answers
103 views

gdalwarp not clipping neatline properly

I have a PDF file (10mb) which contains a NEATLINE polygon in metadata. $ gdalinfo I5601_GeoPDF.pdf Driver: PDF/Geospatial PDF Files: I5601_GeoPDF.pdf Size is 5906, 3543 Coordinate System is: ...
1
vote
0answers
137 views

How can I combine these two commands? (Get a clipped raster ASC-II file from PostGIS)?

My objective is to clip a polygon from PostGIS raster table, and then save this raster (draw polygon shape) as an ASC-II in a folder. What I did now are 2 steps. 1st Step: Insert the clipped raster ...
1
vote
0answers
34 views

GDAL Plugin with a different version of GDAL aready built into QGIS

My qgis windows version is 1.8.0 Lisboa and use the gdaltools in version1.9.2. I want to use the plugin NetCDFBrowser, but this plugin need the gdal 1.10. I install 1.10 stand-alone gdal find here ...
1
vote
0answers
92 views

Python GDAL/OGR: check lines geometries connectivity

I'm working with shapefiles representing roads, so basically it's stored as (poly)lines. I need to perform a road connectivity checking. In other words, I want to ensure that roads meant to ...
1
vote
0answers
127 views

setup Mapserver WMS from tile files

I'm totally new to Mapserver and I've been banging my head trying to figure out how to mount a WMS server that will seed images of a small geographical area to a OpenLayers HTML page. I already ...
1
vote
0answers
507 views

Using gdal to convert from UTM to equidistant lat/lon

I am getting rather confused with converting a UTM gridded dataset into a equidistant/rectangular latlon gridded dataset. Can anybody help me? I have a UTM set with the following projection: ...
1
vote
0answers
132 views

How can I call gdal_translate from Java code?

How can I execute gdal_translate in Java? For example, will something like this work? gdal_translate "PG:host=localhost dbname='xyzy' schema=public' table='xyz'" "d:\xyx\ab.tiff"
1
vote
0answers
166 views

using GDAL/OGR (faster)

Eventually, I was able to get the GDAL installed and compiled correctly using Visual C++ 2010. I even tested some of its methods on some geo-referenced raster files, and it looks like that all is ...
1
vote
0answers
118 views

Normalize with gdal_translate

I do not understand the command line when I try to work with gdal_translate in QGIS. I found this at the forum: Normalize raster values to 0-255 And tried this: *"If not, from within QGIS choose ...
1
vote
0answers
199 views

How can I fix badly specified geographic coordinate system with units in arc-seconds?

Check out this depth-to-rock soil data. It is linked from here. The spatial reference system is a latlon, but the units are arcseconds. I could not reason out how to override the SRS even after ...
1
vote
0answers
174 views

Question about geotiffs

My colleague and I are georeferencing a bunch of .tif images. After each image, we'll click "update georeferencing". The resultant .tif files are different between our computers and I'm trying to ...
1
vote
0answers
122 views

Are there limitations with gdal virtual raster option?

I am new at Postgis and I'm groping... I pretend to rasterize some weighty vector layers in Postgis, do some map algebra and reclassification with them and visualize intermediate as well as final ...
1
vote
0answers
94 views

Coarse heat map placement using gdal

I have generated a coarse heatmap of where data is available by 1 degree x 1 degree cells that I want to overlay on top of a world image. Using a python script I merged two images (heat map and a ...
1
vote
0answers
56 views

How do I link an existing install of GDAL to FileGDB in Linux?

The distribution I use (CentOS) only has GDAL 1.8 available through elgis. I've downloaded the FileGDB API from ESRI and I'd like to be able to link it after an install. How can I do this?
1
vote
0answers
185 views

How to use GDAL to burn a selected value from one raster into another?

I think this should be a relatively simple process, but I'm not quite sure how to approach it. I have a shapefile containing 1000's of polygons. I also have 19 rasters. All of the polygons fall ...
1
vote
0answers
101 views

OSGB36 referenced PNG to Google Earth

I am creating a PNG that is referenced to the OSGB36 datum (1 pixel per 50m easting/northing that is colored according to the result of some calculation code) I want to display this on Google Earth ...
1
vote
0answers
332 views

OGR-Driver GML axis order handling (GDAL, QGIS)

By default the OGR-Driver will swap the axis order to longitude, latitude for geographic coordinate systems. "Since OGR 1.8.0, the GML driver has coordinate system support. This is only reported ...
1
vote
0answers
82 views

Size of output geotiff file not as in the .vrt file

I have the following vrt file (input.vrt) with GCPs in EPSG:3857 <VRTDataset rasterXSize="5409" rasterYSize="5408"> <Metadata /> <Metadata domain="IMAGE_STRUCTURE"> <MDI ...
1
vote
0answers
109 views

VSIMalloc3() and UnifiedSrcDensity mask errors when running gdal_warp in OSGEO4w shell

This is related to gdal.py as part of bat script in windows 7 I am getting ERROR 2: VSIMalloc3(): Out of memory allocating 854522136 bytes. ERROR 2: Out of memory allocating 119813400 bytes ...
1
vote
0answers
267 views

ArcGIS 10.1 crashes when a button in an addin invokes GDAL DLL through SWIG CSharp wrapper

I am developing an addin with a button that tries to do coordinate transformation using GDAL OSR. When I create a new CoordinateTransformation object ArcMap crashes and close. If a debug against a ...
1
vote
0answers
80 views

resampling errors

I have to reproject my images, but it contains some resampling errors. Any one has any idea on significance of resampling error. Are they neglectable? Is there any good documentation to refer to? I'll ...
1
vote
0answers
73 views

Create 6 band IMG file with GDAL out of two JPG files, one with WorldFile coordinates

I have two JPG files 1 ortomap: http://imageshack.us/photo/my-images/703/ortofoto1.jpg/ with WorldFile: 0.683772574816354 0.000000000000000 0.000000000000000 -0.683739485024190 ...
1
vote
0answers
155 views

`gdalbuildvrt` and selecting individual bands

I want to produce a set of VRTs (GDAL's internal virtual rasters) that comprise a re-arrangement of bands of a number of individual GeoTIFF files. Say I have my initial set of GeTIFF files, all having ...
1
vote
0answers
94 views

Calculating tiles extents using some level data

I have some data looking like this: <?xml version="1.0" encoding="utf-8"?> <levels> <level nr="0" tileWidth="512" tileHeight="512" envSizeX="5000" envSizeY="5000" minX="465000" ...
1
vote
0answers
233 views

Homebrew GDAL installation ECW or MrSID

I have experimented with installing gdal 1.9.1 via homebrew. It works wonderfully. Unfortunately, their is no real JP2000 support. I realize that JP2000 is supported by OpenJp2000, but that ...
1
vote
0answers
172 views

gdallocationinfo from C++?

I need a C++ API to convert x,y Image Coordinates in their native projection to output image pixel,line coordnates. I see the gdallocationinfo documentation, but no API for it. Where in online GDAL ...
1
vote
0answers
125 views

Is Django's LayerMapping capable of putting a static value into all features?

I'm using GeoDjango's LayerMapping utility to try and import a set of features. The catch is that aside from the attributes in the source layer, I want all the destination features to have the same ...
1
vote
0answers
309 views

ValueError when reading a raster band as a numpy array

Here's the problem shown in Python shell: In [32]: from osgeo import gdal In [33]: geotiff = '/Users/calvin/work/luc/proj_public/media/GLO ...
1
vote
0answers
143 views

Geoserver Custom CRS IAU_Moon_2000

I am trying to setup a custom basemap using a GTiff in the IAU's Moon2000 projection (30100 or 30101). The basemap image was originally created through gdal_merge and gdalinfo returns the expected, ...
0
votes
0answers
4 views

Python 2.6 Module for GDAL 1.9.2, ImportError

I work on an internal network(IE I cant get the latest 1.11 stuff) and have GDAL 1.9.2 setup on a Red Hat Linux Enterprise Server 6. My linux box has python 2.6.6 however the python binding I received ...
0
votes
0answers
20 views

Batch Process - How to merge images per day using gdal-merge and cmd(.bat file)

There are 4 to 5 projected images in one day. Let's say I have all the projected images for 3 days. For day 1, the file name are T2010232034000.L2.Tera.tif, T2010232035000.L2.Tera.tif, ...
0
votes
0answers
54 views

PALSAR 50m Orthorectified K & C Mosaic it's possible to create interferometry?

JAXA is providing free access to PALSAR 50m Orthorectified K & C Mosaic. According to the website, data are provided in raw (HDR) and png. ALOS PALSAR K&C Mosaic or New global 50m-resolution ...
0
votes
0answers
8 views

Ignore empty tiles in gdal2tiles.py when generating tiles

I am in the process of generating a large set of tiles from a .vrt file composed of many images. The image is following a road section, being very long and small strip of images and due the shape of ...
0
votes
0answers
10 views

gdalbuildvrt not working with star in path on windows

In my directory C:\dev\Jpeg-geo\del1 I have a set of .jpg and .jgw files that I am trying to put together with gdalbuildvrt but its giving me a weird error. ...
0
votes
0answers
15 views

Adding a column to several feature classes in a file gdb using ogr2ogr

I'm trying to add the same column to 25 feature classes in an ESRI FileGDB. I managed to do this for one feature class by using ogr2ogr: ogrinfo C:\temp\test.gdb -sql "ALTER TABLE table01 ADD COLUMN ...
0
votes
0answers
23 views

Precision loss when transforming coordinates systems in GDAL (vector layers)

Anybody knows the order of magnitude of the error introduced when transforming coordinates using GDAL? More specifically, I would like to know if significant error is introduced when the ...
0
votes
0answers
27 views

Dxf to Shp - Coordinates wrong

I've converted a dxf file to a shape file using the folowing ogr2ogr command: ogr2ogr -s_srs EPSG:31464 -where "LAYER='0'" -f "ESRI Shapefile" dgm.shp DGM1.dxf Or ogr2ogr -where "LAYER='0'" -f ...
0
votes
0answers
15 views

Confuse about how to make a tiles to support retina device

I have a map picture(15772*13796) so I break them to tiles using gdal2tiles(512*512px not 256*265px) for support retina device but is this correct? Because I a had read this issue and this. They ...
0
votes
0answers
31 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
0answers
55 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
37 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 ...
0
votes
0answers
51 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 ...
0
votes
0answers
51 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 ...
0
votes
0answers
26 views

Gdal_Calc.py random pixel error

I'm using the gdal_calc.py in python in order to do the same thing as the Con tool in ArcGIS. I'm calling the .py using os.system: os.system('gdal_calc.py -A%s -B%s --outfile=%s ...