- QGIS Python Programming Cookbook
- Table of Contents
- QGIS Python Programming Cookbook
- Credits
- About the Author
- About the Reviewers
- www.PacktPub.com
- Preface
- 1. Automating QGIS
- Introduction
- Installing QGIS for development
- Using the QGIS Python console for interactive control
- Using the Python ScriptRunner plugin
- Setting up your QGIS IDE
- Debugging QGIS Python scripts
- Navigating the PyQGIS API
- Creating a QGIS plugin
- Distributing a plugin
- Creating a standalone application
- Storing and reading global preferences
- Storing and reading project preferences
- Accessing the script path from within your script
- 2. Querying Vector Data
- Introduction
- Loading a vector layer from a file sample
- Loading a vector layer from a spatial database
- Examining vector layer features
- Examining vector layer attributes
- Filtering a layer by geometry
- Filtering a layer by attributes
- Buffering a feature intermediate
- Measuring the distance between two points
- Measuring the distance along a line sample
- Calculating the area of a polygon
- Creating a spatial index
- Calculating the bearing of a line
- Loading data from a spreadsheet
- 3. Editing Vector Data
- Introduction
- Creating a vector layer in memory
- Adding a point feature to a vector layer
- Adding a line feature to a vector layer
- Adding a polygon feature to a vector layer
- Adding a set of attributes to a vector layer
- Adding a field to a vector layer
- Joining a shapefile attribute table to a CSV file
- Moving vector layer geometry
- Changing a vector layer feature's attribute
- Deleting a vector layer feature
- Deleting a vector layer attribute
- Reprojecting a vector layer
- Converting a shapefile to KML
- Merging shapefiles
- Splitting a shapefile
- Generalizing a vector layer
- Dissolving vector shapes
- Performing a union on vector shapes
- Rasterizing a vector layer
- 4. Using Raster Data
- Introduction
- Loading a raster layer
- Getting the cell size of a raster layer
- Obtaining the width and height of a raster
- Counting raster bands
- Swapping raster bands
- Querying the value of a raster at a specified point
- Reprojecting a raster
- Creating an elevation hillshade
- Creating vector contours from elevation data
- Sampling a raster dataset using a regular grid
- Adding elevation data to line vertices using a digital elevation model
- Creating a common extent for rasters
- Resampling raster resolution
- Counting the unique values in a raster
- Mosaicing rasters
- Converting a TIFF image to a JPEG image
- Creating pyramids for a raster
- Converting a pixel location to a map coordinate
- Converting a map coordinate to a pixel location
- Creating a KML image overlay for a raster
- Classifying a raster
- Converting a raster to a vector
- Georeferencing a raster from control points
- Clipping a raster using a shapefile
- 5. Creating Dynamic Maps
- Introduction
- Accessing the map canvas
- Changing the map units
- Iterating over layers
- Symbolizing a vector layer
- Rendering a single band raster using a color ramp algorithm
- Creating a complex vector layer symbol
- Using icons as vector layer symbols
- Creating a graduated vector layer symbol renderer
- Creating a categorized vector layer symbol
- Creating a map bookmark
- Navigating to a map bookmark
- Setting scale-based visibility for a layer
- Using SVG for layer symbols
- Using pie charts for symbols
- Using the OpenStreetMap service
- Using the Bing aerial image service
- Adding real-time weather data from OpenWeatherMap
- Labeling features
- Changing map layer transparency
- Adding standard map tools to the canvas
- Using a map tool to draw points on the canvas
- Using a map tool to draw polygons or lines on the canvas
- Building a custom selection tool
- Creating a mouse coordinate tracking tool
- 6. Composing Static Maps
- Introduction
- Creating the simplest map renderer
- Using the map composer
- Adding labels to a map for printing
- Adding a scale bar to the map
- Adding a north arrow to the map
- Adding a logo to the map
- Adding a legend to the map
- Adding a custom shape to the map
- Adding a grid to the map
- Adding a table to the map
- Adding a world file to a map image
- Saving a map to a project
- Loading a map from a project
- 7. Interacting with the User
- Introduction
- Using log files
- Creating a simple message dialog
- Creating a warning dialog
- Creating an error dialog
- Displaying a progress bar
- Creating a simple text input dialog
- Creating a file input dialog
- Creating a combobox
- Creating radio buttons
- Creating checkboxes
- Creating tabs
- Stepping the user through a wizard
- Keeping dialogs on top
- 8. QGIS Workflows
- Introduction
- Creating an NDVI
- Geocoding addresses
- Creating raster footprints
- Performing network analysis
- Routing along streets
- Tracking a GPS
- Creating a mapbook
- Finding the least cost path
- Performing nearest neighbor analysis
- Creating a heat map
- Creating a dot density map
- Collecting field data
- Computing road slope using elevation data
- Geolocating photos on the map
- Image change detection
- 9. Other Tips and Tricks
- Introduction
- Creating tiles from a QGIS map
- Adding a layer to geojson.io
- Rendering map layers based on rules
- Creating a layer style file
- Using NULL values in PyQGIS
- Using generators for layer queries
- Using alpha values to show data density
- Using the __geo_interface__ protocol
- Generating points along a line
- Using expression-based labels
- Creating dynamic forms in QGIS
- Calculating length for all selected lines
- Using a different status bar CRS than the map
- Creating HTML labels in QGIS
- Using OpenStreetMap's points of interest in QGIS
- Visualizing data in 3D with WebGL
- Visualizing data on a globe
- Index
Expressions are a kind of mini-programming language or SQL-like language found throughout different QGIS functions to select features. One important use of expressions is to control labels. Maps easily become cluttered if you label every single feature. Expressions make it easy to limit labels to important features. You can filter labels using expressions from within Python, as we will do in this recipe.
You will need to download the zipped shapefile and decompress it to a directory named ms in your qgis_data directory from the following:
https://geospatialpython.googlecode.com/files/MS_UrbanAnC10.zip
We'll use the QGIS PAL labeling engine to filter labels based on a field name. After loading the layer, we'll create our PAL settings and write them to the layer. Finally, we'll add the layer to the map. To do this, we need to perform the following steps:
- First, we'll set up the path to our shapefile:
pth = "/Users/joellawhead/qgis_data/ms/MS_UrbanAnC10.shp" - Next, we'll load our layer:
lyr = QgsVectorLayer(pth, "Urban Areas", "ogr") - Now, we create a labeling object and read the layer's current labeling settings:
palyr = QgsPalLayerSettings() palyr.readFromLayer(lyr)
- We create our expression to only label the features whose population field is greater than 50,000:
palyr.fieldName = 'CASE WHEN "POP" > 50000 THEN NAME10 END' - Then, we enable these settings:
palyr.enabled = True - Finally, we apply the labeling filter to the layer and add it to the map:
palyr.writeToLayer(lyr) QgsMapLayerRegistry.instance().addMapLayer(lyr)
-
No Comment
..................Content has been hidden....................
You can't read the all page of ebook, please click here login for view all page.