PS6-DCT-Soln-correction
|
|
- Joan Prudence Fox
- 6 years ago
- Views:
Transcription
1 PS6-DCT-Soln-correction Unknown Author March 18, 2014 Part I DCT: Discrete Cosine Transform DCT is a linear map A R N N such that the N real numbers x 0,..., x N 1 are transformed into the N real numbers X 0,..., X N 1 according to the formula N 1 ( π X k = x x i cos N (k + 1 ) 2 )i. (1) i=1 1 Question 1: Consider the case where N = 2. Write Python function using (1) to compute X 0, X 1 given x 0, x 1. Specifically your function should take in x 0, x 1 as arguments and return X 0, X 1. In [1]: import numpy as np """ Your code goes here """ def f(x0,x1): N=2 X0=x0+2*x1*np.cos(np.pi/N*(1+1/2)*0) X1=x0+2*x1*np.cos(np.pi/N*(1+1/2)*1) return X0,X1 def vald(i,k): N=8 return 2*np.cos((np.pi/N) * i * (k + 0.5) ) #x[n]*cos(pi*(k+0.5)*n/n) Verify that your function works by computing X 0, X 1 given x 0 = 2, x 1 = 4. In [2]: np.set_printoptions(precision=3) """ Check your function here """ y=f(2.0,4.0) print "%.2f "*len(y) % y
2 2 Question 2: JPEG compression uses the DCT with N = 8. Create a matrix called D that is the matrix representation of the discrete cosine transform, i.e. when you multiply D by the vector x 0. x N 1 you should get the vector X 0. X N 1 In [3]: """ your code goes here""" N = 8 D = np.zeros((n, N)) for i in range(n): D[0,i] = 1 for i in xrange(1,n): for k in xrange(n): D[i,k] = 2*np.cos((np.pi/N) * i * (k + 0.5) ) D=D.T print D [[ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ]] It is interesting to note that the rows of the matrix D are basis vectors which oscillate with successively higher spatial frequencies. Below is code to plot them. In [4]: import matplotlib.pyplot as plt %matplotlib inline # Graphing helper function from our last assignment def setup_graph(title=, x_label=, y_label=, fig_size=none): fig = plt.figure() if fig_size!= None: fig.set_size_inches(fig_size[0], fig_size[1]) ax = fig.add_subplot(111) ax.set_title(title) ax.set_xlabel(x_label) ax.set_ylabel(y_label) fig=plt.figure(figsize=(9,12)) for u in xrange(n): setup_graph(title= u= +str(u), x_label=, y_label=, fig_size=(6,3)) _=plt.plot(d[u, :]) _=plt.plot(d[u, :], ro )
3 <matplotlib.figure.figure at 0x10f068490>
4
5
6 3 Question 3: Consider the vector x = Compute the DCT using the the map D you defined above and compare it to the result you get when you use the function dct from scipy.fftpack. Use type=3 in the dct function.
7 In [6]: from scipy.fftpack import dct, idct x=np.array([[1.0, 2, 3, 4, 5, 6, 7, 8]]) """ Your Code Goes Here""" print np.dot(d,x.t) print (dct(x,type=3)).t [[ ] [ ] [ ] [ ] [ 6.549] [ ] [ 2.184] [ ]] [[ ] [ ] [ ] [ ] [ 6.549] [ ] [ 2.184] [ ]] 4 Question 4: Now, we will take the image we have provided you with and apply the 2D DCT to it. In [7]: import matplotlib.image as mpimg img = mpimg.imread( montypython.png ) p=plt.imshow(img, origin= upper )
8 The above picture is of Monty Python (Python s namesake). Below we select a portion of the image to which we will apply the 2-D DCT. Let s just look at values from one tiny 8 x 8 block (which is what s used JPEG compression). Then, we use a false color spectrum to visualize the pixel intensity. In [8]: img.shape b = img[50:58, 50:58, 0] def show_image(img): plt.imshow(img) plt.colorbar() show_image(b)
9 Now, we display the actual values in b. In [9]: Out [9]: b array([[ 0.157, 0.161, 0.165, 0.165, 0.157, 0.153, 0.157, 0.145], [ 0.165, 0.161, 0.165, 0.165, 0.161, 0.161, 0.157, 0.133], [ 0.161, 0.153, 0.149, 0.149, 0.149, 0.153, 0.161, 0.133], [ 0.165, 0.157, 0.149, 0.145, 0.145, 0.153, 0.169, 0.145], [ 0.169, 0.161, 0.149, 0.149, 0.149, 0.157, 0.173, 0.157], [ 0.157, 0.149, 0.141, 0.141, 0.141, 0.149, 0.165, 0.157], [ 0.165, 0.157, 0.141, 0.137, 0.133, 0.141, 0.153, 0.161], [ 0.153, 0.153, 0.153, 0.153, 0.145, 0.145, 0.153, 0.161]], dtype=float32) The 2-D DCT is just the 1-D DCT applied to every column first and then applied to every row. Write a function to compute the 2-D DCT. Compute the 2-D dct of b and display it as an array of numbers as well as using the false color spectrum to visualize as above. In [15]: """ your Code for computing 2-D DCT goes here""" def dodct(grid): return np.dot(np.dot(d, grid), D.T) """ your Code for viualizing""" print dodct(b) show_image(dodct(b))
10 [[ 1.589e e e e e e e e-02] [ e e e e e e e e-02] [ 2.862e e e e e e e e-02] [ e e e e e e e e-02] [ 1.485e e e e e e e e-02] [ e e e e e e e e-03] [ 3.353e e e e e e e e-02] [ e e e e e e e e-02]] Now, write a function to undo the 2-D DCT, that is write a function to compute the inverse of the 2-D DCT. Print out the values of the inverse of the 2-D dct applied to b and use the false color spectrum visualization technique to view the result. Compare to show_image(b).
11 In [16]: """ Your code goes here""" def undodct(grid): Di=np.linalg.inv(D) return np.dot(np.dot(di, grid), Di.T) print undodct(dodct(b)) tiny_do_undo = undodct(dodct(b)) show_image(tiny_do_undo) # Yup, looks the same. [[ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ]] 5 Question 5: Apply the built in scipy functions dct and idct to the image b and compare what you have done above using the show_image function. You should get the same images. Note that the scipy function dct does the 1-D DCT. So, you must figure out how to use it to do the 2-D DCT. In [17]: """Your code for computing the 2-D dct using dct() and for plotting it using show_image g=dct(dct(b, type=3).t, type=3) show_image(g)
12 In [18]: """Your code for computing the inverse 2-D dct using idct() and for plotting it using s ig=idct(idct(g, type=3).t, type=3) show_image(ig) In []:
MATPLOTLIB. Python for computational science November 2012 CINECA.
MATPLOTLIB Python for computational science 19 21 November 2012 CINECA m.cestari@cineca.it Introduction (1) plotting the data gives us visual feedback in the working process Typical workflow: write a python
More informationIntroduction to Matplotlib: 3D Plotting and Animations
1 Introduction to Matplotlib: 3D Plotting and Animations Lab Objective: 3D plots and animations are useful in visualizing solutions to ODEs and PDEs found in many dynamics and control problems. In this
More informationHW0 v3. October 2, CSE 252A Computer Vision I Fall Assignment 0
HW0 v3 October 2, 2018 1 CSE 252A Computer Vision I Fall 2018 - Assignment 0 1.0.1 Instructor: David Kriegman 1.0.2 Assignment Published On: Tuesday, October 2, 2018 1.0.3 Due On: Tuesday, October 9, 2018
More informationImage Compression Techniques
ME 535 FINAL PROJECT Image Compression Techniques Mohammed Abdul Kareem, UWID: 1771823 Sai Krishna Madhavaram, UWID: 1725952 Palash Roychowdhury, UWID:1725115 Department of Mechanical Engineering University
More informationINTRODUCTION TO DATA VISUALIZATION WITH PYTHON. Working with 2D arrays
INTRODUCTION TO DATA VISUALIZATION WITH PYTHON Working with 2D arrays Reminder: NumPy arrays Homogeneous in type Calculations all at once Indexing with brackets: A[index] for 1D array A[index0, index1]
More informationPartial Differential Equations II: 2D Laplace Equation on 5x5 grid
Partial Differential Equations II: 2D Laplace Equation on 5x5 grid Sam Sinayoko Numerical Methods 5 Contents 1 Learning Outcomes 2 2 Introduction 3 3 Laplace equation in 2D 3 4 Discretisation 3 4.1 Meshing:
More informationUser-Defined Function
ENGR 102-213 (Socolofsky) Week 11 Python scripts In the lecture this week, we are continuing to learn powerful things that can be done with userdefined functions. In several of the examples, we consider
More informationInterpolation and curve fitting
CITS2401 Computer Analysis and Visualization School of Computer Science and Software Engineering Lecture 9 Interpolation and curve fitting 1 Summary Interpolation Curve fitting Linear regression (for single
More information1 CSE 252A Computer Vision I Fall 2017
Assignment 1 CSE A Computer Vision I Fall 01 1.1 Assignment This assignment contains theoretical and programming exercises. If you plan to submit hand written answers for theoretical exercises, please
More informationNumerical Calculations
Fundamentals of Programming (Python) Numerical Calculations Sina Sajadmanesh Sharif University of Technology Some slides have been adapted from Scipy Lecture Notes at http://www.scipy-lectures.org/ Outline
More informationHomework 11 - Debugging
1 of 7 5/28/2018, 1:21 PM Homework 11 - Debugging Instructions: Fix the errors in the following problems. Some of the problems are with the code syntax, causing an error message. Other errors are logical
More informationPython for Scientists
High level programming language with an emphasis on easy to read and easy to write code Includes an extensive standard library We use version 3 History: Exists since 1991 Python 3: December 2008 General
More informationThe SciPy Stack. Jay Summet
The SciPy Stack Jay Summet May 1, 2014 Outline Numpy - Arrays, Linear Algebra, Vector Ops MatPlotLib - Data Plotting SciPy - Optimization, Scientific functions TITLE OF PRESENTATION 2 What is Numpy? 3rd
More informationIntro to Research Computing with Python: Visualization
Intro to Research Computing with Python: Visualization Erik Spence SciNet HPC Consortium 20 November 2014 Erik Spence (SciNet HPC Consortium) Visualization 20 November 2014 1 / 29 Today s class Today we
More informationSkills Quiz - Python Edition Solutions
'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ EGR 103L Fall 2017 Skills Quiz - Python Edition Solutions Michael R. Gustafson II Name (please print): NetID (please print): In keeping with the Community
More informationARTIFICIAL INTELLIGENCE AND PYTHON
ARTIFICIAL INTELLIGENCE AND PYTHON DAY 1 STANLEY LIANG, LASSONDE SCHOOL OF ENGINEERING, YORK UNIVERSITY WHAT IS PYTHON An interpreted high-level programming language for general-purpose programming. Python
More informationMS6021 Scientific Computing. TOPICS: Python BASICS, INTRO to PYTHON for Scientific Computing
MS6021 Scientific Computing TOPICS: Python BASICS, INTRO to PYTHON for Scientific Computing Preliminary Notes on Python (v MatLab + other languages) When you enter Spyder (available on installing Anaconda),
More informationLab 9 - Linear Model Selection in Python
Lab 9 - Linear Model Selection in Python March 7, 2016 This lab on Model Validation using Validation and Cross-Validation is a Python adaptation of p. 248-251 of Introduction to Statistical Learning with
More informationSomoclu Python Documentation
Somoclu Python Documentation Release 1.7.5 Peter Wittek, Shi Chao Gao Mar 01, 2018 Contents 1 Introduction 1 1.1 Copyright and License.......................................... 1 1.2 Acknowledgment.............................................
More informationData Science and Machine Learning Essentials
Data Science and Machine Learning Essentials Lab 3B Building Models in Azure ML By Stephen Elston and Graeme Malcolm Overview In this lab, you will learn how to use R or Python to engineer or construct
More informationLatent Semantic Analysis. sci-kit learn. Vectorizing text. Document-term matrix
Latent Semantic Analysis Latent Semantic Analysis (LSA) is a framework for analyzing text using matrices Find relationships between documents and terms within documents Used for document classification,
More informationlof April 23, Improving performance of Local outlier factor with KD-Trees
lof April 23, 2014 1 Improving performance of Local outlier factor with KD-Trees Local outlier factor (LOF) is an outlier detection algorithm, that detects outliers based on comparing local density of
More informationLecture 15: High Dimensional Data Analysis, Numpy Overview
Lecture 15: High Dimensional Data Analysis, Numpy Overview Chris Tralie, Duke University 3/3/2016 Announcements Mini Assignment 3 Out Tomorrow, due next Friday 3/11 11:55PM Rank Top 3 Final Project Choices
More informationSimpleITK Spatial Transformations
SimpleITK Spatial Transformations Summary: 1. Points are represented by vector-like data types: Tuple, Numpy array, List. 2. Matrices are represented by vector-like data types in row major order. 3. Default
More informationPart VI. Scientific Computing in Python. Alfredo Parra : Scripting with Python Compact Max-PlanckMarch 6-10,
Part VI Scientific Computing in Python Compact Course @ Max-PlanckMarch 6-10, 2017 63 Doing maths in Python Standard sequence types (list, tuple,... ) Can be used as arrays Can contain different types
More informationCME 193: Introduction to Scientific Python Lecture 6: Numpy, Scipy, Matplotlib
CME 193: Introduction to Scientific Python Lecture 6: Numpy, Scipy, Matplotlib Nolan Skochdopole stanford.edu/class/cme193 6: Numpy, Scipy, Matplotlib 6-1 Contents Homeworks and Project Numpy Scipy Matplotlib
More informationImage Processing (1) Basic Concepts and Introduction of OpenCV
Intelligent Control Systems Image Processing (1) Basic Concepts and Introduction of OpenCV Shingo Kagami Graduate School of Information Sciences, Tohoku University swk(at)ic.is.tohoku.ac.jp http://www.ic.is.tohoku.ac.jp/ja/swk/
More informationmaking connections general transit feed specification stop names and stop times storing the connections in a dictionary
making connections 1 CTA Tables general transit feed specification stop names and stop times storing the connections in a dictionary 2 CTA Schedules finding connections between stops sparse matrices in
More informationL15. 1 Lecture 15: Data Visualization. July 10, Overview and Objectives. 1.2 Part 1: Introduction to matplotlib
L15 July 10, 2017 1 Lecture 15: Data Visualization CSCI 1360E: Foundations for Informatics and Analytics 1.1 Overview and Objectives Data visualization is one of, if not the, most important method of communicating
More informationCHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT
CHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT 9.1 Introduction In the previous chapters the inpainting was considered as an iterative algorithm. PDE based method uses iterations to converge
More informationBi 1x Spring 2014: Plotting and linear regression
Bi 1x Spring 2014: Plotting and linear regression In this tutorial, we will learn some basics of how to plot experimental data. We will also learn how to perform linear regressions to get parameter estimates.
More informationImage Manipulation in MATLAB Due Monday, July 17 at 5:00 PM
Image Manipulation in MATLAB Due Monday, July 17 at 5:00 PM 1 Instructions Labs may be done in groups of 2 or 3 (i.e., not alone). You may use any programming language you wish but MATLAB is highly suggested.
More informationIntroduction to Python Practical 1
Introduction to Python Practical 1 Daniel Carrera & Brian Thorsbro October 2017 1 Introduction I believe that the best way to learn programming is hands on, and I tried to design this practical that way.
More information1 2 (3 + x 3) x 2 = 1 3 (3 + x 1 2x 3 ) 1. 3 ( 1 x 2) (3 + x(0) 3 ) = 1 2 (3 + 0) = 3. 2 (3 + x(0) 1 2x (0) ( ) = 1 ( 1 x(0) 2 ) = 1 3 ) = 1 3
6 Iterative Solvers Lab Objective: Many real-world problems of the form Ax = b have tens of thousands of parameters Solving such systems with Gaussian elimination or matrix factorizations could require
More informationLecture 11 : Discrete Cosine Transform
Lecture 11 : Discrete Cosine Transform Moving into the Frequency Domain Frequency domains can be obtained through the transformation from one (time or spatial) domain to the other (frequency) via Fourier
More informationPYTHON DATA VISUALIZATIONS
PYTHON DATA VISUALIZATIONS from Learning Python for Data Analysis and Visualization by Jose Portilla https://www.udemy.com/learning-python-for-data-analysis-and-visualization/ Notes by Michael Brothers
More informationLogistic Regression with a Neural Network mindset
Logistic Regression with a Neural Network mindset Welcome to your first (required) programming assignment! You will build a logistic regression classifier to recognize cats. This assignment will step you
More informationROI Based Image Compression in Baseline JPEG
168-173 RESEARCH ARTICLE OPEN ACCESS ROI Based Image Compression in Baseline JPEG M M M Kumar Varma #1, Madhuri. Bagadi #2 Associate professor 1, M.Tech Student 2 Sri Sivani College of Engineering, Department
More informationLab 1 - Basic ipython Tutorial (EE 126 Fall 2014)
Lab 1 - Basic ipython Tutorial (EE 126 Fall 2014) modified from Berkeley Python Bootcamp 2013 https://github.com/profjsb/python-bootcamp and Python for Signal Processing http://link.springer.com/book/10.1007%2f978-3-319-01342-8
More informationINVESTIGATION OF ASSOCIATIVE PROPERTY OF MATRIX MULTIPLICATION IN DISCRETE COSINE TRANSFORMS FOR FIXED-POINT CALCULATION SYSTEMS
УДК 621.3 Ayman M.M. Ahmed Author contact information: a_ymn2002@yahoo.com INVESTIGATION OF ASSOCIATIVE PROPERTY OF MATRIX MULTIPLICATION IN DISCRETE COSINE TRANSFORMS FOR FIXED-POINT CALCULATION SYSTEMS
More informationMANY image and video compression standards such as
696 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS FOR VIDEO TECHNOLOGY, VOL 9, NO 5, AUGUST 1999 An Efficient Method for DCT-Domain Image Resizing with Mixed Field/Frame-Mode Macroblocks Changhoon Yim and
More informationData Science and Machine Learning Essentials
Data Science and Machine Learning Essentials Lab 3C Evaluating Models in Azure ML By Stephen Elston and Graeme Malcolm Overview In this lab, you will learn how to evaluate and improve the performance of
More informationMatPlotTheme Documentation
MatPlotTheme Documentation Release 0.1.2 James Yu July 31, 2014 Contents 1 Contents 3 1.1 Overview................................................. 3 1.2 Gallery..................................................
More informationSession 1 Use test driven development (i.e. write the tests first) to design functions to give the square, cube and an arbitary power N for a number a. In [1]: import unittest def square(a): return a**2
More informationENGR (Socolofsky) Week 07 Python scripts
ENGR 102-213 (Socolofsky) Week 07 Python scripts A couple programming examples for this week are embedded in the lecture notes for Week 7. We repeat these here as brief examples of typical array-like operations
More informationPyPlot. The plotting library must be imported, and we will assume in these examples an import statement similar to those for numpy and math as
Geog 271 Geographic Data Analysis Fall 2015 PyPlot Graphicscanbeproducedin Pythonviaavarietyofpackages. We willuseapythonplotting package that is part of MatPlotLib, for which documentation can be found
More informationTutorial Four: Linear Regression
Tutorial Four: Linear Regression Imad Pasha Chris Agostino February 25, 2015 1 Introduction When looking at the results of experiments, it is critically important to be able to fit curves to scattered
More informationHardware Optimized DCT/IDCT Implementation on Verilog HDL
Hardware Optimized DCT/IDCT Implementation on Verilog HDL ECE 734 In this report, I explore 4 implementations for hardware based pipelined DCT/IDCT in Verilog HDL. Conventional DCT/IDCT implementations
More informationzap Documentation Release 1.0.dev86 Kurt Soto
zap Documentation Release 1.0.dev86 Kurt Soto February 03, 2016 Contents 1 Installation 3 1.1 Requirements............................................... 3 1.2 Steps...................................................
More informationCSE 152 : Introduction to Computer Vision, Spring 2018 Assignment 5
CSE 152 : Introduction to Computer Vision, Spring 2018 Assignment 5 Instructor: Ben Ochoa Assignment Published On: Wednesday, May 23, 2018 Due On: Saturday, June 9, 2018, 11:59 PM Instructions Review the
More informationCSC 1315! Data Science
CSC 1315! Data Science Data Visualization Based on: Python for Data Analysis: http://hamelg.blogspot.com/2015/ Learning IPython for Interactive Computation and Visualization by C. Rossant Plotting with
More informationSearch. The Nearest Neighbor Problem
3 Nearest Neighbor Search Lab Objective: The nearest neighbor problem is an optimization problem that arises in applications such as computer vision, pattern recognition, internet marketing, and data compression.
More informationCME 193: Introduction to Scientific Python Lecture 5: Numpy, Scipy, Matplotlib
CME 193: Introduction to Scientific Python Lecture 5: Numpy, Scipy, Matplotlib Sven Schmit stanford.edu/~schmit/cme193 5: Numpy, Scipy, Matplotlib 5-1 Contents Second part of course Numpy Scipy Matplotlib
More informationSOME CONCEPTS IN DISCRETE COSINE TRANSFORMS ~ Jennie G. Abraham Fall 2009, EE5355
SOME CONCEPTS IN DISCRETE COSINE TRANSFORMS ~ Jennie G. Abraham Fall 009, EE5355 Under Digital Image and Video Processing files by Dr. Min Wu Please see lecture10 - Unitary Transform lecture11 - Transform
More informationPlotting with an introduction to functions
Plotting with CERN@school: an introduction to functions Twitter: @nicoleshearer93 N. Shearer a, T. Whyntie b, c a Durham University, b Langton Star Centre, c Queen Mary University of London Coding with
More informationPractical 06: Plotting and the Verlet integrator Documentation
Practical 06: Plotting and the Verlet integrator Documentation Release 1.0 Oliver Beckstein February 07, 2013 CONTENTS 1 Practical 06 3 1.1 IPython and pylab............................................
More informationCSE237A: Final Project Mid-Report Image Enhancement for portable platforms Rohit Sunkam Ramanujam Soha Dalal
CSE237A: Final Project Mid-Report Image Enhancement for portable platforms Rohit Sunkam Ramanujam (rsunkamr@ucsd.edu) Soha Dalal (sdalal@ucsd.edu) Project Goal The goal of this project is to incorporate
More informationMathematical Programming
Mathematical Programming Example #9-1. Minimize the following problem related to the linear programming. Minimize: f(x, y) = x + 4y Subject to: 3x + y 6 x + 2y 4 y 3 where x TA: Junhee Lee Python Code
More informationTF Mutiple Hidden Layers: Regression on Boston Data
TF Mutiple Hidden Layers: Regression on Boston Data This is adapted from Frossard's tutorial (http://www.cs.toronto.edu/~frossard/post/tensorflow/). This approach is not batched, and the number of layers
More informationAPPM 2360 Project 2 Due Nov. 3 at 5:00 PM in D2L
APPM 2360 Project 2 Due Nov. 3 at 5:00 PM in D2L 1 Introduction Digital images are stored as matrices of pixels. For color images, the matrix contains an ordered triple giving the RGB color values at each
More informationFor examples, documentation, tutorials, etc, see Astropy at ( # For retrieving an image from a URL
Astronomy example 1 Downloading images and writing FITS files For examples, documentation, tutorials, etc, see Astropy at http://www.astropy.org (http://www.astropy.org) In [1]: import scipy as sp import
More informationLab 8 - Subset Selection in Python
Lab 8 - Subset Selection in Python March 2, 2016 This lab on Subset Selection is a Python adaptation of p. 244-247 of Introduction to Statistical Learning with Applications in R by Gareth James, Daniela
More informationIntroduction to Wavelets
Lab 11 Introduction to Wavelets Lab Objective: In the context of Fourier analysis, one seeks to represent a function as a sum of sinusoids. A drawback to this approach is that the Fourier transform only
More informationPython Crash Course Numpy, Scipy, Matplotlib
Python Crash Course Numpy, Scipy, Matplotlib That is what learning is. You suddenly understand something you ve understood all your life, but in a new way. Doris Lessing Steffen Brinkmann Max-Planck-Institut
More informationPython Matplotlib. MACbioIDi February March 2018
Python Matplotlib MACbioIDi February March 2018 Introduction Matplotlib is a Python 2D plotting library Its origins was emulating the MATLAB graphics commands It makes heavy use of NumPy Objective: Create
More informationNonlinear curve-fitting example
Nonlinear curve-fitting example Implementation of curve-fitting in Python. Compare with results of Mathematica for same data sets: see pythontest.nb. In [1]: import scipy as sp from scipy.optimize import
More informationRandom Numbers Random Walk
Random Numbers Random Walk Computational Physics Random Numbers Random Walk Outline Random Systems Random Numbers Monte Carlo Integration Example Random Walk Exercise 7 Introduction Random Systems Deterministic
More informationProgramming for Engineers in Python
Programming for Engineers in Python Autumn 2016-17 Lecture 11: NumPy & SciPy Introduction, Plotting and Data Analysis 1 Today s Plan Introduction to NumPy & SciPy Plotting Data Analysis 2 NumPy and SciPy
More informationVideo Compression MPEG-4. Market s requirements for Video compression standard
Video Compression MPEG-4 Catania 10/04/2008 Arcangelo Bruna Market s requirements for Video compression standard Application s dependent Set Top Boxes (High bit rate) Digital Still Cameras (High / mid
More informationHW2. October 24, CSE 252A Computer Vision I Fall Assignment 2
HW2 October 24, 2018 1 CSE 252A Computer Vision I Fall 2018 - Assignment 2 1.0.1 Instructor: David Kriegman 1.0.2 Assignment Published On: Wednesday, October 24, 2018 1.0.3 Due On: Wednesday, November
More informationStatistical Data Analysis: Python Tutorial
1 October 4, 2017 Statistical Data Analysis: Python Tutorial Dr A. J. Bevan, Contents 1 Getting started 1 2 Basic calculations 2 3 More advanced calculations 4 4 Data sets 5 4.1 CSV file input.............................................
More informationMAS212 Scientific Computing and Simulation
MAS212 Scientific Computing and Simulation Dr. Sam Dolan School of Mathematics and Statistics, University of Sheffield Autumn 2017 http://sam-dolan.staff.shef.ac.uk/mas212/ G18 Hicks Building s.dolan@sheffield.ac.uk
More informationInverse Ray Shooting Tutorial. Jorge Jiménez Vicente Dpto. Física Teórica y del Cosmos Universidad de Granada Spain
Inverse Ray Shooting Tutorial Jorge Jiménez Vicente Dpto. Física Teórica y del Cosmos Universidad de Granada Spain Final goal Session I Introduction to Python Solving the lens equation Ray shooting basics
More informationIntroduction ti to JPEG
Introduction ti to JPEG JPEG: Joint Photographic Expert Group work under 3 standards: ISO, CCITT, IEC Purpose: image compression Compression accuracy Works on full-color or gray-scale image Color Grayscale
More informationMatplotlib Python Plotting
Matplotlib Python Plotting 1 / 6 2 / 6 3 / 6 Matplotlib Python Plotting Matplotlib is a Python 2D plotting library which produces publication quality figures in a variety of hardcopy formats and interactive
More informationUNIVERSITETET I OSLO
UNIVERSITETET I OSLO Det matematisk-naturvitenskapelige fakultet Examination in: IN1900 Introduction to programming with scientific applications Day of examination: Tuesday, October 10, 2017 Examination
More informationConvolutional Neural Networks (CNN)
Convolutional Neural Networks (CNN) By Prof. Seungchul Lee Industrial AI Lab http://isystems.unist.ac.kr/ POSTECH Table of Contents I. 1. Convolution on Image I. 1.1. Convolution in 1D II. 1.2. Convolution
More informationProf. Dr. Rudolf Mathar, Dr. Arash Behboodi, Emilio Balda. Exercise 5. Friday, December 22, 2017
Fundamentals of Big Data Analytics Prof. Dr. Rudolf Mathar, Dr. Arash Behboodi, Emilio Balda Exercise 5 Friday, December 22, 2017 Problem 1. Discriminant Analysis for MNIST dataset (PyTorch) In this script,
More informationFinal Exam, Version 3 CSci 127: Introduction to Computer Science Hunter College, City University of New York
Final Exam, Version 3 CSci 127: Introduction to Computer Science Hunter College, City University of New York 22 May 2018 1. (a) What will the following Python code print: i. a = "one+two+three+four+five+six"
More informationIntroductory Scientific Computing with Python
Introductory Scientific Computing with Python More plotting, lists and FOSSEE Department of Aerospace Engineering IIT Bombay SciPy India, 2015 December, 2015 FOSSEE (FOSSEE IITB) Interactive Plotting 1
More informationInteractive Mode Python Pylab
Short Python Intro Gerald Schuller, Nov. 2016 Python can be very similar to Matlab, very easy to learn if you already know Matlab, it is Open Source (unlike Matlab), it is easy to install, and unlike Matlab
More informationPyPlot. The plotting library must be imported, and we will assume in these examples an import statement similar to those for numpy and math as
Geog 271 Geographic Data Analysis Fall 2017 PyPlot Graphicscanbeproducedin Pythonviaavarietyofpackages. We willuseapythonplotting package that is part of MatPlotLib, for which documentation can be found
More informationFacial Recognition Using Eigenfaces
Lab 11 Facial Recognition Using Eigenfaces Load the Data Lab Objective: Use the singular value decomposition to implement a simple facial recognition system. Suppose we have a large database containing
More informationLECTURE 22. Numerical and Scientific Computing Part 2
LECTURE 22 Numerical and Scientific Computing Part 2 MATPLOTLIB We re going to continue our discussion of scientific computing with matplotlib. Matplotlib is an incredibly powerful (and beautiful!) 2-D
More information4. BASIC PLOTTING. JHU Physics & Astronomy Python Workshop Lecturer: Mubdi Rahman
4. BASIC PLOTTING JHU Physics & Astronomy Python Workshop 2016 Lecturer: Mubdi Rahman INTRODUCING MATPLOTLIB! Very powerful plotting package. The Docs: http://matplotlib.org/api/pyplot_api.html GETTING
More informationNAVIGATING UNIX. Other useful commands, with more extensive documentation, are
1 NAVIGATING UNIX Most scientific computing is done on a Unix based system, whether a Linux distribution such as Ubuntu, or OSX on a Mac. The terminal is the application that you will use to talk to the
More informationVisualisation in python (with Matplotlib)
Visualisation in python (with Matplotlib) Thanks to all contributors: Ag Stephens, Stephen Pascoe. Introducing Matplotlib Matplotlib is a python 2D plotting library which produces publication quality figures
More information5 File I/O, Plotting with Matplotlib
5 File I/O, Plotting with Matplotlib Bálint Aradi Course: Scientific Programming / Wissenchaftliches Programmieren (Python) Installing some SciPy stack components We will need several Scipy components
More informationDigital Image Representation Image Compression
Digital Image Representation Image Compression 1 Image Representation Standards Need for compression Compression types Lossless compression Lossy compression Image Compression Basics Redundancy/redundancy
More information15 Data Compression 2014/9/21. Objectives After studying this chapter, the student should be able to: 15-1 LOSSLESS COMPRESSION
15 Data Compression Data compression implies sending or storing a smaller number of bits. Although many methods are used for this purpose, in general these methods can be divided into two broad categories:
More informationJune 10, 2014 Scientific computing in practice Aalto University
Jussi Enkovaara import sys, os try: from Bio.PDB import PDBParser biopython_installed = True except ImportError: biopython_installed = False Exercises for Python in Scientific Computing June 10, 2014 Scientific
More informationmpl Release latest May 17, 2017
mpl a nimationmanagerdocumentation Release latest May 17, 2017 Contents 1 NOTE: Documentation is curently in development!!! 1 1.1 Matplotlib animation manager (GUI) 1.0a1...............................
More informationDerek Bridge School of Computer Science and Information Technology University College Cork
CS468: Artificial Intelligence I Ordinary Least Squares Regression Derek Bridge School of Computer Science and Information Technology University College Cork Initialization In [4]: %load_ext autoreload
More informationOpenMSI Arrayed Analysis Toolkit: Analyzing spatially defined samples in mass spectrometry imaging
OpenMSI Arrayed Analysis Toolkit: Analyzing spatially defined samples in mass spectrometry imaging Introduction This is the accompanying notebook of the manuscript OpenMSI Arrayed Analysis Toolkit: Analyzing
More information2. For each of the regular expressions, give a string that will matches it:
CMP 464-C401 Sample Final Exam, Spring 2016 1. What will the following code draw: n = 10 X = np.arange(n) Y1 = X/2.0 Y2 = X/4.0 plt.bar(x, +Y1, facecolor= blue ) plt.bar(x, -Y2, facecolor= red ) for x,
More informationAMath 483/583 Lecture 28 June 1, Notes: Notes: Python scripting for Fortran codes. Python scripting for Fortran codes.
AMath 483/583 Lecture 28 June 1, 2011 Today: Python plus Fortran Comments on quadtests.py for project Linear vs. log-log plots Visualization Friday: Animation: plots to movies Binary I/O Parallel IPython
More informationSPEAKER RECOGNITION. 1. Speech Signal
SPEAKER RECOGNITION Speaker Recognition is the problem of identifying a speaker from a recording of their speech. It is an important topic in Speech Signal Processing and has a variety of applications,
More informationAMS 27L LAB #2 Winter 2009
AMS 27L LAB #2 Winter 2009 Plots and Matrix Algebra in MATLAB Objectives: 1. To practice basic display methods 2. To learn how to program loops 3. To learn how to write m-files 1 Vectors Matlab handles
More informationPlotting With matplotlib
Lab Plotting With matplotlib and Mayavi Lab Objective: Introduce some of the basic plotting functions available in matplotlib and Mayavi. -D plotting with matplotlib The Python library matplotlib will
More informationData Science and Machine Learning Essentials
Data Science and Machine Learning Essentials Lab 3A Visualizing Data By Stephen Elston and Graeme Malcolm Overview In this lab, you will learn how to use R or Python to visualize data. If you intend to
More information