if 'int8' in im_data.: datatype = gdal.GDT_Byte elif 'int16' in im_data.: datatype = gdal.GDT_UInt16 else: datatype = gdal.GDT_Float32 if len(im_data.shape) == 3: im_bands, im_height, im_width = im_data.shape elif len(im_data.shape) == 2: im_data = np.array([i...
#gdal .GDT_UInt16, gdal.GDT_Int16, gdal.GDT_UInt32, gdal.GDT_Int32, #gdal.GDT_Float32, gdal.GDT_Float64 #判断栅格数据的数据类型 if 'int8' in im_data.dtype.name: datatype = gdal.GDT_Byte elif 'int16' in im_data.dtype.name: datatype = gdal.GDT_UInt16 else: datatype = gdal...
if 'int8' in im_data.dtype.name: datatype = gdal.GDT_Byte elif 'int16' in im_data.dtype.name: datatype = gdal.GDT_UInt16 else: datatype = gdal.GDT_Float32 if len(im_data.shape) == 3: im_bands, im_height, im_width = im_data.shape elif len(im_data.shape) == 2: im_d...
data = dataset.ReadAsArray() if data.dtype == 'uint16': driver = gdal.GetDriverByName('GTiff') out_dataset = driver.Create(save_file, cols, rows, 1, gdal.GDT_UInt16) out_dataset.SetProjection(projection) out_dataset.SetGeoTransform(trans) out_dataset.GetRasterBand(1).WriteArray(data) ...
datatype=gdal.GDT_Byteelif'int16'inim_data.dtype.name: datatype=gdal.GDT_UInt16else: datatype=gdal.GDT_Float32#判读数组维数iflen(im_data.shape) == 3: im_bands, im_height, im_width=im_data.shapeelse: im_bands, (im_height, im_width)= 1, im_data.shape#创建文件driver = gdal.Get...
###structval=rb.ReadRaster(px,py,1,1,buf_type=gdal.GDT_UInt16) 解释一下,px是算的,见上面公式,是坐标减去栅格最左值,除以像元大小,就是第几个像元了,同理,py;1,1是计算一个像元的意思,横着1,竖着1.。。。后面就是16位,,, 得到的结果可以存到csv里 :...
(0,0,0,0,0,0), im_proj=""):if'int8'indata.dtype.name: datatype = gdal.GDT_Byteelif'int16'indata.dtype.name: datatype = gdal.GDT_UInt16else: datatype = gdal.GDT_Float32iflen(data.shape) ==3: im_bands, im_height, im_width = data.shapeeliflen(data.shape) ==2: data =...
im_geotrans,im_data#写文件,以写成tif为例defwrite_img(self,filename,im_proj,im_geotrans,im_data):#gdal数据类型包括#gdal.GDT_Byte,#gdal .GDT_UInt16, gdal.GDT_Int16, gdal.GDT_UInt32, gdal.GDT_Int32,#gdal.GDT_Float32, gdal.GDT_Float64#判断栅格数据的数据类型if'int8'inim_data.dtype....
Byte elif 'int16' in im_data.dtype.name: datatype = gdal.GDT_UInt16 else: datatype = gdal.GDT_Float32 if len(im_data.shape) == 3: im_bands, im_height, im_width = im_data.shape elif len(im_data.shape) == 2: im_data = np.array([im_data]) im...
(im_data, im_geotrans, im_proj, path): if 'int8' in im_data.dtype.name: datatype = gdal.GDT_Byte elif 'int16' in im_data.dtype.name: datatype = gdal.GDT_UInt16 else: datatype = gdal.GDT_Float32 if len(im_data.shape) == 3: im_bands, im_height, im_width = im_data....