Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  848 1384 (N14E15 L=112)   234     40        2 ( 833/1355- 869/1417)  90.0 E16-E13 N11-N15
      2 1466  792 (S13W31 L=158)   353    100        0 (1440/ 740-1499/ 830)  71.6 W29-W32 S16-S10
      3 1718  643 (S17W56 L=183)   378     62        2 (1694/ 619-1746/ 673) -31.0 W52-W58 S18-S15
      4 1429  530 (S29W36 L=164)   372    111        0 (1396/ 492-1477/ 578)  66.8 W34-W38 S32-S25
      5 1552  513 (S28W47 L=175)   509    147        0 (1494/ 480-1608/ 568)  49.4 W43-W49 S30-S23
      6 1654  522 (S26W56 L=184)   495     97        1 (1608/ 490-1692/ 553)  50.2 W53-W59 S28-S22
      7 1430  464 (S34W40 L=167)   373    105        0 (1380/ 450-1477/ 486)   6.3 W34-W43 S34-S31
      8 1642  422 (S32W63 L=191)   509     94        0 (1586/ 407-1683/ 441)  15.3 W57-W68 S33-S30
      9 1322  171 (S57W63 L=190)   633     58        1 (1300/ 146-1344/ 203)  -9.0 W55-W74 S58-S54
  

    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