- 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
Map bookmarks store important locations on a map, so you can quickly find them later. You can programmatically navigate to bookmarks using the Python sqlite3 library in order to access the bookmarks database table in the QGIS user database and then use the PyQGIS canvas API.
We will use a census tract layer to test out the bookmark navigation. You can download the zipped shapefile from https://geospatialpython.googlecode.com/files/GIS_CensusTract.zip.
Manually load this layer into QGIS after extracting it from the ZIP file. Also, make sure that you have completed the previous recipe, Creating a map bookmark. You will need a bookmark named BSL for an area of interest in this shapefile.
We will retrieve a bookmark from the QGIS user database and then set the map's extents to this bookmark. To do this, perform the following steps:
- First, import the Python
sqlite3library:import sqlite3 - Next, get the location of the user database from the QGIS data:
dbPath = QgsApplication.qgisUserDbFilePath() - Now, we connect to the database:
db = sqlite3.connect(dbPath) - Then, we need a database cursor to run queries:
cursor = db.cursor() - Now, we can get the bookmark information for the bookmark named BSL:
cursor.execute("""SELECT * FROM tbl_bookmarks WHERE name='BSL'""") - Now, we'll get the complete results from the query:
row = cursor.fetchone() - Then, we split the values of the result into multiple variables:
id,mark_name,project,xmin,ymin,xmax,ymax,srid = row - Now, we can use the bookmark to create a QGIS extent rectangle:
rect = QgsRectangle(xmin, ymin, xmax, ymax) - Next, we reference the map canvas:
canvas = qgis.utils.iface.mapCanvas() - Finally, we set the extent of the canvas to the rectangle and then refresh the canvas:
canvas.setExtent(rect) canvas.refresh()
-
No Comment
..................Content has been hidden....................
You can't read the all page of ebook, please click here login for view all page.