Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  948  866 (S05E05 L=224)   312    147        1 ( 888/ 844-1013/ 902) -29.7 E08-E01 S06-S02
      2  473  727 (S12E39 L=191)   306     41        2 ( 450/ 702- 496/ 755) -90.0 E40-E37 S13-S10
      3 1393  663 (S20W23 L=252)   237     83        2 (1372/ 625-1406/ 697) -78.7 W21-W23 S23-S18
      4  451  637 (S18E43 L=187)   457     72        1 ( 423/ 609- 492/ 663)  18.4 E45-E39 S19-S16
      5 1691  624 (S26W49 L=279)   399     61        1 (1675/ 600-1709/ 657) -76.0 W47-W51 S27-S23
      6  395  523 (S26E53 L=177)   698     84        1 ( 363/ 495- 437/ 550) -26.6 E56-E49 S28-S24
      7 1461  464 (S36W31 L=260)   340     61        1 (1442/ 439-1482/ 494)  63.4 W29-W32 S37-S33
      8  404  429 (S33E59 L=171)   782     74        2 ( 391/ 386- 433/ 473)  36.0 E66-E54 S37-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