Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  992 1762 (N38E25 L=271)   748    115        1 ( 951/1738-1049/1797)  -6.3 E29-E20 N35-N39
      2  374 1618 (N15E63 L=233)  2294    243        2 ( 332/1544- 400/1689) -39.8 E68-E56 N10-N20
      3  472 1458 (N06E45 L=250)   790    181        2 ( 414/1416- 527/1505) -14.0 E52-E40 N04-N08
      4  900 1333 (N07E14 L=281)  1241    298        3 ( 866/1209- 932/1471) -61.9 E18-E09 N00-N15
      5 1707 1087 (N14W40 L=336)   891    164        1 (1666/ 997-1740/1165) -41.6 W36-W45 N10-N18
      6  481 1088 (S15E35 L=261)   389    160        0 ( 413/1070- 549/1113)   5.7 E39-E29 S15-S14
      7  717 1065 (S12E19 L=277)   225     80        0 ( 675/1039- 754/1085)  -9.5 E21-E15 S12-S11
      8  237 1018 (S23E54 L=242)  1777    349        1 ( 153/ 928- 340/1083)  38.2 E61-E44 S29-S18
      9 1307  651 (S21W27 L=322)   236     59        0 (1284/ 625-1328/ 692) -76.0 W25-W27 S22-S18
     10 1743  547 (S13W64 L=360)  1551    197        0 (1672/ 512-1814/ 608)  42.0 W58-W69 S16-S07
  

    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