tif栅格转geojson

### 将 TIF 栅格文件转换为 GeoJSON 文件

为了将 TIF 格式的栅格数据转换成 GeoJSON 文件，可以采用 ArcGIS 软件中的工具完成这一过程。具体而言，在 ArcGIS 中先利用 “Raster to Polygon” 工具把 TIF 文件转化为矢量面要素[^1]。

import arcpy

# 设置工作空间
arcpy.env.workspace = "C:/data"

# 输入栅格数据集
in_raster = "input.tif"
# 输出的矢量要素类
out_polygon_features = "output.shp"

# 执行 Raster To Polygon 转换
arcpy.RasterToPolygon_conversion(in_raster, out_polygon_features, "NO_SIMPLIFY", "VALUE")

一旦完成了从栅格到矢量的转变之后，则可以通过使用其他 GIS 库比如 Fiona 和 Shapely 来进一步处理这些矢量数据，并最终将其保存为 GeoJSON 文件格式[^2]。

import fiona
from shapely.geometry import shape, mapping
import json

# 定义要读取的 Shapefile 的路径以及目标 CRS (坐标参照系)
src_file = 'output.shp'
dst_crs = {'init': 'epsg:4326'}

with fiona.open(src_file) as src:
    meta = src.meta.copy()
    meta['crs'] = dst_crs
    
    features = []
    
    for feat in src:
        geom = shape(feat['geometry'])
        
        properties = {
            key: value 
            for key, value in feat["properties"].items()
        }
        
        feature = {
            "type": "Feature",
            "id": feat.get('id'),
            "geometry": mapping(geom),
            "properties": properties,
        }

        features.append(feature)

geojson_data = {"type": "FeatureCollection", "features": features}

with open("result.geojson", "w") as dest:
    json.dump(geojson_data, dest)

上述 Python 代码片段展示了如何加载由 ArcGIS 创建的 Shapefile 并将其导出为 GeoJSON 文件。此过程中涉及到的空间操作依赖于第三方库 `Fiona` 和 `Shapely`，它们提供了便捷的方法来进行地理信息系统的数据分析与处理。