Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1399 1644 (N45W33 L=131)  1022    186        3 (1353/1583-1431/1688)  41.6 W29-W38 N39-N48
      2  578 1537 (N36E36 L= 63)   281     84        0 ( 559/1495- 591/1567) -51.3 E37-E33 N33-N38
      3  919 1473 (N32E08 L= 91)   331     93        2 ( 904/1429- 938/1517) -71.6 E09-E06 N29-N35
      4 1719 1460 (N30W55 L=154)   432     84        2 (1694/1438-1746/1479) -26.6 W53-W57 N28-N31
      5 1544 1312 (N21W36 L=135)   476     97        1 (1496/1274-1590/1346) -29.7 W32-W39 N18-N23
      6  764 1223 (N15E17 L= 82)   568    141        4 ( 716/1151- 811/1297) -66.0 E20-E13 N11-N20
      7 1544 1224 (N15W35 L=133)   539    157        1 (1474/1203-1607/1253)  10.3 W29-W40 N14-N17
      8 1087 1163 (N12W04 L=102)   360    103        0 (1066/1116-1116/1221) -63.4 W02-W05 N09-N15
      9 1548 1096 (N08W34 L=133)   387    130        0 (1491/1075-1601/1128) -20.6 W30-W38 N06-N09
     10  479  993 (N01E36 L= 63)   432    146        1 ( 424/ 945- 537/1034)  32.0 E39-E31 S02-N03
     11 1349  784 (S12W21 L=119)   228     69        0 (1333/ 763-1381/ 815) -36.9 W19-W23 S12-S09
     12  773  741 (S14E16 L= 83)   219     51        2 ( 757/ 715- 798/ 771)  76.0 E17-E14 S16-S12
     13  383  708 (S18E45 L= 53)   400    108        1 ( 366/ 661- 397/ 760)  63.4 E47-E43 S20-S14
     14 1567  627 (S22W39 L=137)   632     95        2 (1530/ 598-1603/ 662)  18.4 W35-W41 S24-S19
     15  703  626 (S22E22 L= 77)   614    127        0 ( 647/ 601- 763/ 656)  -7.1 E25-E17 S23-S20
     16  709  376 (S41E26 L= 73)   411     73        0 ( 682/ 355- 731/ 408)  45.0 E28-E23 S42-S38
     17  982  350 (S43E03 L= 95)  1026    208        4 ( 916/ 316-1053/ 408)  11.3 E08-W02 S45-S37
  

    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