Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  923 1843 (N42E41 L=135)   902    125        1 ( 862/1806- 984/1876)   0.0 E49-E34 N39-N43
      2  452 1666 (N17E63 L=114)  1421    179        0 ( 368/1620- 499/1697)   0.0 E73-E54 N14-N19
      3  791 1259 (N00E19 L=157)   767    229        1 ( 732/1165- 859/1348) -22.6 E25-E13 S03-N04
      4  157 1087 (S23E68 L=109)  4663    284        6 ( 104/ 979- 240/1188)  45.0 E79-E57 S30-S16
      5 1689 1122 (N19W38 L=215)  1047    169        2 (1625/1078-1743/1164)   5.2 W32-W43 N17-N21
      6  291 1128 (S19E52 L=125)   325     73        0 ( 272/1091- 312/1160)  68.2 E53-E49 S21-S16
      7  856 1101 (S07E11 L=165)   290     76        0 ( 827/1065- 890/1133)  68.2 E12-E09 S09-S04
      8  259  979 (S28E53 L=124)   359     53        0 ( 238/ 946- 281/1008)  53.1 E55-E50 S30-S26
      9  960  674 (S27W07 L=184)  1252    455        1 ( 783/ 590-1166/ 746)  49.6 E01-W16 S37-S17
     10 1325  658 (S17W29 L=205)  1113    297        0 (1186/ 628-1471/ 682)  34.5 W21-W37 S22-S11
     11  633  366 (S60E06 L=171)   343     53        0 ( 607/ 344- 655/ 385)  -8.1 E09-E02 S61-S58
  

    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