Re: Function Graphics Map Projection Woes [message #77674 is a reply to message #77673] |
Mon, 19 September 2011 08:33   |
David Fanning
Messages: 11724 Registered: August 2001
|
Senior Member |
|
|
David Fanning writes:
>
> Folks,
>
> Sigh...
>
> OK, I'm on to my next topic on my list in my effort
> to learn about function graphics: map projections.
>
> Can anyone explain to me why this code only shows
> the Northern hemisphere?
>
> data = Dist(200)
> imgObj = Image(data, limit=[-90,-180,90, 180], $
> grid_units=2, $
> map_projection='Equirectangular')
>
> Does *anything* in function graphics work correctly!?
Seriously. I am going to absolutely tear my hair out!
Function graphics are so &*^$#ed up they are impossible
to work with!
I have some globally gridded data.
IDL> help, tnmin
TNMIN FLOAT = Array[96, 73]
I have *COMPLETELY* regularly spaced longitude and latitude
vectors.
IDL> help, lon,lat
LON DOUBLE = Array[96]
LAT DOUBLE = Array[73]
IDL> print, lon
0.00000000 3.7500000 7.5000000 11.250000
15.000000 18.750000 22.500000 26.250000
30.000000...
IDL> print, lat
90.000000 87.500000 85.000000 82.500000
80.000000 77.500000 75.000000 72.500000
70,000000 ...
I try to display this image in a map projection like this:
imgObj = Image(tnmin, lon, lat, limit=[-90,-180,90, 180], $
grid_units=2, map_projection='Cylindrical')
And I get this error:
% IMAGE: X and Y parameters must be evenly spaced
What the @*%$!
Cheers,
David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.idlcoyote.com/
Sepore ma de ni thui. ("Perhaps thou speakest truth.")
|
|
|