Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  386 1613 (N46E63 L= 52)   588     62        1 ( 372/1584- 404/1643) -21.8 E69-E58 N43-N47
      2  884 1544 (N39E07 L=108)   430    133        1 ( 808/1522- 941/1561)  -5.2 E13-E02 N37-N40
      3 1059 1441 (N31W05 L=120)   420    110        0 (1020/1404-1098/1483) -45.0 W02-W07 N28-N33
      4  507 1336 (N26E35 L= 80)   783    257        2 ( 466/1258- 553/1412) -63.4 E39-E30 N21-N31
      5  369 1311 (N25E48 L= 68)   624    112        0 ( 355/1252- 386/1368) -69.4 E49-E45 N20-N28
      6  276 1240 (N20E56 L= 60)   857    119        4 ( 251/1182- 302/1304) -76.0 E58-E52 N16-N24
      7 1715 1212 (N11W50 L=165)   580    116        2 (1684/1157-1756/1253) -63.4 W47-W52 N07-N13
      8  232 1173 (N16E59 L= 56)   501    121        2 ( 218/1124- 251/1227)  81.9 E61-E57 N12-N19
      9 1729  945 (S06W48 L=163)   603    183        2 (1687/ 866-1765/ 995) -76.0 W44-W51 S10-S02
     10 1714  751 (S18W48 L=163)  2956    326       10 (1616/ 548-1815/ 925)  67.4 W40-W57 S31-S07
     11  262  871 (S03E56 L= 60)   472     70        2 ( 229/ 842- 285/ 900) -45.0 E59-E53 S04-N00
     12 1610  856 (S10W38 L=153)   332    136        0 (1564/ 823-1664/ 893)  20.6 W34-W42 S11-S08
     13 1212  800 (S10W10 L=126)   274    124        0 (1178/ 758-1259/ 856)  45.0 W08-W13 S12-S07
     14 1616  776 (S15W39 L=154)   468    116        3 (1571/ 729-1663/ 805)  -8.1 W35-W42 S18-S12
     15  277  672 (S15E59 L= 56)   562    119        1 ( 234/ 622- 307/ 709) -51.3 E63-E55 S18-S12
     16  772  533 (S25E20 L= 95)   929    169        5 ( 732/ 456- 816/ 600)  55.0 E24-E16 S31-S20
     17  652  425 (S33E33 L= 83)   985    175        2 ( 569/ 401- 735/ 447)   0.0 E40-E25 S35-S31
  

    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