Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  217 1293 (S03E65 L=100)   862    216        1 ( 192/1227- 246/1355) -22.6 E72-E59 S06-N00
      2 1401 1255 (N20W18 L=184)   323     73        0 (1371/1213-1422/1290) -53.1 W15-W19 N17-N21
      3  759 1147 (S01E18 L=148)   274     93        1 ( 740/1099- 783/1191) -78.7 E19-E16 S03-N02
      4  203 1045 (S17E59 L=107)   431     62        0 ( 182/1017- 225/1078)  53.1 E61-E56 S19-S15
      5 1884  882 (N09W67 L=233)  2715    345        3 (1854/ 778-1921/1038) -72.6 W61-W74 N02-N18
      6  420  968 (S18E38 L=128)   247    101        0 ( 399/ 929- 460/ 994)  53.1 E39-E34 S20-S16
      7  250  922 (S24E52 L=114)  1210    405        0 ( 175/ 824- 339/ 981)  39.1 E60-E44 S31-S18
      8  891  913 (S11E05 L=160)   255     85        1 ( 844/ 894- 942/ 930)  38.7 E07-E02 S13-S09
      9  732  803 (S21E13 L=153)   491    172        3 ( 648/ 775- 794/ 841)  33.7 E18-E09 S24-S17
     10 1475  696 (S11W35 L=201)   986    203        2 (1431/ 614-1507/ 791) -74.7 W31-W37 S16-S05
     11 1760  596 (S10W65 L=230)  1390    167        4 (1696/ 544-1796/ 665)  58.0 W59-W70 S13-S05
     12 1274  505 (S27W29 L=195)   543    124        0 (1203/ 479-1328/ 544)  49.4 W25-W31 S30-S23
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y