Did you try GDAL? It comes with a python interface and supports a large number of exotic raster formats. I'm not sure if it supports 16bits per channel, but most probably it does. You should check it out.