Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1098 1834 (N52W42 L=263)   434     59        0 (1081/1799-1114/1866)  36.9 W37-W46 N48-N54
      2  500 1531 (N49E26 L=196)   602     73        2 ( 468/1494- 535/1562)   0.0 E29-E22 N46-N51
      3 1026 1506 (N32W14 L=236)   654    147        1 ( 953/1481-1099/1528) -32.0 W09-W17 N28-N34
      4  709 1388 (N33E10 L=212)   348    112        1 ( 660/1355- 756/1427)  14.0 E14-E06 N31-N34
      5  651 1192 (N23E18 L=204)   285    137        0 ( 586/1176- 716/1201) -26.6 E22-E14 N20-N24
      6  165  918 (N17E66 L=156)   857     96        1 ( 143/ 870- 181/ 962)  59.0 E68-E62 N13-N19
      7 1772  853 (S26W46 L=268)   870    198        2 (1727/ 765-1805/ 936)  41.6 W40-W50 S30-S21
      8  719  299 (S32E45 L=177)   355     51        1 ( 701/ 269- 735/ 326)  51.3 E47-E42 S34-S29
      9 1450  284 (S59W11 L=233)  1886    204        4 (1353/ 250-1542/ 320) -14.0 W00-W24 S62-S53
  

    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