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
3answers
146 views

Extract raster values gdal looping rasters

I have this code, most of them getting from this forum. I need to loop over a folder and sub-folders, or just to run a list (putting previously all the raster inside, I don't mind one way or the ...
1
vote
1answer
497 views

Using GDAL to read dted files

How do I get the elevation at a certain point using gdal to read in a dted file? I don't see much useful information out on the web. I'm sorry I should have been more specific. How would I ...
2
votes
2answers
292 views

Driver Error : DTM to Raster using GDAL_GRID

I am looking to interpolate a DTM (eventually over 4,000 DTMs) that is in a personal geodatabase format. I would like to do this in as few step as possible, eventually writing a batch script to do ...
0
votes
1answer
156 views

Combine overlapping rasters in gdal_calc.py, taking the maximum value from each pixel?

I want to combine 3 overlapping rasters to create 1 raster. Each have the same dimensions, resolution and CRS. I want to populate the new raster with the maximum value from each of the layers. ...
2
votes
1answer
35 views

How to define a projection for a file with geo points?

I have a HDF file that I want to convert to tif or ENVI file. I know that the +Proj of the original file is: "+proj=cea +lon_0=0 +lat_ts=30 +ellps=WGS84 +units=m" so: gdal_translate -of "ENVI" ...
1
vote
0answers
19 views

raster tiles from vector data with GDAL - how to avoid aliasing artefacts?

Prequesites I have vector data in a Shapefile N47E007_CONTOURS.shp produced with wget -O N47E007.zip http://e4ftl01.cr.usgs.gov/SRTM/SRTMGL1.003/2000.02.11/N47E007.SRTMGL1.hgt.zip unzip N47E007.zip ...
2
votes
2answers
143 views

gdal_grid unable to get nodata outside of original data area

I have a vector based csv file that I want to convert into a grid using gdal_grid, but I want to force null data values strictly around the area where there is data. I can get gdal_grid to work but ...
1
vote
0answers
17 views

Output formats not showing in heatmap plugin (QGis 2.10.1 in ubuntu 14.04)

I just upgraded to QGIS 2.10 (using the debian repos) and the heatmap plugin doesn't show any formats in the dropdown menu. I'm able to open every raster format I have at hand and, of course, the ...
0
votes
1answer
18 views

Dtm format save only heights

I need dtm format that store only heights and scaling and left up position for saving storage I glad to know if any stadard save heights that way
0
votes
1answer
9 views

C# ogrinfo with virtual vector source doesn't work GDAL 1.10

I am trying to use GDAL 1.10 with a vector virtual source which is nothing more than the example union of two shapefile layers. from the commandline, ogrinfo works fine on the source. But the C# ...
1
vote
2answers
109 views

Using GDAL to get vertices of a multipolygon

I am successfully using GDAL to get the vertices of a polygon from a ShapeFile using .NET code. This link showed me how to do that: ...
1
vote
3answers
320 views

How to subtract rasters using the gdal_calculate command

I just want to know how to subtract two raster using gdal_calculate command? Or what is the proper way of subtracting using this? Any idea would be greatly appreciated. Note: gdal_calc is different ...
6
votes
1answer
748 views

How to read Geotiff tags and pixel per pixel values using GDAL C#

I'm new here and I'm new at programing with GDAL C#. My question is how I can get the value of SRTM GEOTIFF pixel to compare with it's neighbors, and how can I get the values of Geotiff tags.
7
votes
5answers
2k views

Are there GUI resources for gdal and ogr?

For those who are new to the programming side of GIS, I feel that a lot of great tools and resources are inacessible. I am interested in GUI based wizards for handling geoprocessing with libraries ...
0
votes
1answer
35 views

GDAL -SQL or -where options

ogr2ogr -overwrite -t_srs EPSG:3978 -f "ESRI Shapefile" -sql "SELECT * FROM 'Layer Name' WHERE FIELD_NAME = 'A CLASS'" dst.shp src.shp help me out here - looking to 1) reproject a shape and 2) ...
0
votes
1answer
16 views

Why does Spatialite layer prompt for CRS if SRID was not explicitly defined in Layer Options?

In QGIS, if I Save As... a SpatiaLite file and don't enter an SRID in Layer Options, I will be prompted for a CRS upon opening the file. Yet GDAL documentation says that this option is only for ...
2
votes
2answers
43 views

Finding the pixel coordinates using GDAL/numpy

Is there a way to get the XY cordinates of a pixel based on the array index using numpy/gdal? I have a new set of extents for a raster and, based on these newer extents, want to determine the ...
1
vote
1answer
61 views

How to create GTiff from Tiff using GDAL

I have a tiff image and I would like to georeference it, and so create a Geotiff. To do that, I have to use GDAL. I use in particular, the executable store in the QGIS file (gdalwarp.exe, ...
1
vote
1answer
693 views

How to generate JPEG tiles using gdal2tiles.py

I am using gdal2tiles.py to generate tiles and packing it into MBTiles using mb-util. gdal2tiles by default generates PNG tiles. When I convert PNG tiles to JPEG large black area around tiles are left ...
2
votes
1answer
203 views

Extract data from netcdf using a shapefile polygon in python

I need to subset data from a NetCDF using a specific shapefile. The data are sea surface temperature and ocean color at 1/4 degree resolution. I have 4 polygons describing the US. Northeast ...
0
votes
1answer
48 views

Geoserver ecw plugin problem

I am trying to publish some ecw files with geoserver 2.7.2 on ubuntu 15.04 x64. I downloaded gdal plugin, gdal_data and native libraries for image-io 1.1.11 and set environment variables; GDAL_DATA, ...
6
votes
1answer
136 views

Extracting the pixelwise proportion of classes of a finer raster in a coarser one

I have two rasters, one has a 1 km resolution. The other has a spatial resolution of 30 m. This latter pixelwise is a classification. I want to overlay the 1 km raster on the classification raster and ...
0
votes
0answers
15 views

find pairs of near features with geopandas, fiona, gdal etc

I have a feature class of points, a few thousands I have now, but I expect to deal with hundreds of thousands of points. I would like find pair of points that is within a specific distance. The ...
1
vote
1answer
19 views

Image footprint in Landsat mosaic using GDAL

I need to produce a Landsat mosaic over part of Indonesia and an extra shape file containing the dates of each part of the mosaic. The steps I used are: Stack 3 bands to produce RGB image Project ...
-1
votes
1answer
23 views

ASC file to GeoTIFF

I would like to convert using Arcpy an ASC file to GeoTiff but I can't find nothing about it. I tried : arcpy.ProjectRaster_management() changing .tiff to .geotiff but it didn't work.
2
votes
1answer
62 views

Converting Lambert Conformal Conic (LCC) to UTM

I have two numpy meshgrids (X, Y both 2D) in Lambert Conformal conic with the following information: I'd like to directly (or indirectly if not possible) project them to UTM. I have both pyproj and ...
3
votes
1answer
505 views

Latest version of ogr2ogr fails with Layer Creation

I am able to successfully run ogr2ogr to add S57 data to PostgresDB. The version I use is GDAL 1.7.0b2, FWTools 2.4.7 and release date 2010/01/19 I upgraded my version to GDAL 1.9.1 and release date ...
14
votes
4answers
15k views

How can I install gdal 1.9.0 on my linux ubuntu server 10.04

apt-get install python-gdal installs version 1.6.0 and this is already installed. I would like to install the 1.9.0 version. Do I have to de-install 1.6.0 first and then compile manually or is there ...
2
votes
3answers
784 views

How to calculate Image boundary / footprint of satellite images using open source tools?

I'm referring to this question&discussion: How do I create a shapefile showing footprints of Rasters? I need to create the polygon outlines of several single layer raster images, not the ...
1
vote
1answer
35 views

ogr2ogr problem mapping shapefile to geojson ETRS89 / Poland CS92

I'm trying to convert official boundaries of administrative regions of Poland from official governmental site (codgik.gov.pl). I want to use those data because they are official ones, and set contains ...
1
vote
2answers
86 views

Suitable raster format for elevation data with sub-meter accuracy?

I've got some elevation data represented in ascii gridded xyz-files which I want to convert to a conventional raster format. What's a bit special is that the elevation values have sub-meter accuracy, ...
1
vote
1answer
37 views

gdal2tiles georeferenced image as overlay not working

Ive been looking for a solution to my problem for some time now. I'm using gdal2tiles-leaflet to create tiles for my Leaflet application. The images i use are GeoTiff's with the embedded projection ...
1
vote
1answer
22 views

Use gdalwarp to crop area with simple coordinate system

I'm new to GIS and trying to cut image with gdalwarp I need image projected like google map and I need to input coordinates of area of interest in latidute longitute system, like: from 45.3 10.2 to ...
16
votes
1answer
5k 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 ...
4
votes
1answer
256 views

How to change field name in raster file?

I´m having a bit of an understanding problem with a DEM TIF file. I have a DEM TIF file which is served with Geoserver. When I click on the layer the WMSGetFeatureInfo gives me the following output: ...
1
vote
0answers
117 views

How to import ERDAS ECW in QGIS 2.6 using GDAL 1.11? [duplicate]

I am using Lubuntu 14.10, QGIS 2.6, and Ubuntu GIS(unstable). QGIS was used to use GDAL 1.10, and can import ERDAS ecw format. But after GDAL upgraded to 1.11, QGIS couldn't import ecw images. How to ...
1
vote
0answers
14 views

How to Divide geo spatial Raster images into tiles and combine back into a complete image

I am reading a 2048 x 2048 image in the form of (204 x204) tiles and bilinear interpolating them to get (819x819) tiles. Then combining several 819x819 tiles to form a 8190x8190 image. Now as you can ...
1
vote
1answer
63 views

-9999 (no data value) becomes 0 when writing array to GDAL memory file

I am trying to reproject and resample (1000m ->> 30m) a raster image (shape 238X406) using gdal.ReprojectImage() function in Python. The input file has -9999 cells as no data value. The result is an ...
20
votes
3answers
7k views

How to call gdal_translate from Python code?

Is it possible to use some gdal API to call gdal_translate from Python code? I do not mean simply executing the gdal_translate.py from the file system, but rather call it somehow in code so I do not ...
0
votes
0answers
30 views

GeoPandas save to shapefile issue: Unable to open EPSG support file

When I try to use to_file() to write a GeoPandas GeoDataFrame to a shapefile I get the following error. df3.to_file('~/Desktop/testout.shp', driver = 'ESRI Shapefile'); ERROR:Fiona:OGR Error 4: ...
4
votes
2answers
252 views

Image processing using Python, GDAL and Scikit-Image

I'm new here and I'm sorry if I do something wrong. I am struggling with a processing and hopefully I will be able to solve here. I work with Remote Sensing applied to Forestry, especially working ...
2
votes
1answer
135 views

Is it possible to read raster files by block with rasterio?

I find rasterio much easier to use that plain ogr/gdal, so I'd prefer using it. However, I need to process large rasters, so reading rasters by block as explained in the Python GDAL course is ...
0
votes
0answers
43 views

gdal2tiles: how to make filenames of tiles to represent xyz coordinates?

I'm converting ECW file to tiles using gdal2tiles. Right now the tiles and folders generated by gdal2tiles are named as numbers starting from 0 successively. The folder structure looks like: How to ...
0
votes
2answers
28 views

Perform band calc in numpy array

I'd like to perform a band calc in some raster bands. What I am trying to do is import the raster (ENVI format) like array with gdal, then I do the calc but when I try to save the raster I get an ...
0
votes
0answers
16 views

GDAL data - change top left coordinate

I have GDAL raster dataset. I want to extract raw values from it. All is working correctly, but problem is with top left coorinates. If I extract data from GDALRasterBand, point[0,0] is at topLeft. ...
4
votes
1answer
2k 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
1answer
35 views

How to get elevation from dted only USGEO format using gdal apps

I need to get hieght by lat long position using gdal utilities
0
votes
1answer
150 views

Tools/Software for QA/QC of data

I currently have written a few custom scripts that run basic checks (null values, geocoding errors, fullness of metadata) using Arcpy(10.2 with almost all licenses) and gdal. I'm interested to know ...
2
votes
2answers
44 views

How can I convert points in EPSG:2956 to EPSG:4326 with pyproj?

I have some points in EPSG:4326 -110.28879750 64.46310187 405.49938863 -110.28879359 64.46310086 405.45967110 -110.28880175 64.46310446 405.93823604 -110.28880245 64.46310258 405.48060850 ...
0
votes
1answer
251 views

GDAL - projection from pixels to geo coordinates

I am using GDAL library in C++. On my dataset, I have obtained geo transform with GetGeoTransform to GT array (GCP is NULL). Than I use this values for this equation (Source): Xgeo = GT(0) + ...