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

Getting the attributes' values

Let's explore the attributes of the world borders to find out why we were unable to get the names.

Edit the if __name__ == '__main__': block:

if __name__ == '__main__':
    world = BoundaryCollection("../data/world_borders_simple.shp")
    print(world.data[0].attributes.keys())

Run the code and look at the output:
```
File imported: ../data/world_borders_simple.shp
['SUBREGION', 'POP2005', 'REGION', 'ISO3', 'ISO2', 'FIPS', 'UN', 'NAME']

Process finished with exit code 0
```
What we did was we got the first item in world.data and then printed its attribute keys. The list shown in the output has a NAME key, but it is all in the uppercase. This is very common for Shapefiles whose data is contained in the DBF files.
Since we don't want to worry if the attributes' names are in the uppercase or lowercase, we have two possible solutions: convert the names at the moment of the import or convert the names on the fly when the attribute value is requested. Depending on your application...

Search

Your notes and bookmarks