Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  465 1718 (N31E60 L=333)   723     92        1 ( 425/1689- 496/1743)  53.1 E64-E56 N28-N32
      2  512 1619 (N23E47 L=347)   325     73        0 ( 489/1588- 544/1647)  76.0 E48-E45 N21-N25
      3 1856 1133 (N16W58 L= 91)   822     94        1 (1836/1086-1877/1185) -59.0 W55-W60 N13-N19
      4  855  816 (S21E07 L= 27)   413    159        0 ( 789/ 793- 928/ 855)  32.0 E11-E03 S22-S17
      5  295  627 (S40E56 L=338)   775     94        0 ( 279/ 577- 318/ 678)  49.4 E60-E50 S43-S36
      6 1403  499 (S31W37 L= 71)   439    124        0 (1371/ 442-1446/ 548) -82.9 W35-W39 S34-S26
      7  571  357 (S58E32 L=  2)   646     71        2 ( 541/ 332- 602/ 390) -23.2 E36-E28 S60-S54
  

    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