Method and apparatus for compressing Z-stack microscopy images
A method and apparatus for compressing Z-stack microscopy images comprising compressing an all-focus image and depth map representation of an original Z-stack of images, computing a prediction residual as a difference between a simulated Z-stack, generated from the all-focus image and the depth-map, and the original Z-stack and concatenating the prediction residual to the compressed all-focus image and the depth map to produce a final compressed file.
Latest Sony Corporation Patents:
- POROUS CARBON MATERIAL COMPOSITES AND THEIR PRODUCTION PROCESS, ADSORBENTS, COSMETICS, PURIFICATION AGENTS, AND COMPOSITE PHOTOCATALYST MATERIALS
- POSITIONING APPARATUS, POSITIONING METHOD, AND PROGRAM
- Electronic device and method for spatial synchronization of videos
- Surgical support system, data processing apparatus and method
- Information processing apparatus for responding to finger and hand operation inputs
1. Field of the Invention
Embodiments of the present invention generally relate to digital microscopy and, more specifically, to a method and apparatus for compressing Z-stack digital microscopy images.
2. Description of the Related Art
A digital microscope images a specimen at different depths to create a sequence of digital images. Each image represents a portion of the specimen at a particular depth of focus. Thus, at a certain depth, only a fraction of the entire specimen is in focus. The sequence of digital images is stacked along a Z dimension corresponding to depth, referred to as a Z-stack. Each image in a Z-stack is focused at a different depth in the captured specimen. The optics of digital microscopes impose a very narrow depth of field, too small to capture the entire specimen in focus with a single picture. Thus, all images in the Z-stack are of the same specimen, but with different focus for each image.
Z-stacks are generally large sized data sets with high resolution, which impose high bandwidth requirements for the storage or transmission of the Z-stack. For example, in several applications, the number of images is 15, 30, or 50, each with a spatial resolution of approximately 100,000 pixels by 100,000 pixels, which imposes high bandwidth requirements for the storage or transmission of Z-stack.
According to some conventional methods for compressing a Z-stack, the digital images in the Z-stack are compressed individually. While simple to implement, such methods for compressing a Z-stack result in large data files. Due to the large size of data files, there are delays in viewing the Z-stack over a bandwidth -constrained (or bandwidth limited) channel. Video compression techniques are also utilized in the field, however video compression delays viewing of the data due to frame-by-frame representation of the Z-stack images.
Thus, there is a need for method and apparatus for compressing a Z-stack of microscopy images.
SUMMARY OF THE INVENTIONAn apparatus and/or method for compressing Z-stack microscopy images, substantially as shown in and/or described in connection with at least one of the figures, as set forth more completely in the claims.
Various advantages, aspects and novel features of the present disclosure, as well as details of an illustrated embodiment thereof, will be more fully understood from the following description and drawings.
So that the manner in which the above recited features of the present invention can be understood in detail, a more particular description of the invention, briefly summarized above, may be had by reference to embodiments, some of which are illustrated in the appended drawings. It is to be noted, however, that the appended drawings illustrate only typical embodiments of this invention and are therefore not to be considered limiting of its scope, for the invention may admit to other equally effective embodiments.
Disclosed herein is an apparatus and/or a method for compressing microscopy Z-stack images, substantially as shown in and/or described in connection with at least one of the figures, as set forth more completely in the claims. According to an embodiment of the present invention, an all-focus image and depth map of a specimen are compressed according to a method known to those of ordinary skill in the art such as JPEG, JPEG-LS, and JPEG 2000, although numerous other possibilities may work equally well. Depth of field is simulated from a Z-stack to produce a simulated Z-stack, as described in the commonly assigned and related pending U.S. patent application Ser. No. 13/553,434, filed Jul. 19, 2012. The approximated/simulated Z-stack is subtracted from the original Z-stack to obtain a prediction residual and then compressed similarly as the all-focus image. The compressed prediction residual is concatenated to the all-focus image and the depth map as a file to produce the final compressed result.
The digital microscope 102 captures digital images of a specimen at different depths for creating a sequence of images. The digital microscope 102 images the specimen at different depths to create the sequence of digital images. For purposes of illustration, the sequence of digital images is represented by following Expression A: I(x, y, z), where (x, y) is a spatial index with x=0 to X−1 and y=0 to Y−1 and z is a depth index with z=0 to Z−1 [Expression A].
In certain embodiments, the digital microscope 102 is coupled to the computing device 104. In other embodiments, a stack of images is created by the microscope 102 and transferred to the computing device 104, e.g., via memory stick, compact disk, network connection or the like.
The computing device 104 comprises a Central Processing Unit (or CPU) 106, support circuits 108, a memory 110 and an I/O device 112. The CPU 106 may comprise one or more commercially available microprocessors or microcontrollers that facilitate data processing and storage. The various support circuits 108 facilitate the operation of the CPU 106 and include one or more clock circuits, power supplies, cache, input/output circuits, and the like. The memory 110 comprises at least one of Read Only Memory (ROM), Random Access Memory (RAM), disk drive storage, optical storage, removable storage and/or the like.
The memory 110 comprises an Operating System (OS) 114 and an image processing module 116. The image processing module 116 facilitates compression and progressive transmission of the sequence of images of a specimen captured using the digital microscope 102. The image processing module 116 comprises a Depth of Field (DOF) simulation module 118, an all focus-image generator 120, a depth-map generator 122 and a codec (or coder decoder module) 124.
The DOF simulation module 118 simulates the sequence of images of a specimen captured using the digital microscope 102. The DOF simulation module 118 generates a low-bandwidth approximation of the sequence of images by applying Point Spread Function (PSF) operations to each pixel of an all-focus image and a blur function to each pixel according to the depth of the pixel as dictated by the depth map. Specifically, the DOF simulation module 118 approximates the sequence of images by applying depth-dependent out-of-focus blurring to the all-focus image, thereby imparting a perception of depth of specimen to users in browsing through the sequence of images in order to simulate the functionality of scrolling through a Z-stack of images without actually storing the full sequence of images in the Z-stack.
The all focus-image generator 120 generates an all focus-image in accordance with the principles of the invention. The depth-map generator 122 generates a depth map. The codec 124 applies residual coding based on a prediction from the approximated sequence of images thereby facilitating compression.
As used herein, the term “all-focus image” refers to the reduction of a Z-stack of a sequence of images of a specimen on a slide from three dimensions to two dimensions such that the result contains the entire specimen in focus, i.e., out of focus information is removed. The all-focus image is formed by weighting each image according to a particular measure and combining the weighted images accordingly, as described in commonly assigned U.S. patent application Ser. No. 13/553,265, filed Jul. 19, 2012, hereby incorporated by reference in its entirety.
The all-focus image is a composite of portions of the images in the Z-stack that are in focus at each focus depth without the out-of-focus blur. The all-focus image is represented by following Expression B: I∞(x, y) [Expression B]. In order to determine the all-focus image, the depth-map generator 120 generates a depth at each pixel to indicate which image of the Z-stack has the best focus at that pixel, where depths for each pixel form a depth-map. The depth-map is represented by following Expression C: d(x, y) [Expression C].
In some applications, there are a large number of images in the Z-stack, for example 15, 30, or 50. This imposes high bandwidth requirements for the storage or transmission of the z-stack.
In certain embodiments, the all-focus image is a suitable replacement of the original Z-stack data. However, in many scenarios, the user has to interpret the Z-stack data based on out-of-focus and depth-dependent blur. Still however, storage and transmission of Z-stack is problematic. The large size of the Z-stack data often occupies excessive disk space in storage systems, and imposes bandwidth constraints. The large size of the Z-stack data and bandwidth constraints cause noticeable and objectionable delays in viewing the Z-stack data.
The method 500 is an exemplary implementation of the image processing module 116 as executed by the CPU 106. The image processing module 116 comprises a DOF simulation module 118, an all focus-image generator 120, a depth-map generator 122 and a codec (or coder decoder module) 124.
The method 500 begins at step 502 and proceeds to step 504. At step 504, the method 500 compresses an all-focus image and a depth map of images in a Z-stack captured using the microscope 102. The all-focus image and depth map are generated according to the correlated, co-pending 3513 application based on the Z-stack images.
The method then proceeds to step 506, where the codec 124 computes a prediction residual. The prediction residual is computed by subtracting a simulation of the Z-stack from the original Z-stack images. The simulation, or approximation, of the Z-stack is described in the 3832 application, as discussed above. After subtracting the simulation Z-stack from the original Z-stack, the prediction residual remains. The subtraction is performed, according to an exemplary embodiment, on a pixel-by-pixel difference basis.
At step 508, the image processing module 116 combines the compressed all-focus image, depth map and the prediction residual into a final compressed result. The final compressed result contains the required data to recreate the original Z-stack once the result is transmitted to a remote computer system, or the like. The method ends at step 510.
As shown in
According to a common method for compressing images and sequences of images, an efficient prediction is created that approximates the original data, i.e. sequences of images. After subtracting the prediction from the original data, the difference (or prediction residual) has lower entropy and is compressed with fewer bits. In certain scenarios involving use of a low-bandwidth all-focus image and depth map to approximate the Z-stack, the residual between the approximated Z-stack and the original Z-stack has lower entropy than the original Z-stack, and thus is represented with fewer bits. In certain other scenarios, perception of reduced delay occurs in viewing the Z-stack on a remote client due to instant availability and accessibility of the approximated Z-stack, while data for the prediction residual is being transmitted in as a background process to the remote client.
In certain embodiments, an aspect of the present invention is use of the DOF simulation as described in the 3832 application. The method for simulating DOF in microscopic imaging comprises approximating a sequence of images captured by a microscope using a depth map and an all-focus image. The DOF simulation module 118 takes as input an all-focus image and corresponding depth map and produce images with simulated depth of field that approximate the appearance of the original Z-stack to reduce the size of transmitted data files. The output of the DOF simulation module 118, i.e. a simulated Z-stack, is used as a prediction for encoding the Z-stack using the encoder. Simulated Z-stacks act as a low-bandwidth approximation of the original Z-stack that may be used to progressively transmit the Z-stack such that the perceived user delay is minimized.
The step of combining the sequence of images captured from different focal distances to form a single all-focus image comprises computing a focus measure at every pixel, finding multiple candidate values for an all-focus pixel based on the focus measure and blending the candidates together according to the focus measure in order to determine a final all-focus image. The focus measures for each pixel are combined to form a depth map.
The all-focus image and depth map are compressed at block 708 according to methods commonly known to those of ordinary skill in the art, for example, JPEG, JPEG-LS, and JPEG 2000, although numerous other possibilities work equally well. To be used as a prediction, the compressed versions of the all-focus image and the depth-map represented by Expressions D and E: I∞(x, y) [Expression D] and d(x, y) [Expression E] are decompressed at block 710. As shown in
The simulation of the Z-stack is performed by the DOF simulation module 118. The simulated Z-stack 714 is subtracted from the input Z-stack 704 at block 714 to form a prediction residual represented by Expression I: r(x, y, z) [Expression I]. The prediction residual is then compressed at block 716 and combined at block 718 with the compressed all-focus image and depth map, represented by the Expressions D and E, to form the final compressed representation of the Z-stack. According to exemplary embodiments, compression methods include, among others, JPEG, JPEG2000, lossless JPEG (JPEG-LS), H.264, MPEG-4, MPEG-2, and 3-D image coding methods such as 3-D SPIHT wavelet-based coding. To decode the compressed Z-stack, the compressed versions of I∞(x, y) and d(x, y) are first split at block 720 and decompressed at blocks 722 and 724 to yield I′∞(x, y) and d′(x, y), which are used by simulate block 725 to form the Z-stack approximation I′Z(x, y). The approximation is added at block 726 to the decompressed Z-stack residual r′(x, y, z) to form the final decoded Z-stack I′(x, y, z) 728.
Although not explicitly shown in
As shown in
In certain embodiments, images in digital pathology are often very large, and thus their compressed versions use methods such as hierarchical- and/or tile-based representations. In hierarchical- and/or tile-based representations, the high-level layout differs from the representation of
In certain embodiments, the importance of the representation 800 of
In many cases, the all-focus image is all that the viewer requires, in which case there is no need to transmit other parts of the compressed Z-stack. However, in other cases the viewer chooses to browse through the actual Z-stack instead of the all-focus image; in such cases, the viewer expects to scroll quickly through the Z dimension without long delays. The timing is as indicated below in
In certain embodiments, at a time represented by Expression J: t=t1>t0 [Expression J], a user begins to scroll or browse through the Z-stack. Due to the progressive nature of the compression, the user is able to instantly see the approximated Z-stack. The quick availability is due to the computationally efficient nature of the DOF simulation module 118, which provides the approximated Z-stack without requiring transmission of any additional information. The user perceives virtually no delay in browsing the Z-stack. While the user is browsing, the image processing module 116 begins downloading the full decompressed representation of Z-stack in the background. Once full image of the decompressed representation of Z-stack are available, they are used for displaying to the viewer for browsing. In certain scenarios, the viewer selects a particular image for detailed analysis, if the decompressed representation of Z-stack comprising the particular image is fully downloaded in the background, then it is available for viewing. In certain other scenarios, there is additional delay in viewing by the time the downloading of the decompressed representation of Z-stack comprising the particular image is completed during which time the user has access to the decompressed approximated all-focus image for the particular image.
The foregoing description, for purpose of explanation, has been described with reference to specific embodiments. However, the illustrative discussions above are not intended to be exhaustive or to limit the invention to the precise forms disclosed. Many modifications and variations are possible in view of the above teachings. The embodiments were chosen and described in order to best explain the principles of the present disclosure and its practical applications, to thereby enable others skilled in the art to best utilize the invention and various embodiments with various modifications as may be suited to the particular use contemplated.
While the foregoing is directed to embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.
Claims
1. A method for compressing Z-stack microscopy images, the method being executed by a processor, comprising:
- compressing an all-focus image to form a compressed all-focus image and compressing a depth-map representation of an original Z-stack of images to form a compressed depth-map;
- computing a prediction residual as a difference between a simulated Z-stack, generated from the all-focus image and the depth-map representation, and the original Z-stack; and
- concatenating the prediction residual to the compressed all-focus image and the compressed depth-map representation to produce a compressed representation of the original Z-Stack.
2. The method of claim 1 further comprising decoding the compressed representation of the original Z-stack.
3. The method of claim 2, wherein decoding the compressed representation of the original Z-stack further comprises retrieving and decompressing the compressed all-focus image and the compressed depth-map to yield a decompressed all-focus image and a corresponding depth-map.
4. The method of claim 1 wherein compressing the all-focus image and the depth-map is performed using one of JPEG, JPEG-LS and JPEG 2000.
5. The method of claim 3 further comprising simulating a Z-stack using the decompressed all-focus image and the corresponding depth-map for generating an approximation of the original Z-stack.
6. The method of claim 5 further comprising adding the simulated Z-stack to the prediction residual to form a final Z-stack.
7. The method of claim 1 further comprising transmitting information to form the simulated Z-stack to a remote client.
8. The method of claim 7 further comprising transmitting the prediction residual as a background process while a user views the simulated Z-stack on the remote client.
9. The method of claim 1, wherein the simulated Z-stack is formed using parameters comprising a correlation between intensity of blurring and depth at each pixel of the images in the simulated Z-stack.
10. The method of claim 9 further comprising storing the parameters in the compressed representation of the original Z-Stack.
11. An apparatus for compressing Z-stack microscopy images comprising:
- a digital microscope for capturing a sequence of images of a specimen at different focus depths forming an original Z-stack; and
- an image processing module for managing compressing the Z-stack comprising: a depth of field (DOF) simulation module for generating a simulated Z-stack based on the original Z-stack; an all-focus image generator for generating an all-focus image for the original Z-stack; a depth-map generator for generating a depth-map corresponding to the all-focus image; and a codec for compressing the all-focus image to generate a compressed all-focus image, compressing the depth-map to generate a compressed depth-map, computing a prediction residual as a difference between the simulated Z-stack and the original Z-stack and concatenating the prediction residual to the compressed all-focus image and the compressed depth-map to produce a compressed representation of the original Z-stack.
12. The apparatus of claim 11 wherein the codec further decodes the compressed representation of the original Z-stack.
13. The apparatus of claim 12, wherein the codec further retrieves and decompresses the compressed all-focus image and corresponding depth-map and to yield a decompressed all-focus image and a corresponding depth-map.
14. The apparatus of claim 11 wherein compressing the all-focus image and the depth map is performed using one of JPEG, JPEG-LS and JPEG 2000.
15. The apparatus of claim 13 wherein the codec further simulates a Z-stack using the decompressed all-focus image and the corresponding depth-map for generating an approximation of the original Z-stack.
16. The apparatus of claim 15 wherein the codec further adds the simulated Z-stack to the prediction residual to form a final Z-stack.
17. The apparatus of claim 11 wherein the image processing module further transmits information to form the simulated Z-stack to a remote client.
18. The apparatus of claim 17 wherein the image processing module further transmits the prediction residual as a background process while a user views the simulated Z-stack on the remote client.
19. The apparatus of claim 11, wherein the DOF simulation module further simulates the Z-stack using parameters comprising a correlation between intensity of blurring and depth at each pixel of the images in the simulated Z-stack.
20. The apparatus of claim 19, wherein the DOF simulation module further stores the parameters in the compressed representation of the original Z-stack.
6456285 | September 24, 2002 | Hayhurst |
7372985 | May 13, 2008 | So et al. |
7593581 | September 22, 2009 | Hecht et al. |
7724937 | May 25, 2010 | So et al. |
8252829 | August 28, 2012 | Duggan |
8331627 | December 11, 2012 | Xiong et al. |
2005198801 | July 2005 | JP |
2011150400 | August 2011 | JP |
- Won-Ki Keong et al., “Interactive Histology of Large-Scale Biomedical Image Stacks”, manuscript received Mar. 31, 2010; accepted Aug. 1, 2010; posted online Oct. 24, 2010; mailed on Oct. 16, 2010, 10 pages.
- Office Action dated May 26, 2010, Japan Patent Application No. 2013-142662, 11 pages.
Type: Grant
Filed: Jul 19, 2012
Date of Patent: Aug 26, 2014
Patent Publication Number: 20140023283
Assignee: Sony Corporation (Tokyo)
Inventors: Ming-Chang Liu (San Jose, CA), Mark Robertson (Cupertino, CA)
Primary Examiner: Anh Do
Application Number: 13/553,428
International Classification: G06K 9/36 (20060101);