nc_file = r'C:\Users\hp\Desktop\project execution\geopandas\New thematic maps\soil_maps(20_50)cm\soil nitrogen total(20_50)cm\masked_nitrogen_total_2.nc' ds = xr.open ...