Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1282 1681 (N35W41 L=148)  1339    226        1 (1186/1643-1382/1705) -33.7 W33-W47 N29-N38
      2  367 1396 (N45E40 L= 67)  1355    176        2 ( 283/1351- 446/1434) -18.4 E52-E31 N40-N49
      3 1822 1111 (S13W55 L=162)  2173    460        4 (1716/ 942-1880/1265)  25.6 W40-W66 S19-S05
      4 1521 1129 (S02W31 L=138)   323     78        0 (1504/1088-1540/1175)  78.7 W30-W32 S05-N00
      5  148 1055 (N28E67 L= 40)  1648    142        2 (  88/1001- 196/1103) -19.3 E83-E60 N23-N31
      6 1576  651 (S32W25 L=132)  1371    275        1 (1531/ 545-1609/ 778)  71.6 W19-W27 S38-S23
      7  398  744 (N05E46 L= 60)   386     79        0 ( 362/ 716- 434/ 774) -51.3 E49-E43 N02-N07
      8  229  726 (N07E64 L= 43)   429     21        1 ( 208/ 707- 248/ 746)   0.0 E66-E61 N06-N08
      9  889  581 (S16E20 L= 87)   271     58        0 ( 874/ 548- 902/ 616)  63.4 E21-E18 S17-S13
     10  747  445 (S20E34 L= 73)   531    112        0 ( 697/ 404- 791/ 491)  16.7 E39-E29 S22-S18
     11 1366  409 (S42W03 L=110)   316     74        0 (1329/ 392-1408/ 431)  -8.1 W00-W07 S42-S40
  

    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