Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  545 1749 (N29E58 L=143)  1569    229        3 ( 501/1709- 593/1777) -15.3 E63-E52 N25-N30
      2  356 1568 (N13E64 L=137)  1853    242        3 ( 321/1489- 388/1636) -55.0 E68-E59 N07-N17
      3 1847  868 (N10W61 L=262)   654    107        0 (1826/ 838-1869/ 908) -20.6 W57-W65 N07-N11
      4  768  847 (S20E11 L=190)   711    252        2 ( 680/ 792- 870/ 896)  36.0 E16-E04 S24-S15
      5 1102  816 (S13W10 L=210)  1407    375        3 ( 944/ 723-1246/ 890)  43.0 W02-W17 S22-S07
      6 1459  807 (S05W31 L=232)   425     63        1 (1431/ 770-1477/ 849) -63.4 W28-W32 S06-S02
      7  751  750 (S26E10 L=191)   238     48        0 ( 719/ 732- 774/ 769)  18.4 E11-E08 S27-S24
      8  277  660 (S45E49 L=151)   442    102        0 ( 249/ 643- 317/ 686)  -8.1 E52-E44 S45-S43
      9  585  627 (S39E19 L=181)  1256    310        0 ( 456/ 599- 727/ 651)  17.7 E30-E08 S42-S33
     10  288  565 (S52E51 L=150)   642    211        0 ( 247/ 545- 331/ 586)  17.1 E57-E44 S54-S49
     11 1367  462 (S27W39 L=240)  2291    672        0 (1047/ 392-1594/ 507)  14.3 W16-W59 S34-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