Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  966 1744 (N43E10 L= 62)   472     81        0 ( 934/1715-1008/1762)  26.6 E12-E06 N40-N44
      2 1107 1692 (N40W02 L= 74)   273     15        0 (1086/1675-1131/1710) -18.4 W00-W03 N38-N41
      3 1197 1664 (N38W08 L= 81)   293     72        1 (1172/1637-1234/1683)  45.0 W06-W11 N35-N39
      4  717 1619 (N32E26 L= 46)   663    141        0 ( 676/1558- 754/1672) -45.0 E30-E21 N27-N35
      5 1282 1628 (N36W15 L= 87)   337     68        0 (1266/1595-1305/1653)  53.1 W13-W16 N33-N38
      6 1344 1597 (N34W19 L= 92)   400     71        1 (1321/1567-1376/1625)  45.0 W17-W21 N31-N36
      7 1444 1551 (N32W27 L= 99)   606    104        2 (1390/1522-1496/1573)   0.0 W22-W30 N29-N33
      8 1606 1499 (N30W39 L=111)   472    103        0 (1578/1471-1646/1516)  26.6 W36-W42 N27-N31
      9 1216 1319 (N14W10 L= 82)   223     88        0 (1200/1267-1228/1363) -80.5 W08-W10 N11-N17
     10  579 1312 (N11E30 L= 43)   284     80        0 ( 540/1285- 610/1327) -21.8 E32-E27 N09-N11
  

    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