Я пытаюсь создать файл .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
основная причина, по которой я хотел скопировать, заключалась в том, чтобы доказать, что я могу читать, а затем записывать обновленные данные на основе вычислений.
как я мог бы записывать выходные данные с использованием цветовой шкалы вместо оттенки только одного цвета?