| Re: /fill_ocean [message #10405] |
Mon, 24 November 1997 00:00 |
davidf
Messages: 2866
Registered: September 1996
|
Senior Member |
|
|
Kelly Dean (dean@cira.colostate.edu) writes:
> The MAP_CONTINENTS in IDL 5.0 has a /FILL_CONTINENT keyword which
> works out great, but I need a /FILL_OCEAN keyword. I am attempting to
> add the worldelev.dat to a map, but I don't need the ocean depths.
> Anyone have an idea how to do a /FILL_OCEAN on a map projection?
I don't have any great suggestions, but have you tried something
like this:
water = Where(worldelv LT 130)
worldelv(water) = black
I don't know if 130 is the right number or not. I picked it
based on the histogram of the image, and it seems to work
reasonably well. Perhaps with some fooling around it would
work for you.
Cheers,
David
-----------------------------------------------------------
David Fanning, Ph.D.
Fanning Software Consulting
E-Mail: davidf@dfanning.com
Phone: 970-221-0438
Coyote's Guide to IDL Programming: http://www.dfanning.com/
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
