Author Archives: lifrye

Conclusion: Diffraction Patterns of C. elegans

I set out to model the diffraction patterns created by C. elegans nematodes using mathematica, specifically using the diversely useful mathematical tool of the Discrete Fourier Transform. The Transform is a quick way to find the diffraction patterns from Fraunhofer diffraction and interference (far-field diffraction). In simpler language, $\left| FT^{2} \right|$ produces the diffraction pattern, which is an analysis of the Electric field strength across the aperture, in this case, the image. The goal of this project was to produce a “library” of shapes and their corresponding diffraction patterns. However, there were many obstacles that made it necessary for me to do quite a bit of relevant research.

When I started, I did not know that the type of diffraction I was studying was Fraunhofer. Diffraction patterns are a result of the type of wave and the type of experimental setup. In this case, the light source, the aperture, and the screen were are far enough away to classify it as far-field. Mathematically, $L>>\frac{b^{2}}{\lambda}=\frac{area. of. aperture}{\lambda}$ . The fact that this is an example of Fraunhofer Diffraction makes it possible to apply the FT (on Mathematica) to yield the diffraction pattern.

I came in with some preliminary knowledge of what Fourier Transforms were in theory. In general, when applied, it changes a function’s variable dependence: $F(t) \leftrightarrow \Phi (v)$. I did quite a bit of research on the derivations of the transform and understood the matrix calculations necessary to do the transform “by hand” (see Project Plan).

To produce good-quality images, it was also necessary to do quite a bit of exploring in the Mathematica Documentation Center, becoming familiar with a variety of image manipulation commands. Those included: sharpening, brightening, finding the pixel count and dimensions, cropping, changing the color scheme to grayscale, partitioning and reassembling images manually based on the pixel dimensions, and more, depending on what degree of manipulation the image needed.

All in all, the images produced are not only pretty, but informative and now available for reference. Because it is mathematically impossible to go from the diffraction image to the shape, it is necessary to have some type of library, like the one I created, to go in that direction. (It is impossible to do it mathematically because when the Fourier Transform is applied, the phase information of the light is lost as the data is converted to the complex space.) Therefore, my “library” is useful because I have made it possible to guess the approximate shape and orientation of the worm from the diffraction pattern for some basic worm shapes. Another application for this analysis is that a compilation of several consecutive diffraction patterns shows the thrashing frequency of the worm. For example, if you took several pictures within a few seconds and found the diffraction patterns, the amount that the patterns change corresponds to the “thrashing frequency,” or the quantifiable amount that the worm wiggles.

If I were to continue this project, I would continue to build upon the library in the same manner, testing to see if different input image qualities would make a difference in the final product (if different pieces of the worms were fluorescing, if the worm was glowing a different color, etc). Finally, I would construct a real library (instead of this evolving blog site) that was easy to navigate, and generally more accessible for research.

Share

Final Data: Diffraction Patterns of C. elegans (revised)

Listed below are five sets of data (3-7), two images each. On the left is the photograph of the fluorescing worm. To the right is its corresponding diffraction pattern, produced with Mathematica.

First, a few comments:

Notice that each worm shape contains a different symmetry (or lack of symmetry). However, every diffraction pattern is symmetrical. The two images tend to share some symmetrical qualities, as noted below each image.

These diffraction patterns are due to the fluorescing parts of the worm (GFP, or green fluorescing protein). I believe that because of my image manipulation (converting to grayscale and brightening the image significantly), the diffraction pattern is approximately due to the entire worm and not just the fluorescing parts.

On that note, I’d like to discuss the parameters I used for image manipulation and Fourier Transforming.

  • The Fourier Parameters I chose were {0, 1} . These are the default parameters when using the FT command on Mathematica. I chose them among the three most common by ruling out the other two: {-1, 1} is used for data analysis (which I was not doing), and {1, -1} is used for signal processing (also irrelevant to my project). My choice of parameters specified the conventions I wanted the program to use when applying the transform.
  • The command ImageAdjust [image, {a, b}] adjusts the contrast of the image by a and the brightness by b. For most of my images, I adjusted the contrast by the same amount, but changed the brightness differently for each photo depending on my judgement.
  • I divided the images by hand and reassembled them according to the image dimensions, rectangles of dimension {col/2, row/2} (where col = the number of pixels in each column of the image, and row = the same for each row). I could have use the command RotateRight, but I did not figure that out until after I developed this equally effective algorithm.

 

IMAGE 3

Screen Shot 2014-04-24 at 12.22.33 PM Screen Shot 2014-04-24 at 12.50.08 PM   

This image has axial symmetry perpendicular to the axial symmetry of the worm. The worm is mostly linear, and so is the diffraction.

(complete file for image 3: book 3 )

 

IMAGE 4

3975048733_8abd2dce44_z Screen Shot 2014-04-24 at 1.11.07 PM

See comment for image 3.

(complete file for image 4: book 4 )

 

IMAGE 5

Screen Shot 2014-04-24 at 1.19.58 PM Screen Shot 2014-04-24 at 1.22.16 PM

Although this does not have an “inversion center” (the worm does not complete a loop), it is significantly less straight than the first few examples. Correspondingly, the diffraction pattern has an axial symmetry perpendicular to the “axis” of the worm. However, it splays out to a higher degree, indicating that the area of the pattern corresponds with the curvature of the worm.

(full file for image 5: book 5 )

 

IMAGE 6

Screen Shot 2014-04-24 at 4.10.28 PM Screen Shot 2014-04-25 at 8.51.20 PM

This further backs up my hypothesis from image 5. The curvature of the worm is great and the area of the pattern is also large. However, this worm does have a sort of “inversion center.” The diffraction image has an eye in the center, which is different than all of the preceding images. The curvature of the worm greatly effects the produced image, which is not surprising, but the overall shape of the image changes dramatically. With more circular worms, an “eye” is formed in the center of the diffraction pattern, and the pattern here has both some rotational symmetry and axial symmetry. With straighter worms, the diffraction is much more linear, containing only axial symmetry.

(full file for image 6: book 6 )

 

IMAGE 7

Screen Shot 2014-04-26 at 1.43.14 PM Screen Shot 2014-04-26 at 1.45.09 PM

See comment for image 6.

(full file for image 7: book 7 )

Share

Preliminary Results: C. elegans Diffraction Pattern Modeling

Staying true (so far) to my tentative project timeline, I acquired images of the C. elegans in various shapes, I have done quite a bit of research on Fourier Transforms and Fraunhofer Diffraction, and so far I have successfully transformed one image into the corresponding diffraction pattern.

 

IMAGE 1

sampleworm1I took this image and used Screen Shot 2014-04-21 at 5.16.08 PMmathematica to sharpen it –> setting it to grayscale and brighten it –> collect dimensional information –> apply a Fourier Transform, yielding (after some similar image manipulation):

This is a great diffraction pattern, but I had issues with the poor resolution and general image quality. To remedy this, I proceeded with images taken with a higher-resolution camera.

(full file for image 1: book 1 )

IMAGE 2

3975815904_e33a49b65b_z

This beautiful image needed some manipulation, similar to image 1: I converted it to grayscale –> brightened it significantly (to make it a more definite shape, and to get rid of the “holes” in the luminescing nematode)  –> collected image dimensions and data –> applied the FT. Unfortunately, I ran into a problem.

The produced image:

Screen Shot 2014-04-21 at 5.31.56 PM

Obviously this is quite different than the first diffraction image.

I had a few hypotheses:

1. The image was saved as a .jpg, but the same image was produced when I tried again with a .png version of the image.

2. The computer is phase- shifting the image so that instead of the origin lying in the center of the product, it is splitting the right side from the left side and lining them up in the wrong order. How can I rearrange and correct the phase shift in the output?

An analogy to the second hypothesis:

Screen Shot 2014-04-21 at 5.45.39 PM

Screen Shot 2014-04-21 at 5.45.53 PM

 

—>

 

 

It is as if, instead of centering the origin in the center of the produced diffraction pattern, the computer is putting the “origin” in a different place, and splitting the image, similar to the parabola I produced above.

My solution is a little underhanded. I divided the image into four equal rectangles, and manually rearranged them to produce what I knew was the true image:Screen Shot 2014-04-22 at 7.52.00 PM

(full file for image 2: book 2a  book 2b )

—————-

Reflection:

It is important to keep this process grounded: how is this relevant to Electromagnetism? The answer is that this entire process is only viable because of the laws of electromagnetism. I am analyzing the images by taking their Fourier Transforms. The diffraction pattern is the FT of the function that describes the electric field strength across the aperture of diffraction. In other words, I am applying an operation (the FT) to the image, which is a direct indication of the electric field strength across the aperture (the microscope slide) to mathematically find the diffraction pattern produced by the specific electric field array created by the shape of the worm.

Specifically, the diffraction pattern here is the Fraunhofer Diffraction pattern, or “far-field” diffraction, which occurs when the distances between the screen, aperture, and light source are appropriately far $L>>\frac{b^{2}}{\lambda}$. Diffraction effects are an outcome of the type of light wave.

It is also essential to realize what information is lost in the computation of these diffraction patterns. I am taking a real image, applying a FT to it, squaring the absolute value of the result, and arriving in a complex space. This process loses the phase information of the light, and as a result, it is possible to go from the image to the diffraction pattern, but impossible to find the image from the FT diffraction image.

Share

Project Plan (revised)

C. elegans Diffraction Pattern Modeling

Sources & Resources

  1. Introduction to Electrodynamics, D. Griffiths
  2. A Student’s Guide to Fourier Transforms, J. F. James (chapters 1 and 9)
  3. Past work by Professor J. Magnes and her assistants
  4. C. elegans nematodes, property of Professor K. Susman, photos taken by her laboratory assistants
  5. Introduction to Optics, 3rd edition, Pedrotti & Pedrotti & Pedrotti (chapter 11)

Necessary steps & what I plan to model

  • Take photos of the C. elegans with the Insight Camera (directly from the microscope)
  • Manipulate the photos (remove noise) using Mathematica
  • Take the Fourier transform of the photos to discover the diffraction pattern
  • ( |Fourier Transform|2 = the diffraction pattern)

Collaborators

I will be working alone, but with Professor K. Susman’s worms and equipment, and Professor J. Magnes’ old work as a reference. I also will be receiving some training in using the microscope’s camera attachment (most likely from one of Professor Susman’s lab assistants).

Tentative Timeline

Starting with week1= April 6-12 (Sunday to Saturday):

Week1: I plan to take photos (if I get permission from Professor Susman), and work on an algorithm for generating reliably noiseless photos. I plan to become comfortable with Fourier transforms on Mathematica, and take some preliminary transformation models.

Week2: From now on, it is going to progress in a semi-continuous fashion. This week I plan to continue working with the Fourier transforms on Mathematica, making improvements to the algorithm and recording the changes. Keeping in mind the goal of this project, it is important to take note of the models in the context of electromagnetic waves. It is also important to keep careful track of the data, and make sure the images are clear and easy to see.

Week3: By now I expect to have some good images. I will continue to improve on them, and possibly will begin resorting to doing the matrix multiplication “by hand” on Mathematica. I hope to have a very good understanding of Fourier transforms of images by the end of this project.

Week4: Now that I hopefully have good images corresponding to several shapes of C. elegans, I hopefully will be able to begin adding at least a small volume of genuinely good data to the Diffraction Symmetries Library.

Week5: In the final week, I will make some final touches to the library. Even though it is primarily about the data, I hope to polish the presentation and add any necessary comments. 

Some Preliminary Information

It is important to discuss what I already have knowledge of: I already know how to grow worms (how to transfer them to new E. coli food dishes so they can reproduce). I have a general idea how to take Fourier Transforms by hand, and what they generate. I do not know how to take a Discrete Fourier Transform of an image, and I do not know how to do it with matrix multiplication. I also know very little about electromagnetic Fraunhofer diffraction.

Some Essential Relationships

The crux of my thesis:

(1)   \begin{equation*}\left|FT|\right^2=Diffraction Pattern\end{equation*}

Basic Fourier Transform equation:

(2)   \begin{equation*} F(t)= \int_{-\infty}^{\infty}  \Phi (v) {e }^{2 \pi i v t } dv \end{equation*}

(3)   \begin{equation*} F(t) \leftrightarrow \Phi (v) \end{equation*}

Discrete Fourier Transform in matrix form:

(4)   \begin{equation*}\begin{bmatrix} A(0)\\ A(1)\\ A(2)\\ ...\\ A(N-1) \end{bmatrix} = \begin{bmatrix} 1&1&1&...&1\\ 1&e^(\frac{2 \pi i}{n})&e^(\frac{4 \pi i}{n})&...&e^(2(N-1)\frac{\pi i}{n})\\ 1&e^(\frac{4 \pi i}{n}) &e^(\frac{8 \pi i}{n}) &...& e^(4(N-1)\frac{\pi i}{n})\\ ... & ... & ... & .... &...\\ 1 & ... & ... &...& e^((N-1)^2(\frac{\pi i}{n})) \end{bmatrix} \begin{bmatrix} a(0)\\ a(1)\\ a(2)\\ ...\\ a(N-1) \end{bmatrix} \end{equation*}

(It can be noted that this matrix multiplication requires $N^2$ multiplications. The Fast Fourier Transform (FFT) method, which can only be run by computing machines, reduces the number of multiplications from $N^2$ to $2Nlog_2(N)$.

Share

Project Proposal: Diffraction Symmetries of C. elegans

The C. elegans nematode is a common subject of biological studies, and has become more and more popular in physics research. I intend to find the diffraction patterns generated by the worms’ shape and log the findings into a Symmetries Library (with the eventual goal of using Group Theory to get the worm shape directly from a diffraction image).

The shape of the worm (photos to be taken with a microscope) will correspond to a particular diffraction pattern. I will model the Fraunhoufer diffraction patterns (Far-Field diffraction) of the electromagnetic waves (light waves) by generating images with Mathematica using the Fourier Transforms. The idea is that $\left | Fourier Transform | \right ^2 $ = the diffraction pattern. This project is a study of the behaviors of light waves.

I will eventually be keeping a log of my findings on the already existing website, the Diffraction Symmetries Library.

Share