Skip to main content

Activity 5: Enhancement by Histogram Manipulation

IN this activity, we will try to enhance an image by manipulating its histogram. By enhance, I mean the improvement of a dark image such that features disguised by the dark areas may be revealed. The histogram of a graylevel image, when normalized to the number of pixels, gives us the probability distribution function (PDF). 

For example, if we have the following color image:
Figure 1a. Dark color image which will be enhanced
 (taken from  http://twoleggedtripod.net/index.php?postid=3  )  
Then, convert it to a graylevel image using Scilab: 
Figure 1b. Grayscaled version of the image in Figure 1a
We get its PDF by dividing each pixel by the total number of pixels.
Figure 1c. Code for obtaining PDF of image
Figure 1d. PDF of the graylevel image
From the PDF, we see that the distribution is biased towards the low values.

 If the image has graylevels, r, with a PDF of p1(r) and a cumulative distribution function (CDF) given by,

then we can map it to a new set of graylevels, z, described by a different CDF described by:
where p2(z) is the new PDF of the image. This technique is called backprojection and is best illustrated by the diagram found in Figure 2.
Figure 2. Illustration of the backprojection method. Look for the graylevel value, r,  of the pixel, and find its corresponding CDF value T(r). Project the value of T(r) for the graylevel on the desired CDF , G(z). Look for the graylevel value which corresponds to the mapped value for G(z). This graylevel value will replace the original value.



For the image in Figure 1b, we determine the CDF using the cumsum() command in SciLab as shown in the code snippet in Figure 3. The CDF is shown in Figure 4.
Figure 3. Code for obtaining CDF of image
Figure 4. CDF of the graylevel image in Figure 1b.
We then map the y values from the original CDF to a desired CDF. In this case we used a linear CDF and the CDF of a Gaussian PDF as shown in Figure 5.
Figure 5a. Linear CDF
Figure 5b. CDF of Gaussian PDF
We then use the interp() function in SciLab to find the corresponding graylevel value, z, for the y value that we mapped on the desired CDF. We then replace the graylevel value for each pixel of the original image by the new grayvalue obtained from the backprojection method. The code used to do this is shown below.



Shown below are the enhanced images using the linear CDF and the CDF of a Gaussian PDF. The corresponding CDFs of the enhanced images are also shown.
Figure 6a. Image enhanced using a linear CDF
Figure 6b. CDF of the image in Figure 6a.

Figure 6c. Enhanced image using CDF for a Gaussian PDF
Figure 6d. CDF of image in Figure 6c
The enhanced images show some feature that we cannot see just by looking at the original grayscale image such as some of the features of the branches and leaves of the surrounding trees. The image which resulted from the projection on a linear CDF is sharper compared to that obtained for the CDF of a Gaussian PDF. However, the image is too bright making some parts of the image to be indistinguishable due to the brightness. For the enhanced image using the CDF of a Gaussian PDF, the contrast is not as sharp and most parts of the original image are still distinguishable.

For this activity, I would like to thank Ms. Eloisa Ventura, Ms. Maria Isabel Saludares and Mr. Benjur Borja for useful discussions. 

 Finally, I give myself a grade of 10/10 for being able to fulfill the requirements in the activity. 

 References:
1.  M. Soriano, A5: Enhancement by Histogram Manipulation Manual, 2012

Comments

Post a Comment

Popular posts from this blog

Activity 10 Applications of Morphological Operation 3 of 3: Looping through images

When doing image-based measurements, we often want to separate the region of interest (ROI) from the background. One way to do this is by representing the ROIs as blobs. Binarizing the image using the optimum threshold obtained from the image histogram simplifies the task of segmenting the ROI. Usually, we want to examine or process several ROIs in one image. We solve this by looping through the subimages and processing each. The binarized images may be cleaned using morphological operations.  In this activity, we want to be able to distinguish simulated "normal cells" from simulated "cancer cells" by comparing their areas. We do this by taking the best estimate of the area of a "normal cell" and making it our reference.  Figure 1 shows a scanned image of scattered punched papers which we imagine to be cells examined under the microscope. These will be the "normal cells." Figure 1. Scattered punched paper digitized using flatbe...
1 comment
Read more

Activity 7: Morphological Operations

When talking about morphology, what immediately comes to mind are the forms and structures or shapes of objects. Hence, performing morphological operations imply that the shape or form of an object is altered.       In this activity, we will perform morphological operations on binary images. In particular, we make use of erosion and dilation . Erosion and dilation were performed on the following: 1. A 5×5 square 2. A triangle, base = 4 boxes, height = 3 boxes 3. A hollow 10×10 square, 2 boxes thick 4. A plus sign, one box thick, 5 boxes along each line Using each of the structuring elements below: 1. 2×2 ones 2. 2×1 ones 3. 1×2 ones 4. cross, 3 pixels long, one pixel thick. 5. A diagonal line, two boxes long, i.e. [[0 1],[1 0] ].      When performing these operations, it is important to note the “anchor” or “origin” of the structuring element in order to give an accurate prediction of the result. For the 2x2 ones, 2...
Post a Comment
Read more

Activity 12: Basic Video Processing

Hello!  In this activity we will try to process a video of a kinematic event in order to extract information such as constants, frequencies, etc. For our group, we took a video of a 3D spring pendulum which we observed in one plane. We would like to trace its path and then try to determine its phase-space plot. The mass was covered in masking tape with the bottom colored red to facilitate easier segmentation. The video was taken using a Canon D10 camera at frame rate of 30fps.  Media 1. Video of the spring pendulum (first 50 frames only) The frames of the video were then extracted using Avidemux 2.5. The mass was then segmented from each frame using parametric segmentation. The patch of the region of interest (ROI) used for color segmentation is shown in Figure 1. Figure 1. Patch used to segment ROI  Using morphological operations, particularly Open and Close operations, the segmented images were cleaned. The extracted frames for different observation...
Post a Comment
Read more