It is my understanding that you would like to specify image data as the initial condition to a time-dependent PDE.
You can do so by interpolating the image data to the coordinates of the vertices in your mesh and then pass those interpolated values as an initial condition.
Suppose your image data, defined at background vertices "xb" and "yb", is stored in the variable "ub0". We now want to get the values, "u0", of this initial condition at your mesh vertices with coordinates "x" and "y". The interp2 function will interpolate the background image data to the mesh vertices:
>> u0 = interp2(xb,yb,ub0,x,y);
You can then pass "u0" as an initial condition to a PDE solver since it is defined at the mesh vertices. For more information on interp2, please refer to:
