comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Question about projection for Google Earth
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Question about projection for Google Earth [message #85455 is a reply to message #85454] Thu, 08 August 2013 21:50 Go to previous messageGo to previous message
timothyja123 is currently offline  timothyja123
Messages: 57
Registered: February 2013
Member
> All of these details are handled automatically for you in cgImage2KML.
>
> http://www.idlcoyote.com/idldoc/cg/cgimage2kml.html
>
> Cheers,
> David

Hi David,

Thanks for the reply. I have already tried to use your procedure however I ran into a problem I wasnt sure how to solve. Your procdure takes an image as an argument which is fine as I can use the image produced by my existing code.
However I then need to supply a map coordinate object (cgMap) from which map projection information and map boundaries for the image overlay can be obtained.

The problem is that your cgMap function says: "Only GCTP projections are allowed. If you wish to use projections normally set up with Map_Set, use the comparable cgMap_Set command. "
But that command does not produce the cgMap object that I need for cgImage2KML.

I'm most likely missing something obvious (as I said I'm new to projections and mapping).

I tired using a map_set projection that had a matching GCTP projection MILLER_CYLINDRICAL and using the produced image with cgImage2KML but this didnt seem to work.

limit = dblarr(4)
limit[0] = min(lat)
limit[1] = min(lon)
limit[2] = max(lat)
limit[3] = max(lon)

polon = limit[1]+0.5*(limit[3]-limit[1])
map_set, 0,polon,0,/MILLER_CYLINDRICAL,limit=limit, /noborder, xmargin=0,ymargin=0


contour, values*mask,lon,lat, levels=levels,C_color=c_levels,/overplot,c_labels = 0,/cell_fill, min_value = 2


most_image = tvrd()
Device, close = 1


name = 'MOST_max_amplitude.kml'
outputOverlayFile = filepath(name, root_dir=sourcepath(), subdir=['output'])

googleMapCoord = Obj_New('cgMap', 'Miller Cylindrical', limit=limit)
cgImage2KML, most_image, googleMapCoord, filename=outputOverlayFile, LATLONBOX=LATLONBOX
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Keyword UNAME not allowed in call to: CW_FSLIDER
Next Topic: parallel sign in postscript output in device fonts (direct graphics)

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Thu Oct 09 23:37:39 PDT 2025

Total time taken to generate the page: 0.96056 seconds