Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  428 1652 (N16E64 L= 74)  3617    273        5 ( 349/1585- 497/1777) -61.9 E74-E54 N11-N26
      2  769 1637 (N23E35 L=103)   510    138        1 ( 748/1574- 797/1693) -41.2 E38-E30 N19-N27
      3  685 1571 (N17E37 L=101)   955    273        4 ( 566/1523- 806/1623)   9.0 E47-E28 N15-N19
      4  610 1504 (N11E39 L= 99)   415    133        1 ( 577/1482- 674/1543)  15.9 E41-E34 N09-N13
      5  767 1480 (N14E28 L=110)   303     23        0 ( 746/1460- 790/1501) -14.0 E30-E26 N12-N14
      6  804  936 (S16E10 L=128)   309     83        0 ( 774/ 897- 833/ 975)  78.7 E11-E09 S18-S13
      7 1472  558 (S17W42 L=181)   301     50        1 (1439/ 541-1495/ 581)  45.0 W40-W44 S18-S15
      8 1022  506 (S35W17 L=155)   241     56        0 ( 989/ 491-1054/ 518)  33.7 W14-W18 S35-S33
      9  617  458 (S52E11 L=127)   952    152        0 ( 552/ 436- 702/ 492)  10.6 E18-E02 S54-S48
  

    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