北京6环边界Geo范围的外边长与区域面积计算

from shapely.geometry import Polygon
import geopandas as gpd
poly = Polygon([(116.072, 39.714),
        (116.075, 39.705),
        (116.078, 39.695),
        (116.099, 39.688),
        (116.122, 39.688),
        (116.167, 39.685),
        (116.202, 39.684),
        (116.275, 39.682),
        (116.37, 39.701),
        (116.415, 39.711),
        (116.459, 39.715),
        (116.516, 39.729),
        (116.543, 39.734),
        (116.604, 39.749),
        (116.627, 39.769),
        (116.648, 39.788),
        (116.665, 39.806),
        (116.673, 39.815),
        (116.685, 39.824),
        (116.691, 39.832),
        (116.697, 39.841),
        (116.707, 39.864),
        (116.712, 39.874),
        (116.715, 39.884),
        (116.718, 39.904),
        (116.718, 39.922),
        (116.714, 39.938),
        (116.715, 39.952),
        (116.716, 39.965),
        (116.709, 39.993),
        (116.7, 40.005),
        (116.69, 40.016),
        (116.661, 40.057),
        (116.656, 40.064),
        (116.652, 40.069),
        (116.643, 40.078),
        (116.64, 40.082),
        (116.637, 40.086),
        (116.633, 40.091),
        (116.63, 40.103),
        (116.629, 40.107),
        (116.627, 40.112),
        (116.624, 40.121),
        (116.623, 40.132),
        (116.617, 40.142),
        (116.613, 40.145),
        (116.608, 40.146),
        (116.601, 40.147),
        (116.584, 40.149),
        (116.542, 40.152),
        (116.501, 40.161),
        (116.468, 40.164),
        (116.399, 40.164),
        (116.344, 40.178),
        (116.279, 40.181),
        (116.173, 40.175),
        (116.15, 40.161),
        (116.138, 40.126),
        (116.129, 40.086),
        (116.11, 40.066),
        (116.08, 40.029),
        (116.071, 39.98),
        (116.099, 39.94),
        (116.122, 39.909),
        (116.122, 39.885),
        (116.11, 39.859),
        (116.096, 39.814),
        (116.091, 39.797),
        (116.084, 39.787),
        (116.077, 39.75),
        (116.073, 39.734)])

crs = {'init': 'epsg:4326'}
gpoly = gpd.GeoSeries(poly, crs=crs)
print('gpoly.crs:', gpoly.crs)

gpoly = gpoly.to_crs(epsg=3395)# epsg=2436
print('gpoly.to_crs:', gpoly.crs)

poly = gpoly[0]
print('area(km*km):', poly.area/1.0e6)
print('length(km):', poly.length/1.0e3)

  北京6环边界Geo范围的外边长与区域面积计算

上一篇:awk小技巧之执行shell命令


下一篇:【DB笔试面试720】在Oracle中,OCR备份恢复有哪几种?