Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1086 1613 (N31E07 L=284)   537    134        1 (1041/1547-1124/1657)  77.5 E08-E04 N25-N34
      2  459 1554 (N18E48 L=243)  1267    153        1 ( 386/1528- 541/1575)  -4.1 E55-E40 N15-N19
      3 1033 1454 (N19E07 L=284)   239     58        0 (1020/1425-1054/1491) -53.1 E08-E04 N17-N21
      4 1685 1226 (N17W39 L=330)   461    105        0 (1633/1206-1744/1247)  11.3 W34-W44 N16-N18
      5  792 1204 (N00E16 L=274)   419    170        0 ( 701/1173- 858/1238)  -5.2 E22-E11 N00-N01
      6 1820 1101 (N13W52 L=343)   787    148        1 (1798/1047-1854/1178) -76.0 W50-W56 N09-N17
      7 1902 1054 (N13W63 L=353)  3752    263        8 (1850/ 942-1947/1176)  90.0 W55-W70 N06-N20
      8 1650 1073 (N08W38 L=329)   274     59        2 (1629/1034-1666/1103) -63.4 W36-W39 N05-N09
      9 1712 1042 (N07W43 L=334)   416     69        1 (1692/1002-1739/1073)  68.2 W41-W45 N04-N09
     10 1728  974 (N04W46 L=336)   269     43        1 (1714/ 946-1744/1004)  90.0 W44-W46 N02-N05
     11 1831  925 (N03W56 L=347)   349     37        1 (1808/ 899-1847/ 944) -33.7 W53-W58 N02-N04
     12  825  810 (S22E09 L=282)   225     66        0 ( 798/ 781- 851/ 846)  68.2 E10-E07 S24-S19
     13 1007  701 (S26W04 L=295)  1393    422        1 ( 831/ 653-1218/ 746)  18.4 E06-W18 S31-S21
     14 1598  502 (S25W52 L=343)  3637    519        8 (1478/ 415-1745/ 666)  -1.5 W39-W77 S29-S16
     15  914  539 (S37W02 L=292)   265     74        0 ( 878/ 525- 957/ 560)   0.0 E01-W04 S38-S36
     16  833  444 (S45E03 L=288)  1306    343        1 ( 684/ 404- 989/ 509)  -4.6 E15-W10 S47-S43
     17 1192  367 (S43W28 L=319)   973    267        1 (1080/ 323-1315/ 421)  -2.5 W17-W40 S45-S41
  

    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