Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  560 1785 (N31E66 L= 75)  1731    114        1 ( 536/1723- 589/1842) -29.1 E79-E59 N26-N36
      2 1499 1702 (N56W18 L=159)  1636    116        3 (1453/1653-1540/1749) -58.0 W12-W23 N50-N61
      3 1204 1659 (N40E07 L=134)   537    115        1 (1145/1628-1250/1678)   0.0 E11-E03 N38-N41
      4  899 1479 (N19E20 L=122)  1254    270        1 ( 845/1388- 963/1597) -35.5 E27-E13 N14-N25
      5 1614 1572 (N49W31 L=172)   424     33        0 (1593/1552-1641/1587)  33.7 W28-W34 N46-N50
      6  358 1534 (N10E62 L= 79)  1359    169        2 ( 301/1484- 418/1581)  58.0 E66-E57 N06-N14
      7  403 1490 (N09E55 L= 87)  1034    161        2 ( 360/1451- 444/1565)  82.9 E57-E52 N05-N13
      8 1329 1530 (N35W05 L=147)   406     68        1 (1301/1507-1363/1559) -21.8 W03-W08 N32-N36
      9 1658 1167 (N21W37 L=178)   437     78        0 (1612/1148-1700/1186)   7.1 W32-W40 N20-N22
     10  396  832 (S30E38 L=104)   327     68        0 ( 367/ 816- 431/ 855)   0.0 E40-E34 S31-S29
     11  713  810 (S24E13 L=128)   357    147        0 ( 666/ 793- 784/ 836)  23.2 E16-E09 S25-S20
     12  738  762 (S26E11 L=131)   708    318        0 ( 632/ 717- 870/ 797)  31.0 E17-E02 S31-S21
     13 1421  572 (S19W38 L=179)   447     47        0 (1400/ 545-1444/ 605) -36.9 W35-W40 S20-S16
     14 1641  584 (S12W54 L=195)   371     27        0 (1624/ 562-1663/ 603)  45.0 W51-W55 S13-S10
     15  872  532 (S37W05 L=146)   392     49        1 ( 841/ 512- 912/ 556)  31.0 W02-W07 S38-S35
     16 1485  426 (S25W51 L=193)   562     98        0 (1437/ 409-1536/ 447)  11.3 W46-W56 S25-S23
     17  954  290 (S52W24 L=165)   484     70        1 ( 930/ 262- 974/ 327) -39.8 W20-W27 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