comp.lang.idl-pvwave archive: archive » bounding box subsetting for lat/lon

Home » Public Forums » archive » bounding box subsetting for lat/lon

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

bounding box subsetting for lat/lon [message #69856]

Mon, 15 February 2010 22:27

mankoff
Messages: 131
Registered: March 2004

Senior Member

Hi,

I'm wondering if a routine already exists that can subset a data set
based on non-square lat/lon limits. I have large 1D arrays of lat,
lon, and data, and want to do something like this:

roi = WHERE( lat LE a AND lat GE b AND lon LE c AND lon GE d )
lat = lat[roi] & lon=lon[roi] & data = data[roi]

But at a minimum I want to specify 4 lat/lon elements that do not
defined the edges of the box but the four corners, so the roi could be
a diamond shape (when projected on, say, the /cylinder projection).
And ideally I could define a path of n points defining a closed region
and get the indices back that lie within that region.

Anyone know if this exists? I don't see it under the ROI section, and
searching for bounding box brings up too much noise.

Thanks,

-k.

Report message to a moderator

Previous Topic:	Window or Widget to front?
Next Topic:	How to match multiple arrays and plot graph?

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Oct 08 18:41:17 PDT 2025

Total time taken to generate the page: 0.00501 seconds