Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1251 1826 (N57E15 L=353)  1441    180        1 (1173/1786-1331/1854)   6.3 E23-E05 N53-N60
      2  618 1733 (N28E53 L=315)  2582    252        6 ( 581/1624- 662/1837) -26.6 E74-E44 N19-N35
      3 1667 1624 (N57W40 L= 49)   901     75        2 (1643/1580-1686/1664)  38.2 W35-W52 N52-N63
      4  647 1571 (N18E40 L=329)   370     87        0 ( 619/1551- 692/1602)  -7.1 E43-E35 N16-N19
      5 1890  908 (N14W67 L= 75)   702    125        2 (1868/ 860-1916/ 950) -24.0 W62-W72 N12-N16
      6  515  862 (S26E28 L=341)   714    175        2 ( 472/ 776- 545/ 927) -58.0 E32-E24 S31-S22
      7  898  765 (S22E00 L=  8)   242     81        0 ( 881/ 734- 919/ 800) -56.3 E01-W01 S23-S19
      8  685  733 (S30E13 L=355)   295    144        0 ( 627/ 715- 752/ 755)  32.0 E16-E08 S32-S27
      9 1365  589 (S19W34 L= 42)   388    131        0 (1299/ 574-1425/ 606)  35.5 W29-W37 S21-S16
     10  625  335 (S62E09 L=359)   585     92        0 ( 598/ 309- 666/ 365)  26.6 E12-E04 S64-S58
     11 1411  281 (S35W60 L= 68)  1276    138        3 (1351/ 258-1473/ 314)  18.4 W52-W65 S38-S32
  

    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