Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1006 1866 (N51E26 L=293)   723     74        3 ( 968/1846-1051/1892) -15.3 E31-E20 N48-N53
      2  734 1704 (N31E35 L=283)  1123    313        1 ( 682/1555- 795/1816) -69.1 E41-E31 N19-N40
      3 1492 1518 (N33W22 L=341)   366     83        0 (1445/1502-1531/1534)  23.2 W18-W25 N31-N35
      4  591  936 (S19E25 L=293)   722    172        2 ( 537/ 862- 650/ 997)  53.1 E28-E21 S23-S14
      5  408  585 (S44E42 L=276)  1160    159        2 ( 361/ 544- 459/ 638) -35.5 E47-E36 S46-S40
      6 1495  576 (S23W41 L=360)   428    131        1 (1446/ 534-1555/ 615)  54.5 W38-W44 S26-S19
      7 1745  486 (S21W72 L= 30)  2943    228        3 (1662/ 408-1816/ 554)  56.3 W63-W78 S27-S15
      8 1288  492 (S33W30 L=349)   329     81        0 (1249/ 474-1332/ 514)  45.0 W27-W32 S34-S30
      9 1098  218 (S57W36 L=355)   858    121        2 (1041/ 185-1146/ 251)  -6.7 W27-W44 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