根据上面的思路利用PYTHON实现
时间: 2024-03-17 12:44:27 浏览: 186
利用PYTHON编写的小程序
以下是利用 Python 实现批量获取卫星影像实际范围4个角点坐标的示例代码:
```python
# 导入必要的库
import rasterio
from rasterio.plot import show
from rasterio.warp import transform
# 定义函数,输入卫星影像文件路径和目标坐标系,输出实际范围4个角点坐标
def get_image_extent(image_path, dst_crs):
# 打开卫星影像
with rasterio.open(image_path) as src:
# 显示卫星影像
show(src)
# 获取卫星影像的四个角点坐标
lonlat_coords = [(src.bounds.left, src.bounds.bottom),
(src.bounds.right, src.bounds.bottom),
(src.bounds.right, src.bounds.top),
(src.bounds.left, src.bounds.top)]
# 将卫星影像的坐标转换为目标坐标系
dst_coords = transform(src.crs, dst_crs, lonlat_coords)
return dst_coords
# 测试函数,获取卫星影像实际范围坐标(WGS84经纬度)
image_path = "path/to/your/image.tif"
dst_crs = "EPSG:4326"
extent_coords = get_image_extent(image_path, dst_crs)
print(extent_coords)
```
代码中使用了 `rasterio` 库来读取和处理卫星影像,其中 `transform` 函数用于将卫星影像的坐标转换为目标坐标系,`show` 函数用于显示卫星影像。需要注意的是,代码中的 `dst_crs` 参数应该根据目标坐标系进行修改。
阅读全文