Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  520 1499 (N10E43 L=280)  2754    527        2 ( 466/1263- 602/1705) -64.3 E49-E34 S03-N24
      2 1771 1468 (N42W45 L=  8)   494     98        1 (1750/1438-1797/1488)  45.0 W42-W48 N39-N44
      3 1657 1389 (N31W32 L=356)   385     89        0 (1615/1358-1681/1407)  -9.5 W28-W34 N29-N32
      4 1314 1305 (N16W10 L=333)   434    132        0 (1278/1237-1343/1377) -49.4 W06-W12 N12-N19
      5  206 1186 (S13E61 L=263)  2286    526        1 ( 130/1146- 319/1271)  -7.1 E72-E48 S15-S08
      6 1726  969 (N08W45 L=  8)   411     76        0 (1705/ 932-1744/1014) -36.9 W42-W47 N05-N10
      7  722  927 (S19E16 L=308)   223     73        0 ( 687/ 916- 767/ 940)  21.8 E17-E12 S20-S18
      8 1290  841 (S11W20 L=343)   215     67        0 (1263/ 810-1319/ 878) -21.8 W17-W22 S12-S09
      9 1434  701 (S15W33 L=356)   515    126        0 (1378/ 659-1503/ 761)  60.9 W30-W35 S18-S09
     10 1709  593 (S12W58 L= 21)   397     76        0 (1694/ 559-1723/ 624) -63.4 W55-W60 S14-S10
     11 1665  557 (S15W56 L= 19)   559    112        3 (1639/ 516-1694/ 609) -63.4 W53-W59 S18-S12
     12 1571  433 (S25W55 L= 18)  1499    206        3 (1513/ 401-1633/ 480)  32.0 W49-W60 S27-S21
  

    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