Hello everyone,
I recently installed the new IDL 8.0 on my machine and am having difficulties overplotting scattered geographical data onto map projections. I have data in single columns of latitude, longitude and elevation. I have been able to do this with direct graphics easily enough in the past, however with the new object oriented graphics functions, I am having some difficulty. I would like simply to create a map of the Arctic, such as this:
map1 = map('polar stereographic',limit=[60,-180,90,180])
cont = map_continents(fill_color='light gray',/current)
and then using the plot function overplot the scattered lat/lon data measurements on the map (a symbol would suffice). All the IDL help examples either use gridded data, or grid it before plotting, however I would simply like to plot the original locations. I have tried the "convertcoord" method associated with the new map function, to try and get my data from geographic coordinates into an associated x,y space:
xy = map1.convertcoord(lon,lat)
but have been unable to get the transformations to work correctly; my locations end up either off the map or in a small cluster in the middle, depending upon the keywords I use. Does anyone have any advice?
Thanks,
Jeremy
|