Projection of grid data onto world map

Hello all,

I have grid data (spacing of grid data depends on climate model) that I would like to overlay on WorldMap, lat[-90 90] & lon [0 360], using scatterm (proportional circle map). When I do this the projection of the data does not look correct. How can I ensure that the data I am projecting onto the World Map is correct? Thank you and I look forward to learning something new!

[Vesp]

I have figured it out. This is what I did (in case it may be helpful to someone else).

In order to use scatterm(lat,lon,z*10,z,'filled') I had to create a matrix with the specific models coordinates. I then extracted the latitudes and longitudes into their own vectors to use in the scatterm. I also transposed the data,z, into vector to match the size of lat and lon.

I hope this was clear, but the main idea for my case was to have the latitude and longitude as its own vector to supply in the scatterm function for the projection to be correct. Initially I had extracted the position of the data in the matrix so I was supplying rows,columns not lat,lon of my data.

Thank you for the help :)

Answer 1

[KSSV]

You may look into the following package which might be helpful to you:

https://www.eoas.ubc.ca/~rich/map.html

[Vesp]

Thank you for your suggestion. I have looked through the examples on the webpage you provided and didn't see anything regarding proportional circle maps using this. Ultimately my goal is to produce this for the data that I have, which has results mostly along the coastlines of countries: http://evankrueger.blogspot.com/2012/04/proportional-circle-map.html

I have been using scatterm and world('world') function to try to produce this, but the projection is not correct.