Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  473 1659 (N37E59 L=342)  1206    128        1 ( 414/1629- 533/1684)  -3.0 E68-E49 N35-N38
      2  582 1547 (N31E40 L=  1)   722    100        1 ( 556/1491- 628/1595) -40.6 E43-E35 N27-N34
      3  915 1418 (N25E11 L= 30)  1233    332        0 ( 823/1343- 964/1540)  85.6 E16-E07 N20-N33
      4 1763 1353 (N28W57 L= 98)  4611    620        7 (1547/1166-1924/1517) -27.7 W37-W79 N16-N38
      5  347 1375 (N17E54 L=347)   344     68        2 ( 318/1351- 374/1392)   0.0 E56-E50 N15-N18
      6 1237 1310 (N21W11 L= 52)   481    222        3 (1154/1268-1315/1358)  31.0 W06-W16 N17-N24
      7 1781 1192 (N18W55 L= 96)  1349    297        2 (1727/1059-1816/1305) -68.2 W50-W58 N09-N24
      8 1690 1148 (N14W45 L= 86)   467    113        0 (1660/1115-1733/1198) -29.7 W42-W49 N12-N17
      9 1392 1129 (N11W22 L= 63)   229     65        0 (1361/1107-1421/1151)  36.9 W20-W24 N09-N12
     10  800  727 (S19E12 L= 29)   693    229        0 ( 769/ 638- 864/ 787)  78.7 E13-E07 S24-S14
     11 1716  622 (S18W57 L= 98)   583     68        0 (1685/ 598-1741/ 653)  36.9 W54-W59 S19-S16
  

    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