Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1692 1427 (N22W56 L=189)   361     46        2 (1668/1403-1714/1447) -45.0 W53-W58 N20-N23
      2 1803 1383 (N18W67 L=200)  1042    101        0 (1761/1360-1856/1401)  -6.7 W61-W78 N16-N19
      3 1867 1112 (N01W65 L=198)   724     73        2 (1838/1076-1887/1150)  36.9 W61-W69 S01-N02
      4 1165 1068 (N07W09 L=142)   461    105        2 (1139/1026-1202/1124) -90.0 W07-W11 N04-N10
      5 1359  799 (S12W19 L=152)   373    159        1 (1334/ 742-1391/ 852)  66.8 W17-W21 S15-S08
      6 1749  725 (S22W51 L=184)  1950    316        2 (1669/ 610-1832/ 832)  40.6 W43-W59 S28-S16
      7 1168  743 (S13W07 L=140)   362    151        1 (1101/ 720-1219/ 774)   0.0 W02-W09 S14-S11
      8 1560  503 (S35W37 L=170)   582     68        1 (1528/ 465-1588/ 527)  90.0 W33-W39 S38-S33
      9 1377  314 (S49W24 L=157)   397     48        2 (1356/ 295-1406/ 341)  11.3 W21-W27 S51-S45
  

    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