- 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
This recipe simply measures the area of a polygon.
For this recipe, we'll use a single-feature polygon shapefile, which you can download from https://geospatialpython.googlecode.com/files/Mississippi.zip
Unzip the shapefile and put it in a directory named qgis_data/ms within your root or home directory.
Perform the following steps to measure the area of a large polygon:
- First, import the QGIS constants library, as follows:
from qgis.core import QGis - Load the layer:
lyr = QgsVectorLayer("/qgis_data/ms/mississippi.shp", "Mississippi", "ogr") - Access the layer's features:
fts = lyr.getFeatures() - Get the boundary feature:
boundary = fts.next() - Create the measurement object instance:
d = QgsDistanceArea() - Pass the polygon list to the
measureArea()method:m = d.measurePolygon(boundary.geometry().asPolygon()[0]) - Convert the measurement from decimal degrees to miles:
d.convertMeasurement(m, QGis.Degrees, QGis.NauticalMiles, True) - Verify that your output looks similar to the following:
(42955.47889640281, 7)
PyQIS has no measureArea() method, but it has a
measurePolygon() method in the QgsDistanceArea object. The method accepts a list of points. In this case, when we convert the measurement output from decimal degrees to miles, we also specify True in the convertMeasurement() method so that QGIS knows that it is an area calculation. Note that when we get the boundary geometry as a polygon, we use an index of 0, suggesting that there is more than one polygon. A polygon geometry can have inner rings, which are specified as additional polygons. The outermost ring, in this case the only ring, is the first polygon.
-
No Comment
..................Content has been hidden....................
You can't read the all page of ebook, please click here login for view all page.