Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  161 1135 (S12E65 L=260)   441    109        0 ( 130/1121- 192/1151)   8.1 E69-E61 S12-S11
      2 1774 1088 (N13W48 L= 13)   483     80        1 (1751/1041-1801/1114)  68.2 W45-W50 N09-N14
      3  311  987 (S20E47 L=278)  2414    687        0 ( 144/ 837- 528/1078)  24.4 E63-E30 S30-S14
      4 1887  895 (N05W64 L= 30)  1807    315        2 (1856/ 827-1936/ 979) -66.0 W58-W72 N00-N09
      5  805  745 (S27E08 L=317)   253     65        0 ( 774/ 714- 832/ 779)  68.2 E10-E06 S29-S24
      6 1477  657 (S18W36 L=  2)   422    139        0 (1408/ 639-1529/ 696)  35.5 W31-W38 S20-S15
      7  363  674 (S39E43 L=283)   306     42        1 ( 340/ 652- 388/ 692)  14.0 E44-E40 S40-S37
      8 1778  627 (S12W63 L= 28)  1479    360        3 (1701/ 566-1845/ 674)  -3.2 W53-W71 S14-S08
      9  324  600 (S44E48 L=277)   458     76        1 ( 298/ 572- 351/ 627)  45.0 E50-E44 S46-S42
     10  602  594 (S40E22 L=304)   291     66        0 ( 573/ 566- 630/ 622)  38.7 E24-E19 S42-S38
     11 1547  408 (S31W53 L= 18)  2715    405        3 (1412/ 372-1725/ 453)  29.5 W42-W67 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