Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1368 1204 (N16W18 L=284)  3974    938        7 (1059/ 827-1619/1523) -32.9 E07-W41 S01-N30
      2  220 1034 (S15E57 L=209)   689     78        1 ( 201/ 984- 242/1072) -55.0 E59-E55 S17-S12
      3  139  991 (S19E67 L=199)   909    133        1 ( 117/ 954- 158/1020)  78.7 E70-E64 S21-S16
      4  559  974 (S12E28 L=238)   239     58        1 ( 543/ 943- 572/1012) -76.0 E28-E26 S14-S10
      5 1388  772 (S09W27 L=293)   479    116        0 (1336/ 752-1448/ 792)  26.6 W23-W30 S10-S07
      6  593  760 (S25E23 L=243)   461    122        1 ( 534/ 725- 649/ 784)   0.0 E27-E18 S26-S23
      7 1094  689 (S20W11 L=277)   768    187        0 (1011/ 669-1199/ 721)   9.5 W05-W17 S21-S17
      8 1477  304 (S37W55 L=321)  1047    187        1 (1448/ 272-1508/ 361) -50.2 W50-W59 S39-S33
      9 1297  300 (S42W41 L=307)  2085    416        2 (1123/ 272-1444/ 337)  10.9 W26-W53 S45-S37
  

    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