参考:Cutting a polygon from TIFF with Python [closed] import rasterio import rasterio.mask import geopandas as gpd dataset = rasterio.open("wc2.1_10m_elev.tif") gdf_africa = gpd.read_file("africa1_map.gpkg") poly = gdf_africa.loc[0, "geometry"] # out_image is the tiff image cut ...
Lat_Res = (LatMax - LatMin) / (float(N_Lat)-1) # 3.3 构建.tiff文件框架 spei_ds = gdal.GetDriverByName('Gtiff').Create(OutTif,N_Lon,N_Lat,1,gdal.GDT_Float32) # 3.4 设置影像的显示范围 geotransform = (LonMin,Lon_Res, 0, LatMin, 0, Lat_Res) spei_ds.SetGeoTransform(geotransfo...
/** 读取一定范围的tiff影像,将数据保存在Mat结构中 @param filename 遥感影像的路径 @param dstMat 要保存数据的Mat @param bufWidth 保存数据数Mat的列数 @param bufHeight 保存数组Mat的行数 @param startCol 读取的起始列 @param startRow 读取的起始行 @param startCol 读取的列数 @param startRow 读取...
out_meta = ds.meta.copy() out_meta.update({"driver": "GTiff", "height": out_img.shape[1], "width": out_img.shape[2], "transform": out_transform, "crs":ds.crs} ) with rasterio.open('out.tif', "w", **out_meta) as dest: dest.write(out_img) from rasterio.plot import sh...
(window) # 读取栅格数据 grid_data = src.read(window=window) # 创建输出文件路径 output_file = os.path.join(output_dir, f"grid_{i+1}_{j+1}.tif") if os.path.exists(output_file): break # 保存裁剪后的栅格数据,使用指定压缩方式 with rasterio.open( output_file, "w", driver="GTiff"...
例如:`driver = 'GTiff'`表示使用GeoTIFF格式,`driver = 'HDF5'`表示使用HDF5格式。 - `count`: 指定栅格数据集中的波段数量。例如:`count = 3`表示栅格数据集包含三个波段。 2. 数据集属性参数 - `bounds`: 返回栅格数据集的范围,即左下角和右上角的坐标。例如:`bounds = (xmin, ymin, xmax, y...
# 将Shapely几何对象转换为rasterio可以理解的边界clipped = src.read(window=window)# 读取并裁剪TIFF数据out_meta = src.meta.copy()out_meta.update({"driver":"GTiff","compress":'lzw'})# 更新元数据,rle,lzw等delsrcwithrasterio.open(output_file,'w', **out_meta)asdest:# 写入裁剪后的TIFF文件...
driver:文件的格式,一般在r和r+模式会省略,因为函数会自动获取其文件后缀判别。若是创建文件则需要指定该参数。若是创建GeoTIFF文件,则需要指定driver=GTiff。此处的driver与gdal类似,可以参照:Raster drivers — GDAL documentation获取不同文件格式的driver参数; ...
该函数可以处理多种栅格图像格式,例如:GeoTIFF、JPEG、PNG 等。通过使用 gdal rasterio 函数,用户可以方便地对栅格图像数据进行操作,满足各种地理信息处理需求。 二、gdal rasterio 函数的功能 gdal rasterio 函数的主要功能包括: 1.读取栅格图像数据:可以使用 gdal rasterio 函数从文件中读取栅格图像数据,并返回一个 ...
1. 读取栅格数据 读取栅格数据的最常见方式是使用rasterio.open()打开文件,然后通过.read()方法读取图像数据。 示例代码 importrasterio# 打开 GeoTIFF 文件withrasterio.open('example.tif')assrc:# 读取所有波段的栅格数据data=src.read()# 获取栅格的元数据print("CRS:",src.crs)print("Width, Height:",src...