Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  786 1803 (N45E32 L=260)  1088    101        3 ( 736/1778- 845/1829)  -4.8 E37-E25 N43-N48
      2 1126 1703 (N40E00 L=292)   300     63        0 (1090/1680-1163/1721)  45.0 E03-W02 N38-N42
      3 1615 1620 (N42W43 L=336)   469    100        2 (1575/1591-1668/1643) -12.5 W39-W48 N40-N43
      4  449 1483 (N19E45 L=248)  1400    257        0 ( 382/1397- 543/1579) -31.0 E56-E35 N13-N25
      5  719 1526 (N23E25 L=267)   317    115        0 ( 699/1473- 740/1572) -54.5 E28-E22 N19-N26
      6 1703 1290 (N19W44 L=336)   718    196        0 (1639/1242-1788/1345) -13.0 W38-W51 N17-N21
      7 1871  876 (S03W64 L=356)   596     66        1 (1847/ 848-1904/ 907)  23.2 W61-W68 S05-S01
      8 1128  834 (S16W08 L=301)   193     54        0 (1101/ 817-1160/ 851) -14.0 W06-W10 S16-S14
      9  638  694 (S29E23 L=269)  1096    345        0 ( 563/ 555- 741/ 831)  53.6 E30-E15 S38-S19
     10 1766  590 (S22W64 L=356)  1187    152        1 (1702/ 544-1826/ 620)  26.6 W58-W73 S25-S19
  

    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