Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1820 (N55W04 L= 72)   436     68        2 ( 933/1798- 996/1845)  29.7 W01-W08 N52-N57
      2  822 1675 (N43E09 L= 59)   253     33        0 ( 803/1652- 844/1691)  33.7 E11-E07 N41-N44
      3  635 1519 (N33E23 L= 45)   685    109        3 ( 598/1475- 685/1554)   0.0 E26-E19 N30-N35
      4 1392 1453 (N22W28 L= 96)   413    131        1 (1332/1418-1452/1476) -15.9 W23-W31 N19-N23
      5 1186 1354 (N17W13 L= 81)   213     78        1 (1162/1316-1223/1393) -68.2 W11-W14 N14-N19
      6  400 1290 (N20E40 L= 28)   749    248        1 ( 335/1225- 481/1349) -36.9 E45-E33 N15-N24
      7 1813 1287 (N09W59 L=127)   878    148        4 (1771/1226-1840/1331) -74.1 W55-W62 N04-N11
      8  228 1021 (N05E56 L= 12)  1174    288        0 ( 173/ 905- 261/1145)  49.8 E63-E51 S01-N12
      9 1787  900 (S14W52 L=120)   333     62        2 (1768/ 871-1804/ 930) -63.4 W50-W53 S16-S12
     10  966  821 (S13E05 L= 63)   517    127        1 ( 930/ 766-1013/ 886)  50.2 E07-E01 S16-S09
     11  556  331 (S42E47 L= 21)   600    117        1 ( 513/ 307- 603/ 352)   0.0 E51-E42 S43-S40
  

    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