GDAL is an open source translator library & utilities for raster geospatial 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 useful commandline utilities for data translation and processing.
2
votes
1answer
51 views
How to test if two GDAL datasets are in the same projection?
I'm attemping to write a funciton to test if a list of datasets are in the same projection. Naively, I thought I could test the text of WKTs, but that does not work in the following case:
...
1
vote
2answers
98 views
ESRI State Plane Coordinates as EPSG
I am working with some State Plane projections from Colorado. I have looked them up on spatialreference.org and they are the following:
ESRI:102254: NAD 1983 HARN StatePlane Colorado Central FIPS ...
2
votes
1answer
65 views
How to flip X and Y coordinates of a grid by GDAL
I have a grid, which has flipped X and Y coordinates (X is represented by Y values and Y is represented by X values). This results in grid being positioned in completely wrong place.
How can I flip ...
0
votes
1answer
70 views
GeoTIFF: How to remove Internal Overview without recompression
"gdaladdo -clean" does not remove the internal overview from the file see http://trac.osgeo.org/gdal/ticket/4940
How do I remove GeoTIFF overview(s) without having to re-compress the main image.
1
vote
0answers
24 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
49 views
gdal_grid failing to set nodata
Using gdal_grid for spatial interpolation I have run across a problem: the nodata argument is ignored. The only solution I have found for this is to then use gdal_translate but it is rather cumbersome ...
4
votes
1answer
69 views
How to use -lco options in gdal_contour?
I'm trying to create contours from geotiff raster with gdal_contour. My output format is Postgresql and I would like to ovewrite the existing table.
When I do:
gdal_contour -a 1 -f PostgreSQL -i 50 ...
1
vote
0answers
52 views
How I can incorporate results of histogram stretching into georeferenced image?
I have one channel 8 bit raster and when I open it in QGIS with histogram stretching "on the fly" it looks fine. How I can recalculate original raster accordingly changed histogram?
4
votes
1answer
111 views
Add band to existing GeoTiff using GDAL
This seems like it should be straightforward but I can't for the life of me figure it out.
Given an existing GeoTiff containing a single band, how can you add another band to that GeoTiff? Is the ...
1
vote
1answer
36 views
A way to log the output from the Msys terminal in OSGEO4W?
Is there a way to log all the text output from commands in the Msys terminal to a text file when using OSGEO4W in Windows? I'm using GDAL commands and specifically wish to record any error messages ...
1
vote
0answers
103 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 ...
4
votes
1answer
161 views
Can't export a polygon shapefile to MySQL using ogr2ogr?
I am having problems with ogr2ogr and keep getting different versions of the error:
ERROR 1: MySQL error message: The used table type doesn't support SPATIAL indexes.
Description: ALTER TABLE ...
1
vote
0answers
73 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
2answers
360 views
Leaflet with EPSG4326 Tile Layer (using maptiler.org)
I have successfully been using Klokan Petr Pridal's MapTiler (http://maptiler.org ) to take GeoTiff files and create TMS tile pyramids for Leaflet (option #1 which projects to EPSG3857, also called ...
0
votes
0answers
52 views
Converting different zoom-level GeoTiff to tms
I have several GeoTiff images I want to convert to a TMS structure of .png files.
The GeoTiff images are different scale maps of an area (1:10000, 1:5000, 1:2500) with increasing detail level.
How ...
0
votes
0answers
137 views
How use thin plate spline transformer to wrap image in c#
I'm new at programing with GDAL C#.
I need to geotransform raster image (using GCP's) and bring it to the Merkator projection. With FWTools I can do this by using gdalwarp -r cubic -tps after ...
2
votes
1answer
55 views
How to interpret the OSGeo4w software version under Advanced Install option?
When choosing to install QGIS Desktop using the OSGeo4w advanced install option, I noticed that there are what looks like ambiguous software versions to choose.
For example, currently QGIS Desktop ...
3
votes
1answer
212 views
How to read/write ENVI metadata using GDAL?
I am trying to use GDAL to read and write images in ENVI format (for example a .bsq file with an associated .hdr file). The GDAL Formats Info page states "Starting with GDAL 1.10, all ENVI header ...
0
votes
1answer
50 views
Should a GDAL Plugin be the same version as GDAL aready built into QGIS
Should a GDAL Plugin be the same version as the GDAL that's aready built into the QGIS version installed?
Or,can another different version of GDAL be built and it installed as a Plugin for QGIS, ...
2
votes
1answer
136 views
How do I solve a “No PROJ.4 translation” error with gdal?
I've installed gdal using the brew utility on OS X so the proj and geos dependencies were also installed. I'm now trying to warp a GeoTIFF so that I can import it into Google Earth. I am getting a ...
3
votes
4answers
388 views
Batch processing DEM rasters in QGIS/GDAL
I am attempting to add relief shading to map using .asc raster files, following this tutorial. It works well, but I have only tried it on single tile of raster data which is 20x20km. I want to ...
3
votes
0answers
74 views
Build Raster Attribute Table with GDAL Utilities
Is there a way to build a raster attribute table using one of the GDAL utilities? I know you can view it using gdalinfo and I have successfully created one using Python, but it would be so much ...
1
vote
1answer
56 views
SHAPE FILE to TIGER with OGR
How do I convert SHAPE FILE to TIGER using GDAL (OGR). I try but throw following error:(layer number)is not a known TIGER/Line layer.
3
votes
0answers
152 views
Polygonize vegetation area from Landsat7 RGB using QGis / Bilko?
I'm trying to extract the areas with vegetation as a polygon-shape from a landsat7-image, 30m resolution, RGB, using QGis 1.8 or Bilko or (hopefully in the near future) ArcGIS. The area is in Northern ...
3
votes
1answer
132 views
Simplify polygon in OGR
We use ArcObjects pretty heavily and simplify polygons using something similar to the ITopologicalOperation::Simplify(void). No tolerance is passed in as the parameter (it uses the XY tolerance of ...
1
vote
1answer
96 views
calculating terrain elevation line
I found too many hits on the web to even see the right direction to go.
I have to provide a small proxy server serving requests for a line
of terrain elevation data between two coordinates.
It is a ...
3
votes
2answers
337 views
How to add color interpretation for raster bands using GDAL?
I have 3 separate 1-band GeoTIFF files. The color interpretation for band is gray. I want a 3-band RGB file. I've used gdal_merge.py
gdal_merge.py -separate file1.tif file2.tif file3.tif -o ...
0
votes
0answers
93 views
How to re-project the EASE (Equal Area Scalable Earth) grid with a ~25 km cylindrical projection to WGS84 0.25 degree?
I have nc files ,here is one file https://echange-fichiers.inra.fr/get?k=LUfqiyO1Dl7az5bihnS. from the metadata ,the projection is cylindrical and the resolution is 25 km:
...
1
vote
1answer
184 views
Slow tiles generation for leaflet js
I am currently trying to use leaflet to display images (jpeg or png) instead of maps, I have followed the instructions on this blog post to create the tiles using gdal2tiles.py here is the code we ...
0
votes
0answers
81 views
Regional statistics from a dataset (point inside a polygon aka region)
I am relatively new to GIS and I am trying to do some geographical analysis on same data I've collected. The dataset I have consists essentially in a file with the following format:
...
8
votes
1answer
507 views
How do I get the pixel value of a GDAL raster under an OGR point without NumPy?
I am working on a computational model of the abundance of wild pollinators across a landscape. The model itself is complete, and I am now struggling with a post-processing step.
I have my GDAL ...
2
votes
0answers
98 views
USGS GeoTiff extent via GDAL
I have 2006 land cover GeoTiff from the USGS site (for a county in VA, USA). The projection is Albers Conical Equal Area Units. The download came with an HTML metadata file that tells me the geo ...
4
votes
2answers
176 views
Create subsetted multiband image from multiple files using gdal or python?
Given a time series of an area (single band), how you subset them and combine the subsets into a single multiband file (eg. vrt) using gdal or gdals' bindings for python?
22
votes
2answers
254 views
A comprehensive, user-friendly GDAL manual?
I'm becoming a real fan of using GDAL tools in a bash terminal. However, I quite often get stuck trying to solve particular problems, or getting syntaxes right. Apart from the details provided in the ...
2
votes
0answers
80 views
What means the gdal_translate message: “Too many command options”? [closed]
I have a little question, regarding the gdal_translate-utility:
I want to extract raster-files from database-stored gridfiles using the gdal_translate:
gdal_translate -ot Float32 -of SAGA PG: ...
5
votes
1answer
91 views
GDAL/OGR Determine if dataset is raster or vector
I am writing a script using the GDAL/OGR Python bindings and it is possible that the user may provide a raster dataset or a vector dataset as input. If the user provides a vector dataset, the script ...
1
vote
3answers
315 views
How to subset raster image using gdal?
I have read pixel values of an raster image using GDAL libraries in visual studio 2010(vc++).
Now I have to crop the image (subset) according to the grid given in shape file.
I have got a code for ...
3
votes
1answer
456 views
How do you open shapefiles with ogr2ogr
I'm following along the D3 tutorial, and on this line:
ogr2ogr \
-f GeoJSON \
-where "adm0_a3 IN ('GBR', 'IRL')" \
subunits.json \
ne_10m_admin_0_map_subunits.shp
run into the following ...
2
votes
2answers
202 views
gdalwarp cutline along with shapefile
I am trying to clip my Raster file according to a shape file geometry. I am using the following code
gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif
But it is resulting a black color outside the ...
2
votes
0answers
212 views
gdal ImportError in python on Windows
Trying to import gdal for Python 2.7.3 on Windows XP:
>>> import gdal
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File ...
0
votes
0answers
37 views
Use Overviews to Select Area of Interest in Very Large GIS Rasters/Images (e.g., GeoTiff's)
I am working with large images in Python and Gdal on a Windows machine and wondering if it is possible to directly access the specific smaller overviews I built into them using gdladdo. I am viewing ...
1
vote
1answer
288 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
2answers
169 views
How to work with GDAL/OGR in Tomcat?
I make web application in Eclipse using GDAL. In Eclipse in Debug Configurations i add variable path=path/to/gdal/folder. In project i add gdal.jar to build path and copy gdal folder in lib directory. ...
4
votes
2answers
189 views
How do I Capture GDAL projection and create matching output (OGR or GDAL) shape files?
Working in C++, using GDAL and OGR: Depending on what the user has supplied for image input, I need to output point (or other) shape files in the same projection as the input raster images. The ...
3
votes
1answer
307 views
Rasterize shapefile to a new gtiff file with gdal_rasterize utility,however,I cannot open the tif file?
Now I'm tring to rasterize a shapefile with a utility "gdal_rasterize" and want to generate a new gtiff file,but when I open the rasterized gtiff file,it prompts "Windows Photo Viewer can't open this ...
9
votes
2answers
334 views
Gdal Dataset.ReadAsArray() crashes Python
I am using Python 2.6.5 (32bit) with Numpy 1.3 and Gdal 1.9.1 installed on Windows 7 64bit. I am trying to read an 800 MB Imagine (.img) raster dataset into a Numpy array to do some raster algebra, ...
0
votes
1answer
58 views
GeoEye BGRN --> RGB
I got many GeoEye files with BGRN color format. The image looks gray.
-co PHOTOMETRIC=RGB ALPHA=NO gives greenish color. Rearranging bands does not work too.
Could someone help me convert this image ...
5
votes
2answers
242 views
Where can the GDAL executables be found in a QGIS installation?
I run QGIS 1.8.0 Lisboa on a Windows 7 machine. Yesterday I attempted to install the TauDEM plugin using these instructions: http://lists.osgeo.org/pipermail/qgis-developer/2012-June/020807.html
...
7
votes
1answer
130 views
Multiple versions of GDAL (32bit and 64bit) on same machine
I currently have GDAL installed on my Windows 7 64bit workstation. It is the 32bit version. I chose the 32 bit version because I wanted to use the GDAL Python bindings with a 32bit installation of ...
3
votes
1answer
141 views
Python equivalent of gdalbuildvrt
Is there a way to perform the same task as the gdalbuildvrt utility using GDAL Python bindings? So far I have not found any way to do this other than creating a vrt of a single dataset and manually ...
