Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  698 1683 (N37E28 L= 82)  3706    672        8 ( 497/1591- 922/1783) -18.0 E52-E09 N30-N46
      2  588 1341 (N13E29 L= 81)   350     86        1 ( 569/1299- 605/1386) -80.5 E30-E27 N10-N16
      3  951 1024 (S05E04 L=106)   261     75        0 ( 918/ 992- 977/1060)  63.4 E06-E02 S06-S02
      4  776  695 (S25E15 L= 96)   993    268        0 ( 663/ 601- 905/ 793)  37.4 E23-E06 S31-S18
      5  239  681 (S26E58 L= 53)   549    124        1 ( 221/ 634- 260/ 727)  56.3 E60-E54 S29-S23
      6 1094  549 (S33W07 L=118)   450    102        1 (1035/ 532-1148/ 568)   7.1 W03-W11 S34-S32
      7  708  319 (S53E26 L= 84)   303     64        0 ( 696/ 293- 722/ 347)  68.2 E28-E24 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