Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1597 1575 (N19W56 L=223)  3414    720        7 (1371/1536-1781/1604) -22.8 W40-W79 N11-N27
      2 1490 1349 (N11W37 L=204)   316    104        0 (1475/1293-1503/1402)  63.4 W35-W39 N07-N13
      3 1908 1178 (S13W67 L=234)  2971    364        8 (1866/1090-1951/1263) -90.0 W61-W74 S19-S07
      4  579 1171 (N26E23 L=143)   798    243        1 ( 490/1070- 649/1244) -63.4 E28-E19 N18-N32
      5 1828  956 (S22W51 L=218)  3011    358        2 (1774/ 794-1860/1118)  81.9 W45-W54 S31-S10
      6  134  931 (N21E67 L= 99)   708    125        1 ( 124/ 866- 142/ 988)  63.4 E70-E64 N16-N24
      7 1681  778 (S27W35 L=201)  1079    277        0 (1657/ 675-1731/ 886)  81.9 W31-W38 S35-S20
      8 1365  738 (S19W12 L=178)   232     53        0 (1336/ 708-1389/ 764) -68.2 W10-W12 S21-S16
      9 1057  627 (S17E09 L=158)   859    281        1 ( 968/ 502-1132/ 742)  33.7 E18-E02 S21-S10
     10  564  605 (S06E39 L=128)  1082    171        4 ( 512/ 533- 626/ 677) -79.7 E41-E36 S11-N00
     11 1191  448 (S32E07 L=160)   434     58        1 (1160/ 421-1224/ 475) -36.9 E08-E04 S33-S30
     12  996  341 (S34E24 L=143)  1280    241        1 ( 924/ 263-1079/ 406)  19.4 E32-E15 S37-S29
  

    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