Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1273 1802 (N52W14 L=261)   664     77        2 (1250/1756-1293/1842)  83.7 W11-W16 N47-N56
      2  548 1661 (N33E44 L=203)   913    103        1 ( 514/1599- 579/1699) -41.6 E49-E39 N27-N35
      3 1479 1483 (N28W27 L=274)   565    107        1 (1431/1444-1530/1510)  14.0 W23-W31 N25-N30
      4 1727 1183 (N12W46 L=293)   504     64        0 (1694/1148-1753/1210) -18.4 W42-W48 N09-N13
      5 1683 1062 (N04W43 L=290)   323     46        1 (1664/1031-1699/1090)  71.6 W40-W43 N02-N05
      6 1711  645 (S20W54 L=301)   541     51        1 (1682/ 621-1735/ 670)  56.3 W51-W57 S21-S18
      7 1351  649 (S25W25 L=273)   217     28        0 (1332/ 633-1379/ 665) -14.0 W23-W27 S25-S23
      8 1539  360 (S41W56 L=303)  1066     77        1 (1502/ 317-1568/ 412) -28.6 W51-W63 S43-S37
      9  571  309 (S57E43 L=204)  1219    136        1 ( 516/ 277- 627/ 351)  15.6 E57-E32 S60-S53
  

    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