| Re: Help needed in remapping data using IDL [message #70785] |
Wed, 05 May 2010 15:38 |
David Fanning
Messages: 11724
Registered: August 2001
|
Senior Member |
|
|
kaba writes:
> Thanks for the link, I have downloaded both the coyote and catalyst
> libraries.
> I have started writing a new code based on the article you send me,
> and I have a question.
>
> (1) My input file in this case is HDF 4 not netCDF. Hence I get an
> error message when I call the subroutine
> FindMapBoundy(). How do I get around this issue?
You will have to get the corner lat/lon points from
the geolocation information in the HDF 4 file and
convert them to XY points using MAP_PROJ_FORWARD.
This will allow you to find the "boundary" XY values.
You can see how I have done this for netCDF files
in ndcf_coord.pro in the catalyst/source/coordinates
directory. In fact, if you want to write a
hdf_coord function, please do! :-)
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: Help needed in remapping data using IDL [message #70786 is a reply to message #70785] |
Wed, 05 May 2010 14:24  |
kaba
Messages: 2
Registered: May 2010
|
Junior Member |
|
|
On May 5, 7:50 am, David Fanning <n...@dfanning.com> wrote:
> kaba writes:
>
>> Hello Guys, I am relative a new IDL programer who would appreciate
>> some help in how to remap a sample data from satellite projection to
>> Lambert Conformal.
>
>> Here are the detail of What I have done so far:
>> ...
>
>> When I plot this using on a window, the map projection looks right but
>> the data does not seem to be remapped correctly.
>> This is my first time using map_proj_init and map_pro_image. Please
>> help.
>
> It seems natural to want to work with latitude and
> longitude coordinates when you are working with map
> projections, but resist that thought! To do this kind
> of reprojection properly, you want to be in XY space.
> (IDL sometimes calls this UV space, but in any case
> it is a rectangular space where things make a lot of
> sense, as opposed to, say, lat/lon space where they
> don't.)
>
> Have a look at this article for additional information:
>
> http://www.dfanning.com/map_tips/warpimage.html
>
> 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.")
Ph.D Fanning
Thanks for the link, I have downloaded both the coyote and catalyst
libraries.
I have started writing a new code based on the article you send me,
and I have a question.
(1) My input file in this case is HDF 4 not netCDF. Hence I get an
error message when I call the subroutine
FindMapBoundy(). How do I get around this issue?
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
