comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Display radar data?
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Display radar data? [message #62865 is a reply to message #62861] Wed, 08 October 2008 08:33 Go to previous messageGo to previous message
pgrigis is currently offline  pgrigis
Messages: 436
Registered: September 2007
Senior Member
pgri...@gmail.com wrote:
> David Fanning wrote:
>> Dave writes:
>>
>>> I have some weather radar data, whose coordinates are distances
>>> (kilometers) and angles to the radar station. So how can I convert
>>> these local coordinates to longitude/latitude cordinates? Please give
>>> any clue.
>>
>> Set up a map projection of the area of interest with MAP_SET
>> or MAP_PROJ_INIT. Then use CV_COORD to convert your data
>> from polar coordinates to rectangular coordinates (which
>> will be the data or lat/lon coordinates, in this case).
>
> While I am sure that this strategy would work, that will be depriving
> the original
> poster of the joy of figuring out how to convert polar coordinates in
> a plane tangent
> to a sphere at a given point to spherical coordinates. Now, that's a
> beautiful
> geometry problem, isn't it? Furthermore, if high-precision is not
> important,
> the problem may be simplified a bit by making the assumption that the
> earth is
> locally flat at the radar position, thus reducing the problem to
> converting plane
> coordinates from polar to cartesian...

In the latter approximation, a back of the envelope computation
yielded (you may want to double check this before actually using it)

lon=lon0+r/R*sin(alpha)*!radeg
lat=lat0+r/R*cos(alpha)/cos(lat0)*!radeg


Cheers,
Paolo

>
> Ciao,
> Paolo
>
>
>>
>> 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.")
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: segmentation fault with ncdf_open
Next Topic: Open file from full address in IDL Workbench

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Sun Oct 12 04:46:44 PDT 2025

Total time taken to generate the page: 2.31943 seconds