Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1361 1225 (N17W17 L=283)  3517    782        6 (1067/ 932-1590/1519) -30.6 E07-W37 N03-N30
      2  153 1009 (S17E65 L=200)  1324    187        2 ( 121/ 951- 182/1072)  10.1 E69-E61 S21-S13
      3  226 1032 (S15E57 L=209)   633     86        0 ( 209/ 982- 247/1071)  90.0 E58-E54 S17-S12
      4  260  961 (S18E52 L=213)   885    246        2 ( 199/ 873- 299/1040)  54.0 E58-E49 S24-S13
      5 1389  764 (S10W27 L=293)   439    113        1 (1329/ 744-1443/ 787)  26.6 W24-W30 S11-S08
      6  602  757 (S25E22 L=243)   432    120        2 ( 545/ 723- 654/ 780)  -7.1 E26-E18 S26-S23
      7 1103  685 (S20W11 L=277)   732    171        1 (1025/ 663-1198/ 721)   5.2 W06-W17 S21-S18
      8 1349  298 (S40W45 L=310)  3256    552        6 (1129/ 267-1512/ 357)  10.3 W26-W59 S45-S33
  

    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