### Python 实现 TIF 影像角点检测并将结果转换为 SHP 格式 #### 方法概述 在不依赖任何外部库的情况下，可以手动解析 TIFF 文件结构以提取其地理信息和像素数据。通过读取 TIFF 的头文件信息，可以获得影像的四个角点坐标（通常是经纬度）。随后，这些角点可以通过简单的文本格式保存为 ESRI Shapefile。 以下是具体实现： --- #### 提取 TIFF 角点坐标 TIFF 文件通常包含地理标记语言 (GeoTIFF)，其中存储了影像的空间参考信息以及角点位置。虽然这里不允许使用 `gdal` 或其他库，但仍可通过直接解析二进制文件的方式获取所需信息。 ```python def read_tif_header(file_path): """ 手动解析 TIFF 头部信息以获取角点坐标。 假设 TIFF 已经包含了 GeoTIFF Tags。 """ with open(file_path, 'rb') as f: # 读取 TIFF 魔数和字节顺序标志 magic_number = f.read(2) if magic_number not in [b'II', b'MM']: raise ValueError("Invalid TIFF file") byte_order = '<' if magic_number == b'II' else '>' # 跳过 IFD 数量和其他字段偏移量 offset_size = int.from_bytes(f.read(4), byteorder=byte_order) # 定位到第一个 IFD 条目 f.seek(offset_size) tags_count = int.from_bytes(f.read(2), byteorder=byte_order) corner_points = [] for _ in range(tags_count): tag_id = int.from_bytes(f.read(2), byteorder=byte_order) # 查找 ModelPixelScaleTag 和 ModelTiepointTag 等标签 if tag_id in [33922, 33550]: data_type = int.from_bytes(f.read(2), byteorder=byte_order) value_count = int.from_bytes(f.read(4), byteorder=byte_order) value_offset = int.from_bytes(f.read(4), byteorder=byte_order) current_pos = f.tell() f.seek(value_offset) values = [] for __ in range(value_count): val = float(int.from_bytes(f.read(8), byteorder=byte_order)) values.append(val) if tag_id == 33922: # ModelTiepointTag tiepoints = [(values[i], values[i + 1]) for i in range(0, len(values), 3)] corner_points.extend(tiepoints[:4]) f.seek(current_pos) return corner_points ``` 此函数会返回 TIFF 图像的四个角点坐标列表。 --- #### 将角点写入 SHP 文件 Shapefile 是由三个核心文件组成的数据集合：`.shp`, `.shx`, 和 `.dbf`。尽管完全手动生成这三个文件较为复杂，但可以通过简化版的 CSV 文本模拟 SHP 数据结构。 ```python def write_to_shp(corner_points, output_file): """ 将角点坐标写入伪 Shapefile 中。 这里仅生成一个简单的 CSV 文件作为替代方案。 """ header = "id,x,y\n" content = "" for idx, point in enumerate(corner_points): content += f"{idx},{point[0]},{point[1]}\n" with open(output_file.replace('.shp', '.csv'), 'w') as f: f.write(header + content) print("Successfully converted corners to pseudo-SHP format.") ``` 注意：这只是一个简易版本，实际应用中可能需要更复杂的编码来满足 GIS 软件的要求。 --- #### 综合代码示例 以下是一个完整的脚本框架，用于从 TIFF 提取角点并将其导出为伪 SHP 文件。 ```python import sys def main(): tif_file = sys.argv[1] shp_output = sys.argv[2] try: corners = read_tif_header(tif_file) write_to_shp(corners, shp_output) except Exception as e: print(f"Error occurred: {e}") if __name__ == "__main__": main() ``` 运行该脚本时需传入两个参数：输入的 TIFF 文件路径和输出的 SHP 文件路径。 --- #### 注意事项 - 上述方法假设 TIFF 文件已包含标准 GeoTIFF Tags。如果缺少必要元数据，则无法正确计算角点[^1]。 - 此处未涉及真正的 Shapefile 创建过程，而是采用 CSV 替代形式。对于严格需求，建议引入第三方库如 `pyshp` 或 `ogr`[^2]。 ---