Home Accessibility Courses Twitter The Mouth Facebook Resources Site Map About Us Contact
For 2021 - online Python 3 training - see ((here)).

Our plans were to retire in summer 2020 and see the world, but Coronavirus has lead us into a lot of lockdown programming in Python 3 and PHP 7.
We can now offer tailored online training - small groups, real tutors - works really well for groups of 4 to 14 delegates. Anywhere in the world; course language English.

Please ask about private 'maintenance' training for Python 2, Tcl, Perl, PHP, Lua, etc.
3d and contour plots through numpy and matplotlib
numpy, scipy and matplotlib example from a Well House Consultants training course
More on numpy, scipy and matplotlib [link]

This example is described in the following article(s):
   • 3D graphics - web site usage - simple matplotlib and python example - [link]
   • Learning more about our web site - and learning how to learn about yours - [link]

Source code: xyz.py Module: Y118
import numpy
import pylab
import mpl_toolkits.mplot3d.axes3d as p3
import re
import sys

# Extract a count of hour by hour accesses from recent log files

hourfinder = re.compile(r'\d\d/\w\w\w/\d\d\d\d:(\d\d):\d\d:\d\d')
counter = []
x = []
y = []

files = []
for k in range(18,31):
        files.append("ac_201009"+ "%02d" % k)
for k in range(1,11):
        files.append("ac_201010"+ "%02d" % k)

for day in range(len(files)):
        counter.append([0] * 24)
        x.append([day] * 24)
        for hit in open(files[day]).xreadlines():
                hour = int((hourfinder.findall(hit))[0])
                counter[day][hour] += 1

# Apply some gentle smoothing
# (wasteful - only needed if "smooth" parameter given!)

smooch = []
for day in range(len(counter)):
        smooch.append([0] * len(counter[day]))
        for hour in range(len(counter[day])):
                tsf = 0.0
                fsf = 0.0
                for k in range(-1,2):
                        kx = day + k
                        if kx < 0: continue
                        if kx >= len(counter): continue
                        for j in range(-1,2):
                                jx = hour + j
                                if jx < 0: continue
                                if jx >= len(counter[kx]): continue
                                factor = 0.3
                                if j == k and j == 0: factor = 1.0
                                tsf += counter[kx][jx] * factor
                                fsf += factor
                smooch[day][hour] = tsf / fsf

# Convert to numpy arrays for plot functions

x = numpy.array(x)
y = numpy.array(y)

if len(sys.argv) == 3 and sys.argv[2] == "smooth":
        z = numpy.array(smooch)
        z = numpy.array(counter)

# Set up for your plot

if len(sys.argv) < 2 or sys.argv[1] == "0":


        # pylab.contour(x,y,z)

        ax = p3.Axes3D(fig)

        if sys.argv[1] == "1": ax.plot_wireframe(x,y,z)
        elif sys.argv[1] == "2": ax.plot_surface(x,y,z)
        elif sys.argv[1] == "3": ax.contour3D(x,y,z)
        else: ax.contourf3D(x,y,z)

        ax.set_zlabel('Hit Count')

# go plot it!

Learn about this subject
This module and example are covered on the following public courses:
 * Learning to program in Python
 * Python Programming
Also available on on site courses for larger groups

Books covering this topic
Yes. We have over 700 books in our library. Books covering Python are listed here and when you've selected a relevant book we'll link you on to Amazon to order.

Other Examples
This example comes from our "numpy, scipy and matplotlib" training module. You'll find a description of the topic and some other closely related examples on the "numpy, scipy and matplotlib" module index page.

Full description of the source code
You can learn more about this example on the training courses listed on this page, on which you'll be given a full set of training notes.

Many other training modules are available for download (for limited use) from our download centre under an Open Training Notes License.

Other resources
• Our Solutions centre provides a number of longer technical articles.
• Our Opentalk forum archive provides a question and answer centre.
The Horse's mouth provides a daily tip or thought.
• Further resources are available via the resources centre.
• All of these resources can be searched through through our search engine
• And there's a global index here.

Purpose of this website
This is a sample program, class demonstration or answer from a training course. It's main purpose is to provide an after-course service to customers who have attended our public private or on site courses, but the examples are made generally available under conditions described below.

Web site author
This web site is written and maintained by Well House Consultants.

Conditions of use
Past attendees on our training courses are welcome to use individual examples in the course of their programming, but must check the examples they use to ensure that they are suitable for their job. Remember that some of our examples show you how not to do things - check in your notes. Well House Consultants take no responsibility for the suitability of these example programs to customer's needs.

This program is copyright Well House Consultants Ltd. You are forbidden from using it for running your own training courses without our prior written permission. See our page on courseware provision for more details.

Any of our images within this code may NOT be reused on a public URL without our prior permission. For Bona Fide personal use, we will often grant you permission provided that you provide a link back. Commercial use on a website will incur a license fee for each image used - details on request.

You can Add a comment or ranking to this page

© WELL HOUSE CONSULTANTS LTD., 2021: 48 Spa Road • Melksham, Wiltshire • United Kingdom • SN12 7NY
PH: 01144 1225 708225 • EMAIL: info@wellho.net • WEB: http://www.wellho.net • SKYPE: wellho

PAGE: http://www.wellho.net/resources/ex.php4 • PAGE BUILT: Sun Oct 11 14:50:09 2020 • BUILD SYSTEM: JelliaJamb