How to overlay a shape file onto a wind density contour plot

Hi,

I am having problems overlaying a shape file of Northern Ireland onto a wind density colour contour map that I have generated within the latitude and longitude of Northern Ireland.

I have attached the wind density contour map

If any one could help me with this, would be much appreciated.

Thanks

Answer 1

[Chad Greene]

For a single landmass,

S = shaperead('filename.shp','usegeocoords',true); 
plot(S.Lon,S.Lat)

If the shapefile contains multiple landmasses, loop the plot line from 1 to N.

[EL]

Thank you Chad for your response, Much appreciated.

Im afraid this isn't working, the command: plot(S.Lon,S.Lat) is not accepted, even if I had my lat and lon defined. Any other suggestions?

Thanks again

[Chad Greene]

What do you mean whenyou say the plot command isn't accepted? Is there an error message?

[Nabeel]

same case is with me i want to over lay shape file on a surf image which i have produced by following code

base=surf(LON,LAT,ETOPO2); shading interp; axis equal, view(0,90);

i use the hold on command for surf and than try to plot shape file but it is not displayed on surf image (topographic)

[Chad Greene]

The plot command is placing lines at z=0. So presumably, the lines are plotted below the z values of your topographic surface. Use

plot3(S.Lon,S.Lat,max(ETOPO2(:))+ones(size(S.Lon)))