Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  842 1736 (N32E32 L=315)   521     85        2 ( 809/1694- 870/1770) -26.6 E36-E28 N29-N33
      2  992 1581 (N24E16 L=331)   312     79        0 ( 974/1544-1009/1628) -51.3 E17-E13 N22-N27
      3  866 1536 (N19E22 L=325)   474    174        0 ( 787/1515- 953/1557)  16.7 E26-E16 N17-N20
      4 1589 1360 (N27W27 L= 14)   482    138        0 (1557/1297-1635/1396)  71.6 W25-W30 N21-N30
      5 1692 1294 (N26W36 L= 24)   384    102        1 (1657/1249-1719/1333) -33.7 W33-W39 N23-N27
      6 1301 1288 (N14W10 L=357)   329     82        2 (1285/1238-1327/1323) -63.4 W08-W10 N11-N16
      7  955 1228 (N03E08 L=339)   561    235        0 ( 837/1206-1059/1263)  26.6 E15-E03 N00-N06
      8  312 1197 (S11E50 L=297)   324    109        2 ( 293/1149- 341/1250)  81.9 E51-E48 S14-S07
      9 1612 1200 (N17W31 L= 18)   381    129        1 (1573/1164-1668/1233)   8.1 W28-W35 N15-N19
     10 1847  916 (N09W57 L= 45)  1739    341        2 (1817/ 777-1880/1069) -56.3 W53-W64 N01-N17
     11  203  976 (S26E59 L=289)  1547    161        5 ( 128/ 938- 280/1009)  14.0 E70-E49 S28-S23
     12 1466  728 (S12W34 L= 21)   241     44        1 (1439/ 705-1494/ 749)  33.7 W32-W35 S13-S10
     13  199  668 (S45E68 L=280)  1104    122        1 ( 139/ 639- 259/ 691)  -6.6 E84-E57 S46-S43
     14 1159  646 (S25W18 L=  5)   337     85        0 (1120/ 626-1199/ 668)  11.3 W15-W20 S26-S23
     15 1444  631 (S18W36 L= 23)   285    114        0 (1394/ 610-1489/ 657)  38.7 W32-W37 S19-S15
     16 1444  506 (S25W42 L= 29)  1839    493        2 (1253/ 458-1689/ 577)  36.2 W31-W57 S32-S13
     17 1770  486 (S15W74 L= 61)   744     67        1 (1736/ 469-1805/ 508)   8.7 W68-W81 S16-S12
  

    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