Sign In Start Free Trial

Book Overview & Buying
Table Of Contents
Feedback & Rating

Geospatial Development By Example with Python

By : Pablo Carreira

5 (4)

Geospatial Development By Example with Python

5 (4)

By: Pablo Carreira

Overview of this book

From Python programming good practices to the advanced use of analysis packages, this book teaches you how to write applications that will perform complex geoprocessing tasks that can be replicated and reused. Much more than simple scripts, you will write functions to import data, create Python classes that represent your features, and learn how to combine and filter them. With pluggable mechanisms, you will learn how to visualize data and the results of analysis in beautiful maps that can be batch-generated and embedded into documents or web pages. Finally, you will learn how to consume and process an enormous amount of data very efficiently by using advanced tools and modern computers’ parallel processing capabilities.

Preface

Preface

What this book covers

What you need for this book

Who this book is for

Conventions

Reader feedback

Customer support

Free Chapter

1. Preparing the Work Environment

1. Preparing the Work Environment

Installing Python

Python packages and package manager

Installing packages and required software

Installing NumPy

Installing GDAL and OGR

Installing Mapnik

Installing Shapely

Installing other packages directly from pip

Installing an IDE

Creating the book project

Programming and running your first example

Transforming the coordinate system and calculating the area of all countries

Sort the countries by area size

Summary

2. The Geocaching App

2. The Geocaching App

Building the basic application structure

Downloading geocaching data

Opening the file and getting its contents

Combining functions into an application

Setting your current location

Finding the closest point

Summary

3. Combining Multiple Data Sources

3. Combining Multiple Data Sources

Representing geographic data

Making data homogeneous

Importing geocaching data

Integrating new functionality into the application

Summary

4. Improving the App Search Capabilities

4. Improving the App Search Capabilities

Working with polygons

Using Shapely to handle geometries

Importing polygons

Getting the attributes' values

Importing lines

Converting the spatial reference system and units

Geometry relationships

Filtering by attributes and relations

Filtering by multiple attributes

Integrating with the app

Summary

5. Making Maps

5. Making Maps

Knowing Mapnik

Creating utility functions to generate maps

Styling maps

Using Python objects as a source of data

Exporting geo objects

Creating the Map Maker app

Summary

6. Working with Remote Sensing Images

6. Working with Remote Sensing Images

Understanding how images are represented

Processing remote sensing images and data

Building an image processing pipeline

Summary

7. Extract Information from Raster Data

7. Extract Information from Raster Data

Getting the basic statistics

Creating color classified images

Blending images

Showing statistics with colors

Summary

8. Data Miner App

8. Data Miner App

Measuring execution time

Code profiling

Storing information on a database

Importing massive amount of data

Searching for data and crossing information

Summary

9. Processing Big Images

9. Processing Big Images

Working with satellite images

Memory and images

Processing images in chunks

Creating image compositions

Summary

10. Parallel Processing

10. Parallel Processing

Multiprocessing basics

Block iteration

Improving the image resolution

Summary

Index

Index

Customer Reviews

5 (4)

5 star

100%

4 star

0

3 star

0

2 star

0

1 star

0

Blending images

We can make our results even more visually appealing and informative if we can combine the colorized image with the shaded relief image. Again, since we are dealing with arrays, we may deduce that this kind of composition can be achieved by performing an arithmetic operation between the two arrays.

In image processing, this is called alpha blending. Basically, a transparency is applied to both of the images and then they are blended into a new one. In the next steps, we are going to create a function that performs this operation:

First, to avoid generating the shaded relief multiple times, let's save it on the disk and edit the if __name__ == '__main__': block of the raster_data.py file:

if __name__ == '__main__':
    raster_data = RasterData('output/dem.tif')
    raster_data.adjust_values().create_hillshade(
        10, 60).write_image('output/shaded.png')

Run the code and check whether the image was correctly written on the disk.
Now, add the alpha_blend method to the RasterData...

Search

Your notes and bookmarks