import xarray as xr
import geemap
import numpy as np
from mayavi import mlab
# Step 1: Download the NetCDF file
url = 'https://github.com/giswqs/leafmap/raw/master/examples/data/wind_global.nc'
filename = 'wind_global.nc'
data = geemap.read_netcdf(filename)
# Step 2: Convert the NetCDF file to a GeoTIFF
tif = 'wind_global.tif'
geemap.netcdf_to_tif(filename, tif, variables=['u_wind', 'v_wind'], shift_lon=True)
print(data)
# Step 3: Read the GeoTIFF file using rioxarray
data = xr.open_rasterio(tif)
# Extract latitude and longitude values from the data
latitude = data.y.values
longitude = data.x.values
# Convert wind data from U/V components to magnitude and direction
u_wind = data.values[0]
v_wind = data.values[1]
wind_magnitude = np.sqrt(u_wind*2 + v_wind*2)
wind_direction = np.degrees(np.arctan2(u_wind, v_wind))
# Create a meshgrid of latitude and longitude
lon_grid, lat_grid = np.meshgrid(longitude, latitude)
что выдает ошибку.
у модуля «xarray» нет атрибута «open_rasterio»
Я пытаюсь прочитать файл, этот файл netcdf, преобразовать его в файл tif и прочитать с помощью xarray [code]import xarray as xr import geemap import numpy as np from mayavi import mlab
# Step 1: Download the NetCDF file url = 'https://github.com/giswqs/leafmap/raw/master/examples/data/wind_global.nc' filename = 'wind_global.nc' data = geemap.read_netcdf(filename)
# Step 2: Convert the NetCDF file to a GeoTIFF tif = 'wind_global.tif' geemap.netcdf_to_tif(filename, tif, variables=['u_wind', 'v_wind'], shift_lon=True) print(data) # Step 3: Read the GeoTIFF file using rioxarray data = xr.open_rasterio(tif)
# Extract latitude and longitude values from the data latitude = data.y.values longitude = data.x.values
# Convert wind data from U/V components to magnitude and direction u_wind = data.values[0] v_wind = data.values[1] wind_magnitude = np.sqrt(u_wind*2 + v_wind*2) wind_direction = np.degrees(np.arctan2(u_wind, v_wind))
# Create a meshgrid of latitude and longitude lon_grid, lat_grid = np.meshgrid(longitude, latitude) [/code] что выдает ошибку. у модуля «xarray» нет атрибута «open_rasterio»