X-Git-Url: http://users.mur.at/ms/git/gitweb/?a=blobdiff_plain;f=aces_1.0.0%2Fpython%2Faces_ocio%2Fgenerate_lut.py;h=e83a7c5c154ef4defe67e494203411b94b9eb800;hb=d5c4b360c121a6b7ec0f3568049d8bfb9301ce1a;hp=a9b84c3a9c18efa85c940e9af58c29781278b4ec;hpb=89360f70a6c71121580324ed3c125a83c3973887;p=OpenColorIO-Configs.git diff --git a/aces_1.0.0/python/aces_ocio/generate_lut.py b/aces_1.0.0/python/aces_ocio/generate_lut.py index a9b84c3..e83a7c5 100755 --- a/aces_1.0.0/python/aces_ocio/generate_lut.py +++ b/aces_1.0.0/python/aces_ocio/generate_lut.py @@ -10,7 +10,6 @@ from __future__ import division import array import os -import sys import OpenImageIO as oiio @@ -25,6 +24,9 @@ __status__ = 'Production' __all__ = ['generate_1d_LUT_image', 'write_SPI_1d', + 'write_CSP_1d', + 'write_CTL_1d', + 'write_1d', 'generate_1d_LUT_from_image', 'generate_3d_LUT_image', 'generate_3d_LUT_from_image', @@ -78,7 +80,13 @@ def generate_1d_LUT_image(ramp_1d_path, ramp.close() -def write_SPI_1d(filename, from_min, from_max, data, entries, channels): +def write_SPI_1d(filename, + from_min, + from_max, + data, + entries, + channels, + components=3): """ Object description. @@ -96,24 +104,235 @@ def write_SPI_1d(filename, from_min, from_max, data, entries, channels): Return value description. """ + # May want to use fewer components than there are channels in the data + # Most commonly used for single channel LUTs + components = min(3, components, channels) + with open(filename, 'w') as fp: fp.write('Version 1\n') fp.write('From %f %f\n' % (from_min, from_max)) fp.write('Length %d\n' % entries) - fp.write('Components %d\n' % (min(3, channels))) + fp.write('Components %d\n' % components) fp.write('{\n') for i in range(0, entries): entry = '' - for j in range(0, min(3, channels)): + for j in range(0, components): entry = '%s %s' % (entry, data[i * channels + j]) fp.write(' %s\n' % entry) fp.write('}\n') +def write_CSP_1d(filename, + from_min, + from_max, + data, + entries, + channels, + components=3): + """ + Object description. + + Parameters + ---------- + parameter : type + Parameter description. + + Returns + ------- + type + Return value description. + """ + + # May want to use fewer components than there are channels in the data + # Most commonly used for single channel LUTs + components = min(3, components, channels) + + with open(filename, 'w') as fp: + fp.write('CSPLUTV100\n') + fp.write('1D\n') + fp.write('\n') + fp.write('BEGIN METADATA\n') + fp.write('END METADATA\n') + + fp.write('\n') + + fp.write('2\n') + fp.write('%f %f\n' % (from_min, from_max)) + fp.write('0.0 1.0\n') + fp.write('2\n') + fp.write('%f %f\n' % (from_min, from_max)) + fp.write('0.0 1.0\n') + fp.write('2\n') + fp.write('%f %f\n' % (from_min, from_max)) + fp.write('0.0 1.0\n') + + fp.write('\n') + + fp.write('%d\n' % entries) + if components == 1: + for i in range(0, entries): + entry = '' + for j in range(3): + entry = '%s %s' % (entry, data[i * channels]) + fp.write('%s\n' % entry) + else: + for i in range(entries): + entry = '' + for j in range(components): + entry = '%s %s' % (entry, data[i * channels + j]) + fp.write('%s\n' % entry) + fp.write('\n') + + +def write_CTL_1d(filename, + from_min, + from_max, + data, + entries, + channels, + components=3): + """ + Object description. + + Parameters + ---------- + parameter : type + Parameter description. + + Returns + ------- + type + Return value description. + """ + + # May want to use fewer components than there are channels in the data + # Most commonly used for single channel LUTs + components = min(3, components, channels) + + with open(filename, 'w') as fp: + fp.write('// %d x %d LUT generated by "generate_lut"\n' % ( + entries, components)) + fp.write('\n') + fp.write('const float min1d = %3.9f;\n' % from_min) + fp.write('const float max1d = %3.9f;\n' % from_max) + fp.write('\n') + + # Write LUT + if components == 1: + fp.write('const float lut[] = {\n') + for i in range(0, entries): + fp.write('%s' % data[i * channels]) + if i != (entries - 1): + fp.write(',') + fp.write('\n') + fp.write('};\n') + fp.write('\n') + else: + for j in range(components): + fp.write('const float lut%d[] = {\n' % j) + for i in range(0, entries): + fp.write('%s' % data[i * channels]) + if i != (entries - 1): + fp.write(',') + fp.write('\n') + fp.write('};\n') + fp.write('\n') + + fp.write('void main\n') + fp.write('(\n') + fp.write(' input varying float rIn,\n') + fp.write(' input varying float gIn,\n') + fp.write(' input varying float bIn,\n') + fp.write(' input varying float aIn,\n') + fp.write(' output varying float rOut,\n') + fp.write(' output varying float gOut,\n') + fp.write(' output varying float bOut,\n') + fp.write(' output varying float aOut\n') + fp.write(')\n') + fp.write('{\n') + fp.write(' float r = rIn;\n') + fp.write(' float g = gIn;\n') + fp.write(' float b = bIn;\n') + fp.write('\n') + fp.write(' // Apply LUT\n') + if components == 1: + fp.write(' r = lookup1D(lut, min1d, max1d, r);\n') + fp.write(' g = lookup1D(lut, min1d, max1d, g);\n') + fp.write(' b = lookup1D(lut, min1d, max1d, b);\n') + elif components == 3: + fp.write(' r = lookup1D(lut0, min1d, max1d, r);\n') + fp.write(' g = lookup1D(lut1, min1d, max1d, g);\n') + fp.write(' b = lookup1D(lut2, min1d, max1d, b);\n') + fp.write('\n') + fp.write(' rOut = r;\n') + fp.write(' gOut = g;\n') + fp.write(' bOut = b;\n') + fp.write(' aOut = aIn;\n') + fp.write('}\n') + + +def write_1d(filename, + from_min, + from_max, + data, + data_entries, + data_channels, + lut_components=3, + format='spi1d'): + """ + Object description. + + Parameters + ---------- + parameter : type + Parameter description. + + Returns + ------- + type + Return value description. + """ + + ocio_formats_to_extensions = {'cinespace': 'csp', + 'flame': '3dl', + 'icc': 'icc', + 'houdini': 'lut', + 'lustre': '3dl', + 'ctl': 'ctl'} + + if format in ocio_formats_to_extensions: + if ocio_formats_to_extensions[format] == 'csp': + write_CSP_1d(filename, + from_min, + from_max, + data, + data_entries, + data_channels, + lut_components) + elif ocio_formats_to_extensions[format] == 'ctl': + write_CTL_1d(filename, + from_min, + from_max, + data, + data_entries, + data_channels, + lut_components) + else: + write_SPI_1d(filename, + from_min, + from_max, + data, + data_entries, + data_channels, + lut_components) + + def generate_1d_LUT_from_image(ramp_1d_path, output_path=None, min_value=0, - max_value=1): + max_value=1, + channels=3, + format='spi1d'): """ Object description. @@ -133,16 +352,17 @@ def generate_1d_LUT_from_image(ramp_1d_path, ramp = oiio.ImageInput.open(ramp_1d_path) - spec = ramp.spec() - width = spec.width - channels = spec.nchannels + ramp_spec = ramp.spec() + ramp_width = ramp_spec.width + ramp_channels = ramp_spec.nchannels # Forcibly read data as float, the Python API doesn't handle half-float # well yet. type = oiio.FLOAT - data = ramp.read_image(type) + ramp_data = ramp.read_image(type) - write_SPI_1d(output_path, min_value, max_value, data, width, channels) + write_1d(output_path, min_value, max_value, + ramp_data, ramp_width, ramp_channels, channels, format) def generate_3d_LUT_image(ramp_3d_path, resolution=32): @@ -173,7 +393,10 @@ def generate_3d_LUT_image(ramp_3d_path, resolution=32): lut_extract.execute() -def generate_3d_LUT_from_image(ramp_3d_path, output_path=None, resolution=32): +def generate_3d_LUT_from_image(ramp_3d_path, + output_path=None, + resolution=32, + format='spi3d'): """ Object description. @@ -189,21 +412,58 @@ def generate_3d_LUT_from_image(ramp_3d_path, output_path=None, resolution=32): """ if output_path is None: - output_path = '%s.%s' % (ramp_3d_path, 'spi1d') + output_path = '%s.%s' % (ramp_3d_path, 'spi3d') + + ocio_formats_to_extensions = {'cinespace': 'csp', + 'flame': '3dl', + 'icc': 'icc', + 'houdini': 'lut', + 'lustre': '3dl'} + + if format == 'spi3d' or not (format in ocio_formats_to_extensions): + # Extract a spi3d LUT + args = ['--extract', + '--cubesize', + str(resolution), + '--maxwidth', + str(resolution * resolution), + '--input', + ramp_3d_path, + '--output', + output_path] + lut_extract = Process(description='extract a 3d LUT', + cmd='ociolutimage', + args=args) + lut_extract.execute() - args = ['--extract', - '--cubesize', - str(resolution), - '--maxwidth', - str(resolution * resolution), - '--input', - ramp_3d_path, - '--output', - output_path] - lut_extract = Process(description='extract a 3d LUT', - cmd='ociolutimage', - args=args) - lut_extract.execute() + else: + output_path_spi3d = '%s.%s' % (output_path, 'spi3d') + + # Extract a spi3d LUT + args = ['--extract', + '--cubesize', + str(resolution), + '--maxwidth', + str(resolution * resolution), + '--input', + ramp_3d_path, + '--output', + output_path_spi3d] + lut_extract = Process(description='extract a 3d LUT', + cmd='ociolutimage', + args=args) + lut_extract.execute() + + # Convert to a different format + args = ['--lut', + output_path_spi3d, + '--format', + format, + output_path] + lut_convert = Process(description='convert a 3d LUT', + cmd='ociobakelut', + args=args) + lut_convert.execute() def apply_CTL_to_image(input_image, @@ -289,14 +549,16 @@ def convert_bit_depth(input_image, output_image, depth): def generate_1d_LUT_from_CTL(lut_path, ctl_paths, lut_resolution=1024, - identity_LUT_bit_depth='half', + identity_lut_bit_depth='half', input_scale=1, output_scale=1, global_params=None, cleanup=True, aces_ctl_directory=None, min_value=0, - max_value=1): + max_value=1, + channels=3, + format='spi1d'): """ Object description. @@ -316,43 +578,45 @@ def generate_1d_LUT_from_CTL(lut_path, lut_path_base = os.path.splitext(lut_path)[0] - identity_LUT_image_float = '%s.%s.%s' % (lut_path_base, 'float', 'tiff') - generate_1d_LUT_image(identity_LUT_image_float, + identity_lut_image_float = '%s.%s.%s' % (lut_path_base, 'float', 'tiff') + generate_1d_LUT_image(identity_lut_image_float, lut_resolution, min_value, max_value) - if identity_LUT_bit_depth != 'half': - identity_LUT_image = '%s.%s.%s' % (lut_path_base, 'uint16', 'tiff') - convert_bit_depth(identity_LUT_image_float, - identity_LUT_image, - identity_LUT_bit_depth) + if identity_lut_bit_depth not in ['half', 'float']: + identity_lut_image = '%s.%s.%s' % (lut_path_base, 'uint16', 'tiff') + convert_bit_depth(identity_lut_image_float, + identity_lut_image, + identity_lut_bit_depth) else: - identity_LUT_image = identity_LUT_image_float + identity_lut_image = identity_lut_image_float - transformed_LUT_image = '%s.%s.%s' % (lut_path_base, 'transformed', 'exr') - apply_CTL_to_image(identity_LUT_image, - transformed_LUT_image, + transformed_lut_image = '%s.%s.%s' % (lut_path_base, 'transformed', 'exr') + apply_CTL_to_image(identity_lut_image, + transformed_lut_image, ctl_paths, input_scale, output_scale, global_params, aces_ctl_directory) - generate_1d_LUT_from_image(transformed_LUT_image, + generate_1d_LUT_from_image(transformed_lut_image, lut_path, min_value, - max_value) + max_value, + channels, + format) if cleanup: - os.remove(identity_LUT_image) - if identity_LUT_image != identity_LUT_image_float: - os.remove(identity_LUT_image_float) - os.remove(transformed_LUT_image) + os.remove(identity_lut_image) + if identity_lut_image != identity_lut_image_float: + os.remove(identity_lut_image_float) + os.remove(transformed_lut_image) -def correct_LUT_image(transformed_LUT_image, - corrected_LUT_image, +def correct_LUT_image(transformed_lut_image, + corrected_lut_image, lut_resolution): """ Object description. @@ -368,7 +632,7 @@ def correct_LUT_image(transformed_LUT_image, Return value description. """ - transformed = oiio.ImageInput.open(transformed_LUT_image) + transformed = oiio.ImageInput.open(transformed_lut_image) transformed_spec = transformed.spec() width = transformed_spec.width @@ -382,14 +646,14 @@ def correct_LUT_image(transformed_LUT_image, height, lut_resolution * lut_resolution, lut_resolution)) - print('Generating %s' % corrected_LUT_image) + print('Generating %s' % corrected_lut_image) # Forcibly read data as float, the Python API doesn't handle half-float # well yet. type = oiio.FLOAT source_data = transformed.read_image(type) - correct = oiio.ImageOutput.create(corrected_LUT_image) + correct = oiio.ImageOutput.create(corrected_lut_image) correct_spec = oiio.ImageSpec() correct_spec.set_format(oiio.FLOAT) @@ -397,7 +661,7 @@ def correct_LUT_image(transformed_LUT_image, correct_spec.height = width correct_spec.nchannels = channels - correct.open(corrected_LUT_image, correct_spec, oiio.Create) + correct.open(corrected_lut_image, correct_spec, oiio.Create) dest_data = array.array('f', ('\0' * correct_spec.width * @@ -417,22 +681,23 @@ def correct_LUT_image(transformed_LUT_image, correct.close() else: # shutil.copy(transformedLUTImage, correctedLUTImage) - corrected_LUT_image = transformed_LUT_image + corrected_lut_image = transformed_lut_image transformed.close() - return corrected_LUT_image + return corrected_lut_image def generate_3d_LUT_from_CTL(lut_path, ctl_paths, lut_resolution=64, - identity_LUT_bit_depth='half', + identity_lut_bit_depth='half', input_scale=1, output_scale=1, global_params=None, cleanup=True, - aces_ctl_directory=None): + aces_ctl_directory=None, + format='spi3d'): """ Object description. @@ -452,42 +717,48 @@ def generate_3d_LUT_from_CTL(lut_path, lut_path_base = os.path.splitext(lut_path)[0] - identity_LUT_image_float = '%s.%s.%s' % (lut_path_base, 'float', 'tiff') - generate_3d_LUT_image(identity_LUT_image_float, lut_resolution) + identity_lut_image_float = '%s.%s.%s' % (lut_path_base, 'float', 'tiff') + generate_3d_LUT_image(identity_lut_image_float, lut_resolution) - if identity_LUT_bit_depth != 'half': - identity_LUT_image = '%s.%s.%s' % (lut_path_base, - identity_LUT_bit_depth, + if identity_lut_bit_depth not in ['half', 'float']: + identity_lut_image = '%s.%s.%s' % (lut_path_base, + identity_lut_bit_depth, 'tiff') - convert_bit_depth(identity_LUT_image_float, - identity_LUT_image, - identity_LUT_bit_depth) + convert_bit_depth(identity_lut_image_float, + identity_lut_image, + identity_lut_bit_depth) else: - identity_LUT_image = identity_LUT_image_float + identity_lut_image = identity_lut_image_float - transformed_LUT_image = '%s.%s.%s' % (lut_path_base, 'transformed', 'exr') - apply_CTL_to_image(identity_LUT_image, - transformed_LUT_image, + transformed_lut_image = '%s.%s.%s' % (lut_path_base, 'transformed', 'exr') + apply_CTL_to_image(identity_lut_image, + transformed_lut_image, ctl_paths, input_scale, output_scale, global_params, aces_ctl_directory) - corrected_LUT_image = '%s.%s.%s' % (lut_path_base, 'correct', 'exr') - corrected_LUT_image = correct_LUT_image(transformed_LUT_image, - corrected_LUT_image, + corrected_lut_image = '%s.%s.%s' % (lut_path_base, 'correct', 'exr') + corrected_lut_image = correct_LUT_image(transformed_lut_image, + corrected_lut_image, lut_resolution) - generate_3d_LUT_from_image(corrected_LUT_image, lut_path, lut_resolution) + generate_3d_LUT_from_image(corrected_lut_image, + lut_path, + lut_resolution, + format) if cleanup: - os.remove(identity_LUT_image) - if identity_LUT_image != identity_LUT_image_float: - os.remove(identity_LUT_image_float) - os.remove(transformed_LUT_image) - if corrected_LUT_image != transformed_LUT_image: - os.remove(corrected_LUT_image) + os.remove(identity_lut_image) + if identity_lut_image != identity_lut_image_float: + os.remove(identity_lut_image_float) + os.remove(transformed_lut_image) + if corrected_lut_image != transformed_lut_image: + os.remove(corrected_lut_image) + if format != 'spi3d': + lut_path_spi3d = '%s.%s' % (lut_path, 'spi3d') + os.remove(lut_path_spi3d) def main(): @@ -514,6 +785,7 @@ def main(): usage='%prog [options]') p.add_option('--lut', '-l', type='string', default='') + p.add_option('--format', '-f', type='string', default='') p.add_option('--ctl', '-c', type='string', action='append') p.add_option('--lutResolution1d', '', type='int', default=1024) p.add_option('--lutResolution3d', '', type='int', default=33) @@ -533,9 +805,10 @@ def main(): options, arguments = p.parse_args() lut = options.lut + format = options.format ctls = options.ctl - lut_resolution_1d = options.lut_resolution_1d - lut_resolution_3d = options.lut_resolution_3d + lut_resolution_1d = options.lutResolution1d + lut_resolution_3d = options.lutResolution3d min_value = options.minValue max_value = options.maxValue input_scale = options.inputScale @@ -551,19 +824,13 @@ def main(): for param in options.ctlRenderParam: params[param[0]] = float(param[1]) - try: - args_start = sys.argv.index('--') + 1 - args = sys.argv[args_start:] - except: - args_start = len(sys.argv) + 1 - args = [] - if generate_1d: print('1D LUT generation options') else: print('3D LUT generation options') print('lut : %s' % lut) + print('format : %s' % format) print('ctls : %s' % ctls) print('lut res 1d : %s' % lut_resolution_1d) print('lut res 3d : %s' % lut_resolution_3d) @@ -587,7 +854,8 @@ def main(): cleanup, ctl_release_path, min_value, - max_value) + max_value, + format=format) elif generate_3d: generate_3d_LUT_from_CTL(lut, @@ -598,7 +866,8 @@ def main(): output_scale, params, cleanup, - ctl_release_path) + ctl_release_path, + format=format) else: print(('\n\nNo LUT generated. ' 'You must choose either 1D or 3D LUT generation\n\n')) @@ -606,4 +875,3 @@ def main(): if __name__ == '__main__': main() -