day1/session5.tex
author Nishanth Amuluru <nishanth@fossee.in>
Tue, 14 Dec 2010 23:15:36 +0530
branchscipyin2010
changeset 456 a27ccfc118fb
parent 448 dc8d666a594a
permissions -rw-r--r--
removed sslc1.txt from circulate since it is not required

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Tutorial slides on Python.
%
% Author: Prabhu Ramachandran <prabhu at aero.iitb.ac.in>
% Copyright (c) 2005-2009, Prabhu Ramachandran
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\documentclass[14pt,compress]{beamer}
%\documentclass[draft]{beamer}
%\documentclass[compress,handout]{beamer}
%\usepackage{pgfpages} 
%\pgfpagesuselayout{2 on 1}[a4paper,border shrink=5mm]

% Modified from: generic-ornate-15min-45min.de.tex
\mode<presentation>
{
  \usetheme{Warsaw}
  \useoutertheme{infolines}
  \setbeamercovered{transparent}
}

\usepackage[english]{babel}
\usepackage[latin1]{inputenc}
%\usepackage{times}
\usepackage[T1]{fontenc}

% Taken from Fernando's slides.
\usepackage{ae,aecompl}
\usepackage{mathpazo,courier,euler}
\usepackage[scaled=.95]{helvet}

\definecolor{darkgreen}{rgb}{0,0.5,0}

\usepackage{listings}
\lstset{language=Python,
    basicstyle=\ttfamily\bfseries,
    commentstyle=\color{red}\itshape,
  stringstyle=\color{darkgreen},
  showstringspaces=false,
  keywordstyle=\color{blue}\bfseries}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Macros
\setbeamercolor{emphbar}{bg=blue!20, fg=black}
\newcommand{\emphbar}[1]
{\begin{beamercolorbox}[rounded=true]{emphbar} 
      {#1}
 \end{beamercolorbox}
}
\newcounter{time}
\setcounter{time}{0}
\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}}

\newcommand{\typ}[1]{\lstinline{#1}}

\newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}}  }

%%% This is from Fernando's setup.
% \usepackage{color}
% \definecolor{orange}{cmyk}{0,0.4,0.8,0.2}
% % Use and configure listings package for nicely formatted code
% \usepackage{listings}
% \lstset{
%    language=Python,
%    basicstyle=\small\ttfamily,
%    commentstyle=\ttfamily\color{blue},
%    stringstyle=\ttfamily\color{orange},
%    showstringspaces=false,
%    breaklines=true,
%    postbreak = \space\dots
% }


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Title page
\title[Exercises]{Exercises}

\author[FOSSEE] {FOSSEE}

\institute[IIT Bombay] {Department of Aerospace Engineering\\IIT Bombay}
\date[] {SciPy 2010, Tutorials}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\pgfdeclareimage[height=0.75cm]{iitmlogo}{iitmlogo}
%\logo{\pgfuseimage{iitmlogo}}


%% Delete this, if you do not want the table of contents to pop up at
%% the beginning of each subsection:
\AtBeginSubsection[]
{
  \begin{frame}<beamer>
    \frametitle{Outline}
    \tableofcontents[currentsection,currentsubsection]
  \end{frame}
}


% If you wish to uncover everything in a step-wise fashion, uncomment
% the following command: 
%\beamerdefaultoverlayspecification{<+->}

%\includeonlyframes{current,current1,current2,current3,current4,current5,current6}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% DOCUMENT STARTS
\begin{document}

\begin{frame}
  \titlepage
\end{frame}


\begin{frame}[fragile]
  \frametitle{Problem 1}
  \begin{columns}
    \column{0.5\textwidth}
    \hspace*{-0.5in}
    \includegraphics[height=2in, interpolate=true]{data/L-Tsq.png}
    \column{0.45\textwidth}
    \begin{block}{Example code}
    \tiny
    \begin{lstlisting}
l = []
t = []
for line in open('pendulum.txt'):
    point = line.split()
    l.append(float(point[0]))
    t.append(float(point[1]))
plot(l, t, '.')
    \end{lstlisting}
    \end{block}
  \end{columns}
  \begin{block}{Problem Statement}
      Tweak above code to plot data in file \typ{pos.txt}.
  \end{block}
\end{frame}

\begin{frame}
  \frametitle{Problem 1 cont...}
  \begin{itemize}
  \item Label both the axes.
  \item What kind of motion is this?
  \item Title the graph accordingly.
  \item Annotate the position where vertical velocity is zero.
  \end{itemize}
\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 2}
  \begin{columns}
    \column{0.5\textwidth}
    \hspace*{-0.5in}
    \includegraphics[height=2in, interpolate=true]{data/triangle}
    \column{0.45\textwidth}
    \begin{block}{Plot points given x and y coordinates}
    \tiny
    \begin{lstlisting}
In []: x = [3, 2, -2, 3]
In []: y = [1, -3, 4, 1]
In []: plot(x, y)
    \end{lstlisting}
    \end{block}
  \end{columns}
  Line can be plotted using arrays of coordinates.
  \pause
  \begin{block}{Problem statement}
    Write a Program that plots a regular n-gon(Let n = 5).
  \end{block}  
\end{frame}


\begin{frame}[fragile]
  \frametitle{Problem 3}
  \begin{columns}
    \column{0.5\textwidth}
    \hspace*{-0.5in}
    \includegraphics[height=2in, interpolate=true]{data/damp}
    \column{0.45\textwidth}
    \begin{block}{Damped Oscillation}
    \tiny
    \begin{lstlisting}
In []: t = linspace(0, 4*pi)
In []: plot(t, exp(t/10)*sin(t))
    \end{lstlisting}
    \end{block}
  \end{columns}
\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 3 cont...}
Create a sequence of images in which the damped
oscillator($e^{-t/10}sin(t)$) slowly evolves over time $t$.
\begin{columns}
\column{0.35\textwidth}
\includegraphics[width=1.5in,height=1.5in, interpolate=true]{data/plot2}
\column{0.35\textwidth}
\includegraphics[width=1.5in,height=1.5in, interpolate=true]{data/plot4}
\column{0.35\textwidth}
\includegraphics[width=1.5in,height=1.5in, interpolate=true]{data/plot6}
\end{columns}
\begin{block}{Hint}
\small
  \begin{lstlisting}
savefig('plot'+str(i)+'.png') #i is some int 
  \end{lstlisting}  
\end{block}
\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 4}
  \begin{lstlisting}
In []: x = imread('smoothing.gif')
In []: x.shape
Out[]: (256, 256)
In []: imshow(x,cmap=cm.gray)
In []: colorbar()
\end{lstlisting}
\emphbar{Replace each pixel with mean of neighboring pixels}
  \begin{center}
  \includegraphics[height=1in, interpolate=true]{data/neighbour}
  \end{center}
\end{frame}

\begin{frame}
  \begin{center}
    \includegraphics[height=3in, interpolate=true]{data/smoothing}    
  \end{center}
\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 4: Approach}
  For \typ{y} being resultant image:
  \begin{lstlisting}
y[1, 1] = x[0, 1]/4 + x[1, 0]/4 
          + x[2, 1]/4 + x[1, 2]/4    
  \end{lstlisting}
   \begin{columns}
    \column{0.45\textwidth}
    \hspace*{-0.5in}
    \includegraphics[height=1.5in, interpolate=true]{data/smoothing}
    \column{0.45\textwidth}
    \hspace*{-0.5in}
    \includegraphics[height=1.5in, interpolate=true]{data/after-filter}
  \end{columns}
   \begin{block}{Hint:}
     Use array Slicing.
   \end{block}
\end{frame}

\begin{frame}[fragile]
  \frametitle{Solution}
  \begin{lstlisting}
In []: y = zeros_like(x)
In []: y[1:-1,1:-1] = x[:-2,1:-1]/4 +
                      x[2:,1:-1]/4 +
                      x[1:-1,2:]/4 +
                      x[1:-1,:-2]/4
In []: imshow(y,cmap=cm.gray)
  \end{lstlisting}
\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 4 cont\ldots}
  \begin{itemize}
      \item Apply the smoothing operation repeatedly to the original
          image

      \item Subtract the smoothed image from the original to obtain
          the edges
  \end{itemize}
\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 5}

  What if you did the following in problem 4?
  \begin{lstlisting}
In []: y1[1:-1,1:-1] = (x[:-2,1:-1] +
                        x[2:,1:-1] +
                        x[1:-1,2:] +
                        x[1:-1,:-2])/4
  \end{lstlisting}

    Are the answers different? 

\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 5 cont\ldots}
  Why? The answer lies in the following:
\begin{lstlisting}
In []: x.dtype
Out[]: dtype('uint8')

In []: print x.itemsize
1

In []: z = x/4.0

In []: print z.dtype
float64

\end{lstlisting}
\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 5 cont\ldots}
What if you did this?

\begin{lstlisting}
x = imread('smoothing.gif')
y2 = zeros_like(x)
y2[1:-1,1:-1] = x[:-2,1:-1]/4. + 
                x[2:,1:-1]/4.  +
                x[1:-1,2:]/4.  +
                x[1:-1,:-2]/4.
\end{lstlisting}
\begin{itemize}
    \item Will the answer be any different from \typ{y}?
    \item What will the dtype of \typ{y2} be?
    \item Discuss what is going on!
\end{itemize}

\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 5 cont\ldots}
  Did you do the right thing to find the edges earlier in problem 4? Fix it! 
  
  Note that:
\begin{lstlisting}
In []: print x.dtype
uint8
In []: x1 = x.astype('float64')
In []: print x1.dtype
float64
In []: print x.dtype.char
d
In []: x.dtype.<TAB> # Explore!
\end{lstlisting}
\end{frame}

\begin{frame}[fragile]
  \frametitle{Problem 6}

  Edge detection looks much nicer with \typ{lena.png}, try it!  The
  gotcha is that it is a 4 component RGBA image with elements in the
  range $[0.0, 1.0]$.

\begin{lstlisting}
In []: x = imread('lena.png')

In []: print x.shape
(512, 512, 4)

In []: print x.min(), x.max()
0.0 1.0
\end{lstlisting}

 Repeat the edge detection with this image.

\end{frame}


\end{document}

%% \begin{frame}
%%   \frametitle{Problem 4}
%%   Legendre polynomials $P_n(x)$ are defined by the following recurrence relation

%% \center{$(n+1)P_{n+1}(x) - (2n+1)xP_n(x) + nP_{n-1}(x) = 0$}\\

%% with $P_0(x) = 1$, $P_1(x) = x$ and $P_2(x) = (3x^2 - 1)/2$. Compute the next three 
%%    Legendre polynomials and plot all 6 over the interval [-1,1].
%% \end{frame}

%% \begin{frame}[fragile] 
%% \frametitle{Problem Set 5}
%%   \begin{columns}
%%     \column{0.6\textwidth}
%%     \small{
%%     \begin{itemize}
%%       \item[3] Consider the iteration $x_{n+1} = f(x_n)$ where $f(x) = kx(1-x)$.  Plot the successive iterates of this process as explained below. 
%%     \end{itemize}}
%%     \column{0.35\textwidth}
%%     \hspace*{-0.5in}
%%   \includegraphics[height=1.6in, interpolate=true]{data/cobweb}  
%% \end{columns}
%% \end{frame}

%% \begin{frame}
%%   \frametitle{Problem Set 5.3}
%%   Plot the cobweb plot as follows:
%%   \begin{enumerate}
%%     \item Start at $(x_0, 0)$ ($\implies$ i=0)
%%     \item Draw a line to $(x_i, f(x_i))$
%%     \item Set $x_{i+1} = f(x_i)$
%%     \item Draw a line to $(x_{i+1}, x_{i+1})$
%%     \item $(i\implies i+1)$ 
%%     \item Repeat from 2 for as long as you want 
%%   \end{enumerate}
%% \inctime{20}
%% \end{frame}