Reprojecting a raster
A core requirement for all geospatial analysis is the ability to change the map projection of data in order to allow different layers to be displayed on the same map. Reprojection can be challenging, but QGIS makes it a piece of cake. Starting with this recipe, we will begin using the powerful QGIS Processing Toolbox. The Processing Toolbox wraps over 600 algorithms into a highly consistent API available to Python and also as interactive tools. This toolbox was originally a third-party plugin named SEXTANTE but is now a standard plugin distributed with QGIS.
Getting ready
As with many recipes in this chapter, we will use the SatImage raster available at https://github.com/GeospatialPython/Learn/raw/master/SatImage.zip.
Place this raster in your /qgis_data/rasters directory.
How to do it...
In this recipe, we will use the gdalwarp algorithm of the processing module to reproject our image from EPSG 4326 to 3722. To do this, we need to perform the following steps:
Start QGIS...
