获取影像footbound 范围

引入库 1

import gdal
import fiona
from fiona.crs import from_epsg
import numpy as np
from osgeo import ogr
from osgeo import gdal, ogr,osr
import numpy as np

保留有效范围的图形是raster  value 为1 的范围

def reservedata(inputshp='footprint.shp',outputshp='yxfootprint.shp'):
    import geopandas as gpd
 
    # 读取 Shapefile
    gdf = gpd.read_file(inputshp)
     
    # 根据 "value" 字段值为 1 来过滤数据
    filtered_gdf = gdf.loc[gdf['value'] == 1]
     
    # 保存过滤后的数据到新的 Shapefile
    filtered_gdf.to_file(outputshp)

栅格转换矢量

def raster_to_vector(input_raster, output_vector):
    raster_ds = gdal.Open(input_raster)
    #logging.info(input_raster)
    if raster_ds is None:
        raise ValueError("无法打开输入栅格文件。")

    driver = ogr.GetDriverByName("ESRI Shapefile")