search

mdsumner / ghrsst.coop

0 stars 0 forks source link

cf xarray #1

Open mdsumner opened 7 months ago

mdsumner commented 7 months ago

this round trips fine

import xarray as xr
from odc.geo.xr import assign_crs, write_cog
#from xarray import Dataset
data_file = "/rdsi/PUBLIC/raad/data/podaac-opendap.jpl.nasa.gov/opendap/allData/ghrsst/data/GDS2/L4/GLOB/JPL/MUR/v4.1/2002/152/20020601090000-JPL-L4_GHRSST-SSTfnd-MUR-GLOB-v02.0-fv04.1.nc"

data = xr.open_dataset(data_file, mask_and_scale=False)
data = assign_crs(dadata = assign_crs(data, crs="OGC:CRS84")
var = "analysed_sst"
cog_file = "/tmp/write_cog.tif"
write_cog(data[var], str(cog_file), overwrite=True)
ta, crs="OGC:CRS84")
var = "analysed_sst"
cog_file = "/tmp/write_cog.tif"
write_cog(data[var], str(cog_file), overwrite=True)