Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  555 1615 (N36E45 L= 66)   436     91        0 ( 516/1575- 583/1651)  71.6 E46-E43 N32-N38
      2 1335 1544 (N38W19 L=130)   803    235        1 (1239/1465-1408/1597) -31.6 W11-W24 N32-N41
      3  814 1465 (N28E18 L= 93)  1092    190        5 ( 754/1392- 885/1538) -33.7 E24-E12 N23-N32
      4  288 1392 (N18E61 L= 50)   477     58        2 ( 270/1365- 304/1429) -76.0 E62-E58 N16-N20
      5 1113 1356 (N23W03 L=114)   309    126        0 (1090/1298-1135/1410)  90.0 W01-W04 N19-N26
      6 1414 1331 (N24W24 L=135)   359    129        1 (1391/1278-1435/1399) -82.9 W22-W25 N20-N28
      7  258 1074 (S02E55 L= 56)   520    101        0 ( 242/1032- 282/1114)  51.3 E56-E52 S04-N01
      8 1717 1025 (N06W48 L=159)   372     67        2 (1693/ 997-1748/1058) -38.7 W45-W50 N04-N08
      9 1026  698 (S18W03 L=114)  2836    726        7 ( 744/ 538-1247/ 933)  44.0 E15-W14 S31-S02
     10  559  829 (S14E29 L= 82)   249     71        0 ( 526/ 814- 597/ 843) -14.0 E31-E26 S14-S12
     11 1666  732 (S12W47 L=158)   335     59        1 (1643/ 707-1690/ 756)  56.3 W45-W49 S13-S10
     12  229  653 (S29E65 L= 46)  1358    123        1 ( 199/ 595- 265/ 700)  20.9 E79-E57 S33-S25
     13 1622  582 (S22W48 L=159)  1089    179        4 (1549/ 549-1691/ 631)  28.6 W41-W53 S24-S18
     14  481  458 (S41E41 L= 70)  1240    147        1 ( 420/ 411- 538/ 510)  25.3 E52-E33 S45-S35
     15 1700  441 (S31W68 L=179)  1026     56        3 (1676/ 414-1725/ 475) -45.0 W63-W75 S33-S28
     16 1608  389 (S36W60 L=171)  1317    150        5 (1582/ 322-1635/ 465) -22.2 W51-W78 S41-S30
     17  683  314 (S51E26 L= 85)   506     93        0 ( 655/ 290- 720/ 337)  29.1 E29-E20 S53-S48
  

    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