1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
#{ NEDOVERLAY -- Overlay NED catalog sources on an image display.
procedure nedoverlay (image)
string image { prompt = "Input image" }
bool append = yes { prompt = "Append display?" }
real size = 0.25 { prompt = "Field size" }
bool galaxies = no { prompt = "Mark Galaxies?" }
bool radios = yes { prompt = "Mark Radio sources?" }
bool xrays = yes { prompt = "Mark X-Ray sources?" }
int mkcolor = 208 { prompt = "Marker color" }
int galcolor = 207 { prompt = "Galaxy marker color" }
bool verbose = no { prompt = "Verbose output?" }
int status = 0 { prompt = "Service status code" }
begin
string img, ned, coords, tvcoords, url, query
string fields, pos_all, pos_l, pos_w, base
real ra, dec, sz
bool verb, app, gal, rad, xry
int mcol, gcol, nobjs
img = image # Get params to local variables
verb = verbose
app = append
mcol = mkcolor
gal = galaxies
rad = radios
xry = xrays
gcol = galcolor
sz = size
base = "http://ned.ipac.caltech.edu/cgi-bin/nph-NEDobjsearch"
base = base // "?search_type=Near+Position+Search&of=xml_main&"
if (imaccess (img)) {
iferr { wcsinfo (img) } then {
error (0, "Cannot determine image coords for `"//img//"'")
} else {
ra = wcsinfo.midx
dec = wcsinfo.midy
sz = max (wcsinfo.width, wcsinfo.height) / 60.0
if (app == no)
display (img, 1, >& "dev$null")
}
} else {
sesame (img, verbose-)
if (app == no)
dss (img, size=sz, use_display+)
ra = sesame.ra
dec = sesame.dec
img = "cache$" // img // ".fits"
}
# Get NED sources
printf ("RA=%.4f&DEC=%.4f&SR=%.4f\n", ra, dec, sz) | scan (query)
url = base // query // "&RUNID=" // vo.runid
ned = mktemp ("/tmp/ned")
# Download and convert the file.
votcopy (url, ned, format="ascii")
# Create temp files for the output
coords = mktemp ("tmp$ned_coords_")
tvcoords = mktemp ("tmp$ned_tvcoords_")
pos_all = mktemp ("tmp$pos_all_")
pos_l = mktemp ("tmp$pos_l_")
pos_w = mktemp ("tmp$pos_w_")
# Select the RA,Dec from the output NED table.
fields (ned, "3,4,5", >& pos_all)
# Expand the list of object types we want to mark.
fields = mktemp ("tmp$ned_fields_")
print ("G_Lens\nGClstr\nGPair\nQSO\n", > fields)
if (rad) print ("Radio\n", >> fields)
if (xry) print ("Xray\n", >> fields)
# Mark the Galaxies on the display.
if (galaxies) {
match (" G ", ned) | fields("STDIN","3,4,5", >& pos_all)
count (pos_all) | scan (nobjs)
wcsctran (pos_all, "c1", img, verb-,
inwcs="world", outwcs="logical", units="n n")
tvmark (frame=1, coords="c1", mark="plus", radii=10, color=gcol,
txsize=1, nxoffset=5, nyoffset=-10)
delete ("c1", verify-, >& "dev$null")
}
if (verbose)
printf ("%d objects found\n", nobjs)
# Mark the requested objects on the display.
list = fields
while (fscan (list, s1) != EOF) {
match (s1, ned) | fields("STDIN", "3,4,5", >& pos_w)
count(pos_w) | scan(nobjs)
if (nobjs > 0) {
if (verbose)
printf ("%-10s : %d found\n", s1, nobjs)
wcsctran (pos_w,pos_l,img,"world","logical",verb-,units="n n") | \
tvmark (frame=1, coords=pos_l, mark="circle", radii=10,
color=mcol, txsize=1, lab+, nxoffset=5, nyoffset=-10)
}
}
list = ""
# Clean up.
delete ("tmp$ned*,tmp$pos_*", verify-, >& "dev$null")
end
|
