Hi
Thanks for any suggestions.
Below is a procedure for reading a structures in the arguments passed to
procedure.
This is my first attempt at writing anything like this and it took me a long
time to write and plot it the way i want.
I have four nights of data on some magnititudes of stars (v1, v2, .., c1,
c2, ...) each of these are fltarr, the star1, 2,.. are string arrays of
names of stars. my goal is to match the star1 to star2 and then where it
matches the star name, put those mags into fltarr. then again match star1 to
star3 and put mags into fltarr where the starname matches and the same for
star1 to star4. then plot those arrays for each night.
What i am doing is trying to find out what kind of differences there are in
magnitudes of stars through a fish eye lens and see if there differences are
consistent on images that i take.
Ok, i can plot those differences that match up to each night, but now i want
to plot a line to each point for each star. that way i can tell if there is
a difference is all over the plot or linear or ?
any suggestions? if you have and my explanation is confusing please let me
know.
thanks alot,
dennis lamenti
pro starplot2, star1, star2, star3,star4, v1, v2, v3, v4, c1, c2, c3, c4,
cv1, cv2, cv3, cv4
cv1 = fltarr(n_elements(star1))
cv2 = fltarr(n_elements(star1))
cv3 = fltarr(n_elements(star3))
cv4 = fltarr(n_elements(star1))
x = intarr(n_elements(star1) + 1)
y = intarr(n_elements(star1) + 1)
z = intarr(n_elements(star1) +1)
h = intarr(n_elements(star1) + 1)
;pline = fltarr(n_elements(star1))
for i=0, n_elements(star1)-1 do begin
cv1[i] = c1[i] - v1[i]
x[i] = 1.0
endfor
for i=0, n_elements(star1)-1 do begin
t=where(star1[i] EQ star2, count)
if count gt 0 then begin
cv2[i] = c2[t] - v2[t]
y[i] = 2.0
endif
endfor
for i=0, n_elements(star1)-1 do begin
t=where(star1[i] EQ star3, count)
if count gt 0 then begin
cv3[i] = c3[t] - v3[t]
z[i] = 3.0
endif
endfor
for i=0, n_elements(star1)-1 do begin
t=where(star1[i] EQ star4, count)
if count gt 0 then begin
cv4[i] = c4[t] - v4[t]
h[i] = 4.0
endif
endfor
plot, x, cv1, psym=3,$
TITLE = 'test1',$
XTITLE = 'night', $
YTITLE = 'c-v', XRANGE= [0,5]
oplot, y, cv2, psym=4
oplot, z, cv3, psym=5
oplot, h, cv4, psym=2
;for i=0, n_elements(star1)-1 do begin
; pline[
end
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
