Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  871 1634 (N34E25 L=195)   760    150        1 ( 819/1589- 917/1664)  45.0 E28-E20 N30-N36
      2 1261 1438 (N29W07 L=228)  2024    605        6 (1086/1189-1421/1590) -33.2 E08-W21 N16-N35
      3 1143 1213 (N12W03 L=224)   659    233        0 (1038/1182-1240/1253)   0.0 E03-W09 N10-N13
      4 1753 1083 (N16W49 L=270)   701     98        2 (1733/1035-1778/1126)  56.3 W47-W51 N13-N19
      5  243  946 (S20E56 L=165)  3282    423        2 ( 123/ 806- 354/1083)  29.4 E84-E45 S31-S09
      6 1763  872 (N04W54 L=274)   706    138        2 (1737/ 799-1797/ 923)  81.9 W51-W56 N00-N07
      7 1080  728 (S18W09 L=230)   303     46        1 (1053/ 704-1104/ 752) -33.7 W07-W10 S19-S16
      8 1630  667 (S11W48 L=269)  1526    251        3 (1517/ 623-1720/ 709)  12.5 W40-W58 S14-S07
      9 1631  568 (S16W53 L=274)   385    121        0 (1576/ 540-1655/ 609)  59.0 W49-W55 S19-S13
     10  772  295 (S57E02 L=218)  1357    270        2 ( 637/ 272- 883/ 321)  19.0 E20-W09 S64-S52
  

    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