Displaying all 60 zones of the UTM projection#

This example displays all 60 zones of the Universal Transverse Mercator projection next to each other in a figure.

First we create a figure with 60 subplots in one row. Next we set the projection of each axis in the figure to a specific UTM zone. Then we add coastlines, gridlines and the number of the zone. Finally we add a supertitle and display the figure.

UTM Projection - Zones, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60
import matplotlib.pyplot as plt

import cartopy.crs as ccrs


def main():
    # Create a list of integers from 1 - 60
    zones = range(1, 61)

    # Create a figure
    fig = plt.figure(figsize=(18, 6))

    # Loop through each zone in the list
    for zone in zones:

        # Add GeoAxes object with specific UTM zone projection to the figure
        ax = fig.add_subplot(1, len(zones), zone,
                             projection=ccrs.UTM(zone=zone,
                                                 southern_hemisphere=True))

        # Add coastlines, gridlines and zone number for the subplot
        ax.coastlines(resolution='110m')
        ax.gridlines()
        ax.set_title(zone)

    # Add a supertitle for the figure
    fig.suptitle("UTM Projection - Zones")

    # Display the figure
    plt.show()


if __name__ == '__main__':
    main()

Total running time of the script: ( 0 minutes 36.923 seconds)

Gallery generated by Sphinx-Gallery