examples/plotting-data.lisp

   1 ;;; -*- mode: lisp -*-
   2
   3 ;;; Time-stamp: <2009-07-06 18:45:27 tony>
   4 ;;; Creation:   <2009-03-12 17:14:56 tony>
   5 ;;; File:       plotting-data.lisp
   6 ;;; Author:     AJ Rossini <blindglobe@gmail.com>
   7 ;;; Copyright:  (c)2009--, AJ Rossini.  BSD, LLGPL, or GPLv2, depending
   8 ;;;             on how it arrives.
   9 ;;; Purpose:    Example of generating plots.
  10
  11 ;;; What is this talk of 'release'? Klingons do not make software
  12 ;;; 'releases'.  Our software 'escapes', leaving a bloody trail of
  13 ;;; designers and quality assurance people in its wake.
  14
  15 ;;; This organization and structure is new to the 21st Century
  16 ;;; version..   Think, "21st Century Schizoid Man".
  17
  18 ;; SETUP FOR PLOT EXAMPLE:
  19
  20 (defpackage :cl-2d-user-x11
  21   (:use :cl :cl-cairo2  :cl-2d :cl-numlib :cl-colors :bind :cls))
  22
  23 (in-package :cl-2d-user-x11)
  24
  25 ;; PLOT EXAMPLE
  26 #+nil
  27 (progn
  28
  29
  30   ;; this is how you create an X11 frame.  If you supply a
  31   ;; background-color to as-frame, each plot will clear the frame with
  32   ;; this color.
  33
  34   (defparameter *frame1* (as-frame (create-xlib-image-context 300 300)
  35                                   :background-color +white+))
  36
  37   ;; or netbook size, picture is similar but on a lower-res display window.
  38   (defparameter *frame2* (as-frame (create-xlib-image-context 200 200)
  39                                   :background-color +white+))
  40
  41   (plot-function *frame1*
  42                  #'exp (interval-of 0 2)
  43                  :x-title "x"
  44                  :y-title "exp(x)")
  45
  46   ;; split the frame, and you can draw on the subframes independently.
  47   ;; I do this a lot.
  48
  49   (bind ((#2A((f1 f2) (f3 f4))
  50              (split-frame *frame2* (percent 50) (percent 50))))
  51         (defparameter *f1* f1)
  52         (defparameter *f2* f2)
  53         (defparameter *f3* f3)
  54         (defparameter *f4* f4))
  55
  56
  57   (bind ((#2A((f1 f2) (f3 f4))
  58              (split-frame *frame2* (percent 75) (percent 25))))
  59         (defparameter *f1* f1)
  60         (defparameter *f2* f2)
  61         (defparameter *f3* f3)
  62         (defparameter *f4* f4))
  63
  64   (plot-function *f1* #'sin (interval-of 0 2) :x-title "x" :y-title "sin(x)")
  65   (plot-function *f2* #'cos (interval-of 0 2) :x-title "x" :y-title "cos(x)")
  66   (plot-function *f3* #'tan (interval-of 0 2) :x-title "x" :y-title "tan(x)")
  67   (plot-function *f4* #'/ (interval-of 0 2) :x-title "x" :y-title "1/x")
  68
  69   (clear *frame1*)
  70
  71
  72   (let* ((n 500)
  73          (xs (num-sequence :from 0 :to 10 :length n))
  74          (ys (map 'vector #'(lambda (x) (+ x 8 (random 4.0))) xs))
  75          (weights (replicate #'(lambda () (1+ (random 10))) n 'fixnum))
  76          (da (plot-simple *frame1* (interval-of 0 10) (interval-of 10 20)
  77                           :x-title "x" :y-title "y")))
  78     (draw-symbols da xs ys :weights weights))
  79
  80   (xlib-image-context-to-png (context *frame1*) "/home/tony/test1.png")
  81   (xlib-image-context-to-png (context *frame2*) "/home/tony/test2.png")
  82   )
  83
  84
  85 ;;; EXAMPLE FOR DSC2009
  86   (defparameter *frame2* (as-frame (create-xlib-image-context 400 400)
  87                                   :background-color +white+))
  88
  89   (bind ((#2A((f1 f2) (f3 f4))
  90              (split-frame *frame2* (percent 50) (percent 50))))
  91         (defparameter *f1* f1)
  92         (defparameter *f2* f2)
  93         (defparameter *f3* f3)
  94         (defparameter *f4* f4))
  95   (plot-function *f1* #'sin (interval-of 0 2) :x-title "x" :y-title "sin(x)")
  96   (plot-function *f2* #'cos (interval-of 0 2) :x-title "x" :y-title "cos(x)")
  97   (plot-function *f3* #'tan (interval-of 0 2) :x-title "x" :y-title "tan(x)")
  98
  99  (num-sequence :from 0 :to 10 :length 30)
 100
 101   (let* ((n 500)
 102          (xs (num-sequence :from 0 :to 10 :length n))
 103          (ys (map 'vector #'(lambda (x) (+ x 8 (random 4.0))) xs))
 104          (weights (replicate #'(lambda () (1+ (random 10))) n 'fixnum))
 105          (da (plot-simple *f4* (interval-of 0 10) (interval-of 10 20)
 106                           :x-title "x" :y-title "y")))
 107     (draw-symbols da xs ys :weights weights))
 108   (xlib-image-context-to-png (context *f1*) "/home/tony/test1.png")
 109   (xlib-image-context-to-png (context *frame2*) "/home/tony/test2.png")
 110   (destroy (context  *frame2*))
 111
 112