geopandas 生成 geojson 文件

1. 创建 GeoDataFrame  输出 geojson 文件

   import geopandas
   
   nan_index = np.where(~np.isnan(data.flatten()))
   # print(data[nan_index].shape)
   
   ss = np.stack((lon.flatten(), lat.flatten()), 1)
   
   ss1 = [Point(i.tolist()) for i in ss[nan_index]]
   print(ss)
   # print(ss1)
   cq = geopandas.GeoDataFrame({'dbz': data.flatten()[nan_index], 'geometry': ss1}, crs='EPSG:4326')
   print(cq)
   cq.to_file('output1.geojson',
              driver='GeoJSON',  # 默认 geojson  可以换成 ESRI Shapefile
              encoding='utf-8')