EDIT: I guess I’d use qGIS python bindings: http://www.qgis.org/wiki/Python_Bindings That’s the easiest way I can think of. I remember hand rolling something before, but it’s ugly. qGIS would be easier, even if you had to make a separate Windows installation (to get python to work with it) then set up an XML-RPC server to run … Read more
scikit-image has implemented a working version of downsampling here, although they shy away from calling it downsampling for it not being a downsampling in terms of DSP, if I understand correctly: http://scikit-image.org/docs/dev/api/skimage.measure.html#skimage.measure.block_reduce but it works very well, and it is the only downsampler that I found in Python that can deal with np.nan in the … Read more
Check that you installed GDAL using this command gdal-config –version Then run this commands: pip download=”some_path” GDAL cd some_path tar -xvzf GDAL-<version>.tar.gz cd GDAL-<version> python setup.py build_ext –include-dirs=/usr/include/gdal/ python setup.py install
After two dozens of comments to understand the situation, it was found that the libhdf5.so.7 was actually a symlink (with several levels of indirection) to a file that was not shared between the queued processes and the interactive processes. This means even though the symlink itself lies on a shared filesystem, the contents of the … Read more
What worked for me is this: https://gis.stackexchange.com/a/193828/66527 Below, I copy that answer: You can download GDAL 2.1 for Windows from GIS Internals. There is an installer and a portable version that doesn’t require installation. GDAL 2.1 is available for Ubuntu 16.04 from the UbuntuGIS-Stable PPA sudo add-apt-repository -y ppa:ubuntugis/ppa sudo apt update sudo apt upgrade … Read more
A[A==NDV]=numpy.nan A==NDV will produce a boolean array that can be used as an index for A
To get the coordinates of the corners of your geotiff do the following: from osgeo import gdal ds = gdal.Open(‘path/to/file’) width = ds.RasterXSize height = ds.RasterYSize gt = ds.GetGeoTransform() minx = gt[0] miny = gt[3] + width*gt[4] + height*gt[5] maxx = gt[0] + width*gt[1] + height*gt[2] maxy = gt[3] However, these might not be in … Read more
In Ubuntu 18.04 I fixed this by sudo apt install libgdal-dev Hope someone find this helpful. Some above answers seems to be outdated and lengthy. In earlier versions (which had apt-get) sudo apt-get install libgdal-dev