Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  935 1426 (N33W04 L=258)  1414    395        5 ( 813/1285-1024/1531) -72.5 E02-W07 N22-N41
      2 1667 1395 (N12W52 L=306)   410     99        1 (1625/1364-1708/1421) -45.0 W49-W54 N09-N14
      3 1348 1256 (N13W25 L=279)   276    121        1 (1287/1230-1389/1281) -33.7 W20-W27 N10-N14
      4  372 1138 (N26E42 L=212)  2670    568        7 ( 215/1022- 503/1249) -35.3 E56-E32 N17-N34
      5 1719  884 (S18W42 L=296)  2043    357        3 (1633/ 746-1767/1066)  48.8 W33-W49 S24-S08
      6  702  986 (N11E20 L=234)   340     65        2 ( 671/ 958- 736/1018) -45.0 E22-E17 N09-N13
      7  828  939 (N06E13 L=241)   347     72        2 ( 795/ 910- 861/ 972)   0.0 E15-E10 N04-N08
      8  226  817 (N08E62 L=192)   455    154        2 ( 196/ 794- 268/ 851) -14.0 E65-E56 N06-N10
      9  813  741 (S05E18 L=236)   480    144        1 ( 733/ 707- 884/ 772)   0.0 E23-E13 S05-S04
     10 1204  593 (S23W01 L=255)   273     90        0 (1178/ 544-1219/ 627)  78.7 W00-W02 S25-S20
     11  563  389 (S23E48 L=206)   433     86        1 ( 523/ 349- 606/ 423) -80.5 E50-E46 S26-S20
  

    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