Physicist's algorithm simplifies biological imaging -- and also solves Sudoku puzzles

Mar 02, 2006
Sudoku puzzle

Cornell physicist Veit Elser has been engrossed recently in resolving a pivotal question in biological imaging. So he hasn't had much time for brainteasers and number games.

But in discovering an algorithm critical for X-ray diffraction microscopy, Elser and colleagues solved two problems. First, they gave researchers a new tool for imaging the tiniest and most delicate of biological specimens. And second, they discovered that the same algorithm also solves the internationally popular numbers puzzle Sudoku.

Not just one puzzle. All of them.

The Sudoku discovery appeals to Elser's whimsical side. But the algorithm, he notes, has potential for all kinds of other endeavors. "There are a lot of problems that you can represent in terms of this language," he says. "We're providing the technique. Whatever people use it for, it's great for us."

The so-called difference-map algorithm, which Elser says could have applications from productivity optimization to nanofabrication, tackles problems for which the solution must meet two independent constraints. In the case of Sudoku, the constraints are simple: Each of nine numbers, considered alone, appears nine times in the grid so that there is only one per row and column. And all nine numbers appear within each of the nine blocks.

In X-ray diffraction microscopy, the constraints are more complex. But the beauty of the algorithm, as Elser demonstrates, is that complexity doesn't matter. By applying the algorithm to the jumble of raw data from such an experiment, researchers can now reconstruct from it a clear, detailed image.

The result, shown in images published with a recent article in the Proceedings of the National Academy of Sciences (PNAS) with lead author David Shapiro of the State University of New York at Stony Brook and other colleagues, is a richly detailed image of a specimen as small as a single yeast cell -- taken without staining, sectioning or otherwise damaging the specimen.

Unlike optical microscopes, which use a lens to focus light on a target (and are therefore only useful for specimens larger than the wavelength of visible light), imaging methods based on X-rays and electrons take advantage of the finer wavelengths provided by these forms of illumination.

But some of these methods can damage the specimen with harmful radiation, require that a specimen be stained or otherwise altered or lack the penetrating power necessary for three-dimensional reconstructions. X-ray diffraction microscopy, which uses "soft" X-rays and measures the resulting diffraction pattern, is often the method of choice because it gives a detail-rich image and leaves the specimen relatively unscathed.

The tricky part comes in using the diffraction pattern to reconstruct an image of the specimen. Instead of directly measuring the pixel-by-pixel contrast within the specimen, researchers are left with data that represents the object broken down into its constituent waves: a vast number of them, of different frequencies and amplitudes, waiting to be added up in a process called a Fourier synthesis to reconstruct the image.

"But if it were just combining waves with a definite oscillation and adding them up, that would be a piece of cake for a computer to do," says Elser.

The challenge is in the waves' phases, which are critical in reconstructing an image from the diffraction data. Without attention to that piece of the puzzle, the resulting image is reduced to noise.

"People used to say that's an impossible problem," Elser says. "Then people got to thinking about the fact that there are going to be constraints coming from some rather mundane facts -- that in principle will make the problem of deducing the phase like the solution of a solvable puzzle."

The mundane fact, in this case, was a basic premise of the X-ray diffraction experiment: that the object in view have a clearly defined boundary (i.e., that all pixel values outside that boundary be set to zero).

"If I set up waves with known amplitudes and synthesize them, I'll find that for essentially all random combinations of phases, the thing I get is an object not confined," Elser says. "It takes a very clever combination of the phases of all those waves to add up to something in only one region; and cancel out everywhere else."

For the second constraint, the researchers required that the wave amplitudes used in the Fourier synthesis matched those measured by the experiment.

With the two constraints in place, the difference-map algorithm completes the job.

Source: Cornell University, By Lauren Gold

Explore further: Neutron tomography technique reveals phase fractions of crystalline materials in 3-dimensions

add to favorites email to friend print save as pdf

Related Stories

Huge waves measured for first time in Arctic Ocean

8 hours ago

As the climate warms and sea ice retreats, the North is changing. An ice-covered expanse now has a season of increasingly open water which is predicted to extend across the whole Arctic Ocean before the middle ...

Underwater elephants

8 hours ago

In the high-tech world of science, researchers sometimes need to get back to basics. UC Santa Barbara's Douglas McCauley did just that to study the impacts of the bumphead parrotfish (Bolbometopon muricatum) on cor ...

Recommended for you

Refocusing research into high-temperature superconductors

10 hours ago

Below a specific transition temperature superconductors transmit electrical current nearly loss-free. For the best of the so-called high-temperature superconductors, this temperature lies around -180 °C – a temperature ...

MRI for a quantum simulation

15 hours ago

Magnetic resonance imaging (MRI), which is the medical application of nuclear magnetic resonance spectroscopy, is a powerful diagnostic tool. MRI works by resonantly exciting hydrogen atoms and measuring ...

50-foot-wide Muon g-2 electromagnet installed at Fermilab

15 hours ago

One year ago, the 50-foot-wide Muon g-2 electromagnet arrived at the U.S. Department of Energy's Fermi National Accelerator Laboratory in Illinois after traveling 3,200 miles over land and sea from Long Island, ...

User comments : 0