|  | 
|  | prs | 
|  | 
|  | dest | 
|  | 
|  | action | 
|  | 
|  | default | 
|  | 
|  | help | 
|  | 
|  | opts | 
|  | 
|  | args | 
|  | 
|  | VERBOSE = opts.verbose | 
|  | 
|  | rt = rtreader.RasterReader(opts.db, opts.table, opts.column, opts.where) | 
|  | 
|  | logging | 
|  | 
| string | out_format = "GTiff" | 
|  | 
|  | out_driver = gdal.GetDriverByName(out_format) | 
|  | 
| def | out_data_type = pt2gdt(rt.pixel_types[0]) | 
|  | 
|  | out_ds = out_driver.Create(opts.output, rt.width, rt.height, rt.num_bands, out_data_type) | 
|  | 
|  | raster = numpy.zeros((rt.height, rt.width), pt2numpy(out_data_type)) | 
|  | Be careful!! Zeros function's input parameter can be a (height x width) array, not (width x height): http://docs.scipy.org/doc/numpy/reference/generated/numpy.zeros.html?highlight=zeros#numpy.zeros.  More... 
 | 
|  | 
|  | pixel = rt.get_value(b, width_index + 1, height_index + 1) | 
|  | 
|  | band = out_ds.GetRasterBand(b) | 
|  |