Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1638 1579 (N29W62 L= 61)   922    196        2 (1595/1514-1694/1611)   6.3 W53-W71 N25-N30
      2 1520 1513 (N28W45 L= 45)   314     77        0 (1482/1494-1559/1528) -18.4 W42-W48 N25-N28
      3 1737 1444 (N19W63 L= 63)  2267    206        2 (1673/1368-1794/1528) -76.0 W58-W69 N13-N25
      4 1341 1456 (N28W29 L= 29)   381    171        0 (1320/1396-1393/1527)  83.7 W27-W31 N23-N32
      5  587 1336 (N31E27 L=333)  1335    286        6 ( 547/1212- 650/1455) -67.6 E30-E21 N22-N39
      6 1265 1411 (N26W22 L= 22)   258     41        1 (1238/1390-1287/1431) -33.7 W20-W23 N24-N27
      7  855 1181 (N18E09 L=351)   208     43        1 ( 829/1164- 881/1203) -14.0 E10-E06 N16-N19
      8 1535 1103 (N03W34 L= 33)   444     80        1 (1499/1075-1562/1138)  11.3 W30-W35 N01-N04
      9 1712 1031 (S05W46 L= 46)   354     99        0 (1700/ 978-1725/1077) -90.0 W45-W47 S07-S01
     10  269  715 (S06E60 L=300)   578     82        1 ( 248/ 674- 289/ 758)  40.6 E64-E56 S08-S02
     11  947  619 (S18E11 L=349)   227     89        0 ( 906/ 593- 983/ 642)  18.4 E13-E07 S19-S16
  

    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