Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1530 1725 (N30W61 L= 71)   974    119        1 (1493/1678-1577/1759)   3.8 W55-W70 N27-N31
      2 1898 1219 (S07W67 L= 78)  1532    148        5 (1864/1155-1928/1282)  60.3 W63-W74 S10-S03
      3 1804 1178 (S07W55 L= 65)  1208    254        3 (1774/1083-1836/1271)  84.8 W50-W59 S12-S01
      4  418 1215 (N27E35 L=335)   296     88        0 ( 394/1184- 451/1254) -59.0 E36-E32 N24-N29
      5 1565 1190 (S00W36 L= 46)   227    132        0 (1526/1142-1613/1216)   7.1 W31-W39 S01-N00
      6  203 1145 (N27E57 L=314)   488     64        1 ( 169/1117- 240/1169)   6.3 E61-E52 N26-N28
      7 1303 1012 (S04W15 L= 26)  1959    587        3 (1098/ 835-1442/1149) -49.9 W04-W20 S16-N05
      8 1580 1081 (S07W34 L= 44)   587    144        2 (1537/1013-1615/1135) -76.0 W31-W36 S10-S02
      9  407 1112 (N21E37 L=333)   300     85        1 ( 375/1080- 437/1134) -45.0 E39-E34 N18-N22
     10 1689 1063 (S10W42 L= 52)   346     87        1 (1660/1026-1720/1093)  18.4 W38-W44 S11-S09
     11  744 1035 (N10E16 L=355)   602    173        2 ( 652/1011- 823/1062) -29.1 E20-E10 N07-N13
     12 1647  978 (S14W37 L= 47)   615    184        3 (1581/ 931-1706/1031)  10.3 W31-W42 S16-S12
     13 1731  800 (S27W42 L= 52)  1501    275        1 (1655/ 765-1841/ 876) -21.3 W34-W53 S30-S20
     14 1542  732 (S26W25 L= 36)   469    115        1 (1526/ 679-1559/ 803) -90.0 W23-W26 S29-S21
     15  888  707 (S12E15 L=355)   808    157        3 ( 818/ 663- 958/ 751)   0.0 E20-E10 S13-S08
     16  749  466 (S23E31 L=339)  3433    811        7 ( 506/ 244-1000/ 690)  15.7 E66-E09 S33-S14
     17 1530  537 (S39W23 L= 33)  2834    432        8 (1380/ 443-1629/ 660)  18.4 W09-W31 S45-S31
     18  470  579 (S10E46 L=325)   388    154        0 ( 403/ 560- 529/ 616)   0.0 E51-E40 S11-S08
  

    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