|
|
|
|
|
|
|
|
| Re: Overlay Point Sources on Maps [message #51683 is a reply to message #51456] |
Tue, 28 November 2006 11:45  |
David Fanning
Messages: 11724
Registered: August 2001
|
Senior Member |
|
|
greg michael writes:
> I spent a while messing around with griddata() before I realised you
> could do this:
>
> restore,"ptsource_carbon.sav"
> device,decomposed=0
> loadct,39
> grid=fltarr(360,180)
> grid[lon+180,lat+90]=soilc
> tvscl, grid
>
> This particular dataset is simpler than I first thought - the points
> are already gridded at 1x1 degree. The odd thing is that many points
> are repeated (up to 8 times). The values, though, seem to be the same
> in these cases. In general I think a gridding approach is more robust
> for this type of problem.
Humm. Guess I should have realized this is possible, too.
It makes using the land mask a bit different.
Well, there are often several ways to skin this cat.
Cheers,
David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
Sepore ma de ni thui. ("Perhaps thou speakest truth.")
|
|
|
|
| Re: Overlay Point Sources on Maps [message #51685 is a reply to message #51456] |
Tue, 28 November 2006 11:37 |
greg michael
Messages: 163
Registered: January 2006
|
Senior Member |
|
|
I spent a while messing around with griddata() before I realised you
could do this:
restore,"ptsource_carbon.sav"
device,decomposed=0
loadct,39
grid=fltarr(360,180)
grid[lon+180,lat+90]=soilc
tvscl, grid
This particular dataset is simpler than I first thought - the points
are already gridded at 1x1 degree. The odd thing is that many points
are repeated (up to 8 times). The values, though, seem to be the same
in these cases. In general I think a gridding approach is more robust
for this type of problem.
regards,
Greg
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
