Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  851 1810 (N63W24 L=195)  1303    238        2 ( 762/1751- 952/1852)  -2.0 W09-W38 N59-N65
      2  599 1657 (N60E14 L=157)   643    106        0 ( 577/1614- 622/1707)  54.5 E17-E10 N55-N63
      3  673 1083 (N19E19 L=152)   383    152        0 ( 630/1036- 743/1112)   7.1 E22-E14 N16-N19
      4 1611  883 (S18W32 L=203)   295     50        1 (1590/ 858-1636/ 905)  18.4 W29-W33 S19-S16
      5 1455  709 (S23W17 L=189)   694    322        2 (1389/ 577-1485/ 844)  51.8 W10-W21 S29-S15
      6  817  584 (S13E24 L=148)   403     79        1 ( 782/ 569- 868/ 605) -36.9 E26-E21 S14-S11
      7  699  531 (S13E33 L=138)   693    138        2 ( 631/ 493- 767/ 568)   0.0 E38-E28 S15-S11
      8 1270  422 (S37E02 L=169)   376     80        0 (1230/ 404-1313/ 441) -11.3 E05-W00 S37-S34
  

    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