shp2tif

:: Ice free
gdal_rasterize -of GTiff -ot Int16 -a_nodata 255 -a COLORCT -where "COLORCT='01'" -burn 000 -burn 100 -burn 255 -ts 1000 1000 -l aari_kar_20140107_pl_a D:\data\proton\aari_kar_20140107_pl_a.shp D:\data\proton\aari_kar_20140107_pl_a.tif

:: Open water
gdal_rasterize -a COLORCT -where "COLORCT='02'" -b 1 -b 2 -b 3 -burn 150 -burn 200 -burn 255 -l aari_kar_20140107_pl_a D:\data\proton\aari_kar_20140107_pl_a.shp D:\data\proton\aari_kar_20140107_pl_a.tif

:: Very open ice
gdal_rasterize -a COLORCT -where "COLORCT='04'" -b 1 -b 2 -b 3 -burn 140 -burn 255 -burn 160 -l aari_kar_20140107_pl_a D:\data\proton\aari_kar_20140107_pl_a.shp D:\data\proton\aari_kar_20140107_pl_a.tif

:: Open ice
gdal_rasterize -a COLORCT -where "COLORCT='05'" -b 1 -b 2 -b 3 -burn 255 -burn 255 -burn 000 -l aari_kar_20140107_pl_a D:\data\proton\aari_kar_20140107_pl_a.shp D:\data\proton\aari_kar_20140107_pl_a.tif

:: Close ice
gdal_rasterize -a COLORCT -where "COLORCT='06'" -b 1 -b 2 -b 3 -burn 255 -burn 125 -burn 007 -l aari_kar_20140107_pl_a D:\data\proton\aari_kar_20140107_pl_a.shp D:\data\proton\aari_kar_20140107_pl_a.tif

:: Very close ice
gdal_rasterize -a COLORCT -where "COLORCT='07'" -b 1 -b 2 -b 3 -burn 255 -burn 000 -burn 000 -l aari_kar_20140107_pl_a D:\data\proton\aari_kar_20140107_pl_a.shp D:\data\proton\aari_kar_20140107_pl_a.tif

:: Fast ice
gdal_rasterize -a COLORSA -where "COLORSA='19'" -b 1 -b 2 -b 3 -burn 150 -burn 150 -burn 150 -l aari_kar_20140107_pl_a D:\data\proton\aari_kar_20140107_pl_a.shp D:\data\proton\aari_kar_20140107_pl_a.tif

:: Optional

:: New ice
::gdal_rasterize -a COLORCT -where "COLORCT='11'" -b 1 -b 2 -b 3 -burn 255 -burn 175 -burn 255 -l aari_kar_20140107_pl_a D:\data\proton\aari_kar_20140107_pl_a.shp D:\data\proton\aari_kar_20140107_pl_a.tif

:: Nilas
::gdal_rasterize -a COLORCT -where "COLORCT='12'" -b 1 -b 2 -b 3 -burn 255 -burn 100 -burn 255 -l aari_kar_20140107_pl_a D:\data\proton\aari_kar_20140107_pl_a.shp D:\data\proton\aari_kar_20140107_pl_a.tif


import os

# The shapefile to be rasterized:
shapefile = 'aari_kar_20140107_pl_a.shp'
print 'Rasterize ' + shapefile
 #get path and filename seperately
(shapefilefilepath, shapefilename) = os.path.split(shapefile)
#get file name without extension
(shapefileshortname, extension) = os.path.splitext(shapefilename)

# The land area to be masked out, also being a shapefile to be rasterized
SvalbardCoast = 's100-landp_3575.shp'

# The raster file to be created and receive the rasterized shapefile
outrastername = shapefileshortname + '.tif'
outraster = '' + outrastername


ss = 'gdal_rasterize -a COLORCT -where \"COLORCT=\'02\'\" -burn 2 -l ' + shapefileshortname +' -tr 1000 -1000 ' +  shapefile + ' ' + outraster

# Rasterize first Ice Type and at same time create file -- call gdal_rasterize commandline
print '\n Open Water'
os.system('gdal_rasterize -a COLORCT -where \"COLORCT=\'02\'\" -burn 2 -l ' + shapefileshortname +' -tr 1000 -1000 ' +  shapefile + ' ' + outraster)

# Rasterize the other Ice types, adding them to the already created file

print '\nVery Open Drift Ice'
os.system('gdal_rasterize -a COLORCT -where \"COLORCT=\'04\'\" -b 1 -burn 3 -l ' + shapefileshortname +' ' +  shapefile + ' ' + outraster)

print '\n Open Drift Ice'
os.system('gdal_rasterize -a COLORCT -where \"COLORCT=\'05\'\" -b 1 -burn 4 -l ' + shapefileshortname +' ' +  shapefile + ' ' + outraster)

print '\n Close Drift Ice'
os.system('gdal_rasterize -a COLORCT -where \"COLORCT=\'06\'\" -b 1 -burn 5 -l ' + shapefileshortname +' ' +  shapefile + ' ' + outraster)

print '\n Very Close Drift Ice'
os.system('gdal_rasterize -a COLORCT -where \"COLORCT=\'07\'\" -b 1 -burn 6 -l ' + shapefileshortname +' ' +  shapefile + ' ' + outraster)

print '\n Fast Ice'
os.system('gdal_rasterize -a COLORCT -where \"COLORCT=\'19\'\" -b 1 -burn 1 -l ' + shapefileshortname +' ' +  shapefile + ' ' + outraster)

# Rasterize Spitsbergen land area on top
#print '\n SvalbardRaster'
#os.system('gdal_rasterize  -b 1 -burn 8 -l s100-landp_3575 '  +  SvalbardCoast + ' ' + outraster)

print "\n Done"