Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  482 1420 (N37E37 L=129)   329     75        1 ( 447/1403- 521/1435)   8.1 E40-E33 N35-N38
      2  867 1365 (N29E06 L=160)   257     61        1 ( 840/1351- 902/1387) -26.6 E07-E03 N27-N30
      3 1654 1208 (N07W45 L=211)   392     87        3 (1635/1167-1674/1252)  45.0 W42-W47 N04-N10
      4  514 1142 (N19E33 L=134)   358    137        1 ( 460/1096- 545/1185) -45.0 E36-E30 N15-N21
      5  442 1081 (N15E38 L=128)   569    203        1 ( 406/1000- 498/1165) -63.4 E41-E34 N10-N20
      6  299  907 (N06E52 L=114)   968    147        3 ( 258/ 838- 340/ 959)  37.9 E57-E47 N01-N08
      7 1594  859 (S12W35 L=201)   289     61        0 (1582/ 823-1612/ 890)  76.0 W33-W36 S14-S10
      8 1589  394 (S47W42 L=208)  1379    154        3 (1519/ 357-1649/ 436)   7.6 W33-W50 S50-S44
      9  502  350 (S33E56 L=110)   613    105        1 ( 483/ 304- 517/ 402)  32.7 E65-E51 S37-S28
     10  610  287 (S39E50 L=116)  1166    118        2 ( 558/ 255- 661/ 319)  14.7 E61-E42 S42-S36
  

    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