Karta Python 地理空间分析包

Karta 是一个开源的 Python 地理空间分析包，可简化处理矢量和栅格数据。

创建矢量地理空间：

point = Point((-130.0, 52.0), crs=LonLatWGS84)

line = read_geojson(\"linedata.json\")

polygon = Polygon([(-515005.78, -1301130.53),
                   (-579174.89, -1282271.94),
                   (-542977.83, -1221147.82),
                   (-437864.05, -1251641.55),
                   (-438160.72, -1252421.48),
                   (-437961.28, -1285314.00)],
                   crs=NSIDCNorth)

执行简单查询：

point2 = Point((-25.0, 48.0), crs=LonLatWGS84)
point.distance(point2)          # Distance in geographical units

line.intersects(polygon)        # True or False

ch = polygon.convex_hull()      # Returns a new polygon
ch.to_shapefile(\"poly.shp\")

加载可操作的格栅数据：

grid = read_gtiff(\"landsat_scene.tif\")  # Leverages GDAL

grid.profile(line)              # Collect data along a line

grid.resample(500.0, 500.0)     # Return a grid resampled at a new resolution