Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  568 1691 (N28E47 L=309)   365     69        0 ( 531/1665- 598/1714)  53.1 E49-E45 N26-N30
      2  784 1487 (N17E23 L=333)   238     57        0 ( 768/1457- 800/1515) -45.0 E24-E21 N15-N19
      3  190 1035 (S16E61 L=295)   830    116        1 ( 167/ 984- 230/1071)  35.5 E64-E56 S19-S14
      4  287  903 (S24E50 L=306)   284     93        0 ( 267/ 861- 314/ 943)  50.2 E52-E47 S26-S20
      5 1714  644 (S14W56 L= 52)   413     41        1 (1685/ 621-1734/ 670)  26.6 W53-W58 S15-S12
      6  549  624 (S38E28 L=329)   345     52        0 ( 521/ 598- 582/ 642)  26.6 E30-E24 S40-S36
      7  842  305 (S57W02 L=358)   576     72        1 ( 798/ 283- 879/ 340)   7.1 E02-W06 S59-S53
  

    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