xarray 生成 tiff 文件

  1. 生成 tiff 文件
    复制代码
    from eccodes import *
    import numpy as np
    import xarray as xr 
    import rioxarray as rio 
    t2 = xr.Dataset(data_vars={'REF': (['longitude', 'latitude'], DATA)},
                    coords={
                        'longitude': (['longitude'], new_longrid[:, 0]),
                        'latitude': (['latitude'], new_latgrid[0]),
                    })
    print(t2)
    t2.rio.set_spatial_dims('latitude', 'longitude')
    t2.rio.set_crs("epsg:4326")
    t2.rio.to_raster(r"t2_raster.tiff")
    复制代码
  2. 可视化 tiff  
    复制代码
    import matplotlib.pyplot as plt
    from pylab import *
    import xarray as xr
    from matplotlib import rcParams
    
    dem = xr.open_rasterio('t2_raster.tiff')
    dem = dem[0] #getting the first band
    dem.plot(cmap='terrain')
    plt.show()
    复制代码

     

  

  参考 :https://www.heywhale.com/mw/project/622b2d808a84f900178933ac

                  https://www.heywhale.com/mw/project/628507931f8f7f5e2b45b284

                  

posted on   闹不机米  阅读(396)  评论(0编辑  收藏  举报

相关博文:
阅读排行:
· 百万级群聊的设计实践
· 永远不要相信用户的输入:从 SQL 注入攻防看输入验证的重要性
· 全网最简单!3分钟用满血DeepSeek R1开发一款AI智能客服,零代码轻松接入微信、公众号、小程
· .NET 10 首个预览版发布,跨平台开发与性能全面提升
· 《HelloGitHub》第 107 期

导航

< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5
点击右上角即可分享
微信分享提示