Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1559 1722 (N58W54 L=163)   709     84        1 (1526/1694-1586/1744)  13.5 W44-W69 N54-N60
      2 1624 1594 (N46W50 L=158)   401     56        0 (1588/1578-1649/1608)  -8.1 W45-W53 N44-N47
      3  996 1384 (N23E06 L=102)   854    407        0 ( 904/1267-1110/1497) -36.9 E14-W02 N16-N29
      4 1496 1427 (N32W30 L=138)   443    102        0 (1456/1380-1523/1470) -56.3 W26-W32 N28-N34
      5  260 1350 (N12E63 L= 46)   432     95        2 ( 242/1322- 290/1387) -38.7 E66-E60 N10-N14
      6  894 1310 (N17E12 L= 97)   597    167        4 ( 863/1218- 919/1379) -84.3 E13-E09 N11-N21
      7 1108 1316 (N20W02 L=110)   648    196        1 (1044/1258-1185/1353) -18.4 E02-W07 N17-N21
      8  810  762 (S17E11 L= 98)   486    122        1 ( 764/ 730- 877/ 790)  35.5 E13-E06 S19-S14
      9 1366  633 (S19W27 L=136)   249     48        1 (1344/ 606-1390/ 656)  90.0 W25-W28 S21-S17
     10 1424  388 (S36W42 L=150)   555    127        1 (1382/ 349-1472/ 424)  60.3 W39-W44 S39-S32
  

    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