Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1757 1466 (N33W50 L=100)   454     53        1 (1737/1435-1772/1495)  76.0 W47-W51 N31-N35
      2 1084 1470 (N22E02 L= 49)   260     34        1 (1055/1454-1110/1491)  18.4 E03-W00 N20-N23
      3  156 1349 (N08E73 L=338)  1019    188        2 ( 116/1288- 188/1398)  50.2 E80-E68 N04-N10
      4 1077 1232 (N07W01 L= 51)   198     36        0 (1058/1206-1100/1251)  71.6 W00-W02 N05-N08
      5  741  852 (S19E16 L= 35)   201     59        0 ( 725/ 827- 765/ 886) -53.1 E16-E13 S20-S16
      6 1798  800 (S07W57 L=108)  1604    250        2 (1721/ 758-1866/ 845)   7.6 W49-W64 S09-S04
      7  540  780 (S25E29 L= 21)   513    142        0 ( 516/ 720- 589/ 826)  63.4 E31-E25 S29-S21
      8  599  721 (S28E25 L= 26)   223    103        0 ( 555/ 710- 649/ 733)   8.1 E28-E21 S29-S27
      9 1619  542 (S26W50 L=100)   580     70        0 (1587/ 504-1647/ 589) -80.5 W47-W51 S28-S22
  

    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