Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  290 1570 (N12E73 L=277)  1662    170        1 ( 251/1527- 323/1616)  60.3 E78-E70 N08-N15
      2  226 1196 (S12E60 L=291)  1421    254        4 ( 176/1100- 282/1271)  76.0 E64-E54 S18-S06
      3 1828 1064 (N20W53 L= 44)  1198    248        2 (1765/1038-1901/1108)  11.3 W47-W62 N18-N22
      4 1412  977 (N03W23 L= 14)   236     49        0 (1395/ 948-1427/1009) -76.0 W22-W24 N01-N05
      5 1687  883 (N06W45 L= 35)   411     73        0 (1656/ 836-1715/ 921)  71.6 W43-W46 N02-N08
      6  900  668 (S27W02 L=352)   451    154        1 ( 806/ 630- 950/ 695)   5.7 E04-W06 S28-S26
      7 1727  575 (S10W62 L= 52)  1565    269        7 (1696/ 512-1766/ 629)  90.0 W57-W70 S13-S05
      8 1614  533 (S15W54 L= 44)   592    102        3 (1578/ 501-1647/ 557)  45.0 W50-W56 S17-S13
      9 1082  386 (S40W25 L= 15)  1939    386        2 ( 956/ 330-1252/ 490)   0.0 W12-W38 S43-S36
     10 1594  367 (S25W63 L= 53)  1049    248        0 (1539/ 352-1663/ 385)  18.4 W57-W70 S26-S22
  

    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