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

2
votes
2answers
51 views

Memory and Processing Efficient Approach to Read a Timeseries in Python

I am trying to store a timeseries of high resolution raster images containing NDVI values into a dictionary for calculations. I am curious as to which approach may the most efficient when loading ...
2
votes
2answers
83 views

Invert gdaldem aspect

I am trying to identify lakes and other water bodies in a DEM file. The source file is located here. I am using gdaldem aspect to identify contiguous pixels with zero slope change. gdaldem aspect ...
2
votes
2answers
113 views

OSM linestring shown as point

I have loaded the following file http://openstreetmapdata.com/data/coastlines to postgreSQL. In visual studio with gdal/ogr and geos enabled: SELECT geom FROM waterpoly WHERE geom IS NOT NULL LIMIT ...
2
votes
2answers
153 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 ...
1
vote
2answers
97 views

GDAL Conversion of US Topo GeoPDFs

I've used GDAL_translate with "--config GDAL_PDF_LAYERS" option to convert US Topo GeoPDFs to GeoTIFF and it works great for creating collarless topos based on the new format. However, I can't get it ...
1
vote
2answers
69 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. ...
1
vote
2answers
18 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
239 views

GDAL/proj4 equirectangular projection coordinate system units assumed by OGRSpatialReference and GDALDataset::GetGeoTransform()?

How do I configure my OGRSpatialReference object to specify an equirectangular projection system http://en.wikipedia.org/wiki/Equirectangular_projection where the affine geo-transform matrix (provided ...
1
vote
2answers
104 views

display GeoReferencing image through MapServer

I have a tiff image and got the world file which contains the needed information for GeoReferenceing, how to start add this image on Mapserver ? if I use C# and use GMAP API to call the google server ...
1
vote
2answers
240 views

Installing python GDAL with pip in local (user) space on RHEL6 server: failure due to missing headers

There are similar themed questions for other platforms, but the solutions out there do not work for me. (e.g. ...
1
vote
1answer
122 views

combine raster with another raster mask file - gdal

I have 2 rasters in GeoTiff format. One has all of the land and water and the other is a mask file that shows either 0 for water or 1 for land. I basically want my output to just have the land values ...
1
vote
1answer
98 views

gdal2tiles doesn't use resampling method correct?

I'm using gdal2tiles to create tiles from MODIS satellite imagery (250 m resolution) and I have noticed that I get the same bad result if I use -r near or -r bilinear (or lanczos for that matter). See ...
1
vote
1answer
52 views

OGRFeature::DestroyFeature fails

I am using OGR to read features from a shape file; when the feature is destroyed I am getting an error 0xC0000005: Access violation reading location 0xffffffffffffffff. This secton of the code looks ...
1
vote
1answer
87 views

SetColorInterpretation() for a GeoTiff using GDAL not working

I have a single band image which I would like to copy to a 3 bands of a new image and set the ColorInterpretation to red, green and blue respectively. I have tried two approaches, one using the GTiff ...
1
vote
1answer
41 views

Can I create a virtual mosaic using other virtual mosaics?

I am following this example - http://linfiniti.com/2009/09/image-mosaicking-with-gdal/ I have 6TB of aerial imagery in 6 folders 1TB each and each folder has many subdirectories including the ...
1
vote
1answer
400 views

GDAL loads in Python, but not IPython (Linux)

I am trying my hardest to utilize the work of the good folks at the GeoDa center, but I cannot seem to get all of the components to load. GDAL, in particular, has been giving me some difficulty. The ...
1
vote
1answer
268 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
322 views

How can I add all the 4 corners as GCP in an image that I want to reproject?

I am working on a Modis - oceancolor image and I want to reproject it to wgs84. I have the Upper left/right and Low left/right values from the metadata of the image and I want to add them all in the ...
1
vote
1answer
235 views

Is there any support for Rasterlite in the Android QGIS?

Is there any support for Rasterlite databases in the Android port of QGIS? I have tried adding a raster layer as I would on desktop QGIS but there are fewer GDAL file formats listed. Is this ...
0
votes
1answer
18 views

how to calculate the average temperature of each band of a .grd file?

I am having a GRD file in which the temperature at different longitude and latitude are given in a 32x35 array for one day(band),i am having 365 such bands and want to calculate the average of all ...
0
votes
1answer
34 views

GDal problem “no suitable image found” when I run “ogr2ogr” command

I am trying to run "ogr2ogr" command macbook$ sudo ogr2ogr -f "PostgreSQL" PG:"dbname=test user=postgres" /Users/macbook/planning.gdb But I am getting this error: dyld: Library not loaded: ...
0
votes
1answer
38 views

Is it possible to get all projections form OGR?

I'm working on a project where the user can give my program so GIS data which may or may not have which projection it is using. So is it possible to get a list of all the projections from OGR or ...
0
votes
1answer
24 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 ...
0
votes
1answer
74 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 ...
0
votes
1answer
22 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 ...
0
votes
1answer
27 views

How to import r.in.gdal from WCS (only region extent)

Can I import WCS-Data only for the current region by r.in.gdal? I tried to import raster data from my WCS into GRASS GIS, using r.in.gdal as proposed in ...
0
votes
1answer
47 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) + ...
0
votes
1answer
76 views

Can I edit pixels values of an image with GDAL?

There is a method with GDAL (used in C/C++) to edit pixel values of a georeferenced image, maintaining the original file? Or I have to create a new dataset? My goal is to maintain ALL originals ...
0
votes
1answer
69 views

Need assistance setting center coordinate with GDALWARP

I'm having some issues converting my project into an MBTiles format. I'm trying to use gdalwarp to convert my main file to geotiff format, but I can't seem to get the center coordinate to how it is ...
0
votes
1answer
48 views

USGS FLT Hdr file interpretation

I have a legacy program that attempts to read USGS 1/3 arcsec data. When I compare results to GDAL, I find that the lat/lon to pixel conversion is off by 1 line. The reason (I think) is because the ...
0
votes
1answer
147 views

Polygonize.py: AttributeError: 'NoneType' object has no attribute 'CreateDataSource'

Given a topographic raster tif of an european country where elevations are simplyfied to either 0, either 200. How to create a polygon ? Comment: With the following polygonize.py command: ...
0
votes
1answer
82 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 ...
0
votes
1answer
426 views

How to clip raster by multiple polygons in multi rasters?

Using QGIS, I want to clip a raster by multiple polygons contained in a single shapefile. I read that gdalwarp can clip my raster by multiple polygons BUT the output is a single merged clipped raster. ...
-1
votes
1answer
147 views

How to Geocode use GeoPy API combined with GDAL Python bindings?

My question is related to the 'How to Geocode 300,000 addresses on the fly?' topic that has been asked by Mehul Nahar. I'm interested in running the processes explained by capooti in response to mehul ...