Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1594  781 (S20W37 L=305)   413     90        3 (1554/ 749-1626/ 807)  18.4 W33-W39 S21-S18
      2 1670  734 (S23W44 L=312)   312     67        1 (1646/ 701-1699/ 760)  56.3 W41-W46 S25-S21
      3  192  652 (S18E71 L=197)  1439    106        2 ( 166/ 604- 217/ 711)  56.3 E77-E66 S20-S14
      4 1707  640 (S29W50 L=318)   497     41        1 (1677/ 618-1733/ 663)  -9.5 W46-W52 S30-S27
      5 1166  532 (S34W06 L=275)   535    116        1 (1117/ 505-1216/ 576)  23.2 W02-W10 S35-S30
      6 1389  499 (S37W24 L=292)  1091    278        0 (1285/ 440-1517/ 573)  24.0 W15-W33 S41-S32
      7 1020  508 (S34E04 L=264)   255     60        0 ( 987/ 490-1052/ 529) -11.3 E06-E01 S35-S32
      8  867  203 (S58E28 L=240)   486     57        1 ( 837/ 182- 895/ 234)   0.0 E32-E25 S60-S54
  

    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