comp.lang.idl-pvwave archive: archive » Re: Newbie needs help...

Home » Public Forums » archive » Re: Newbie needs help...

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: Newbie needs help... [message #23135 is a reply to message #23134]

Wed, 10 January 2001 12:49

Liam E. Gumley
Messages: 378
Registered: January 2000

Senior Member

David Fanning wrote:
> (Where *is* Liam, anyway!?)

I'm right here, trying to get my semi-annual research report written...

> I guess if I were in your shoes, I might try to
> reformat the part of the image I was interested
> in onto a regular lat/lon grid. I'd probably try
> to use TRIANGULATE and TRIGRID first, to see if
> that would work.

In my experience, TRIANGULATE and TRIGRID work just fine when the input
arrays are small (say a few hundred rows and columns). But when you have
satellite images with thousands of rows and columns, TRIANGULATE and
TRIGRID just don't work.

> There are probably other ways to warp the image
> onto the map, but there is no getting around the
> need to do it. What you would prefer, I'm sure, is
> to have a map warped to the image. But this is NOT
> possible in IDL.

Displaying a satellite image on a map projection is one problem, and my
imagemap procedure provides one solution.

However resampling a satellite image to a regular lat/lon grid in such a
way that data values are preserved is a different problem, and there is
no easy solution in IDL. Typically, the approach is to define your
lat/lon grid (in some map projection) and then for each map grid cell,
find the closest matching location in the satellite image (nearest
neighbor). Refinements include bilinear or cubic spline interpolation.
However this algorithm is not easy to implement efficiently in IDL,
because of the penaly associated with loops. Most people implement a
resampling algorithm in some other language (e.g. C or FORTRAN), and
then read the results in IDL.

You mentioned that you are using MODIS data at 1 km resolution (1354 x
2030 pixels pre granule). One approach we use is to define a global
equal area grid; we happen to use 25 km x 25 km grid cells. We loop over
each 1 km pixel, and based on it's lat/lon, we accumulate the following
statistics in each grid cell:

Number of observations,
Sum of observations,
Minimum observation,
Maximum observation.

From these, we can compute mean and standard deviation. it is quite
straightforward to then resample the equal area grid to an equal angle
grid which can be visualized in IDL.

If you are interested in more details, contact me directly.

Now back to my semi-annual report (sigh).

Cheers,
Liam.
http://cimss.ssec.wisc.edu/~gumley

PS: Here's a couple of MODIS images created in IDL:
http://earthobservatory.nasa.gov/Newsroom/NewImages/images.p hp3?img_id=4549
http://earthobservatory.nasa.gov/Newsroom/NewImages/images.p hp3?img_id=4383

Report message to a moderator

[Message index]

		Re: Newbie needs help... By: Bruce Bowler on Thu, 11 January 2001 05:29
		Re: Newbie needs help... By: James Kuyper on Wed, 10 January 2001 14:55
		Re: Newbie needs help... By: Craig Markwardt on Wed, 10 January 2001 14:46
		Re: Newbie needs help... By: Liam E. Gumley on Wed, 10 January 2001 12:49
		Re: Newbie needs help... By: davidf on Wed, 10 January 2001 12:18
		Re: Newbie needs help... By: Bruce Bowler on Wed, 10 January 2001 11:44
		Re: Newbie needs help... By: davidf on Wed, 10 January 2001 09:25
		Re: Newbie needs help... By: Bruce Bowler on Wed, 10 January 2001 08:33
		Re: Newbie needs help... By: davidf on Wed, 10 January 2001 08:08
		Re: Newbie needs help... By: Liam E. Gumley on Thu, 11 January 2001 07:53
		Re: Newbie needs help... By: Paul van Delst on Thu, 11 January 2001 06:45
		Re: Newbie needs help... By: davidf on Thu, 11 January 2001 06:59

Previous Topic:	Re: IDL compatibility
Next Topic:	Tickmarks and ticklabels

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Oct 11 08:17:19 PDT 2025

Total time taken to generate the page: 0.24246 seconds