match_2d [message #70781] |
Thu, 06 May 2010 01:29 |
d.poreh
Messages: 406 Registered: October 2007
|
Senior Member |
|
|
Hello everyone,
I have a problem that I can’t fix.
I am using match_2d from J.D.Smith's library to match between two
catalogues. Below is what i am doing:
IDL> print, laser
13.51132 78.08015 2.244529
13.5041 78.07921 2.416421
13.49686 78.07828 2.420123
13.48957 78.07736 2.668057
13.48225 78.07644 2.652986
13.4749 78.07554 2.604654
13.46749 78.07465 2.373039
13.46009 78.07375 2.699916
…………………
IDL> print, sat
13.50834 78.08333 5.965842
13.50001 78.08333 3.459386
13.50001 78.075 3.092191
13.49168 78.075 3.444704
13.48334 78.075 3.364458
13.47501 78.075 3.108061
13.46668 78.075 3.121016
13.45834 78.075 3.017546
13.45001 78.075 3.15583
………….
What I want to do is find the closest data from array laser to each of
array sat and then average them (of course for some points of sat I do
not have close laser data). At the end I want to plot sat data and
averaged laser to comparison.
x1=laser[0,*]
y1=laser[1,*]
x2=sat[0,*]
y2=sat[1,*]
match=match_2d(x1,y1,x2,y2,0.008,MATCH_DISTANCE=md)
Any help highly appreciated
Cheers
Dave
|
|
|