Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1477 1634 (N29W47 L=358)  2469    416        3 (1273/1564-1630/1696) -32.5 W33-W60 N21-N36
      2 1729 1510 (N16W62 L= 14)   619     59        2 (1691/1493-1762/1531)  -7.1 W58-W67 N14-N17
      3 1428 1460 (N19W35 L=346)   315     65        0 (1391/1442-1466/1475) -21.8 W32-W37 N17-N20
      4  536 1276 (N25E27 L=285)  2296    560        1 ( 315/1214- 805/1345)  -7.5 E45-E07 N22-N29
      5 1839 1263 (S00W62 L= 13)   728     88        0 (1798/1231-1868/1304)  10.3 W56-W67 S01-N01
      6 1638 1102 (S05W40 L=351)   259     26        1 (1611/1082-1663/1118) -14.0 W37-W41 S06-S04
      7 1871  866 (S24W60 L= 11)  1373    178        2 (1841/ 769-1901/ 942)  74.7 W56-W64 S30-S19
      8 1103  781 (S14W00 L=311)   257     81        0 (1086/ 735-1126/ 818)  63.4 E01-W02 S16-S11
      9 1272  748 (S19W10 L=321)   298     85        0 (1241/ 710-1292/ 795)  45.0 W07-W11 S21-S16
     10  454  420 (S24E56 L=255)   654    104        0 ( 436/ 370- 475/ 481)  28.6 E62-E51 S27-S19
  

    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