Write a Matlab code for the following:

1. Load image Balloon.tif and convert it into a grayscale image, denoted with img1

2. Decompose img1 with wavelet transform using function wavedec2 (you can decide the number of decomposition level). Display the wavelet coefficients as an image. Note: you will need to rescale the values to be in the range of [0 255]; otherwise, the image appears black.

3. Reconstruct the edge map from the detail subbands using function waverec2. Hint: set the approximation subband to zero. Save the reconstructed edge map into an image.

4. Using morphology operator to remove the minor edge segments and reduce the width of the edges to one pixel size. Save the result image.

5. Compute the energy, entropy, and homogeneity features for both processed edge maps and calculated the differences.