Re: Matching Lats and Lons from two arrays [message #62142 is a reply to message #62141] |
Tue, 26 August 2008 10:54   |
pgrigis
Messages: 436 Registered: September 2007
|
Senior Member |
|
|
Jean H wrote:
> pgrigis@gmail.com wrote:
>> Would something like this work for sorted x (plus some fix for first
>> and last element)?
>> There's going to be an overhead for sorting x if not already sorted
>> however.
>>
>> x=[3,3.5,4,6.5,7]
>> y=[3.4, 3.0, 6.8, 6.3]
>>
>> a=value_locate(x,y)
>> result=a+( (y-x[a]) GT (x[a+1]-y))
>> print,result
>>
>>
>> Ciao,
>> Paolo
>
> this would work if you only have 1 coordinate (latitude), not with 2
> (lat,long)...
> Jean
Yes, this is not going to help the original poster.
But I think Brian's routine below is 1-dimensional.
Ciao,
Paolo
>
>>
>> Brian Larsen wrote:
>>> I have tried this on several occasions (for a little different
>>> application but I think its the same) and have had no luck eliminating
>>> the for loop, so I just wrote it in a function to hide it from
>>> myself. This is my try at this based on value locate:
>>> http://people.bu.edu/balarsen/Home/IDL/Entries/2008/1/7_roun d2array_(updated_7Jan2008).html
>>>
>>> If others know how to eliminate the for loop that would be fantastic.
>>>
>>>
>>> Cheers,
>>>
>>> Brian
>>>
>>> ------------------------------------------------------------ --------------
>>> Brian Larsen
>>> Boston University
>>> Center for Space Physics
>>> http://people.bu.edu/balarsen/Home/IDL
|
|
|