Bamfield Resources

Online Resources


Charts & Maps

For a line coastline:

load BMS_coastline
plot(BC(:,1),BC(:,2));
set(gca,'dataaspectratio',[1 cosd(49) 1]); % so land doesn't appear squashed
set(gca,'tickdir','out','linewidth',2);
axis([-125.25 -125 48.78 48.94]);
xlabel('Longitude');
ylabel('Latitude');
title('Barkley Sound Area');

For a filled coastline:

load BMS_coastline
ii=find(isnan(BC(:,1))); % fine separators between lines
landcolorRGB=[.8 .8 .8];
for k=1:length(ii)-1,
  patch(BC(ii(k)+1:ii(k+1)-1,1),BC(ii(k)+1:ii(k+1)-1,2),landcolorRGB);
end;
set(gca,'dataaspectratio',[1 cosd(49) 1]); % so land doesn't appear squashed
set(gca,'tickdir','out','linewidth',2);
axis([-125.25 -125 48.78 48.94]);
xlabel('Longitude');
ylabel('Latitude');
title('Barkley Sound Area');