Alan D. Snow
Alan D. Snow
> The background here is that I want to use rioxarray's clip() functionality on this and other netCDF files. But perhaps rioxarray is not the best way to achieve this....
Looks like your x and y coordinates are in some kind of projected coordinate system: ```python >>> xds.X10_153 array([-760., -750., -740., -730., -720., -710., -700., -690., -680., -670., -660., -650.,...
There are several options for the case where the data is unevently spaced and you want to convert it to a regular grid; - `pyresample`: https://pyresample.readthedocs.io/en/latest/swath.html (maybe @djhoese could assist...
`geocube`: Step 1: https://gis.stackexchange.com/questions/384581/raster-to-geopandas/384691#384691 ```python rds = xarray.open_dataset("path_to_file.nc") df = rds.squeeze().to_dataframe().reset_index() geometry = gpd.points_from_xy(df.x, df.y) gdf = gpd.GeoDataFrame(df, crs=rds.rio.crs, geometry=geometry) ``` Step 2: https://corteva.github.io/geocube/stable/examples/rasterize_point_data.html ```python from geocube.api.core import make_geocube from...
Not much point in adding nodatavals: https://github.com/mapbox/rasterio/issues/1859
Yes, it is also related to #296
What about this: https://xarray.pydata.org/en/stable/generated/xarray.Dataset.to_zarr.html GDAL 3.4 added support for Zarr. Or, are you specifically needing GeoTIff?
The xcog implementation looks pretty neat :+1:. My initial thoughts: - Would be fun to call the multi-file COG output format `czar` or `czarr` :smile: - Having it as its...
Related: - #525 - #573