I am trying to plot geospatial sea surface temperature correlation data on the world map (only Northern hemisphere). I am using the image function (plot correlation data) along with the map function (plot coastlines). On their own, I get them to work. My problem is I am not able to overlay the map on the image, the map just shows up as a small point at the left corner at the bottom of the plot.
lon has the structure
0.0000 0.9375 ... 358.1250 359.0625 (384 entries) .
lat has the structure
0.4675309 1.4025927 ... 88.3570035 89.2842275 (96 entries).
str(cor_nao_map)
num [1:384, 1:96] 0.0376 0.0443 0.0476 0.0485 0.0444 ...
image(x=lon,y=lat,z=cor_nao_map,col=diverge_hcl(100),zlim=c(-0.43,0.43),xlab="Longitude in °",ylab="Latitude in °N")
map('world',projection="gilbert",add=T)
It seems as a new user I am not able to post a screenshot, sorry
EDIT: The packages I am using are graphics for the image function and maps for the map function.