python地理空间分析_Karta

软件简介

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