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

UnIntersected part of gdal-geometry?

i have a rectangle and polygon geometries . polygon.Instersection(rectangle) gives me the polygon as it is inside the rectangle. But , what i want is the unintersected(non-intersecting) part of both ...
1
vote
2answers
97 views

How to batch “gdal_merge” on Win7 in one command line command?

I want to create a mosaic of all tif files in a directory on windows 7. How could I create a one-line command in the windows command line to loop through all tifs and merge them all together using ...
1
vote
2answers
179 views

I need to extract a specific layer from a Geospatial PDF using GDAL/OGR

I have GDAL/OGR 1.10 installed on a 64bit Windows 7 machine. I have some Geospatial PDFs that I want to crop to the outline of the map and convert to Rasters. However I cannot clip to the neatline. ...
2
votes
1answer
49 views

GDAL not accepting *.tif to list through files in directory

I am using gdal_merge.py to mosaic tif. If I explicitly write the file names it works, but if I use *.tif it fails with the following error. Traceback (most recent call last): File "C:\Program ...
1
vote
1answer
108 views

install gdal extention in GeoServer

I am trying to install the GDAL extension (I want to use ecw files in GeoServer). So here is what I did so far: 1) downloaded the GDAL extension and extracted to: C:\Program Files ...
0
votes
0answers
73 views

command line alternates for softwares like tilemill, maptiler

tilemill works just fine for my needs but , i want to do the whole process on command line . since new geotiff files are going to be added in my project, and i need to automate this whole process . ...
0
votes
0answers
52 views

gdal2tiles colorpallete error

merged few files using gdal_merge gdal_merge.py -pct -o 713.tif sa713.tif na713.tif wa713.tif ea713.tif When i try to make tiles it says gdal2tiles.py: error: Please convert this file to RGB/RGBA ...
0
votes
1answer
118 views

How to get the pixels from a geotiff file in gdal python for a given bbox?

I have a GeoTIFF file in WSG 84 Projection reference (it's a GeoTIFF of SRTM). I can use python & gdal to read the whole lot of data into an array (with ReadAsArray()), but I want to only have a ...
0
votes
1answer
70 views

Convert a KML ImageOverlay to World Files with GDAL

I have many raster outputs from an Open Source radio transmission topographic path tool called SPLAT, georeferenced in KML image overlay format. I have figured out an easy way to use tools to ...
1
vote
1answer
120 views

PostGIS Raster in MapServer

I'm having trouble getting a PostGIS raster to load and display in MapServer. I'm not sure, but suspect it's a version problem. My mapfile layer and the error message are shown below. LAYER ...
2
votes
2answers
695 views

How to convert a .tif with .tfw to a GeoTiff?

I have a tif image file, and a .tfw world file (created from OSM data with Mapertive). How can I (using linux command line) convert that to a GeoTIFF file? (i.e. I want to have one file, not 2). I'm ...
0
votes
1answer
131 views

Create Tiles for a raster(georeferenced) to use in leaflet

I have raster, with each pixel containing some data. What i want is to split a large image with each pixel's color is determined by its value, to show on leaflet layer.How to go about this ? raster ...
1
vote
0answers
61 views

Extra row and column in a raster after GDAL grid import

I'm exporting an ESRI grid to TIFF format and after that I'm reading this exported raster file into SAGA GIS. However when Saga reads this file in there appears to be an additional row and column in ...
0
votes
0answers
78 views

GEOS assertion failed

The following script fails, if I import the 'rasterstats' library. import rasterstats import ogr geom1 = ogr.CreateGeometryFromWkt('POLYGON ((415924.806155724043492 ...
2
votes
1answer
68 views

How to use GDAL to import lat/lon points from ESRI File Geodatabas into MySQL

I have a file geodatabase from here: ftp://ftp.csc.noaa.gov/temp/MarineCadastre/AIS.SampleData.zip Using ogrinfo, I can see that it contains a layer called "Broadcast_" that has lat/lon points: ...
3
votes
3answers
134 views

how to digitize areas of a raster and snap to the raster cells?

I need to digitize the extent of an ecw raster image where the RGB value is not black.(see attatchment)ie. not the envelope, just the area where the image is shown without the border areas. The black ...
1
vote
1answer
80 views

WKT for local mine grid

I have been trying to create a WKT entry to use in gdalwarp etc. I need an affine transform from coordinates in MGA94 Zone 54 to our local mine grid. All efforts have so far failed, including using ...
0
votes
0answers
28 views

How do i install liblas using OSGeoW4

I have downloaded and installed the OSGeo4W program on my system. I recently downloaded gdal file and liblas. When i try opening the osgeow4 shell, it keeps telling me the gdal file is missing. when i ...
0
votes
1answer
40 views

Missing Georeferencer

I installed Qgis Master 64 bits on win7 64 16gb of ram and I am unable to locate the georeferencer. All other Gdal application are there, Gdal is active, but no georeferencer. I tried yesterday, ...
1
vote
1answer
443 views

Clip raster layer with *raster* mask layer in QGIS

Given two input layers 1. The clip layer - an arbitrary layer (though in my case it's a float 1 band geotiff) 2. The mask layer - a paletted raster in which every pixel is either NO DATA or a single ...
1
vote
2answers
32 views

Comparing projections of two datasets when SRID not available

Similar to this question: How to test if two GDAL datasets are in the same projection? However, the above does not cover cases where an SRID is not available. It can often be the case that an SRID is ...
1
vote
2answers
482 views

os.system /subprocess.call error handling when error in command line function?

In my Python script I call a GDAL utility using: subprocess.call('gdal_rasterize -a ICE_TYPE -where \"ICE_TYPE=\'Open Water\'\" -b 1 -burn 5 -l ' + shapefileshortname +' ' + shapefile + ' ' + ...
0
votes
1answer
207 views

GDAL WMS Question

I'm new to GDAL and am currently trying to hit a WMS using gdal. I need a way to save the image to my desktop. I currently have an xml file with the WMS description syntax that gdal wants -> ...
3
votes
2answers
155 views

How to check which version of GDAL do I have and how many?

I have installed in my computer and running the following software: PosgresSQL 9.3.2.3 PostGIS 2.1.1.1 QGIS 2.0.1 Dufour (this one installed through the OSGeo4w package) QGIS 1.7.0 Wroclaw I am ...
0
votes
0answers
40 views

GDAL - least cost path analysis height data

I would like to use Pythons GDAL bindings to find the least cost path through an array and then show the route on top of the original array. I am using DEM data and would like the slope to be taken ...
1
vote
2answers
251 views

Processing HDF4 data in R

I am trying to process MODIS files in R, which are in HDF4 format (link to data). I have experience reading in HDF5 and netCDF files in R, but HDF4 seems to be a common problem in R. I am using a MAC ...
3
votes
2answers
677 views

How to convert a NetCDF file to GeoTIFF using GDAL

I have a NetCDF file which contains data for the number of rain days in the UK. I'd like to convert it to a GeoTIFF, but can't get GDAL correctly georeference it. See dumps from gdalinfo below. I ...
4
votes
2answers
245 views

how would one get the end points of a polyline?

I am attempting to write a script to get the end points of a series polylines, and then use those end points to create a polygon. I'm not sure which tools to begin to use on this project. I would ...
2
votes
3answers
110 views

Geotiff reprojection

I am new for GIS programming. I try to change geotiff coordinate system from epsg 32635 to epsg 4326. I used the following gdal command: gdalwarp -rc -s_srs epsg:32635 -t_srs epsg:4326 -dstnodata 0 ...
3
votes
2answers
135 views

Python, GDAL, VRT: help wth code efficiency and processing time

I'm going to attempt to make this question as succinct as possible while still including relevant information, as I have trouble with being concise and am trying to avoid writing unnecessary stuff. ...
1
vote
1answer
83 views

How can I use the proj string in java code?

I want to do projection from UTM/oblique mercator CRS to WGS84. And I know the solution to use a proj.4 string in commandline like+proj=omerc +lat_0=51.4 +lonc=7 +alpha=-10 +k=1 +x_0=0 +y_0=0 +gamma=0 ...
0
votes
0answers
30 views

Installing GarminCustomMap: “missing _gdal_array module” error

I've used the GarminCustomMap plugin in QGIS 1.8 without troubles. But now, with QGIS 2.0, when I try to install it, I get the following message "missing _gdal_array module". Any idea how to fix this? ...
1
vote
1answer
145 views

how to create CF-compliant NetCDF from ESRI Grid?

been trying to make a CF-compliant netCDF file from an ESRI Grid. I want to serve it w/the THREDDS server. Been using gdal_translate. I know the netcdf output option is definitely not CF compliant, so ...
1
vote
1answer
95 views

Trying to build GDAL with ECW SDK 5.0 on Windows

I'm using Visual Studio 2010 to open and build the GDAL 1.10 source (using the "makegdal10.sln"). This builds without an issue. Now I want to include the ECW support by including the ECW v5.0 SDK ...
0
votes
0answers
60 views

Update a NetCDF file using GDAL

I have to edit a NetCDF file using GDAL in C++ code. When I try to update a subset I have this error message: ERROR 6: The NETCDF driver does not support update access to existing datasets. How can ...
3
votes
2answers
168 views

Merge two rasters with no data value

I have two rasters: Raster A and B. Raster A has has only values of 0. The white area in Raster B are No Data values. I want to merge them to one raster (Raster C). So basically I want all pixels in ...
1
vote
1answer
162 views

Why does OGR driver.Open() for File Geodatabase give AttributeError: NoneType object has no attribute Open?

There are a number of existing Questions that hint at a possible Answer to this Question, but none seem to be providing me with the solution to my problem: How to access feature classes in file ...
2
votes
1answer
190 views

“ImportError: No module named gdalsrsinfo” What's the issue?

in the python console I can call import gdal, and gdal is recognized as a module; however, most gdal programs, such as gdalinfo and gdalsrsinfo are not recognized. I checked the GDAL source files to ...
2
votes
0answers
90 views

Zonal Statistics using Gdal and Python

I recently came across a loop_zonal_stats function by @ustroetz found here: Issue Trying to create Zonal Statistics using Gdal and Python.. I am using a (single) polygon shapefile and a raster file to ...
1
vote
1answer
66 views

gdal_grid has some problem with reading vrt file

I am using gdal_grid module(from gdal 1.10) and it has some problem with reading the vrt file. I checked this link but it does not work. gdal_grid problem reading vrt file The command is ...
0
votes
0answers
137 views

Add labels for printed map (GeoTools)

I have some problems with labels for printing map with GeoTools. Next code (with some cuts, code on Scala language) can help me for printing map: val map = new MapContent() val parser = new ...
2
votes
2answers
90 views

How to show only raster attribute table using GDAL?

I have a binary ArcGRID raster (AIG/Arc/Info Binary Grid) with an attribute table. I would like to get only this attribute table but gdalinfo outputs also projection information and bounds and then ...
1
vote
0answers
62 views

Write raster with its attribute table and symbology

I am classifying lots of rasters and I'd like to automate the legend creation and symbology displaying for easy manipulation. Here is the kind of raster I'm dealing with: library(raster) r <- ...
0
votes
0answers
21 views

Is there a way of converting polygon vertices to points using GDAL python? [duplicate]

I have an ArcGIS standard license and I want to convert polygon vertexes to points in a script. Since I don't have a good enough ArcGIS license I've been using ET Geowizards to get the points. The ...
0
votes
2answers
70 views

Is it possible to alter raster values without Arcpy? [closed]

I'd like to be able to alter raster values, through fairly simple equations, without using the module Arcpy. Is that possible? It seems like gdal, numpy, and rasterlang might make this possible but ...
1
vote
2answers
192 views

how to georeference my png using gdal

Hello I have created my own fantasy map, I have the map in a full size png and the map broken down into tiles ( with no geo data ) I installed MS4W and tried to use it to get the georeference data ...
1
vote
1answer
66 views

gdal_merge Sum/Agregation files into single output band

I need to create a single tiff files from several others (each have the same extent and resolution) and the output file need to have a single band containing the aggregated/sum of the values in the ...
0
votes
0answers
46 views

How to assign no value to pixels with (255,255,255)

I have Imagery where I want to replace white edges (255,255,255) with no value. I have used gdal translate within QGIS and put in no value 255 but this give no value to other pixels eg 255,240,200. ...
0
votes
1answer
24 views

Break Mosaic into Large Tiles

I have 607 aerial images of Massachusetts in .jp2 format. They are in a mosaic however I need to break down the 607 images into groupings of 16(4x4) spatially so that I can more easily process ...
4
votes
1answer
196 views

Why is tiff different in epsg 3857 and epsg 900913 in gdal

When I transform tiff from 4326 to 3857 and from 4326 to 900913 in gdal, I get different result(when display in Qgis) eventhough these two projections(3857 and 900913) should be the same. What is ...