Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1349 1691 (N35W48 L=149)  3370    525        6 (1210/1536-1551/1768)  -7.1 W31-W71 N28-N39
      2  608 1142 (N24E22 L= 79)   395     87        1 ( 580/1104- 638/1192) -80.5 E23-E20 N21-N27
      3  200 1011 (N24E61 L= 40)   518    126        0 ( 151/ 968- 245/1036) -22.6 E68-E55 N20-N25
      4 1542  959 (S11W29 L=130)   232     59        1 (1518/ 937-1571/ 989) -53.1 W27-W30 S12-S08
      5 1735  892 (S21W42 L=144)   817    192        0 (1653/ 865-1800/ 936) -28.3 W35-W49 S23-S15
      6 1813  769 (S32W50 L=152)  3222    406        6 (1711/ 630-1936/ 912)  15.3 W38-W71 S38-S26
      7  243  874 (N15E58 L= 44)   440     92        1 ( 223/ 829- 261/ 908)  59.0 E60-E55 N12-N17
      8  175  795 (N11E70 L= 32)  2537    148        4 ( 127/ 743- 226/ 844) -32.5 E79-E63 N07-N14
      9 1433  758 (S19W17 L=119)   359    102        0 (1418/ 700-1454/ 805)  63.4 W14-W19 S22-S16
     10  384  571 (S05E55 L= 47)   507     89        0 ( 339/ 558- 423/ 591)   0.0 E59-E50 S06-S03
     11  771  359 (S26E37 L= 65)   411     53        0 ( 738/ 330- 797/ 388)  23.2 E40-E33 S28-S24
  

    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