Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  441 1666 (N25E59 L=148)   484     58        0 ( 420/1638- 459/1691)  90.0 E61-E56 N23-N27
      2 1588 1570 (N40W32 L=239)   330     45        1 (1565/1548-1613/1589)   0.0 W30-W34 N38-N41
      3 1418 1429 (N25W18 L=225)   294     82        1 (1377/1405-1459/1455)  -8.1 W14-W21 N24-N26
      4 1617 1343 (N24W33 L=240)   292    116        0 (1571/1297-1658/1381) -23.2 W29-W37 N21-N25
      5  228 1124 (S11E58 L=149)   402     90        2 ( 207/1082- 252/1166)  90.0 E60-E55 S13-S08
      6  136 1102 (S13E70 L=137)  2072    135        6 ( 100/1039- 176/1163)  41.2 E78-E63 S16-S08
      7  335  808 (S29E46 L=160)  1389    267        3 ( 227/ 773- 473/ 837)   6.6 E59-E33 S32-S27
      8 1518  538 (S25W45 L=251)  3509    728        6 (1255/ 454-1862/ 697)  31.5 W27-W71 S35-S07
      9 1714  642 (S14W56 L=263)   355     56        1 (1694/ 612-1738/ 663) -36.9 W53-W58 S16-S13
  

    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