Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  741 1600 (N19E35 L=134)  1085    219        2 ( 692/1546- 807/1660) -10.3 E40-E29 N15-N22
      2 1843 1123 (N25W52 L=221)   507     85        1 (1815/1088-1866/1153) -31.0 W48-W54 N23-N26
      3 1733  858 (N06W49 L=218)   296     73        0 (1696/ 846-1761/ 872)  18.4 W45-W51 N05-N07
      4 1050  802 (S17W08 L=177)   663    192        0 ( 963/ 753-1151/ 863)  57.5 W04-W11 S22-S11
      5 1739  716 (S01W55 L=224)   940    211        2 (1710/ 619-1764/ 818) -61.4 W51-W59 S06-N05
      6 1247  620 (S21W25 L=194)   716    290        1 (1144/ 563-1356/ 714)  -3.4 W16-W33 S22-S18
      7 1704  582 (S09W59 L=228)   509     67        1 (1673/ 557-1727/ 607)  18.4 W55-W61 S09-S06
      8  783  500 (S43W02 L=170)   827    236        2 ( 673/ 474- 895/ 525)  26.6 E07-W09 S47-S39
      9 1464  323 (S30W56 L=225)  1441    346        4 (1378/ 292-1570/ 365)  19.4 W47-W64 S32-S25
  

    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