Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  943 1785 (N49E27 L=253)   465     69        0 ( 915/1760- 978/1806)  80.5 E29-E25 N46-N52
      2  764 1684 (N37E35 L=245)   858    167        0 ( 725/1620- 821/1751) -26.6 E43-E27 N33-N41
      3 1046 1697 (N44E14 L=266)   547     94        0 (1000/1663-1084/1725) -12.5 E18-E09 N41-N45
      4  527 1626 (N28E51 L=229)   785    153        0 ( 464/1593- 585/1657)   0.0 E58-E43 N26-N30
      5 1209 1357 (N23W06 L=286)   605    263        1 (1109/1286-1295/1434) -21.0 E01-W12 N19-N26
      6 1662 1269 (N25W41 L=321)   285     60        0 (1635/1251-1694/1291) -21.8 W38-W43 N24-N26
      7  441 1154 (S02E40 L=240)   376     62        1 ( 415/1126- 467/1180) -21.8 E42-E37 S03-N00
      8 1701  643 (S12W56 L=335)   575     85        1 (1671/ 607-1720/ 682)  78.7 W53-W58 S14-S09
      9 1102  522 (S30W15 L=295)   316     72        1 (1071/ 504-1142/ 536)  26.6 W12-W17 S31-S29
  

    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