number 09, September 1, 1996
![]()
a forum for virtual meetings
of scientists and students involved in a research activity on
CONTEMPORARY PHYSICS
Editors:
Copyright (C) 1996 by Zbigniew Koziol.
IN THIS ISSUE:
Letter from the Editor
Multifractal Scaling With Wavelet Packets, by Cameron L. Jones
Animated GIFS - A Simple Tutorial Example, by Zbigniew Koziol
"Lawrence and His Laboratory" - a history of physics site
![]()
LETTER FROM THE EDITOR
![]()
Multifractal Scaling With Wavelet Packets
Cameron L. Jones
Centre for Applied Colloid and BioColloid Science
Swinburne University of Technology, School of Chemical Sciences
Hawthorn, Victoria, 3122, Australia
CJONES@swin.edu.au
Abstract
| A multifractal Wavelet Packet approach is described for characterising one dimensional signals using the Hurst exponent. The Hurst exponent quantifies scaling behaviour between data points, and classifies these trends as persistent, antipersistent, or random walk. A multifractal approach provides additional information regarding local scaling behaviour. This algorithm uses commercially available software to perform the Wavelet Packet Transform, with additional manipulation carried out using an iterative spreadsheet function to determine the regression slope in lagged windows. The method is demonstrated in a biophysical reaction- diffusion system where image-analysis was used to collect raw data concerning the spatial distribution of local chemical (enzyme) concentrations. These one dimensional signals (Enzyme Walks) display antipersistent scaling at small lags, but random walk and persistent scaling at longer lags. This emphasises the importance of local coordination during enzyme secretion, which breaks down at longer lags as a function of the advancing diffusion front. |
Introduction
This paper uses WPA to estimate the Hurst exponent. In addition, we demonstrate the method by examining a complex biophysical system which displays fractal-like scaling trends during enzyme secretion1. This gives rise to one dimensional signals called 'Enzyme Walks', which can be characterised using time series tools. Enzyme Walks are an example of a diffusion controlled spatial deposition process. From a practical viewpoint, we also extend the global method to estimate multifractal, or local scaling trends. This therefore extends the utility of this method to heterogeneous systems which display complex behaviours at multiple length or time scales. It is hoped that this method will find widespread application for the rapid analysis of complex, multifractal type signals encountered in physics, chemistry, biology or in other signal analysis situations.
Global versus Local Scaling
Self-affine Scaling
Biological Background
|
| Figure 1. Three day old colony of Pycnoporus cinnabarinus grown on a microporous membrane overlaid onto Malt Extract Agar, and stained to reveal the exo-enzyme laccase. |
| Figure 2. Single digital line profile overlaid on the image [ALL]. |
| Figure 3. Pseodocolour representation to reveal the radial double banding pattern of enzyme activity concentration. The arrow indicates the centre-point inoculum, and the letter a and b indicate the high laccase expression region. Scale bar == 10mm. |
Numerical Approach
|
| Figure 4. (a) Single laccase enzyme walk for ALL. (b) Frequency partitioning for the Wavelet Packet Transform for ALL shown in (a), plotting the best basis (x-axis) versus scale level (y-axis). (c) Enzyme walk for HALF. (d) Wavelet Packet Transform frequency partition with the best basis coefficients identified as black blocks. |
To examine sensitivity to resolution scale changes, one can double the magnification and take a line profile from the centre-point to the right hand edge [HALF]. Figure 4c and 4d shows the respective 1-D profile and frequency partition with the d.4 wavelet packet.
Computer Software and Methodology
Mathematical Notation
![]()
(1)
The slope,![]()
(2)
Note that in 2-D, the following relationship holds between the Fractal Dimension, D and the Hurst exponent, H: D=2-H.Multifractal Data Manipulation
| Figure 5. (a) Spreadsheet calculation of the nth data set for specified lags. (b) x versus y plotting scheme between coefficient energy magnitude and index position for different lags. |
Results and Discussion
| Figure 6. Sample plot for HALF showing the 12 regressions for each lag position. |
By taking the average slope for all data sets for the two treatments, it was possible to obtain mean statistics for a multifractal description of local enzyme activity concentration. This provides an overview of the scaling behaviour between nearest neighbour regions separated by a defined lag position (measured in pixel units). This data is plotted for both ALL and HALF in Figure 7.
|
| Figure 7. Mean values for Hlocal for different lagged positions at each magnification scale (ALL or HALF). |
Biological Implications
Conclusions
References
- Jones, C.L., Lonergan, G.T. and Mainwaring, D.E. (1996). Wavelet Packet Computation of the Hurst Exponent. J. Phys. A: Math. Gen. 29: 2509-2527.
- Peters, E.E. (1991). Chaos and Order in the Capital Markets - A New View of Cycles, Prices, and Market Volatility. John Wiley & Sons, New York. pp. 61-80.
- Jones, C.L. (1996). Fractal Dimension Estimation with Wavelet packets. Virtual Physics. August 1 ,No. 7. URL: http://www.isisnet.com/MAX/science/physics/vp/vp07.html
- Lonergan, G.T., Panow, A., Jones, C.L., Schliephake, K., Ni, C.J. and Mainwaring, D.E. (1995). Physiological and Biodegradative Behaviour of the White-rot Fungus, Pycnoporus cinnabarinus in a 200 litre Packed Bed Bioreactor. Australasian Biotechnology. 5(2), March-April: 107-111.
- Field, J.A., de Jong, E., Feijoo-Costa, G. and de Bont, J.A.M. (1993). Screening for Ligninolytic Fungi Applicable to the Biodegradation of Xenobiotics. Tibtech. 11: 44-49.
- Jurasek, L. and Paice, M.G. (1988). Biological Treatment of Pulps. Biomass. 15: 103-108.
- Coifman, R.R. and Wickerhauser, M.V. (1993). Wavelets and Adapted Waveform Analysis - A Toolkit for Signal Processing and Numerical Analysis. [Computer program and Book] -Wavelet Packet Laboratory for Windows (Digital Diagnostics Corporation - Yale University; Email to: dcc@cs.yale.edu). A. K. Peters, Massachusetts.
- Hastings, H.M and Sugihara, G. (1993). Fractals - A User's Guide for the Natural Sciences. Oxford University Press, New York.
- Jones, C.L. (1996) 2-D Wavelet Packet Analysis of Structural Self-Organization and Morphogenic Regulation in Filamentous Fungal Colonies. Complex Systems Conference - "From Local Interactions to Global Phenomena" July 14-17, 1996 Charles Sturt University; Albury - Wodonga, N.S.W, Australia. Also in: Complex Systems 96 - From Local Interactions to Global Phenomena. (Stocker, R., Jelinek, H., Durnota, B. and Bossomaier, T. eds.) IOS Press, Amsterdam, 1996. pp. 12-23. Available online: http://www.swin.edu.au/chem/bio/cs96/camjones.htm
- Jones, C.L. (1996) Wavelet Packet Fractal Analysis - Software Operating Instructions. [A series of 7 software functions which calculate the mean or global Fractal Dimension of 2-D objects in digital images. Requirements include: S-Plus and S+Wavelets for Windows. Programs available for the following image sizes: 64x64, 128x128, 192x192, 256x256, 384x384, 512x512, and 300x200.] Available online: http://www.swin.edu.au/chem/bio/s+code/wpafrac1.htm
![]()
Animated GIFS - a simple tutorial example
Zbigniew Koziol, WebEx@ra.isisnet.com
CONTENTS
- Nonlinear Diffusion - a One-Dimensional Case
- Animated demonstration
- Numerical Modelling of Nonlinear Diffusion
- References
- EVERYTHING ABOUT GIF89a-based Animation for the WWW
- SELECTED TOOLS
- GIFTOOL
- GIF Construction Set for Windows
- MORE EXAMPLES
![]()
![]()
Nonlinear Diffusion - a One-Dimensional Case
Animated demonstration
![[image: a nonlinear diffusion equation ]](beta.gif)
Here, the following cases can be viewed, corresponding to various time-dependences of the external magnetic field, H(t), at the surface of a flat superconductor (see also the figures):
- MODE 0: External field H is a step function of time
- MODE 1: External field has a triangle-dependence on time
- MODE 2: External field is a periodic function of time
Magnetic field penetration into an infinite slab of a superconductor,
for kappa = 0 (classical, linear diffusion limit)
![]() |
![]() |
time from 0 to 2500 images taken every 100 steps of iteration |
![]() |
![]() |
time from 0 to 2500 images taken every 100 steps of iteration |
![]() |
![]() |
time from 1000 to 2000 images taken every 20 steps of iteration |
Numerical Modelling of Nonlinear Diffusion
A listing of the essential part of a Pascal program to perform numerical modelling of the nonlinear diffusion equation. This is an example only, not a complete source code of a ready-to-use program.
CONST maxSTEP=5000; { number of loops, 50 iterations each }
kappa=3; { parameter determining nonlinearity }
PERIOD = 25000 { time-period of external AC field }
a = 0.5; { 'a' must be less or equal than about}
{ 0.5, to have calculations stable }
VAR f: array[0..200,0..50] of real; xs, ts, step : integer;
BEGIN { magnetic field =0 inside of the sample }
for xs:=1 to 199 do f[xs,0]:=0;
step:=0;
REPEAT
FOR ts:=0 to 50 do
BEGIN { field = sin(2*(*time/PERIOD) on both surfaces }
{ any other function could be placed there }
f[0,ts]:= sin(2*3.14159*(50*step+ts)/PERIOD);
f[200,ts]:=f[0,ts];
END;
FOR ts:=0 to 49 do
FOR xs:=1 to 199 DO
f[xs,ts+1]:=f[xs,ts] +a *( pwr(f[xs+1,ts]-f[xs,ts],kappa+1)
-pwr(f[xs,ts] - f[xs-1,ts], kappa+1) );
{The function pwr(x,z) is x to power z;
this function is defined elsewhere }
{Here, the data can be stored to the disk. Now,the magnetic}
{field distribution at t=50 is assumed to be the starting }
{condition for the next 50 steps of iteration }
FOR xs:=0 to 200 do f[xs,0]:=f[xs,50];
inc(step);
UNTIL step=maxSTEP;
END.
References
- Z. Koziol, Ph.D. Thesis, University of Amsterdam, 1994.
- Z. Koziol and P.F. de Chatel, IEEE Trans. Magn. 30 (1994) 1169.
- Nonlinear diffusion in superconductors
- DIFUSION.EXE , A WINDOWS 3.1 DEMO which should be used together with a WINDOWS HELP file DYFUZJA.HLP .
![]()
EVERYTHING ABOUT GIF89a-based Animation for the WWW
The site http://members.aol.com/royalef/gifanim.htm, maintained by Royal E. Frazier is the largest and the best source of information about all the aspects related to the implementation of animated gifs on www pages. Here is a few short fragments from this page.- What is GIF89a. Well, its the technical rules from 1989 that explain how GIFs can be put together. You see most GIFs over the years have only one image per file. According to the technical specifications from 1987, a GIF could have had more than one image per file, making it like a slide show presentation and not a single image. However, most programs that work with GIF are designed around the idea of one image per file. So the multi-image aspect of GIFs was forgotten. In 1989, they added timing and various other abilities to the GIF format, including transparency. Nobody used these new additions either. Then the Web took off. Transparency and interlacing became features people wanted to use and software companies began supporting those features.
- [ The site ] contains a tutorial of over 35 printed pages, 80K of readable text, animated illustrations, 1 megabyte of images and data, bugs reports and more. Please check back at the Overview and look to the table of contents for more info.
![]()
SELECTED TOOLS
GIFTOOL
It is available from http://www.homepages.com/tools/giftool/.-
This tool does a variety of things, add/remove comments and interlace
GIF images. It can do batch convertion of large sets of images.
- Binaries are available for all the major platforms known: Sun Sparc Solaris 2.X (SunOS 5.X) and 1.X (SunOS 4.X), Dec Alpha OSF 2.0, Hewlett Packard HPUX 9.05, IBM RS6000, SGI Irix 4.0.5 and 5.2, x86 Linux, PC MS-DOS. The source code is also available.
A source code of a perl script for animation is also available at the same site.
-
This tool does a variety of things, add/remove comments and interlace
GIF images. It can do batch convertion of large sets of images.
GIF Construction Set for Windows
It is available from http://www.mindworkshop.com/alchemy/gifcon.html, a site of Alchemy Mindworks- GIF Construction Set for Windows is a powerful collection of tools to work with multiple-block GIF files. It will allow you to assemble GIF files containing image blocks, plain text blocks, comment blocks and control blocks. It includes facilities to manage palettes and merge multiple GIF files together.
![]()
MORE EXAMPLES
- The 1st Internet Gallery of GIF Animation
- Victorian Animated Gifs
- Animated .Gifs, Egor and Java Applets
- Music from Poland - Muzyczna Jedynka
![]()
"Lawrence and His Laboratory" - a history of physics site
Received: August 29, 1996
http://www.lbl.gov/Science-Articles/Research-Review/Magazine/1981/
chronicles Berkeley Lab's "Lawrence Years," beginning with the Ernest Lawrence's founding of the Lab in 1931. It features hundreds of scientific images from the Lab's historical photo archive. Images include:
-
- photos of the first cyclotron and later particle-accelerating machines
- - handwritten pages from lab notebooks describing Nobel Prize-winning discoveries
- - engineer's sketches for devices that made possible probing the secrets of the atomic nucleus
- - handwritten pages from lab notebooks describing Nobel Prize-winning discoveries
If you have further questions, please contact Jeffery Kahn at jbkahn@lbl.gov or (510) 486-4019.
Berkeley Lab is a U.S. Department of Energy national laboratory located in Berkeley, California. It conducts unclassified research and is managed by the University of California.
Editors:
Virtual Physics URL addresses:
Université de Liège, Sart Tilman, B-4000 Liège, Belgium, tel. (+32 41) 66 37 52
Umeå University, S-907 42 Umeå, Sweden, tel. +46-(0)90-167717
Swinburne University of Technology, P.O. Box 218 Hawthorn, Victoria, 3122 Australia, tel. +613 9214 8935, fax +613 9819 0834
2-6032 Compton Ave., Halifax, Nova Scotia, B3H 1E7 Canada, tel. (902) 423 2149
Warsaw University, Hoza 69, 00-681 Warsaw, Poland, tel. (+48 2) 628 3031
University of Pennsylvania, Rm. C-501 Richards Bldg., Philadelphia, PA 19104-6089, U.S.A., tel. (215) 898-6396
To subscribe a F R E E e-mail
version or submit materials for publication, write to Zbigniew Koziol.
Copyright (C) 1996 by Zbigniew Koziol.
this copyright notice concerns the whole of the Virtual Physics
edition but not specific articles published there which are
property of their respective copyright owners![]()
![[Virtual Physics]](../images/vp.gif)
![[ animation for kappa=0, mode No 0 ]](k0m0.gif)
![[.]](m0.gif)
![[ animation for kappa=0, mode No 1 ]](k0m1.gif)
![[.]](m1.gif)
![[ animation for kappa=0, mode No 2 ]](k0m2.gif)
![[.]](m2.gif)
![[ WebExperts ]](../images/we5.gif)