Cycles Volume Scattering with Voxel Data
For those wondering since there are not really any supported volume formats (as of spring 2016) supported in Blender Cycles hacks must be used. I settled on a solution that used the answer posted in a related question. It is important to note this solution is a general method to create 3d image from a 2d image and is gpu accelerated. I was able to use this method to create 3d volume scattering.



Here is an ugly but fully working one that I produced with the following python code to create the input image. I found this approach to also be very memory intensive as expected since each image was 200x200x200 pixels.

import numpy as np
import png # using pypng

n = 200
x = np.linspace(0, 10, n)
array = np.zeros(n**3).reshape(n, n**2)
for i in range(n):
    for j in range(n):
        for k in range(n):
            array[i, j+k*n] = (np.sin(x[i]) * np.sin(x[j]) * 
np.sin(x[k]))**2 * 2**16
array = np.floor(array).astype(np.uint16)
image = png.from_array(array, 'L')
image.save('hello.png')