Browse Library
Sign In Start Free Trial
Close icon Search icon
Account
Browse Library Advanced Search Sign In Start Free Trial

Add to playlist

Create a Playlist

Modal Close icon
You need to login to use this feature.
  • QGIS Python Programming Cookbook, Second Edition
  • Toc
  • feedback
QGIS Python Programming Cookbook, Second Edition

QGIS Python Programming Cookbook, Second Edition

By : Joel Lawhead
1.5 (2)
close
QGIS Python Programming Cookbook, Second Edition

QGIS Python Programming Cookbook, Second Edition

1.5 (2)
By: Joel Lawhead

Overview of this book

QGIS is a desktop geographic information system that facilitates data viewing, editing, and analysis. Paired with the most efficient scripting language—Python, we can write effective scripts that extend the core functionality of QGIS. Based on version QGIS 2.18, this book will teach you how to write Python code that works with spatial data to automate geoprocessing tasks in QGIS. It will cover topics such as querying and editing vector data and using raster data. You will also learn to create, edit, and optimize a vector layer for faster queries, reproject a vector layer, reduce the number of vertices in a vector layer without losing critical data, and convert a raster to a vector. Following this, you will work through recipes that will help you compose static maps, create heavily customized maps, and add specialized labels and annotations. As well as this, we’ll also share a few tips and tricks based on different aspects of QGIS.
Table of Contents (10 chapters)
close
close
Preface
Preface
What this book covers
What you need for this book
Who this book is for
Sections
Conventions
Reader feedback
Customer support
Free Chapter
1
1. Automating QGIS
chevron up
1. Automating QGIS
Introduction
Installing QGIS 2.18 for development
Using the QGIS Python console for interactive control
Automatically starting the Python console
Using the Python's Script Runner plugin
Setting up your QGIS IDE
Debugging QGIS Python scripts
Navigating the PyQGIS API
Creating a traditional QGIS plugin
Creating a Processing Toolbox plugin
Distributing a plugin
Building a standalone application
Storing and reading global preferences
Storing and reading project preferences
Accessing the script path from within your script
2
2. Querying Vector Data
2. Querying Vector Data
Introduction
Loading a vector layer from a file sample
Loading a vector layer from a geodatabase
Examining vector layer features
Examining vector layer attributes
Filtering a layer by geometry
Filtering a layer by attributes
Buffering a feature
Measuring the distance between two points
Measuring distance along a line
Calculating the area of a polygon
Creating a spatial index
Calculating the bearing of a line
Loading data from a spreadsheet
Accessing layer metadata
3
3. Editing Vector Data
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
Changing vector layer geometry
Changing a vector layer feature's attribute
Removing data from a vector layer
Deleting a vector layer feature's attribute
Reprojecting a vector layer
Converting a shapefile to KML or GeoJSON
Merging shapefiles
Splitting a shapefile
Generalizing a vector layer
Dissolving vector features
Performing a union on vector shapes
Rasterizing a vector layer
Exporting a layer to the GeoPackage format
4
4. Using Raster Data
4. Using Raster Data
Introduction
Loading a raster layer
Loading a NetCDF file
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 a 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
5. Creating Dynamic Maps
5. Creating Dynamic Maps
Introduction
Accessing the map canvas
Changing the map units
Iterating over layers
Symbolizing a vector layer
Setting a transparent layer fill
Using a filled marker symbol
Rendering a single band raster using a color ramp algorithm
Setting a feature's color using a column in a CSV file
Creating a complex vector layer symbol
Using icons as vector layer symbols
Using an outline for font markers
Using arrow symbols
Creating a graduated vector layer symbol renderer
Creating a categorized vector layer symbol
Using live layer effects
Creating inverted polygon shapeburst fills
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 2.5D renderer
Using XYZ tiled map services
Using Esri ArcGIS map services
Labeling a feature
Using null symbols
Change 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
6. Composing Static Maps
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 a map
Adding a North arrow to the map
Adding a logo to the map
Adding a vertical legend to the map
Adding a horizontal legend
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
7. Interacting with the User
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 a dock widget
Displaying a message in the status bar
Pushing messages to the message bar
Pushing widgets to the message bar
Creating tabs
Stepping the user through a wizard
Keeping dialogs on top
8
8. QGIS Workflows
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 DEM from LIDAR
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
Adjusting imprecise building footprints
Visualizing multitemporal data
9
9. Other Tips and Tricks
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-definition file
Using NULL values in PyQGIS
Using generators for layer queries
Using alpha values to show data density
Using the __geo_interface__ protocol
Getting the output file names from processing algorithms
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 points of interest in QGIS
Changing the QGIS web proxy
Running QGIS scripts as scheduled tasks
Visualizing data in 3D with WebGL
Visualizing data on a globe
Make a globe-like azimuthal orthographic projection
Animating a layer

Using the QGIS Python console for interactive control

The Python console allows you to interactively control QGIS; you can test out ideas or just do some quick automation. The console is the simplest way to use the API.

How to do it...

In the following steps, we'll open the console, create a vector layer in the memory, and display it on the map:

  1. Start QGIS.
  2. From the Plugins menu, select Python Console.
  3. The following code will create a point on the map canvas:
            layer =  QgsVectorLayer('Point?crs=epsg:4326','MyPoint',"memory") 
        pr = layer.dataProvider() 
        pt = QgsFeature() 
        point1 = QgsPoint(20,20) 
        pt.setGeometry(QgsGeometry.fromPoint(point1)) 
        pr.addFeatures([pt]) 
        layer.updateExtents() 
        QgsMapLayerRegistry.instance().addMapLayers([layer])

How it works...

This example uses a memory layer to avoid interacting with any data on either a disk or a network so as to keep things simple. Notice that when we declare the layer type, we add the parameter for the coordinate reference system (CRS) as EPSG:4326. Without this declaration, QGIS will prompt you to choose one. There are three parts, or levels, of abstraction even to create a single point on the map canvas, as shown here:

  1. First, create a layer that is of the type geometry. Next, reference a data provider to accept the data source.
  2. Then, create a generic feature object, followed by the point geometry.
  3. Next, we add the feature to the layer via the underlying data provider, then add the layer to the map canvas.

The layer type is memory, meaning that you can define the geometry and attributes in the code itself rather than an external data source. In this recipe, we just define the geometry and skip defining any attributes.

End of Section 4
Next Section
bookmark search playlist download
font-size

Change the font size

margin-width

Change margin width

day-mode

Change background colour

Close icon Search
Country selected
Close icon Your notes and bookmarks

Delete Bookmark

Modal Close icon
Are you sure you want to delete it?
Cancel
Yes, Delete