-
创建 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')