使用python读取mid/mif文件

1.使用geopandas读取:

   df = geopandas.read_file(path, encoding='UTF-8')
   gdf = geopandas.GeoDataFrame(df, geometry='geometry', crs='+init=epsg:4326')

此方法简单。但是缺点较为明显,当文件较大时,会比较吃内存。

2.使用fiona读取:

  下面这个方法直接适应fiona去读取mid/mif,此方法读取速度与1相同。但是占内存小很多。

def parse_mif(path, columns=None, **config):
    """
    读取mid/mif,shp文件,该函数占用内存相较于read_mif少
    :param path: 文件所在路径
    :param columns: mif文件所需要的列
    :param config: 配置,例如 编码
    :return:
    """

    import fiona
    from shapely.geometry import shape

    with fiona.open(path, **config) as features:
        crs = features.crs
        record_list = []
        logging.info("start")
        for line in features:
            if columns is None:
                d = {'geometry': shape(line['geometry']) if line['geometry'] else None}
                d.update(line['properties'])
                record_list.append(d)
            else:
                d = {'geometry': shape(line['geometry']) if line['geometry'] else None}
                for col in columns:
                    d[col] = line['properties'][col]
                record_list.append(d)

        logging.info("end")
        columns = columns + ["geometry"] if columns is not None else list(features.meta["schema"]["properties"]) + [
            "geometry"]
        gdf = gpn.GeoDataFrame(record_list, columns=columns)
        gdf.crs = crs
        return gdf

    return gpn.GeoDataFrame()

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值