Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1291 1658 (N39W07 L=238)  1488    402        0 (1134/1568-1458/1747) -15.1 E06-W20 N35-N43
      2  396 1633 (N23E61 L=170)  1265    120        5 ( 359/1587- 428/1677) -49.4 E67-E56 N19-N26
      3  771 1654 (N29E29 L=202)   286     51        0 ( 741/1631- 795/1676)  56.3 E30-E28 N26-N30
      4  540  896 (S21E29 L=202)   536    177        0 ( 493/ 834- 567/ 968)  83.7 E32-E26 S25-S16
      5  922  878 (S16E04 L=227)   201     95        0 ( 900/ 828- 940/ 920)  80.5 E04-E03 S19-S13
      6  150  833 (S28E70 L=161)   807     69        0 ( 115/ 813- 186/ 851)   3.6 E79-E63 S29-S27
      7 1313  788 (S15W22 L=252)   531    148        1 (1245/ 748-1369/ 832)  49.4 W18-W24 S18-S11
      8  767  767 (S26E12 L=219)   207     93        0 ( 723/ 739- 800/ 805)  51.3 E14-E10 S27-S22
      9 1591  761 (S11W41 L=272)   730    192        0 (1512/ 721-1678/ 799)  28.6 W36-W47 S14-S07
     10 1718  636 (S15W57 L=287)   811     99        1 (1685/ 590-1746/ 677)  80.5 W53-W58 S18-S12
     11 1150  607 (S29W16 L=247)  1645    444        1 ( 970/ 541-1360/ 646)  29.5 W05-W28 S35-S22
     12  841  577 (S36E04 L=227)   372    100        0 ( 799/ 540- 881/ 617) -23.2 E07-W00 S38-S34
  

    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