Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  449 1644 (N22E54 L=270)   392     75        1 ( 424/1626- 470/1664)  71.6 E56-E51 N20-N23
      2  186 1160 (S11E62 L=262)  1108    237        2 ( 125/1088- 238/1247)  45.0 E68-E57 S14-S05
      3 1765 1050 (N10W47 L= 12)   661    183        0 (1721/ 960-1781/1128) -53.1 W42-W49 N05-N14
      4 1875  944 (N08W61 L= 25)  2073    538        7 (1816/ 825-1940/1128) -55.5 W51-W73 N00-N17
      5  498 1078 (S12E33 L=291)   220     70        1 ( 473/1059- 524/1111) -33.7 E35-E31 S13-S10
      6  322  983 (S20E46 L=278)  2114    689        0 ( 147/ 834- 545/1079)  25.9 E62-E29 S30-S13
      7  824  746 (S27E07 L=317)   200     47        0 ( 799/ 720- 848/ 777)  63.4 E08-E05 S28-S24
      8  356  643 (S41E44 L=281)   970    256        3 ( 305/ 572- 411/ 719)  40.2 E50-E37 S46-S35
      9 1489  652 (S18W37 L=  2)   424    133        0 (1420/ 633-1544/ 691)  35.5 W32-W40 S20-S14
     10 1778  623 (S12W63 L= 27)  1216    318        2 (1716/ 570-1833/ 668)   0.0 W54-W71 S14-S10
     11 1554  404 (S31W54 L= 18)  2958    454        7 (1424/ 365-1730/ 448)  27.6 W43-W68 S35-S22
  

    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