Stats zonales

100%

✓ Published0🌍 Public

Last edited Nov 13, 2012

Created on Nov 13, 2012

This example uses QGIS’s `QgsZonalStatistics` API to compute zonal statistics from a population raster (`popgrid.tif`) over vector polygons representing urban morphological zones (`UMZ_2000.shp`). It loads the vector layer via `QgsVectorLayer`, prompts the user for confirmation, then runs the calculation with a `ZS_` prefix for new attribute columns. The script prints the layer’s name, feature count, geometry type, and EPSG code before timing the operation, demonstrating a straightforward batch workflow for enriching polygon attributes with raster data.

AI-generated description

Similar vizzes

Loading thumbnail…

Gist f3b3c7d224f1701a9d7f

This Shiny application lets a user upload a CSV file of XY coordinates, which are then read with the built-in `read.table` function and converted into a spatial points object using `SpatialPointsDataFrame` from the `sp` package. The app renders a basic scatterplot of these points via `plot()` and offers a download button that exports the coordinates as an ESRI Shapefile using `writeOGR` from `rgdal`, zipping the resulting files for delivery. The interface is built with `bootstrapPage` and `fileInput`, and it validates that data has been selected before plotting.

RRCura

55% match

Loading thumbnail…

Gist 745c7c03b5bf6b6f59d9

This example splits a single ESRI Shapefile into multiple new Shapefiles, one for each unique value found in a specified attribute field. It demonstrates how to use the OGR library from GDAL to read the original dataset, extract all distinct attribute values, and then create filtered copies. The process first enumerates unique values, then iterates through them to create new Shapefile layers, applies an SQL attribute filter to exclude features not matching the current value, and deletes the remaining unwanted features.

RRCura

54% match

Loading thumbnail…

Gist 5053ebcce793d11bc765

This example maps all Dutch municipalities in 2014, highlighting those with fewer than 20,000 residents in red on top of the full country outline. It uses the `rgdal` and `rgeos` libraries to read shapefile data from a ZIP downloaded from twiav.nl, then computes population density via `gArea()` and identifies the least dense municipality. The visualization shows a static choropleth-style base map with an overlaid selection, while also printing tabular summaries like counts per safety region and the names of the least dense area. It does not animate or use interactive frameworks, relying solely on base R plotting.

FFrieseWoudloper

54% match

Loading thumbnail…

Funcionalidades básicas de QGIS (III)

This example demonstrates how to create a point layer from earthquake data and work with raster files in QGIS. It shows the process of downloading earthquake records from the USGS catalog as a CSV file, importing this delimited text layer into QGIS, and specifying the coordinate reference system (WGS84, EPSG:4326) for geographic coordinates. The tutorial then applies a size assistant to scale point symbols based on earthquake magnitude, filters the data using expressions like "mag" > 5, and adds supplementary layers such as Natural Earth administrative boundaries and tectonic plate shapefiles with categorized styling and labels. The workflow covers loading the CSV via QGIS’s delimited-text layer dialog, configuring X/Y coordinate fields, and managing layer properties for visual customization.

LLuisSevillano

53% match

Loading thumbnail…

Raster & Vector III

This example shows how vector geographic features (U.S. counties) can be overlaid on raster map tiles from OpenStreetMap, with both layers staying aligned during pan and zoom interactions. The vector layer is rendered as a single SVG path using `d3.geoPath` and a `d3.geoMercator` projection, while raster tiles are fetched and positioned via `d3.tile` and `d3-zoom`. Instead of reprojecting, the code updates the SVG transform and stroke-width for efficiency, using `d3.zoomIdentity` and the `on("zoom")` event.

Mmbostock

53% match

Loading thumbnail…

Gist 9587685

This example maps user-uploaded shapefiles by plotting their polygon boundaries in a main panel. It uses R’s `shiny` framework with `server.R` and `ui.R` files, relying on the `maptools` library’s `readShapePoly` function to load the spatial data. A reactive `uploadShpfile` expression handles file renaming and directory switching for multiple uploaded components, while `renderPlot` and `plotOutput` display the resulting map. The interface provides a sidebar `fileInput` for selecting shapefiles and a central plot area for visualization.

RRCura

52% match

Loading thumbnail…

Scoring buildings based on proximity to neighborhood features

This example visualizes building scores in Portland based on proximity to nearby amenities, assigning higher ranks to structures near supermarkets, rail stops, and parks. It shows a ranked map of buildings, with scores aggregated from overlapping amenity zones. The workflow uses PostGIS spatial SQL functions like `st_expand` and `st_intersects` to generate buffered zones, then sums scores per building via a `union` query and stores results in a new table. The rendering likely relies on a mapping library such as Leaflet or Mapbox, though the code excerpt focuses on the spatial data preparation rather than the front-end visualization.

CCaged

52% match

Loading thumbnail…

Funcionalidades básicas de QGIS (II)

This example shows how QGIS vector geoprocessing tools create new spatial datasets, such as intersecting layers, dissolving shared boundaries, and generating fixed-distance buffers around features. The code demonstrates a Vue workflow that filters municipal polygons, applies a dissolve to derive an island contour, and uses a buffer to define a protected perimeter. It relies on QGIS’s native vector analysis APIs for operations like intersection, asymmetric difference, union, and clip, and references Open Street Map extracts from Geofabrik alongside Red Natura 2000 data.

LLuisSevillano

52% match