Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1200 1831 (N54W08 L=282)   580     50        1 (1176/1805-1221/1861)  60.3 W05-W09 N51-N58
      2 1350 1833 (N57W26 L=300)   417     65        0 (1321/1811-1381/1848)  24.4 W20-W31 N54-N59
      3 1603 1538 (N34W40 L=314)   391     84        0 (1588/1496-1624/1584)  81.9 W38-W42 N30-N37
      4  808 1476 (N21E18 L=256)   287     95        0 ( 788/1433- 823/1525)  90.0 E18-E16 N18-N24
      5 1323 1255 (N11W17 L=291)   434    112        0 (1300/1207-1352/1316) -66.8 W14-W18 N08-N15
      6 1299  800 (S16W19 L=293)   217     86        0 (1274/ 757-1326/ 827)  63.4 W17-W20 S18-S14
      7 1423  711 (S20W29 L=303)   479    158        0 (1381/ 664-1474/ 762)  60.3 W26-W31 S23-S16
      8  979  461 (S41W02 L=276)   537    105        0 ( 927/ 443-1031/ 482)  15.9 E02-W05 S42-S38
      9 1294  433 (S40W27 L=301)   605     62        1 (1271/ 399-1315/ 479) -80.5 W25-W30 S42-S36
  

    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