Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1579 1517 (N34W45 L= 89)   589     69        0 (1549/1488-1613/1552) -68.2 W42-W47 N31-N36
      2 1754 1513 (N32W62 L=106)   796    203        1 (1717/1496-1798/1535)  -6.3 W57-W66 N30-N33
      3  404 1490 (N34E49 L=355)   960    217        0 ( 348/1441- 467/1530) -23.2 E55-E41 N30-N36
      4  746 1387 (N27E19 L= 25)   236     73        1 ( 711/1371- 792/1399)   0.0 E21-E15 N26-N27
      5  274 1219 (N16E55 L=349)  1144    213        2 ( 246/1122- 303/1286) -55.0 E58-E51 N09-N19
      6 1655 1220 (N14W44 L= 88)   585    153        1 (1612/1142-1692/1286) -71.6 W41-W46 N08-N18
      7 1350  975 (N00W20 L= 64)   232     83        0 (1308/ 954-1396/ 996)   9.5 W17-W23 S01-N01
      8  241  665 (S20E62 L=342)   724    152        0 ( 229/ 611- 258/ 704)  54.5 E64-E59 S23-S16
      9 1485  655 (S21W31 L= 75)   274    114        0 (1434/ 644-1532/ 666)   0.0 W27-W34 S21-S19
     10 1228  572 (S25W13 L= 57)   333     97        0 (1177/ 559-1282/ 584)  -8.1 W09-W16 S26-S24
     11  431  545 (S27E46 L=358)   735    156        0 ( 392/ 514- 478/ 568)  16.7 E50-E40 S29-S25
  

    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