{"nbformat":4,"cells":[{"execution_count":null,"source":"import os\nimport numpy as np\nimport tifffile as tiff\nimport cv2\nimport matplotlib.pyplot as plt\n\n\ndef stretch_8bit(bands, lower_percent=2, higher_percent=98):\n    out = np.zeros_like(bands)\n    for i in range(3):\n        a = 0 \n        b = 255 \n        c = np.percentile(bands[:,:,i], lower_percent)\n        d = np.percentile(bands[:,:,i], higher_percent)        \n        t = a + (bands[:,:,i] - c) * (b - a) / (d - c)    \n        t[t<a] = a\n        t[t>b] = b\n        out[:,:,i] =t\n    return out.astype(np.uint8)    \n    \ndef M(image_id):\n    filename = os.path.join('..', 'input', 'sixteen_band', '{}_M.tif'.format(image_id))\n    img = tiff.imread(filename)    \n    img = np.rollaxis(img, 0, 3)\n    return img\n\nimage_id = '6120_2_2'\nm = M(image_id)\nimg = np.zeros((837,851,3))\nimg[:,:,0] = m[:,:,0] #red\nimg[:,:,1] = m[:,:,4] #green\nimg[:,:,2] = m[:,:,2] #blue\nplt.imshow(stretch_8bit(img))\n","outputs":[],"metadata":{"_cell_guid":"0f245602-1b68-4fb9-f25b-959673f80fa6","collapsed":true,"_uuid":"f2c9f6108ff43fcc157ab3a0e7da1eb7ccdb1c19"},"cell_type":"code"}],"metadata":{"_change_revision":0,"kernelspec":{"name":"python3","display_name":"Python 3","language":"python"},"_is_fork":false,"language_info":{"version":"3.6.4","mimetype":"text/x-python","pygments_lexer":"ipython3","file_extension":".py","codemirror_mode":{"version":3,"name":"ipython"},"nbconvert_exporter":"python","name":"python"}},"nbformat_minor":1}