Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  729 1447 (N32E12 L= 77)   352    105        0 ( 676/1408- 774/1470)  15.9 E15-E08 N30-N33
      2  531 1401 (N33E27 L= 62)   414     90        0 ( 485/1365- 573/1431)  23.2 E31-E23 N31-N34
      3 1692 1052 (S09W44 L=132)   373     69        0 (1671/1019-1705/1091)  63.4 W41-W45 S10-S06
      4 1887 1029 (S14W64 L=153)  1418    196        4 (1848/ 968-1930/1082) -66.8 W58-W71 S17-S10
      5  218  986 (N12E57 L= 31)   405     56        0 ( 201/ 960- 235/1023) -76.0 E59-E55 N10-N14
      6 1274  890 (S11W12 L=101)   887    295        1 (1205/ 767-1336/1009)  45.0 W06-W18 S16-S04
      7  449  909 (N04E38 L= 51)   905    183        0 ( 366/ 871- 514/ 970)  14.0 E44-E32 N02-N06
      8  681  818 (S05E24 L= 65)   279     52        1 ( 651/ 799- 714/ 836)   0.0 E26-E21 S05-S03
      9  678  671 (S14E28 L= 61)   575    108        4 ( 651/ 637- 721/ 717)  45.0 E30-E25 S15-S10
  

    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