Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  354 1598 (N12E64 L=266)   839    110        1 ( 327/1571- 379/1622) -26.6 E68-E60 N10-N13
      2  278 1111 (S19E51 L=279)  1555    596        1 ( 115/1027- 415/1193)  -8.4 E72-E38 S21-S15
      3  170 1021 (S25E60 L=270)   763    179        0 ( 131/ 954- 210/1079)  58.0 E64-E57 S29-S21
      4  481 1046 (S18E34 L=297)   214     66        1 ( 455/1025- 510/1076)  63.4 E35-E32 S20-S16
      5 1806  973 (N13W51 L= 22)   379    114        1 (1793/ 918-1815/1039) -60.3 W48-W53 N09-N16
      6  546  987 (S21E28 L=302)   644    181        1 ( 460/ 953- 620/1015)  24.4 E33-E22 S24-S18
      7  840  927 (S17E08 L=322)   304    117        0 ( 778/ 879- 879/ 978) -20.6 E12-E04 S19-S15
      8 1332  792 (S12W24 L=354)   217     60        2 (1302/ 772-1363/ 811)  45.0 W22-W25 S13-S09
      9  438  773 (S36E33 L=297)   413    147        1 ( 376/ 743- 507/ 800)   0.0 E38-E27 S36-S34
     10  456  588 (S48E31 L=299)   634    183        1 ( 407/ 537- 507/ 631)  39.3 E36-E25 S53-S44
     11 1379  418 (S30W43 L= 13)  2778    587        3 (1158/ 364-1617/ 467)  18.9 W26-W61 S35-S22
  

    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