Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1288 1684 (N35W36 L=177)  1665    325        5 (1189/1616-1406/1720)  -3.2 W27-W45 N32-N37
      2  440 1581 (N51E35 L=106)   328     63        1 ( 425/1556- 460/1610) -76.0 E36-E32 N49-N53
      3  576 1547 (N45E20 L=120)   360    103        1 ( 524/1529- 629/1565)  -5.7 E25-E15 N43-N45
      4 1493 1037 (S07W28 L=168)  1162    373        3 (1422/ 872-1562/1189)  43.0 W20-W35 S14-N00
      5 1733  841 (S24W42 L=183)   279    116        1 (1688/ 802-1771/ 867)  15.9 W38-W45 S25-S23
      6 1509  781 (S22W24 L=165)   335     65        2 (1474/ 753-1538/ 802)  11.3 W21-W26 S23-S20
      7 1088  481 (S31E09 L=131)   850    216        2 (1024/ 407-1113/ 571)  48.4 E15-E05 S35-S25
      8  664  434 (S24E38 L=103)   952    347        2 ( 639/ 312- 689/ 566)  50.9 E44-E30 S32-S16
      9 1514  535 (S39W22 L=163)   412     68        1 (1491/ 499-1543/ 565)  45.0 W20-W24 S40-S36
     10 1388  455 (S41W11 L=152)   254     48        1 (1363/ 437-1420/ 474)   0.0 W09-W14 S41-S39
  

    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