Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  971 1853 (N57E02 L= 36)  1258    123        3 ( 920/1812-1028/1889)  16.7 E06-W04 N52-N61
      2 1691 1547 (N29W55 L= 94)   391     50        2 (1661/1532-1716/1562)   8.1 W51-W58 N28-N30
      3 1794 1443 (N22W63 L=101)  4261    240        8 (1707/1348-1877/1516) -42.3 W54-W76 N16-N27
      4  654 1357 (N19E23 L= 15)   419    123        0 ( 606/1303- 679/1398) -56.3 E26-E21 N15-N21
      5  414 1343 (N19E41 L=357)   274    104        0 ( 388/1297- 440/1386)  59.0 E42-E39 N16-N21
      6  226 1096 (N05E56 L=342)  1693    259        4 ( 145/ 992- 271/1209)  45.0 E66-E51 N00-N11
      7  225  666 (S21E65 L=333)  2848    272        8 ( 171/ 560- 273/ 755) -81.3 E72-E60 S27-S14
      8 1833  695 (S24W64 L=102)   723    116        0 (1805/ 642-1858/ 752) -81.9 W62-W66 S27-S20
      9 1542  621 (S29W36 L= 74)   319    104        0 (1493/ 591-1589/ 646) -16.7 W31-W41 S30-S26
     10  929  630 (S27E08 L= 31)   229     93        0 ( 883/ 616- 973/ 642)  -9.5 E10-E04 S27-S25
     11 1398  525 (S35W26 L= 64)  1013    296        0 (1340/ 409-1482/ 627)  65.0 W23-W31 S43-S28
     12  784  286 (S51E27 L= 11)   812    117        1 ( 730/ 257- 828/ 330)  17.1 E34-E21 S53-S47
  

    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