Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  730 1708 (N58W03 L=338)  2140    328        4 ( 600/1591- 847/1801)  11.6 E13-W21 N52-N61
      2 1501 1728 (N30W63 L= 38)  1102    123        3 (1449/1695-1561/1757) -50.2 W57-W67 N26-N33
      3 1589 1702 (N25W69 L= 44)   764     42        1 (1571/1675-1614/1725)   5.7 W64-W74 N23-N26
      4 1609 1633 (N21W63 L= 38)  1768    291        5 (1464/1605-1719/1655) -17.2 W48-W78 N15-N25
      5 1663 1564 (N16W63 L= 38)   854    210        0 (1570/1522-1753/1597) -37.6 W56-W70 N10-N20
      6  750 1327 (N30E08 L=327)   421    101        0 ( 695/1292- 789/1365)   7.1 E12-E04 N29-N31
      7  973 1224 (N18W03 L=337)   189     62        0 ( 953/1200-1002/1245)  14.0 W00-W04 N16-N19
      8 1824  842 (S29W51 L= 26)   421     36        0 (1801/ 826-1853/ 860) -33.7 W48-W54 S30-S26
      9 1421  758 (S20W16 L=351)   453    123        0 (1404/ 701-1437/ 827)  60.3 W13-W18 S23-S16
     10 1246  631 (S22W02 L=337)   311    125        1 (1213/ 582-1272/ 693)  39.8 E01-W05 S24-S19
     11 1773  617 (S43W46 L= 21)   428     45        0 (1752/ 588-1792/ 638)  36.9 W43-W49 S44-S40
     12 1356  539 (S32W06 L=341)   336     54        0 (1336/ 511-1374/ 573)  53.1 W04-W08 S33-S29
     13 1340  409 (S41W02 L=337)   699    140        0 (1319/ 365-1381/ 476) -84.3 W00-W04 S45-S35
     14  912  401 (S28E26 L=309)   649    173        1 ( 841/ 361- 989/ 441) -60.9 E28-E23 S32-S23
     15  542  314 (S24E59 L=276)   661    112        0 ( 496/ 279- 578/ 345) -80.5 E61-E56 S27-S21
     16  692  286 (S30E49 L=286)   433     92        1 ( 672/ 253- 708/ 325)  35.5 E52-E45 S32-S26
     17  694  227 (S34E55 L=280)   502     83        1 ( 647/ 214- 736/ 238) -20.6 E58-E50 S35-S32
  

    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