Changeset 4145 for spike/wktraster/scripts/rtrowdump.py
- Timestamp:
- 06/05/09 10:35:02 (3 years ago)
- Files:
-
- 1 modified
-
spike/wktraster/scripts/rtrowdump.py (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
spike/wktraster/scripts/rtrowdump.py
r4144 r4145 109 109 logit("Target GeoTIFF: %s" % opts.output) 110 110 111 raster_size = ( rt.width, rt.height ) 112 111 113 out_format = "GTiff" 112 114 out_driver = gdal.GetDriverByName(out_format) 113 115 out_data_type = pt2gdt(rt.pixel_types[0]) 114 out_ds = out_driver.Create(opts.output, r t.width, rt.height, rt.num_bands, out_data_type)116 out_ds = out_driver.Create(opts.output, raster_size[0], raster_size[1], rt.num_bands, out_data_type) 115 117 116 118 … … 118 120 logit("--- BAND %d ---------------------------------" % b) 119 121 120 raster = numpy.zeros( (rt.height, rt.width), pt2numpy(out_data_type))121 for y in range(0, r t.height):122 for x in range(0, r t.width):122 raster = numpy.zeros(raster_size, pt2numpy(out_data_type)) 123 for y in range(0, raster_size[1]): 124 for x in range(0, raster_size[0]): 123 125 pixel = rt.get_value(b, x + 1, y + 1) 124 126 raster[y, x] = pixel
