Spring 2018 CS543/ECE549
Assignment 1: Shape from Shading
Due date: Monday, February 19, 11:59:59 PM
MATLAB Instructions  see here for Python
The goal of this assignment is to implement shape from shading as described in Lecture 4
(see also Section 2.2.4 of Forsyth & Ponce 2nd edition).
 Download the sample code and data. The data, in the croppedyale directory,
consists of 64 images each of four subjects from the Yale Face database.
The light source directions are encoded in the file names. The code consists of several
.m functions. Your task will be to add some code to the toplevel script, run_me.m,
and to fill in the code in photometric_stereo.m and get_surface.m, as
explained below. The remaining files
are utilities to load the input data and display the output.
 For each subject (subdirectory in croppedyale), read in the images and
light source directions. This is accomplished by the function LoadFaceImages.m
provided in the zip file (which, in turn, calls getpgmraw.m
to read the PGM files). LoadFaceImages returns the images for the 64 light source
directions and an ambient image (i.e., image taken with all the light sources
turned off).
 Preprocess the data: subtract the ambient image from each image in the light source stack,
set any negative values to zero, rescale the resulting intensities to between 0 and 1
(they are originally between 0 and 255).
Hint: these operations can be done without
using any loops whatsoever. You may want to look into MATLAB's bsxfun function.
 Estimate the albedo and surface normals. For this, you need to fill in
code in photometric_stereo.m, which is a function taking as input
the image stack corresponding to the different light source directions and the matrix
of of the light source directions,
and returning an albedo image and surface normal estimates. The latter should be stored in
a threedimensional matrix. That is, if your original image dimensions are h x w,
the surface normal matrix should be h x w x 3, where the third dimension corresponds
to the x, y, and zcomponents of the normals. To solve for the albedo and the normals,
you will need to set up a linear system as shown in slide 20 of Lecture 4.
Hints:
 To get the leastsquares solution of a linear system, use MATLAB's backslash operator.
That is, the solution to Ax = b is given by x = A\b.
 If you directly implement the formulation of slide 20 of the lecture, you will
have to loop over every image pixel and separately solve a linear system in each iteration.
There is a way to get all the solutions at once by stacking the unknown g vectors for
every pixel into a 3 x npix matrix and getting all the solutions with a single
application of the backslash operator.
 You will most likely need to reshape your data in various ways before and after
solving the linear system. Useful MATLAB functions for this include reshape
and cat.
 You may also need to use elementwise operations. For example, for two equalsize
matrices X and Y, X .* Y multiplies corresponding elements, and X.^2
squares every element. As before, bsxfun can also be a very useful function
here.
 Compute the surface height map by integration. The method is shown in slide 23
of Lecture 4, except that instead of continuous integration of the partial derivatives
over a path, you will simply be summing their discrete values. Your code implementing
the integration should go in the get_surface.m file.
As stated in the slide, to get the best results, you should compute integrals
over multiple paths and average the results. You should implement the following variants of
integration:
 Integrating first the rows, then the columns. That is, your path first goes along
the same row as the pixel along the top, and then goes vertically down to the pixel.
It is possible to implement this without nested loops using the cumsum function.
 Integrating first along the columns, then the rows.
 Average of the first two options.
 Average of multiple random paths. For this, it is fine to use
nested loops. You should determine the number of paths experimentally.
 Display the results using functions display_output and plot_surface_normals
included in the zip archive.
Extra Credit
On this assignment, there are not too many opportunities for "easy" extra credit. This said, here are some ideas for exploration:
 Generate synthetic input data using a 3D model and a graphics renderer and run your method on this data. Do you get better results than on the face data?
How close do you get to the ground truth (i.e., the true surface shape and albedo)?
 Investigate more advanced methods for shape from shading or surface reconstruction from normal fields.
 Try to detect and/or correct misalignment problems in the initial images and see if you can improve the solution.
 Using your initial solution, try to detect areas of the original images that do not meet the assumptions of the method (shadows, specularities, etc.).
Then try to recompute the solution without that data and see if you can improve the quality of the solution.
If you complete any work for extra credit, be sure to clearly mark that work in your report.
Grading Checklist
You should turn in both your code and a report
discussing your solution and results. For full credit, your report should include a section for each of
the following questions:
 Briefly describe your implemented solution, focusing especially on the more
"nontrivial" or interesting parts of the solution. What implementation choices did you
make, and how did they affect the quality of the result and the speed of computation?
What are some artifacts and/or limitations of your implementation, and what are
possible reasons for them?
 Discuss the differences between the different integration methods
you have implemented for #5 above. Specifically, you should choose one subject,
display the outputs for all of ad (be sure to choose viewpoints that make the differences
especially visible), and discuss which method produces the best results and why.
You should also compare the running times of the different approaches. For timing, you
can use tic and toc functions. For the remaining subjects (see below),
it is sufficient to simply show the output of your best method, and it is not necessary
to give running times.
 For every subject, display your estimated albedo maps and screenshots of height maps
(use display_output and plot_surface_normals).
When inserting results images into your report, you
should resize/compress them appropriately to keep the file size manageable  but make sure
that the correctness and quality of your output can be clearly and easily judged.
For the 3D screenshots, be sure to choose a viewpoint that makes the structure as clear
as possible (and/or feel free to include screenshots from multiple viewpoints).
You will not receive credit for any results you have obtained, but failed to include
directly in the report PDF file.
 Discuss how the Yale Face data violate the assumptions of the shapefromshading method
covered in the slides. What features of the data can contribute to errors in the results?
Feel free to include specific input images to illustrate your points. Choose one subject and
attempt to select a subset of all viewpoints that better match the assumptions of the method.
Show your results for that subset and discuss whether you were able to get any improvement
over a reconstruction computed from all the viewpoints.
Instructions for Submitting the Assignment
Your submission should consist of the following:
 All your code in a single zip file.
The filename should be lastname_firstname_a1.zip.
 Your report in a single PDF file with all your results (including images) and discussion.
The filename should be lastname_firstname_a1.pdf.
The files will be submitted through Compass 2g. Upload instructions:
 Log into https://compass2g.illinois.edu and find this class and assignment.
 Upload (1) your PDF report and (2) the zip file containing your code as two attachments.
 Hit Submit.
Multiple attempts will be allowed but only your last submission before the deadline will be graded.
We reserve the right to take off points for not following directions.
Late policy: You lose 25% of the points for every day the assignment is late.
If you have a compelling reason for not being able to submit the assignment on time
and would like to make a special arrangement, you must send me email at least a week before the due date
(any genuine emergency situations will be handled on an individual basis).
Academic integrity: Feel free to discuss the assignment with each other in general terms,
and to search the Web for general guidance (not for complete solutions). Coding should be done
individually. If you make substantial use of some code snippets or information from outside sources,
be sure to acknowledge the sources in your report. At the first instance of cheating (copying from
other students or unacknowledged sources on the Web), a grade of zero will be given for the assignment.
At the second instance, you will automatically receive an F for the entire course.
