Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  689 1680 (N37E29 L= 82)  3824    769       13 ( 485/1587- 916/1777) -16.8 E53-E10 N30-N45
      2  571 1344 (N14E30 L= 81)   402    104        0 ( 554/1289- 589/1391)  90.0 E31-E28 N10-N16
      3 1379 1233 (N10W21 L=132)   628    157        4 (1319/1156-1419/1297) -58.0 W17-W23 N05-N13
      4 1857 1036 (N02W60 L=171)   492     74        1 (1840/ 991-1879/1072)  76.0 W57-W62 N00-N04
      5  934 1029 (S04E05 L=106)   299     71        0 ( 900/ 996- 959/1063)  53.1 E07-E03 S06-S02
      6  757  695 (S25E16 L= 95)  1059    267        1 ( 652/ 603- 882/ 800)  39.1 E24-E07 S31-S18
      7 1436  678 (S23W29 L=141)   247     41        0 (1409/ 661-1464/ 695)  18.4 W27-W31 S24-S21
      8 1086  551 (S33W07 L=118)   436     97        1 (1025/ 536-1137/ 569)   0.0 W02-W10 S34-S32
      9  699  318 (S53E28 L= 84)   362     71        0 ( 682/ 289- 713/ 349)  51.3 E30-E25 S55-S50
  

    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