Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  641 1820 (N33E54 L=121)  1592    263        3 ( 598/1778- 673/1875) -36.9 E60-E48 N29-N38
      2 1321 1654 (N40W00 L=176)   382     46        0 (1293/1624-1348/1683) -21.8 E02-W03 N37-N41
      3  667 1627 (N19E39 L=136)   527     84        1 ( 631/1600- 713/1654)  31.0 E41-E36 N16-N21
      4  758 1516 (N14E29 L=146)   513    176        0 ( 730/1445- 787/1593) -45.0 E33-E25 N10-N18
      5  917  978 (S12E05 L=170)   177     98        0 ( 867/ 965- 961/ 998)  31.0 E07-E02 S13-S10
      6 1533  980 (N04W30 L=206)   346    104        0 (1487/ 961-1585/ 997)  33.7 W27-W33 N01-N05
      7  756  812 (S26E10 L=165)   595    198        0 ( 654/ 783- 836/ 853)  38.7 E16-E05 S29-S21
      8 1628  735 (S06W45 L=220)   995    145        2 (1576/ 669-1654/ 812) -69.4 W41-W47 S09-S01
      9 1416  666 (S16W33 L=209)   755    154        2 (1341/ 630-1490/ 725)  60.9 W29-W36 S19-S10
     10 1295  587 (S24W29 L=204)   775    169        0 (1221/ 552-1370/ 628)   0.0 W23-W35 S24-S22
     11 1039  435 (S40W21 L=196)  1720    417        3 ( 861/ 404-1245/ 477)  15.4 W06-W35 S45-S35
     12 1541  376 (S27W57 L=232)  2267    390        3 (1375/ 345-1659/ 408)  31.0 W45-W67 S32-S20
  

    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