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
1answer
19 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" ...
0
votes
1answer
17 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
34 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
42 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 ...
0
votes
0answers
14 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
18 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 ...
0
votes
1answer
8 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
1answer
60 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, ...
0
votes
0answers
22 views

Unexpected No values when reading grib2 file with gdal [closed]

I am trying to read this .grib2 file using gdal. Here's the code I use in Python: from osgeo import gdal from osgeo.gdalconst import GA_ReadOnly fname = 'ten_meter_20150816_00z_01.grib' dataset = ...
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, ...
1
vote
1answer
35 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 ...
2
votes
1answer
60 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 ...
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
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 ...
1
vote
1answer
34 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 ...
0
votes
0answers
26 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: ...
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, ...
2
votes
1answer
134 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 ...
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 ...
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. ...
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 ...
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 ...
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 ...
1
vote
0answers
21 views

Python-gdal create geotiff from array with colormapping

I am creating grayscale, single channel geotiff files from numpy arrays with the below module, but I would like to apply a fairly standard colormap to the data. I've tried three channels of R,B and G ...
2
votes
1answer
38 views

How can I export a shapefile from PostGIS with a GeoJSON query?

I tried to export a shapefile from my PostGIS installation using pgsql2shp as follow pgsql2shp -f tiles.shp -h <host> -u <user> -P <pass> <database> "SELECT id, the_geom FROM ...
0
votes
0answers
28 views

How to transform ENVI Hdr to PNG using GDAL?

I am trying to transform a ENVI Hdr image into a PNG (or other RGB image format) in order to visualize this image. I am selecting the three default bands (55, 41, and 12). I am using GDAL by using ...
0
votes
0answers
26 views

Aggregate/merge data upon insert into PostGIS (via gdal)

I am working on a project where I take a large raster file (1.5GB+), put it through a number of image-processing operations, and finally convert it to vector data. Being that the raster is large, I am ...
2
votes
1answer
49 views

python crashes with returned gdal band object from a function

I wrote these two functions in a module when I want to read the band in NetCDF or HDF files and return the band outside the function. Each HDF and NetCDF has many subdatasets. Each subdataset in a ...
1
vote
0answers
32 views

gdal_merge shifts rasters after processing, layer stack is slightly off

I have experienced that gdal shifts rasters. The outputs of gdal_warp are often skewed, have different pixel sizes and so forth. A known error which has also been discussed here. But now, I had to ...
0
votes
0answers
24 views

Python crashes when I try to import GDAL after re-installing it

A while back I went through the installation process for GDAL by installing the core components and then the Python/GDAL bindings. Everything was fine when I would type in from osgeo import gdal. ...
3
votes
1answer
45 views

How to re-compile ogrinfo against gdal used by QGIS?

I just solved a problem with missing libraries blocking use of ogrinfo and ogr2ogr, only to find that these utilities are garbled in one more way: whenever I run any of these, I get following error ...
0
votes
2answers
33 views

How do I run gdalwarp on a number of csv of coordinates?

I have a CSV of Ground Control Points in NAD83 UTM 12N and I need to convert them to WSG84. Is it possible to do this with gdalwarp or something similar?
1
vote
1answer
56 views

Erase shapfile (polygon) area from GeoTIFF

I've been reading several posts to do so and one option would be to use gdalwarp -cutline, but I would need to keep what it is NOT inside the polygon area. GDAL or any command-line tool does that? ...
0
votes
1answer
62 views

Problem doing a coordinate transform with GDAL python

I'm trying to do a concordant transform from the format of the shapefile (3309) to google maps(4326). This code works on one Centos 6.6 box with the SCL python 2.7 and a custom gdal 2.0 build. ...
-1
votes
1answer
19 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.
1
vote
0answers
31 views

Serving GeoTIFFs of different resolutions as single layer

I'm trying to serve a 30 arc second DEM with global coverage combined with high resolution 1/3 arc second DEM for a very small area. I believed at first I could combine them into a single file, ...
0
votes
0answers
27 views

WCS elevation with Cesium + Geoserver

I'm trying to develop an offline demonstrator showing hi-resolution bathymetry data using Cesium 1.11 and Geoserver. I've merged a number of 1/3 arc second areas around Puerto-Rico from NOAA ...
2
votes
2answers
11 views

Fusion table connexion issue

I'm trying to import some data from a SHP to GFT (Google Fusion Table). I get a 401 error. I'm using an Ubuntu 14.04 LTS Linux 64 bits with GDAL 1.11.0, released 2014/04/16. GFT is available (tested ...
0
votes
1answer
27 views

OSGeo4W GDAL Polygonize Batch Convert?

I'm trying to use GDAL polygonize to run a batch convert but I cannot for the life of my figure out a way to get it to read a text file as input that contains all of the commands. Do any of you know ...
0
votes
2answers
26 views

I cannot find osr after installing OSGeo4W package

I wanted to be able to import OSR from GDAL on my Python 2.7 (Anaconda) Windows 7 64-bit. I installed OSGeo4W after reading this post and eventually I'd like to import osr from osgeo import osr ...
1
vote
1answer
63 views

How do I convert WGS84 to Lat, Long to Google Maps in R?

I have this coordinate which I believe it's in WGS84 format x=1005740.867110, y=243957.356623 I want to plot that in Google Maps which I need to convert that to Lon, Lat. I've found this script ...
3
votes
2answers
84 views

How to have GDAL print layers of GeoPDF AND say which are raster vs vector

My Objective: I would like to use GDAL to convert a GeoPDF. I want the vector layers as shp files and the raster layers as tif files. I want to do this in a programmatic way. Edit: In reality, I ...
0
votes
0answers
16 views

Using GDAL Convert .SHP to GEO TIFF in Command Line Python [duplicate]

I have dozen Shape files (.SHP) with Georeferenced. I need to convert to GeoTiff or other format that GDAL primarly gdal2tiles.py to convert to MBTiles format. All I want is convert the .shp files ...
0
votes
0answers
23 views

Installing GDAL Python bindings on Windows [duplicate]

This has been addressed before (How to install GDAL with Python on windows?), but I wanted to bring the question up again, to see if anybody has come up with a better method in the last 5 years, and ...
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
0answers
26 views

What are the units in GMTED2010 elevation data?

I have downloaded the GMTED2010 elevation data from http://topotools.cr.usgs.gov/gmted_viewer/viewer.htm I am merging all the files/tiles (to create a global dataset) with: gdal_merge.py There is ...
1
vote
0answers
22 views

Non-Square Pixel on georeferencing of images using GDAL

I have a image in .jpg format. I have some GCP and also know the pixels values of these GCp in that image. Now i trying to convert this image into geotiff. for that I have coded the following lines. ...
2
votes
0answers
38 views

Multiple layer intersection

I am trying to find a way with qgis/python to get the intersection area of more than 2 layers. All the tools in Qgis take just the 2 input layers. ='( For example in the image below, I have 4 ...