Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  881 1796 (N43E28 L= 92)  1689    189        4 ( 810/1705- 935/1858) -28.8 E39-E19 N35-N47
      2  622 1634 (N26E39 L= 81)   706    156        1 ( 556/1579- 680/1672) -13.0 E45-E32 N22-N28
      3 1171 1629 (N34E01 L=120)   347     81        1 (1149/1602-1197/1667) -53.1 E02-W01 N32-N36
      4 1664 1592 (N44W43 L=164)   609     58        1 (1640/1559-1693/1617)  31.0 W39-W49 N40-N46
      5 1134 1492 (N24E01 L=120)   288     93        0 (1124/1438-1143/1542) -81.9 E02-W00 N20-N27
      6 1761 1321 (N25W48 L=168)  1729    372        3 (1646/1171-1844/1446) -39.5 W37-W56 N17-N31
      7  159 1051 (S14E66 L= 55)  8915    616       16 ( 103/ 808- 223/1315)  76.4 E77-E57 S28-N01
      8  264  887 (S24E53 L= 67)  2389    464        5 ( 235/ 659- 288/1082)  81.3 E58-E50 S38-S11
      9 1714  945 (N01W46 L=166)   577    203        0 (1643/ 912-1789/ 988)  22.6 W39-W51 N00-N05
     10  705  841 (S22E17 L=103)   224     53        2 ( 683/ 817- 732/ 874) -45.0 E18-E15 S23-S19
     11 1694  835 (S05W47 L=167)   644    196        0 (1614/ 777-1758/ 864)   0.0 W40-W54 S07-S03
     12 1658  630 (S18W51 L=172)   501     80        1 (1620/ 600-1689/ 665)  68.2 W48-W53 S19-S14
     13 1657  536 (S23W56 L=177)   625     77        1 (1617/ 511-1697/ 566)   0.0 W51-W61 S23-S21
     14 1054  484 (S39W13 L=133)   298     49        0 (1035/ 453-1072/ 516) -53.1 W11-W14 S40-S36
     15 1521  432 (S32W51 L=171)   350     83        0 (1486/ 418-1563/ 452)  33.7 W47-W53 S33-S29
     16 1629  408 (S31W64 L=184)   447     56        2 (1603/ 383-1650/ 429) -56.3 W59-W67 S32-S29
     17 1494  360 (S37W54 L=175)   357     43        1 (1475/ 340-1514/ 379)  90.0 W51-W56 S38-S35
  

    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