Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1605 1707 (N25W73 L= 40)   944     79        2 (1569/1678-1633/1730) -45.0 W68-W80 N23-N26
      2 1087 1658 (N41W28 L=355)   252     50        0 (1058/1643-1117/1675)   0.0 W24-W30 N40-N42
      3 1439 1576 (N25W45 L= 13)   657    206        1 (1353/1509-1491/1611) -66.0 W41-W48 N19-N28
      4 1346 1488 (N23W34 L=  2)   774    288        0 (1226/1438-1439/1527) -33.7 W25-W38 N17-N26
      5 1130 1329 (N20W15 L=342)  1065    465        4 (1023/1176-1284/1468) -77.7 W11-W18 N08-N31
      6  306 1433 (N50E52 L=276)   557     66        0 ( 283/1405- 334/1461)  36.9 E56-E47 N47-N52
      7  544 1314 (N36E24 L=304)  1121    204        3 ( 462/1268- 641/1358) -29.7 E30-E16 N31-N40
      8 1631  987 (S13W35 L=  3)  1052    271        0 (1570/ 874-1682/1108)  37.6 W28-W42 S17-S07
      9  282  967 (N20E52 L=276)  1616    317        6 ( 227/ 873- 341/1084) -90.0 E57-E46 N14-N27
     10  881 1007 (N09E08 L=319)   604    213        0 ( 809/ 942- 963/1069) -59.0 E11-E04 N04-N14
     11  580  947 (N13E28 L=299)   310     69        1 ( 550/ 923- 609/ 971) -45.0 E29-E26 N11-N14
     12 1476  844 (S16W21 L=349)   213    104        0 (1423/ 827-1526/ 862) -33.7 W18-W24 S18-S13
     13 1440  728 (S22W16 L=344)   315    132        0 (1384/ 701-1493/ 768)   0.0 W12-W20 S22-S20
     14 1635  709 (S30W31 L=358)   815    102        3 (1592/ 664-1680/ 746)  23.2 W26-W34 S32-S26
     15  265  729 (N07E60 L=267)   379     70        0 ( 231/ 716- 301/ 742)  -6.3 E64-E55 N06-N07
     16 1192  582 (S23E03 L=325)   219    120        0 (1151/ 548-1218/ 631)  38.7 E05-W00 S24-S20
     17  839  333 (S30E34 L=293)   629     90        2 ( 807/ 300- 885/ 360)  12.5 E38-E29 S31-S27
     18 1089  285 (S41E21 L=306)   573    122        0 (1057/ 246-1113/ 345)  37.9 E26-E17 S44-S36
  

    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