Thread Subject: Interpolate Data in Stream Banks

I have a shoreline data from NOAA's coastline extractor and plotted sampling sites within the creeks represented on this shoreline. I then have data files with data readings at 15 minute intervals for each sampling site. I would like to interpolate the data from each sampling site through the creek between each sampling point. In otherwords the ultimate goal is to have, say water temperature, represented by cool colors for cold water and warm colors for warm water, interpolated from the sampling sites through the creeks and contained within the shoreline.dat area.
Any suggestions for helpful functions are appreciated.

Thanks

load ircoast.dat % loads the lat/long data file into Matlab
lat = ircoast (:,2); % defines all the rows (:) in the second column (2) as the variable lat
lon = ircoast (:,1); % defines all the rows (:) in the first column (1) as the variable lon
plot (lon, lat, 'k'); % cross-plot of variables defined in lon and lat in the color black

% Reads the stations.txt file % Reads the first column as a string and stores it to as station_name. % Reads the second column as a floating point number and stores it as lat. % Reads the third column as a floating point number and stores it as lon. station_name,lat,lon? = textread('stations.txt','%s %f %f', 'delimiter',',');

% Loops through all the lines in the file
%load 'station_name.txt'
[station_name,lat,lon] = textread('station_name.txt','%s %f %f', 'delimiter',',');
max(size(station_name));% = the number of station_names and thus the number of lines in the file
for i = 1:max(size(station_name));

% Writes the station_name to the screen at the lat and lon cordinates. % The font size is small to fit all the stations to the screen.
text(lat(i),lon(i),station_name(i),'Color',[0 0 1], 'FontSize', 8);
end

Most specifically I am interested in using a kriging interpolation. I believe the hard part is identifying the shoreline information as the barrier to the interpolation display.

Thanks again