Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1579 1581 (N19W54 L=328)   459     93        0 (1553/1541-1601/1613)  53.1 W51-W56 N17-N21
      2 1736 1566 (N13W69 L=343)   586    173        0 (1714/1541-1770/1585)  -5.7 W65-W75 N11-N13
      3  715 1466 (N39E07 L=266)   372    116        1 ( 658/1436- 780/1487) -12.5 E11-E02 N36-N40
      4 1429 1452 (N17W37 L=310)   526    105        2 (1388/1425-1488/1481) -51.3 W34-W39 N14-N19
      5 1304 1241 (N10W22 L=295)  1896    693        4 ( 970/1131-1575/1290) -11.3 W00-W40 N02-N14
      6  773 1104 (N16E12 L=261)   315    131        0 ( 719/1083- 841/1147) -45.0 E14-E08 N13-N19
      7  290  941 (N18E50 L=224)  3092    654        6 ( 114/ 887- 491/1013) -18.4 E67-E34 N11-N24
      8 1555  746 (S25W25 L=298)  2226    437        7 (1426/ 582-1638/ 913)  31.8 W12-W33 S33-S17
      9  827  488 (S20E27 L=246)   550    117        2 ( 768/ 452- 886/ 510) -23.2 E30-E23 S23-S17
     10  353  466 (S11E64 L=209)  1055    137        3 ( 307/ 422- 396/ 502) -50.2 E68-E60 S14-S08
     11  764  398 (S24E35 L=238)   658    111        2 ( 726/ 362- 807/ 448) -81.9 E37-E33 S27-S20
     12  532  303 (S25E58 L=215)  2606    406        4 ( 469/ 188- 604/ 416)  68.2 E64-E54 S36-S16
  

    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