import xarray as xr
import geopandas as gpd
import rioxarray
import pymannkendall as mk
import numpy as np
netcdf = r"C:\Users\paeli\Documents\Studio\netCDF\2025\water_balance.nc"
shape = r"C:\Users\paeli\Documents\Daten\Kazakhstan\shape\Kazakhstan.geojson"
water_balance = xr.open_dataset(netcdf)
region = gpd.read_file(shape)
# Reproject the shapefile to match the NetCDF file's CRS
water_balance = water_balance.rio.write_crs("EPSG:4326") # adjust EPSG as needed
region = region.to_crs(water_balance.rio.crs) # Reproject the shapefile
# Mask the NetCDF data using the shapefile
masked_water_balance = water_balance.rio.clip(region.geometry, region.crs, drop=True)
# Compute mean water balance over the masked region
water_balance_time = masked_water_balance["water_balance"].mean(dim=["lat", "lon"])
# Perform Mann-Kendall Trend Test
mk_result = mk.original_test(water_balance_time.values)
# Print results
print("Mann-Kendall Test Results for Selected Region:")
print(f"Trend: {mk_result.trend}")
print(f"P-value: {mk_result.p}")
print(f"Tau (correlation coefficient): {mk_result.Tau}")
< /code>
При этом я получаю эту ошибку: attributeError: у модуля «fiona» нет атрибута «Путь». Вы имели в виду: 'path'?
он запускается в строке 15, какая строка, которая пытается загрузить Geojson.>
Я пытаюсь запустить следующий код для вычисления теста на тренд Манн-Кендалла. < /p> [code]import xarray as xr import geopandas as gpd import rioxarray import pymannkendall as mk import numpy as np
# Reproject the shapefile to match the NetCDF file's CRS water_balance = water_balance.rio.write_crs("EPSG:4326") # adjust EPSG as needed region = region.to_crs(water_balance.rio.crs) # Reproject the shapefile
# Mask the NetCDF data using the shapefile masked_water_balance = water_balance.rio.clip(region.geometry, region.crs, drop=True)
# Compute mean water balance over the masked region water_balance_time = masked_water_balance["water_balance"].mean(dim=["lat", "lon"])
# Perform Mann-Kendall Trend Test mk_result = mk.original_test(water_balance_time.values)
# Print results print("Mann-Kendall Test Results for Selected Region:") print(f"Trend: {mk_result.trend}") print(f"P-value: {mk_result.p}") print(f"Tau (correlation coefficient): {mk_result.Tau}") < /code> При этом я получаю эту ошибку: attributeError: у модуля «fiona» нет атрибута «Путь». Вы имели в виду: 'path'? [/code] он запускается в строке 15, какая строка, которая пытается загрузить Geojson.>
У меня есть фрагмент кода, который работал нормально до прошлой недели, но теперь он выходит из строя из-за следующей ошибки:
AttributeError: у модуля «fiona» нет атрибута «path»
Я убедился, что все необходимые библиотеки установлены и...
Недавно я занимался некоторой статистикой на основе OpenStreetMaps. Когда я использовал эстакад для экспорта немецкого Autobahnen и Geopandas для расчета их длины, я обнаружил, что длины не соответствуют официальным цифрам. Перепускной запрос уже...
Недавно я занимался некоторой статистикой на основе OpenStreetMaps. Когда я использовал эстакад для экспорта немецкого Autobahnen и Geopandas для расчета их длины, я обнаружил, что длины не соответствуют официальным цифрам. Перепускной запрос уже...
Недавно я занимался некоторой статистикой на основе OpenStreetMaps. Когда я использовал эстакад для экспорта немецкого Autobahnen и Geopandas для расчета их длины, я обнаружил, что длины не соответствуют официальным цифрам. Перепускной запрос уже...
Недавно я занимался некоторой статистикой на основе OpenStreetMaps. Когда я использовал эстакад для экспорта немецкого Autobahnen и Geopandas для расчета их длины, я обнаружил, что длины не соответствуют официальным цифрам. Перепускной запрос уже...