Search Results

I read your question before being edited and reopened and it could be preferable to calculate medians for each particular raster instead to use QgsZonalStatistics. So, you can re-write your code as fo …

I tried out the next code with a polygon shapefile (as active layer) where, first, it is produced a "negative" buffer (-10000 m) as memory layer. Afterward, this memory layer is used as difference lay …

For this situation: complete code for getting only standard deviation is: from qgis.analysis import QgsZonalStatistics mapcanvas = iface.mapCanvas() layers = mapcanvas.layers() provider = laye …

Your code works perfectly but, you cannot visualize changes immediately because it's necessary to add polygonLayer to registry. However, it can be loaded manually. A simplified version of your code ( …

I can handle partially overlapping pixels adapting code (for PyQGIS in QGIS 3) in the answer to this question: Robust Zonal stats in R and QGIS. import processing registry = QgsProject.instance() …

I found out a solution mapping your original SE_FireHistory Image Collection as list. As your Fedlands_SE Feature Collection has 316 elements, most of them outside the area of interest, for greatly re …

A Robust Zonal stats in QGIS can be implemented with PyQGIS. Following code was run with your layers by using a filter to select only ID_sub.pat between 19 and 27 (showed in your image). import proce …

For some reason your code is not producing a QgsZonalStatistics object. Print 'zoneStats' to corroborate this. It should be observed some similar to: <qgis._analysis.QgsZonalStatistics object at 0x9f …

For computing zonal statistics in PyQGIS with overlapping features you can use 'unary_union' and 'polygonize' methods from shapely python module. For trying out my approach I used next shapefile (with …

The rasterstats python module is based in rasterio, however, it has crop = True option activated by default (when I use raster_out = True for getting all unmasked array values) and I couldn't find an …

QGIS calculates NULL because both layers have not the same projection. Layer named isos has EPSG:4326 and BY-AT raster has EPSG:54009 (World_Mollweide-Projected). I got your same result by using these …

I followed your code line by line and you don't need to create an image that contains one band with the pixel's area and a second band with the class information. You only need to retrieve both landCo …

As it was pointed out by Steven Kay, there's some sample code on the GDAL/OGR cookbook about 'zonal statistics', so it's some extensive when it is compared with the mdsumner's R code. I know that R po …

Zonal Statistics, which include partial pixels, can be performed by using PyQGIS. I tried out PyQGIS code with a small raster (20 rows x 20 columns, pixel size 30 m) and a shapefile with two features …