Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1654 1521 (N21W59 L=195)   554    101        0 (1604/1501-1700/1541) -32.0 W55-W63 N18-N23
      2 1849 1310 (N04W69 L=205)   514     54        1 (1826/1288-1867/1340)  31.0 W64-W72 N02-N05
      3 1891 1121 (S08W67 L=203)   869    140        1 (1866/1070-1912/1184) -63.4 W64-W70 S11-S03
      4 1617  949 (S10W36 L=172)   326    106        1 (1596/ 903-1645/ 993)  38.7 W33-W38 S12-S07
      5 1803  860 (S21W52 L=188)   540     82        1 (1772/ 828-1829/ 904)  38.7 W48-W55 S22-S18
      6 1637  831 (S18W36 L=172)   612    176        0 (1590/ 797-1693/ 871)   0.0 W32-W40 S19-S14
      7  900  611 (S16E16 L=120)   693    243        0 ( 858/ 528- 990/ 659)  29.1 E19-E09 S20-S13
      8 1428  452 (S39W17 L=153)   456     93        0 (1389/ 421-1461/ 480)  26.6 W13-W19 S40-S37
      9 1074  206 (S51E20 L=116)   442     52        0 (1055/ 188-1100/ 235)  80.5 E23-E17 S54-S48
  

    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