- added computation of normal vector
- bilinear interpolation: treat missing value as zero