Python извлекает несколько lat / long из файлов NETCDF с помощью xarray

У меня есть файл ЧПУ (время, широта, долгота), Загрузите отсюда и Я пытаюсь извлечь временные ряды нескольких станций (точки широты и долготы Загрузить отсюда ). Итак, я попробовал прочитать координаты и извлечь ближайшие значения из файла ЧПУ следующим образом:

import pandas as pd
import xarray as xr
nc_file = r"C:\Users\lab\Desktop\harvey\example.nc"
NC = xr.open_dataset(nc_file)
csv = r"C:\Users\lab\Desktop\harvey\stations.csv"
df = pd.read_csv(csv,delimiter=',')
Newdf = pd.DataFrame([])
# grid point lists
lat = df["Lat"]
lon = df["Lon"]
point_list = zip(lat,lon)
for i, j in point_list:
    dsloc = NC.sel(lat=i,lon=j,method='nearest')
    DT=dsloc.to_dataframe()
    Newdf=Newdf.append(DT,sort=True)

Код работает нормально и возвращает это:

                        EVP     lat      lon
time                                        
2019-01-01 19:00:00  0.0546  40.063  -88.313
2019-01-01 23:00:00  0.0049  40.063  -88.313
2019-01-01 19:00:00  0.0052  41.938  -93.688
2019-01-01 23:00:00  0.0029  41.938  -93.688
2019-01-01 19:00:00  0.0101  52.938 -124.938
2019-01-01 23:00:00  0.0200  52.938 -124.938
2019-01-01 19:00:00  0.1644  39.063  -79.438
2019-01-01 23:00:00 -0.0027  39.063  -79.438

Однако мне нужно связать идентификатор станции (из моего исходного файла широты / долготы) для каждой из координат следующим образом:

  Station-ID       Lat        Lon            time     EVP     lat      lon
0        Bo1  40.00620  -88.29040  1/1/2019 19:00  0.0546  40.063  -88.313
1                                  1/1/2019 23:00  0.0049  40.063  -88.313
2        Br1  41.97490  -93.69060  1/1/2019 19:00  0.0052  41.938  -93.688
3                                  1/1/2019 23:00  0.0029  41.938  -93.688
4        Brw  71.32250 -156.60917  1/1/2019 19:00  0.0101  52.938 -124.938
5                                  1/1/2019 23:00  0.0200  52.938 -124.938
6        CaV  39.06333  -79.42083  1/1/2019 19:00  0.1644  39.063  -79.438
7                                  1/1/2019 23:00 -0.0027  39.063  -79.438

Есть мысли, как объединить мои фреймы данных, как в приведенном примере?


python pandas data-extraction netcdf python-xarray
person Seji    schedule 31.10.2019    source источник

Ответы (1)


arrow_upward
4
arrow_downward

Как насчет того, если вы включите имя станции в свою команду zip, а затем вставите идентификатор в строку фрейма данных pandas следующим образом (кстати, я не мог получить доступ к вашему CSV-файлу, поэтому я немного упростил пример с фиктивным списком) .

import pandas as pd
import xarray as xr
nc_file = "example.nc"
NC = xr.open_dataset(nc_file)

#dummy locations and station id as I can't access the CSV
lat=[40,42,41]
lon=[-100,-105,-99]
name=["a","b","c"]

Newdf = pd.DataFrame([])

for i,j,id in zip(lat,lon,name):
    dsloc = NC.sel(lat=i,lon=j,method='nearest')
    DT=dsloc.to_dataframe()

    # insert the name with your preferred column title:
    DT.insert(loc=0,column="station",value=id)
    Newdf=Newdf.append(DT,sort=True)

print(Newdf)

Это дает мне:

                        EVP     lat      lon station
time                                                
2019-01-01 19:00:00  0.0527  39.938  -99.938       a
2019-01-01 23:00:00  0.0232  39.938  -99.938       a
2019-01-01 19:00:00  0.0125  41.938 -104.938       b
2019-01-01 23:00:00  0.0055  41.938 -104.938       b
2019-01-01 19:00:00  0.0527  40.938  -98.938       c
2019-01-01 23:00:00  0.0184  40.938  -98.938       c
person Adrian Tompkins    schedule 31.10.2019
comment
Хорошая идея, которая работает! да, архивирование с идентификатором и вставка в фрейм данных работают отлично. Спасибо. - person Seji; 31.10.2019
comment
да, однажды я получу 15 репутаций :) сейчас я голосую, но это не отображается до моей 15 репутации! - person Seji; 31.10.2019
comment
Не знал этого! Кстати, добро пожаловать в stackexchange ???? - person Adrian Tompkins; 02.11.2019