Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1349 1602 (N36W09 L=245)   952    227        2 (1285/1508-1405/1680) -39.8 W02-W14 N30-N40
      2  733 1583 (N21E31 L=205)   833    169        3 ( 713/1500- 754/1663) -57.5 E34-E27 N15-N26
      3 1879 1100 (N18W59 L=295)   482     60        1 (1850/1074-1909/1115)   6.3 W55-W64 N16-N19
      4  224  941 (S25E57 L=179)  1663    312        3 ( 173/ 853- 334/1071) -38.2 E63-E45 S28-S17
      5  358  762 (S34E44 L=192)   274     62        1 ( 322/ 750- 387/ 777)  -8.1 E47-E40 S35-S32
      6  370  493 (S53E51 L=185)  1104    136        2 ( 330/ 464- 442/ 543)  10.3 E61-E39 S55-S49
      7 1226  323 (S44W38 L=273)   863    201        1 (1172/ 260-1298/ 385) -12.5 W31-W49 S46-S39
      8 1520  329 (S35W62 L=298)   899    148        2 (1468/ 306-1585/ 354)  29.1 W57-W68 S37-S31
  

    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