Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  685 1745 (N29E49 L=235)  1428    142        2 ( 638/1703- 744/1812) -16.4 E58-E41 N25-N33
      2  409 1623 (N15E64 L=220)  2681    211        6 ( 373/1525- 464/1701) -47.5 E73-E56 N08-N20
      3 1690 1566 (N52W40 L=324)  1593    201        2 (1607/1491-1762/1613)  -8.5 W28-W49 N48-N56
      4  561 1537 (N13E44 L=240)   336     89        0 ( 540/1497- 589/1576) -78.7 E45-E43 N10-N15
      5 1743 1214 (N27W43 L=327)   457     54        0 (1721/1176-1762/1244) -63.4 W40-W44 N25-N29
      6  289 1088 (S19E51 L=233)   347     51        1 ( 266/1060- 310/1111) -45.0 E53-E48 S20-S17
      7 1757  912 (N10W50 L=334)   483    116        0 (1743/ 858-1769/ 973) -56.3 W47-W52 N07-N13
      8  433  815 (S32E34 L=250)   574    163        1 ( 363/ 776- 500/ 858)  36.9 E40-E28 S36-S27
      9  832  599 (S34W01 L=285)   331     58        2 ( 800/ 576- 868/ 619)  63.4 W00-W02 S36-S32
     10 1592  493 (S17W56 L=340)  3371    417        3 (1503/ 326-1735/ 606) -76.6 W50-W67 S27-S06
     11 1333  419 (S30W41 L=325)   354     85        1 (1288/ 402-1377/ 433)  38.7 W38-W44 S31-S27
  

    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