Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  398 1531 (N31E48 L=160)   385     42        2 ( 372/1511- 424/1551)  -9.5 E51-E45 N29-N31
      2 1611 1390 (N20W41 L=249)  1053    294        1 (1485/1309-1719/1469) -32.0 W32-W49 N14-N24
      3  188 1206 (N10E62 L=146)  5076    428        9 ( 132/1077- 253/1405) -64.5 E70-E54 N02-N23
      4  501 1292 (N14E34 L=175)   720    216        1 ( 409/1245- 598/1315)  -8.1 E40-E26 N10-N15
      5  689  831 (S14E21 L=187)   893    284        0 ( 643/ 712- 719/ 967)  71.6 E25-E19 S21-S06
      6 1278  865 (S13W16 L=224)   196     56        1 (1246/ 845-1303/ 893) -26.6 W13-W17 S13-S11
      7 1255  794 (S17W14 L=223)   321    117        2 (1204/ 779-1302/ 825)  -9.5 W11-W17 S17-S15
      8 1652  448 (S39W57 L=265)   878    109        2 (1594/ 431-1710/ 470)   3.4 W48-W65 S40-S37
      9  688  457 (S39E27 L=181)   315     76        0 ( 648/ 445- 737/ 469)   0.0 E30-E23 S39-S37
  

    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