Вывод растра GDAL

Я пытаюсь создать файл .tif с помощью GDAL в Python. Он создает файл, но при каждом просмотре говорит: «Предварительный просмотр недоступен». Прямо сейчас я просто пытаюсь заставить его сделать копию входного файла. Вот мой код:



gdal.AllRegister()

inDs = gdal.Open("C:\\Documents and Settings\\patrick\\Desktop\\tiff elevation\\EBK1KM\\color_a2.tif")
if inDs is None:
  print 'Could not open image file'
  sys.exit(1)
else:
    print "successfully opened input file"

rows = inDs.RasterYSize
cols = inDs.RasterXSize
myband = inDs.GetRasterBand(1)
elev_data = myband.ReadAsArray(0,0,cols,rows)
driver = inDs.GetDriver()
outDs = driver.Create('C:\\Documents and Settings\\patrick\\Desktop\\tiff elevation\\EBK1KM\\new.tif', cols, rows, 1, GDT_Int32)

if outDs is None:
    print "couldn't open output file"
    sys.exit(1)

outBand = outDs.GetRasterBand(1)
outData = numpy.zeros((rows,cols),numpy.int16)
outBand.WriteArray(elev_data)
outBand.FlushCache()
outBand.SetNoDataValue(-99)
outDs.SetGeoTransform(inDs.GetGeoTransform())
outDs.SetProjection(inDs.GetProjection())
del outData

============================ обновление ============== =========================== сделал несколько открытий ... Я изучил способы преобразования из одного числового формата в другой с помощью статистической нормализации. Я обработал входные данные и преобразовал их в uint8, используя следующий алгоритм:
std = elev_data.std() #standard dev
avg = elev_data.mean()
arr = numpy.zeros((rows,cols),numpy.uint8)
for _i_ in _range_(_rows_):
    for _j_ in _range_(_cols_):
        arr[i,j] = (((out_elev[i,j]-avg)/std)*127)+128 #normalization formula
        #this puts all vals in range 1 to 255 (uint8)
dr = gdal.GetDriverByName("GTiff")
outDs = dr.Create("name",cols,rows,3,GDT_Byte) 
#creates and RGB file, accepts uint8 for input
outDs.GetRasterBand(1).WriteArray(arr) #write the output as shades of red
#this writes out a format viewable by microsoft products

основная причина, по которой я хотел скопировать, заключалась в том, чтобы доказать, что я могу читать, а затем записывать обновленные данные на основе вычислений.

как я мог бы записывать выходные данные с использованием цветовой шкалы вместо оттенки только одного цвета?

7
задан Pat 8 August 2011 в 16:35
поделиться