Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  611 1521 (N20E38 L= 43)   313     38        0 ( 585/1502- 644/1538)  33.7 E39-E35 N18-N21
      2 1423 1142 (N13W22 L=103)   552    183        1 (1353/1087-1490/1211) -20.0 W16-W27 N10-N16
      3  802  924 (S12E11 L= 70)   277    102        0 ( 790/ 875- 815/ 981) -71.6 E12-E10 S14-S08
      4 1568  834 (S02W38 L=120)   452    105        1 (1536/ 784-1593/ 891)  90.0 W36-W39 S04-N01
      5 1151  615 (S23W17 L= 98)   230     98        0 (1098/ 601-1196/ 624)   9.5 W13-W19 S24-S22
      6 1653  530 (S18W56 L=138)  1428    256        2 (1566/ 474-1713/ 591)  45.0 W50-W60 S22-S13
      7 1521  541 (S20W44 L=126)   317     65        0 (1497/ 515-1546/ 570)  68.2 W42-W45 S22-S17
      8 1712  507 (S18W64 L=145)   728    200        2 (1686/ 474-1743/ 537)  -6.3 W59-W68 S19-S15
  

    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