Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1598 1648 (N22W64 L=210)  1132    130        2 (1553/1602-1655/1693)   0.0 W55-W76 N20-N22
      2  360 1539 (N57E46 L=100)   454     45        0 ( 334/1520- 397/1554)  -5.2 E51-E39 N55-N58
      3  779 1455 (N37E03 L=144)   352    131        0 ( 706/1438- 833/1470) -29.7 E07-W00 N34-N39
      4 1825 1397 (N01W68 L=215)  1116    118        1 (1799/1333-1849/1460)  50.2 W63-W72 S02-N05
      5  155  866 (N17E69 L= 77)  1394    117        4 ( 126/ 806- 184/ 936)  81.9 E74-E64 N13-N20
      6 1359  781 (S17W12 L=158)   650    294        0 (1245/ 713-1464/ 877)   6.7 W04-W21 S18-S14
      7  999  736 (S09E09 L=137)   656    251        1 ( 921/ 630-1053/ 835) -85.9 E11-E07 S16-S02
      8  956  580 (S17E16 L=130)   221    108        0 ( 935/ 548- 970/ 636)  51.3 E18-E14 S18-S13
      9 1449  619 (S30W14 L=161)   318     72        0 (1409/ 601-1485/ 635) -18.4 W11-W17 S31-S27
     10 1104  413 (S32E14 L=132)   973    186        1 (1054/ 329-1143/ 476)  67.4 E18-E11 S38-S26
  

    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