Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1346 1425 (N19W31 L= 52)   645    211        0 (1254/1403-1460/1442) -36.0 W25-W37 N14-N22
      2 1754 1362 (N02W58 L= 79)   517    153        0 (1720/1327-1792/1393)   5.7 W53-W63 N01-N03
      3 1683 1311 (N02W49 L= 71)   354     67        0 (1666/1272-1704/1337)  45.0 W46-W51 N00-N03
      4 1473 1037 (S06W26 L= 47)   862    361        2 (1392/ 974-1623/1151) -37.6 W21-W35 S11-N01
      5 1767  960 (S20W45 L= 67)   531     98        2 (1734/ 921-1792/1001)  15.9 W41-W48 S20-S17
      6 1885  868 (S30W55 L= 77)  1724    307        1 (1813/ 807-1942/ 955)  17.1 W47-W63 S32-S26
      7 1585  815 (S22W28 L= 49)   747    216        1 (1558/ 724-1616/ 942)  59.7 W24-W32 S27-S15
      8  397  892 (N14E42 L=339)   282    135        0 ( 354/ 876- 454/ 926) -29.7 E44-E37 N12-N16
      9  516  743 (N03E38 L=344)   383    105        1 ( 478/ 718- 566/ 778) -51.3 E40-E34 N00-N05
     10 1249  719 (S17W04 L= 26)   197     50        1 (1224/ 692-1268/ 746)  71.6 W02-W05 S18-S15
     11  367  479 (S10E62 L=320)  1049    196        5 ( 339/ 435- 398/ 532)  71.6 E66-E57 S12-S06
  

    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