Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1468 1754 (N47W44 L=317)   807    117        0 (1412/1738-1533/1768)   0.0 W37-W53 N45-N48
      2 1476  722 (S21W30 L=303)   318     76        0 (1441/ 695-1506/ 756)  36.9 W27-W31 S23-S19
      3  283  668 (S23E57 L=216)   839    130        0 ( 264/ 613- 319/ 729)  82.9 E59-E54 S26-S18
      4  976  581 (S30E04 L=269)   660    173        1 ( 906/ 524-1029/ 650)  48.8 E08-W00 S34-S25
      5 1444  518 (S35W31 L=305)  1460    346        5 (1379/ 357-1514/ 648)  84.3 W29-W34 S47-S26
      6 1109  381 (S45W06 L=280)   531     65        2 (1075/ 346-1139/ 412)  45.0 W03-W08 S48-S42
      7  936  321 (S50E09 L=264)   814    158        2 ( 861/ 292-1009/ 367) -24.8 E15-E02 S52-S45
      8  567  288 (S51E52 L=222)  1087    117        0 ( 510/ 266- 632/ 314)   5.2 E63-E41 S52-S49
  

    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