Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  530 1558 (N31E45 L=180)   958    130        2 ( 472/1505- 586/1599)  -8.7 E51-E38 N27-N34
      2  286 1527 (N26E73 L=152)  1060     69        1 ( 263/1499- 314/1558)  -5.2 E88-E65 N24-N27
      3 1772 1459 (N35W65 L=290)   617     78        1 (1748/1436-1793/1483) -21.8 W59-W70 N33-N36
      4  193  843 (S17E64 L=161)  2320    156        2 ( 156/ 777- 246/ 913) -82.9 E70-E56 S21-S13
      5 1309  879 (S05W19 L=244)   198     57        0 (1279/ 871-1350/ 893)  11.3 W16-W21 S05-S04
      6  409  820 (S16E41 L=184)   331     41        0 ( 386/ 803- 443/ 844)   0.0 E42-E37 S17-S14
      7  774  792 (S15E14 L=211)   564    195        1 ( 681/ 761- 854/ 812)   0.0 E20-E08 S15-S13
      8 1314  673 (S18W22 L=247)   376    118        0 (1293/ 621-1343/ 719)  90.0 W20-W23 S21-S14
      9 1433  663 (S18W30 L=255)   268     58        0 (1404/ 643-1461/ 680) -14.0 W28-W32 S18-S16
     10 1528  494 (S28W44 L=269)   470     65        2 (1502/ 472-1558/ 522)  36.9 W40-W46 S30-S26
     11  499  422 (S44E42 L=183)   654     79        1 ( 474/ 386- 530/ 468)  41.2 E47-E37 S47-S40
     12 1414  431 (S34W37 L=262)   429     42        1 (1388/ 409-1439/ 453)  33.7 W34-W39 S35-S32
  

    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