Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1190 1466 (N24W08 L=218)   438     88        0 (1143/1444-1237/1484)   8.1 W04-W11 N22-N25
      2  287 1352 (N13E54 L=156)   456     64        1 ( 265/1317- 308/1385)  76.0 E56-E52 N10-N14
      3 1628 1149 (N07W38 L=248)   999    211        0 (1557/1073-1709/1237) -39.3 W32-W45 N03-N12
      4 1742 1211 (N12W48 L=258)   294     61        1 (1707/1189-1761/1231) -11.3 W44-W49 N11-N13
      5  822 1070 (S03E12 L=198)   515    213        2 ( 803/1006- 842/1140) -81.9 E13-E10 S06-N01
      6 1278 1005 (S04W15 L=226)   185     30        0 (1260/ 979-1294/1030) -56.3 W14-W16 S05-S02
      7 1373  917 (S09W22 L=232)   194     56        0 (1345/ 886-1405/ 946)  53.1 W20-W24 S10-S06
      8 1850  854 (S08W62 L=272)   771     97        1 (1819/ 815-1869/ 901)  63.4 W58-W64 S10-S04
      9 1263  860 (S13W16 L=226)   337    102        0 (1210/ 837-1311/ 876)  -9.5 W12-W18 S13-S11
     10 1549  731 (S18W37 L=247)  2092    355        0 (1386/ 621-1706/ 792)  29.7 W27-W49 S26-S13
     11  336  659 (S30E49 L=161)   703    145        0 ( 305/ 610- 367/ 714)  36.9 E53-E45 S32-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