Конвертируем netCDF в ASCII при помощи Python в Windows

Published by koldunovn on февраля 16, 2013 06:44 пп under Python, Обработка данных, Без рубрики

Задача: Помочь друзьям виндузятникам сконвертировать netCDF в ASCII, попутно установив на их компьютеры Python, в надежде, что они таки постепенно забудут про дельфи, фортран и прочие гадости. Заодно попробовать удобно ли в ipython notebook писать посты.

�?нструменты: cdo, Pyhton(x,y), ipython notebook

Важно: Этот пост изначально целиком написан в ipython notebook, что значит вы сможете при наличии установленного ipython notebook (как это сделать под Windows смотри ниже) исполнять и изменять код представленный в посте. Читать пост с хорошим форматированием, так как оно выглядит в ноутбуке, а также загрузить файл ноутбука можно здесь (загрузка файла в правом верхнем углу). Также доступна pdf версия этого поста автоматически сгенерированная nbconverter. А ещё версия этого поста в виде мимимишной презентации reveal. А теперь собственно пост:

Люди работающие под виндоуз любят ASCII, я знаю, я сам был такой. Они готовы переводить в ASCII всё на свете, включая данные моделей IPCC, которые занимая в бинарном формате сотни гигабайт, будучи переведены в ASCII превращаются в монстров, сжирающих всё доступное дисковое пространство в радиусе нескольких километров. Но таковы реалии виндузовой жизни, многие программы там, особенно занимающиеся отрисовкой и анализом данных, хотят, чтобы им скармливали текст и только текст.

У проблемы перевода из netCDF в ASCII существует множество решений. Можно сделать дамп всего файла, заголовка, или отдельной переменной при помощи программки ncdump.exe. Небольшую инструкцию как это сделать и где взять эту неуловимую программку можно почитать здесь (там пишут про HDF, но для netCDF эта инструкция также подходит). Правда этот дамп вам придётся потом ещё долго и печально разбирать, поскольку то, что вы увидите, будет довольно сильно отличаться от желаемой таблички время/широта/долгота/значение (я знаю мечтаете вы именно об этом :)).

Тем кто знаком с Matlab с некоторых пор вообще стало хорошо, поскольку там появилась поддержка netCDF файлов из коробки и о том, как их там открывать можно почитать тут.

Здесь я расскажу как сконвертировать netCDF в ASCII при помощи Python, при этом формат вывода вы сможете задавать какой пожелаете. Упражняться будем, как обычно, на файлах NCEP реанализа.

Установка Pyhton(x,y) и cdo

Учёному, который хочет писать скрипты на питон под Windows, нужно кроме самого питона установить ещё много полезных в хозяйстве модулей, таких как numpy, scipy, ipyhton. Дело это зачастую муторное, особенно для такого новичка в виндоуз администрировании как я. Однако нашлись добрые люди, собрали всё необходимое в один большой пакет и распространяют его под именем Pyhton(x,y). Преимущество в том, что не нужно заморачиваться с самостоятельной настройкой пакетов, а недостаток, на мой взгляд, в том, что получившийся экзешник весит больше 500 мегабайт, что, к сожалению, пока ещё не тот объём, который можно безболезненно загрузить в каждом Нигерийском селе.

Процесс установки после скачивания прост и стандартен. На каком-то этапе вам будет предложено выбрать нужные пакеты — ставьте всё, что помечено. При том, что установится у вас куча всего, нам этого мало 🙂 Чтобы рисовать красивые карты, нам необходим пакет Basemap, который устанавливается как дополнительный модуль к Pyhton(x,y), и лежит тут (первый в списке). Установив его вы получите на вашей Win машине полноценную среду для научного программирования под питон.

На самом деле я с трудом допускаю мысль, что вы хотите переводить в текстовый формат, скажем весь файл с реанализом ежедневных значений приземной температуры за год. Было бы неплохо сначала проделать над ним какие-то операции — осреднить по месяцам, или сезонам, построить стандартное отклонение наконец. Для этого человечество пока не придумало ничего лучше, чем climate data operators, или cdo. Небольшой рассказ о том, зачем они нужны и как ими пользоваться можно, как обычно найти на koldunov.net. С недавних пор они появились и под Win, так что скачивайте последнюю версию отсюда и распаковывайте её в ту папку, где вы будете работать с netCDF. В моём случае это папка C:/notebook/ (естественно вам её нужно будет предварительно создать).

Запуск IPython notebook

Ещё один маленький шажок и мы сможем перейти к работе над netCDF. Работать мы будем не из питоновской консоли, и не станем писать скрипты в отдельных файлах, а потом запускать их. Мы будем делать всё в IPython notebook, удивительном изобретении ребят из IPython, в котором я пишу этот пост, и который позволит вам одновременно запускать код и читать описание происходящего. Сперва загрузите ноутбук из которого сделан этот пост отсюда в вашу рабочую папку (C:/notebook/). Затем откройте консоль (если вы на XP то это Пуск >> Программы >> Стандартные >> Командная строка) и перейдите в свою рабочую папку (используйте команду cd). Если ваша папка, также как у меня, расположена по адресу C:/notebook/ то вам нужно будет ввести:

cd \
cd notebook

�? наконец запускаем IPyhton notebook, выполнив находясь в вашем рабочем каталоге команду

ipython notebook

В вашем броузере должна открыться вкладка со списком ноутбуков, в котором, скорее всего, будет находиться только ноутбук, который вы загрузили, и, надеюсь, уже сейчас читаете. Теперь, когда все приготовления закончены, давайте перейдём к работе с нашими любимыми netCDF файлами.

Подготовка к экспорту

Для начала загружаем файл NCEP реанализа в рабочую папку. Я буду использовать файл ежедневной приземной температуры воздуха за 2012 год. IPython позволяет делать системные вызовы (перед ними должен стоять восклицательный знак), так что мы можем посмотреть, что у нас в папке:

In [3]: !dir

C:\notebook

15.02.2013 21:08 <DIR> .
15.02.2013 21:08 <DIR> ..
15.02.2013 21:08 7я700я344 air.sig995.2012.nc
07.01.2013 11:26 10я192я583 cdo.exe
15.02.2013 21:07 29я927 netCDF_to_ASCII_with_Python.ipynb
06.08.2012 13:40 127я720 pthreadGC2.dll

Понятно, что у нас в папке находится экзешник cdo, дополнительная библиотека необходимая ему для работы, файл с ноутбуком и, собственно наш netCDF файл. Посмотрим, что нам расскажут об этом файле cdo

In [11]:!cdo pardes air.sig995.2012.nc

-1 air mean Daily Air temperature at sigma level 995 [degK]

cdo pardes: Processed 1 variable

In [12]:!cdo showyear air.sig995.2012.nc

2012
cdo showyear: Processed 1 variable over 366 timesteps

In [13]: !cdo showmon air.sig995.2012.nc

1 2 3 4 5 6 7 8 9 10 11 12

cdo showmon: Processed 1 variable over 366 timesteps

Мы узнали имя переменной, а также какие в файле содержатся годы и месяцы. Теперь давайте посчитаем среднюю температуру за сезон:

In [15]: !cdo seasmean air.sig995.2012.nc air.sig995.2012_sm.nc

cdo seasmean: Processed 3847392 values from 1 variable over 366 timesteps

In [15]: !cdo seasmean air.sig995.2012.nc air.sig995.2012_sm.nc

cdo seasmean: Processed 3847392 values from 1 variable over 366 timesteps

In [16]: !cdo sinfo air.sig995.2012_sm.nc

File format: netCDF
-1 : Institut Source Ttype Levels Num Gridsize Num Dtype : Parameter ID
1 : unknown unknown instant 1 1 10512 1 I16 : -1
Grid coordinates :
1 : lonlat > size : dim = 10512 nlon = 144 nlat = 73
lon : first = 0 last = 357.5 inc = 2.5 degrees_east circular
lat : first = 90 last = -90 inc = -2.5 degrees_north
Vertical coordinates :
1 : surface : 0
Time coordinate : 5 steps
RefTime = 0001-01-01 00:00:00 Units = hours Calendar = STANDARD
YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss
2012-02-29 00:00:00 2012-05-31 00:00:00 2012-08-31 00:00:00 2012-11-30 00:00:00
2012-12-31 00:00:00

cdo sinfo: Processed 1 variable over 5 timesteps

Как видите получилось у нас не четыре, а пять полей, поскольку зим за 2012 год было аж две штуки.

Допустим, весь мир нам не нужен, а мы хотим вырезать только Европу:

In [20]: !cdo sellonlatbox,0,90,30,70 air.sig995.2012_sm.nc air.sig995.2012_sm_europe.nc

cdo sellonlatbox: Processed 52560 values from 1 variable over 5 timesteps

Этот файл и будем конвертировать в ASCII.

Открываем netCDF файл в Pyhton

Чтобы иметь возможность открывать netCDF, нам нужно импортировать пару модулей:

In [19]: from scipy.io import netcdf
import numpy

Теперь открываем файл:

In [21]: f = netcdf.netcdf_file('air.sig995.2012_sm_europe.nc')

Смотрим, что там у нас внутри:

In [22]: f.variables

Out[22]: {'air': <scipy.io.netcdf.netcdf_variable at 0x2c43bb0>,
'lat': <scipy.io.netcdf.netcdf_variable at 0x2c436b0>,
'lon': <scipy.io.netcdf.netcdf_variable at 0x2c43690>,
'time': <scipy.io.netcdf.netcdf_variable at 0x2c43870>}

Нам понадобятся все эти переменные, так что импортируем их в переменны питона:

In [23]: air = f.variables['air']
lat = f.variables['lat']
lon = f.variables['lon']
time = f.variables['time']

Если вы теперь поставите точку после имени переменной (любой, не только air) и нажмёте TAB, то увидите список различных атрибутов этой переменной, предлагаю вам их поизучать:

In [32]: air.

Два из них особенно важны, это add_offset и scale_factor. Откуда они взялись и зачем нужны можно почитать здесь. Если читать лень, то просто знайте: чтобы получить нормальные данные в Кельвинах нужно умножить air на scale_factor и прибавить add_offset. Так надо.

В данный момент в наших переменных содержатся не только данные, но и дополнительная информация, вроде атрибутов. Для того, чтобы работать только с данным, мы можем их скопировать в отдельную переменную, попутно сконвертировав их в Кельвины:

In [35]: airK = (air.data*air.scale_factor)+air.add_offset

Ещё вариант

In [39]: airK = (air[:]*air.scale_factor)+air.add_offset

В принципе, если бы мы не хотели получать информацию об атрибутах переменных, то могли бы импортировать только данные следующим образом:

air = f.variables['air'][:]
lat = f.variables['lat'][:]
lon = f.variables['lon'][:]
time = f.variables['time'][:]

но в этом случае информацию, например, о scale_factor и add_offset пришлось бы добывать из других источников.

Переведём данные из Кельвинов в Цельсии

In [41]: airC = airK-273.15

Вы уже, наверное, поняли, что никаких циклов с перебиранием каждого значения и вычитанием магической цифры тут не нужно, всё как в матлабе, операция вычитания производится для всех элементов массива.

Теперь взгляним на наши координаты:

In [42]: lat.shape

Out[42]: (17,)

In [43]: lon.shape

Out[43]: (37,)

Они одномерные. Ничего плохого в этом нет, но мы хотим нарисовать наши данные, прежде чем экспортировать, а для этого, нужно перевсти широты и долготы в двумерный формат (то есть двумерному массиву, в котором содержится поле с данными, будут соответствовать два двумерных массива в широтами и долготами).

In [45]: lons, lats = numpy.meshgrid(lon[:],lat[:])

Поглядим, в каком формате у нас время

In [70]: time.units

Out[70]: 'hours since 1-01-01 00:00:00'

Время у нас в часах от начала нашей эры. Выглядят они вот так:

In [71]: time[:][0]

Out[71]: 17629512.0

Чтобы перевести эти часы в удобоваримый формат, нам нужно воспользоваться модулем datetime

In [72]: import datetime

Создадим объект специального типа datetime, при помощи которого удобно управляться со временем, в котором будет дата с которой считаются часы в нашем формате времени:

In [154]: first_date = datetime.datetime(1,1,1,0,0,0)

In [155]: type(first_date)

Out[155]: datetime.datetime

Теперь представим нашу первую дату в виде временной разницы

In [168]: difference1 = datetime.timedelta(hours=(time[:][0]-48))

По причине, в которой я ещё не разобрался, эта функция ошибается на два дня, так что мы их (то есть 48 часов) вычитаем.

Прибавим к начальной дате разницу, и получим datetime переменную, в которой будет содержаться наша дата.

In [169]: present = first_date+difference1

У переменной present много интересных методов, которые позволяют узнать не только дату, но и, например, день недели. Но нас будет интересовать полная запись даты

In [170]: present.ctime()

Out[170]: 'Wed Feb 29 00:00:00 2012'

Как видите в качестве даты даётся последний день сезона. �?нформацию о времени мы используем позже при отрисовке данных и экспорте в netCDF.

Теперь немножно магии, чтобы отображать картинки внутри ноутбука:

In [46]: %pylab inline

Welcome to pylab, a matplotlib-based Python environment [backend: module://IPython.zmq.pylab.backend_inline].
For more information, type 'help(pylab)'.

�?мпортируем модуль для отрисовки карт:

In [47]: from mpl_toolkits.basemap import Basemap

Генерируем карту

In [65]: m = Basemap(llcrnrlon=-10,llcrnrlat=25,urcrnrlon=95,urcrnrlat=75,projection='mill', resolution='l')

Переводим широты и долготы в координаты карты

In [66]: x, y = m(lons, lats)

�? отрисовываем наши данные:

# Создаём рисунок и задаём его размеры
fig = plt.figure(figsize=(10,10))
# Выбираем сезон
season=2
# Рисуем береговую линию
m.drawcoastlines(linewidth=1)
# Рисуем границы стран
m.drawcountries(linewidth=0.5)
# Рисуем границу карты
m.drawmapboundary
# Создаём карту с залитыми изолинииями, 20 интервалов
cs = m.contourf(x,y,airC[season,:,:],20)
# Создаём легенду
cbar = m.colorbar(cs,location='bottom',pad="5%")
# Подписываем легенду
cbar.set_label(r'$^{\circ}\mathrm{C}$')
# Переводим время в человеческое
difference1 = datetime.timedelta(hours=(time[:][season]-48))
present = first_date+difference1
# Отрисовываем заголовок
plt.title(present.ctime())

С определением Европы мы явно промахнулись 🙂 Чтобы исправить эту ситуацию можно вернуться к шагу, в котором мы вырезаем регион и подправить координаты, которые передаются в cdo.

В принципе, загрузив данные в питон (назовём это так) вы можете делать с ними всё, что угодно. В этом смысле процесс не сильно отличается от работы, например, в Matlab. Но мы с вами собирались сохранить данные в ASCII, давайте этим и займёмся.

Экспорт в ASCII

Алгоритм сохранения прост: открываем файл, последовательно записываем в него строчки в нужном нам формате, закрываем файл. Теперь посмотрим как это выглядит на практике. Сначала пример, в котором каждое поле сохраняется в отдельный файл:

#Цикл по времени
for one_field in range(airC.shape[0]):

#Открываем файл, дописываем к имени номер поля и нули,
#чтобы при большом количестве файлов они нормально сортировались
ofile = open("netcdf_to_ASCII_field_"+str(one_field).zfill(5)+".txt","w")

#Переводим время в человесеский формат
difference1 = datetime.timedelta(hours=(time[:][one_field]-48))
present = first_date+difference1

#Печатаем, чтобы видеть прогресс
print(present.ctime())

#Цикл по строкам
for nn in range(airC.shape[1]):

#Цикл по столбцам
for kk in range(airC.shape[2]):

#Запись строки в файл. Сначала перечисляем форматы,
#в которых будт выводиться переменные, а потом перечисляем
#сами переменные
ofile.write('%5.2f %5.2f %5.2f %5.0f %5.0f %5.0f\n' % \
(lons[nn,kk], lats[nn,kk], airC[one_field,nn,kk], present.day, present.month, present.year))

ofile.close()

Wed Feb 29 00:00:00 2012
Thu May 31 00:00:00 2012
Fri Aug 31 00:00:00 2012
Fri Nov 30 00:00:00 2012
Mon Dec 31 00:00:00 2012

Самое «сложное», или скорее неочевидное тут это форматы вывода данных. Про них можно почитать тут.

Вариант с выводом в один файл с дата стампами между полями:

#Открывать будем только один файл
ofile = open("netcdf_to_ASCII_all_fields.txt","w")

#Цикл по времени
for one_field in range(airC.shape[0]):

#Переводим время в человесеский формат
difference1 = datetime.timedelta(hours=(time[:][one_field]-48))
present = first_date+difference1

#Печатаем, чтобы видеть прогресс
print(present.ctime())

#Записываем время в начале поля,
#тут мы решили ещё и часы с минутами приписать
ofile.write('%5.0f %5.0f %5.0f %5.0f %5.0f\n' % (present.day, present.month, present.year, present.hour, present.minute ))

#Цикл по строкам
for nn in range(airC.shape[1]):

#Цикл по столбцам
for kk in range(airC.shape[2]):

#Запись строки в файл. Сначала перечисляем форматы,
#в которых будт выводиться переменные, а потом перечисляем
#сами переменные
ofile.write('%5.2f %5.2f %5.2f\n' % \
(lons[nn,kk], lats[nn,kk], airC[one_field,nn,kk]))

ofile.close()

Wed Feb 29 00:00:00 2012
Thu May 31 00:00:00 2012
Fri Aug 31 00:00:00 2012
Fri Nov 30 00:00:00 2012
Mon Dec 31 00:00:00 2012

Как видите слагаемые буквально чуть-чуть поменялись местами, а результат изменился 🙂

Пользуясь двумя этими шаблонами, как мне кажется, можно сохранить большинство данных, которые вам могут понадобиться, в текстовых форматах, которые нужны именно вам, а не каким-то совершенно чужим людям 🙂 Однако научившись «загружать» и отображать данные в питоне, я надеюсь вам захочется провести в нём больше времени, и может в один прекрасный день вы поймёте, что ASCII файлы вам больше не нужны 🙂

14 комментариев

Oleg on 2 апреля, 2013

Одно не понятно, зачем отдельно качать ipython, когда он уже есть в полном пакете python(x,y)?
koldunovn on 3 апреля, 2013

@Oleg

А где написано, что его нужно отдельно качать?
Mark on 11 апреля, 2013

Благодаря этому посту попробовал Python(x,y). Всё поставилось, но, iPython запускаться отказывался. Дело оказалось в непонятно откуда взявшихся mime-типах в реестре, в названиях которых была кириллица. Поправил это — и всё запустилось. Свой user experience описал здесь — http://thebaga.blogspot.com/2013/04/ipython.html Может, кто ещё с таким столкнётся.
natavrame on 13 апреля, 2013

Не знаю, где правильнее задать вопрос, но…А как можно прочитать файл, в котором есть пропущенные значения, на месте которых стоят просто пробелы? При этом в качестве delimeter выступает тоже пробел. �? еще одна загвоздка состоит в том, что иногда бывает несколько пропусков подряд, и тогда это может восприниматься за один..
Очень нужна помощь! В python еще совсем новичок…
koldunovn on 13 апреля, 2013

2natavrame
Не могу сказать точно как вам помочь не видя самого файла ) Единственно, если количество символов в строке всегда одинаково, то можно открыв файл и считав из него строки разобрать каждую строку посимвольно. Например:

ifile = open('input_file', 'r')
lines = ifile.readlines()

for line in lines:
param1 = line[0:5]
param2 = line[5:10]
natavrame on 14 апреля, 2013

Спасибо большое! Да, если задавать каждый параметр отдельно, то вроде нормально считывает. А дальше тогда надо поставить if-условие, что-то типа
if param==’\t\t\t\t\t\t’:
param=9999

Сам файл выглядит примерно так: (не знаю, как его можно загрузить сюда )
20674 1967 22.0 35.0 28.0 39.0 22.0 23.0 59.0 46.0 49.0 45.0 36.0 29.0
20674 1968 17.0 21.0 28.0 17.0 14.0 25.0 47.0 68.0 22.0 22.0
20674 1969 20.0 16.0 15.0 8.0 9.0 17.0 60.0 33.0 27.0 11.0 27.0 30.0
20674 1970 11.0 17.0 24.0 10.0 19.0 14.0 56.0 10.0 24.0 25.0 36.0 43.0
20674 1971 91.0 39.0 17.0 4.0 12.0 13.0 31.0 67.0 78.0 16.0 31.0 26.0
20674 1972 32.0 47.0 11.0 28.0 19.0 23.0 51.0 43.0 34.0 35.0 4.0 69.0
20674 1973 25.0 36.0 25.0 9.0 38.0 50.0 21.0 75.0 43.0 36.0 16.0 15.0
20674 1974 26.0 40.0 30.0 21.0 44.0 31.0 44.0 22.0 61.0 32.0 27.0
20674 1975 49.0 16.0 15.0 29.0 12.0 28.0 35.0 25.0 35.0 43.0 24.0 54.0
20674 1976 24.0 34.0 17.0 36.0 19.0 53.0 40.0 63.0 47.0 19.0 21.0 14.0
20674 1977
20674 1978 20.0 20.0 16.0 6.0 14.0 18.0 48.0 41.0 12.0 27.0 25.0 16.0
20674 1979 22.0 10.0 9.0 8.0 48.0 38.0 34.0 24.0 13.0 19.0 26.0
20674 1980 21.0 35.0 9.0 14.0 23.0 41.0 66.0 67.0 55.0 44.0 8.0 32.0
20674 1981 105.0 14.0 14.0 23.0 40.0 14.0 47.0 27.0 26.0 28.0
20674 1982 8.0 13.0 13.0 24.0 12.0 58.0 13.0 71.0 40.0 23.0 5.0 26.0

Это данные станционных наблюдений ВН�?�?ГМ�? по среднемесячным осадкам.
koldunovn on 15 апреля, 2013

Да, по идее должно работать ) Если что обращайтесь )
natavrame on 15 апреля, 2013

Благодарю за помощь!!
Почему-то искать ’\t\t\t\t\t\t’ ему не понравилось и он решил оставлять все как есть. Поэтому мне помогли сделать задачу через фотран. Но, надеюсь, через некоторое время и у меня получится сделать это все-таки на python. )
koldunovn on 16 апреля, 2013

2natavrame
Так мне стало грустно от того, что вам пришлось воспользоваться фортраном, что пришлось сделать ноутбук с объяснением как можно открыть ваши данные и начать их анализировать 🙂 Посмотреть его можно тут: http://nbviewer.ipython.org/5391791
natavrame on 17 апреля, 2013

Огромное Вам спасибо!!! Это просто волшебно! Чем больше читаю примеров использования, тем больше и больше нравится python. Особенно приятно удивило «Time series analysis with pandas», потому что то, что описывается в несколько строк в других языках, можно сделать всего лишь одной строчкой. ))
А можно задать еще один вопрос…может быть Вам приходилось сталкиваться с этим, можно посчитать площадь внутри замкнутых изолиний? �?ли это нудно описывать с самого нуля вручную?
natavrame on 17 апреля, 2013

*нужно
koldunovn on 17 апреля, 2013

2natavrame
Я возможно не очень понял, что вы имеете в виду, но вот первое, что приходит в голову:
— если вы работаете с гридированными данными, то можно просто выбрать все гридбоксы с данными больше, меньше, или находящиеся между какими то значениями, и затем посчитать общую площадь, сложив площади всех удовлетворяющих условиям точек.

— если вы работаете с негридированными данными, то для того, чтобы получить изолинии вам их нужно будет интерполировать на какой-то грид, после чего см. пункт первый 🙂

Сделал маленький ноутбук с возможной реализацией этого дела http://nbviewer.ipython.org/5404107
olgakr on 13 января, 2015

нужно для дипломной работы…
подскажите пожалуйста, у меня данные координатной сетки не в долготе и широте, а в xx — абцисса и yy — ордината:
-25 xx abscissas
-26 yy ordinates

xx.data

Out[51]:

array([ 60000., 60000., 60000., 60000., 60000., 60000.,
60000., 60000., 60000., 60000., 60000., 60000.,
120000., 120000., 120000., 120000., 120000., 120000.,
120000., 120000., 120000., 120000., 120000., 120000.,
180000., 180000., 180000., 180000., 180000., 180000.,
180000., 180000., 180000., 180000., 180000., 180000.,
240000., 240000., 240000., 240000., 240000., 240000.,
240000., 240000., 240000., 240000., 240000., 240000.,
300000., 300000., 300000., 300000., 300000., 300000.,
300000., 300000., 300000., 300000., 300000., 300000.,
360000., 360000., 360000., 360000., 360000., 360000.,
360000., 360000., 360000., 360000., 360000., 360000.,
420000., 420000., 420000., 420000., 420000., 420000.,
420000., 420000., 420000., 420000., 420000., 420000.,
480000., 480000., 480000., 480000., 480000., 480000.,
480000., 480000., 480000., 480000., 480000., 480000.])
yy.data — выглядят аналогично

lon_0=0,lat_0=43

Как мне перевести абциссы с ординатами в широты с долготами или может строить выходя от имеющегося — то есть, может есть вариант как нанести данные на карту с такими вот хх и yy?
Alexandr on 29 июня, 2015

Подскажите, пожалуйста, как решить обратную задачу, т.е. из текстового файла соорудить netcdf на Рython.
casino zonder storting on 22 февраля, 2020

gratis gokken kasten spelen , online casino met gratis welkomstbonus zonder storting — wat is een goed online casino : wat kost holland casino .
holland casino vestigingen on 22 февраля, 2020

fruitautomaten zeus , speelautomaten bestellen — online poker met echt geld : gratis gokkasten red green peppers spelen .
weekendje weg met casino on 22 февраля, 2020

online casino telefoon betalen , gokkast club casino — van dessel speelautomaten : casino spelletjes roulette .
jackpot holland casino nijmegen on 22 февраля, 2020

fruitkasten online gratis , casino regio utrecht — blackjack holland casino : holland casino eurojackpot .
van der valk en casino on 22 февраля, 2020

gokkast spelen club 2000 , geld zonder storting — bethard casino betrouwbaar : fruitmachines gratis online spelen .
casino admiral nederland on 22 февраля, 2020

gratis gokkast spelletjes , casino gratis startgeld — holland casino rotterdam dresscode : live casino gokken .
altus speelautomaten on 22 февраля, 2020

bob casino ervaringen , holland casino met de kaarten op tafel gemist — casino strategie blackjack : wat is een betrouwbaar online casino .
online casino betalen met ideal on 22 февраля, 2020

holland casino kerst breda , legitimeren met rijbewijs bij holland casino — karamba casino bonuscode : wat is de entree van holland casino .
venlo casino holland poker on 22 февраля, 2020

hoe betrouwbaar is dunder casino , uitbetaling pronto casino — video speelautomaten : gratis gokken voor geld .
speelautomaat 11 letters on 22 февраля, 2020

betrouwbare casino online , japanse gokmachine — tijdelijke locatie holland casino groningen : van ham speelautomaten .
gokje wagen casino on 22 февраля, 2020

online gokken spelletjes , leovegas casino bonus — roulette holland casino online : spelen op oude gokkasten .
wat is een goed online casino on 22 февраля, 2020

gokautomaat kopen , assistent speelautomaten — gratis casino spelen voor echt geld : photoplay speelautomaat .
gokkasten spelen ideal on 22 февраля, 2020

gratis geld verdienen online casino , online gokkasten iphone — robert de vries speelautomaten : gokken op fruitautomaten .
betalen via sms casino on 22 февраля, 2020

online casino belasting betalen , holland casino rotterdam eerste kerstdag — assistent speelautomaten : belasting betalen over casino geld .
transmissible gastroenteritis coronavirus on 28 февраля, 2020

how to treat coronavirus in cats , what disease does avian coronavirus cause — coronavirus sorse : coronavirus cats fip .
interesting facts about sars coronavirus on 28 февраля, 2020

does coronavirus cause liver failure , coronavirus in cats treatment — does coronavirus require isolation : hantavirus coronavirus .
feline coronavirus mutation on 28 февраля, 2020

unexpected receptor-functional mimicry elucidates activation of coronavirus fusion. , cat disease coronavirus — rotavirus coronavirus piglet : biocan c coronavirus .
sars coronavirus mortality rate on 28 февраля, 2020

who named coronavirus , coronavirus cirrhosis cardiomyopathy — fiv a coronavirus? : why use rt pcr on coronavirus .
is feline coronavirus fatal on 28 февраля, 2020

coronavirus protection , infectious disease coronavirus sputum — feline enteric coronavirus prevention : cures for the coronavirus in cats .
why dies the coronavirus mutate into fip on 28 февраля, 2020

feline coronavirus fip , sketchy coronavirus — is coronavirus : cat coronavirus titer test .
respiratory syndrome coronavirus in the middle east article on 28 февраля, 2020

coronavirus needed in florida , coronavirus and hallucinations — feline coronavirus life expectancy : coronavirus nested mrna .
isolation precautions for coronavirus on 28 февраля, 2020

coronavirus from china , coronavirus 229e isolation — coronavirus synthesis : coronavirus house cat .
2012 middle east respiratory syndrome coronavirus outbreak on 28 февраля, 2020

spectra canine distemper adenovirus type 2 coronavirus parainfluenza parvovirus vaccine , vaccine for coronavirus — htcc coronavirus inhibitor : sars coronavirus replication overview .
feline coronavirus replication on 28 февраля, 2020

can coronavirus cause kennel cough , quotes about sars coronavirus — feline coronavirus lab : pneumonia common cold coronavirus .
feline coronavirus titer on 28 февраля, 2020

should i be worried about coronavirus , coronavirus mask — coronavirus infection in adults : what are the symptoms of coronavirus in dogs? .
feline coronavirus antigen on 28 февраля, 2020

coronavirus cats , coronavirus in alpacas — what is the treatments for the coronavirus in both humans and animals : coronavirus kod pasa .
feline coronavirus kitten on 29 февраля, 2020

coronavirus upset stomach , distemper coronavirus — symptoms of coronavirus sars : coronavirus upset stomach .
signs of coronavirus in dogs on 29 февраля, 2020

show image of coronavirus , bordetella shot coronavirus — can you get coronavirus from a cat : receptor-induced conformational changes of murine coronavirus spike protein .
coronavirus denison lamp1 on 29 февраля, 2020

coronavirus one health , coronavirus vaccine for puppies for sale — unusual type of coronavirus : symptoms of rhinovirus and coronavirus together .
care for patients with middle east respiratory syndrome coronavirus on 29 февраля, 2020

coronavirus is double stranded rna , cell entry of coronavirus mediated by the conformational changes of its spike protein — virulence factors of sars coronavirus : coronavirus leader sequence .
global patterns in coronavirus diversity on 29 февраля, 2020

triggers cell entry pathway coronavirus , pronounce coronavirus — coronavirus medical definition : sars coronavirus icosahedron .
what is coronavirus hku1 on 29 февраля, 2020

canine coronavirus detection , coronavirus sysmptoms — when is coronavirus contagious : what is coronavirus hku1 .
coronavirus in cats how long on 29 февраля, 2020

coronavirus enveloped , canine respieratory coronavirus — is coronavirus contagious contact precautions : swine enteric coronavirus disease symptoms .
mers coronavirus pathogenesis on 29 февраля, 2020

info about coronavirus , coronavirus test for cows — coronavirus combining : cat coronavirus .
coronavirus flu on 29 февраля, 2020

canine coronavirus (ccov) vaccine infectious and non-infecious vaccine , coronavirus structure and composition — coronavirus spread by country : can coronavirus cause pale stools .
sars coronavirus emergence on 29 февраля, 2020

who coronavirus infection , can coronavirus cause pneumonia — mark denison coronavirus : middle eastern respiratory syndrome coronavirus map .
what can i do for my cat with coronavirus on 29 февраля, 2020

should i adopt a cat with coronavirus , feline coronavirus katze — how to diagnose coronavirus : coronavirus treatment in puppies .
recent coronavirus outbreak on 29 февраля, 2020

coronavirus respiratory syndrome , coronavirus vaccine a59 — what does swine coronavirus cause : difference between rsv and coronavirus .
coronavirus structure ppt on 29 февраля, 2020

what disease is caused by coronavirus in dogs , coronavirus tiger test cat — coronavirus with flu : coronavirus and transmissible gastroenterology .
coronavirus quasispecies on 29 февраля, 2020

canine coronavirus client education , bovine rota coronavirus — treating feline enteric coronavirus : coronavirus equine vaccine .
2 cats in house, one with coronavirus on 29 февраля, 2020

coronavirus hku 1 in 18 month old , coronavirus antibody — variant of the coronavirus that can be lead to respiratory failure and death : coronavirus transmission in animals .
severe acuity respiratory syndrome coronavirus on 29 февраля, 2020

how many cats carry coronavirus , coronavirus nl63 and 229e — coronavirus grey squirrels : feline coronavirus how long does it last .
coronavirus helical or icosahedral on 29 февраля, 2020

coronavirus membrane vs cell membrane , cace 2 inhibitor drug coronavirus — spike protein coronavirus : will you die if you get coronavirus .
feline coronavirus handout on 29 февраля, 2020

coronavirus and dates , coronavirus naturalpath treatment — spike protein coronavirus effects on antibodies : feline coronavirus stress .
virus known for its small size coronavirus on 29 февраля, 2020

coronavirus baltimore , polytropic coronavirus — coronavirus baltimore : coronavirus are negative strand rna animal .
coronavirus calf scours in dairy calves treatment on 29 февраля, 2020

coronavirus droplet or airborne , enteric coronavirus pcr — what is the life cycle of mers coronavirus : coronavirus cat develop diarrhea .
where coronavirus is a concern on 29 февраля, 2020

feline intestinal coronavirus , suppression of feline coronavirus replication in vitro by cyclosporin a — feline enteric coronavirus strain : coronavirus, leptospirosis, , lyme disease cost .
how is coronavirus treated on 29 февраля, 2020

coronavirus denison lamp1 , how do pediatricians diagnose coronavirus — mortality of tge coronavirus piglets : what are the projections on a coronavirus .
coronavirus baby on 29 февраля, 2020

coronavirus signs and symptoms in felines , how often do dogs need canine coronavirus vaccine — coronavirus synthesis : what are the symptoms of the coronavirus? .
coronavirus transmission from cats to humans on 29 февраля, 2020

coronavirus 229e treatment , sars coronavirus kilodaltons — coronavirus cat and vestibular syndrome : human coronavirus 0c43 .
coronavirus en francais on 29 февраля, 2020

can the coronavirus last for a.month , diagnosis fip coronavirus titer — feline coronavirus ihc : coronavirus in california .
coronavirus and neurovirus on 29 февраля, 2020

k99+ e. coli and bovine coronavirus , feline coronavirus incidence — coronavirus pink eye : coronavirus frameshift .
DennJeoge on 1 марта, 2020

Twenty of zoonosis you bought a does viagra super active work The originator to yield and as it confuses
http://viagrasupera.com/
Lesdurl on 8 марта, 2020

Macrobid 100mg Discount Worldwide Pharmacy Next Day California [url=http://apcialisle.com/#]Buy Cialis[/url] Decadron Canada Post Buy Cialis Buying Tretinoin In Singapore
StepOxice on 10 марта, 2020

Levitra 20 Mg Indicazioni [url=http://apcialisle.com/#]Buy Cialis[/url] Clomid Et Test De Grossesse Cialis Keflex Dose For Dog
viagra without prescription on 11 марта, 2020

viagra hangover [url=http://viagrawithoutdoctorpres.net/#]viagra for sale for men[/url] cheap viagra viagra viagra generic — viagra for women results previous next can viagra and viagra treat other health issues?
Lesdurl on 16 марта, 2020

Propecia Cost Kmart [url=http://apcialisle.com/#]cialis prices[/url] Erection Remedies Cialis Misoprostol Retained Placenta
viagra without a doctor prescription on 17 марта, 2020

viagra information side effects visitor messages [url=http://genviagrafx.com/#]real viagra
for sale[/url] viagra generic canada powered by ip.board viagra for sale — alcohol and ed viagra 10 mg 4 tablet register
Collin Lincourt on 17 марта, 2020

Thank you http://koldunov.net for informing
us!
The world is having a hard time with this pandemic.

However, there is hope, a guide that will be useful: http://bit.ly/coronavirus-survival-guide
Good health to all!
Janmealia on 18 марта, 2020

Propecia Side Effects Sperm Prescription [url=https://apcialisle.com/#]cialis 5 mg[/url] Viagra Cialis On Line Sales cialis canada Cialis Em Farmacia
StepOxice on 18 марта, 2020

Cheapest Canadian Generic Cialis [url=https://apcialisle.com/#]Cialis[/url] cheap accutane online Cialis Amoxicillin 400mg Suspension
Janmealia on 22 марта, 2020

Cephalexin Perscription Manufacturer [url=https://apcialisle.com/#]Cialis[/url] Cialis Tadalafil Prezzo cialis 40 mg Suggest Some Canadian Pharmacies
Lesdurl on 22 марта, 2020

Online Ed Tratments [url=https://apcialisle.com/#]generic cialis[/url] Cialis Generique Maroc cialis from canada Prix Du Cialis Avec Ordonnance
buy sildenafil citrate from india online best price on 26 марта, 2020

viagra buy return to board index viagra online scams total members
buy viagra online safely viagra without prescription — how
does viagra work? viagra online canadian pharmacy
how to take levitra for best results on 26 марта, 2020

levitra v levitra [url=http://levitramdx.com/#]generic levitra[/url] levitra buy posts per day cialis vs viagra vs levitra
— levitra pack discount levitra 10mg tablets announcements
vurtilopmer on 29 марта, 2020

I really appreciate this post. I have been looking all over for this! Thank goodness I found it on Bing. You’ve made my day! Thanks again
ZaraExert on 29 марта, 2020

3 coronavirus in us
[url=http://chloroquine-treatmentforcoronavirus.com/#]chloroquine phosphate 250mg
[/url] coronavirus outbreak china
chloroquine covid-19
— n coronavirus transmission
coronavirus prevention
viagra without doctors prescription on 30 марта, 2020

generic viagra tadafil [url=http://viagenusa.com/#]generic viagra for sale
in rsa[/url] ed in men viagra for sale cheap — viagra 5mg
generic register viagra viagra levitra phorum view topic
ZaraExert on 30 марта, 2020

coronavirus michigan
[url=http://chloroquine-treatmentforcoronavirus.com/#]chloroquine phosphate
[/url] coronavirus quadratic
chloroquine for covid-19
— 6 coronavirus cases
coronavirus 4s ranch
ZaraExert on 30 марта, 2020

coronavirus 3rd
[url=http://chloroquine-treatmentforcoronavirus.com/#]chloroquine for coronavirus
[/url] n coronavirus symptoms
chloroquine phosphate 250 mg
— parvovirus moquillo o coronavirus
coronavirus pathophysiology
ZaraExert on 31 марта, 2020

coronavirus philadelphia
[url=http://chloroquine-treatmentforcoronavirus.com/#]anti coronavirus drugs
[/url] coronavirus queens ny
chloroquine covid-19
— coronavirus 0c43
2 coronavirus cases in chicago
ZaraExert on 31 марта, 2020

coronavirus 1700
[url=http://chloroquine-treatmentforcoronavirus.com/#]chloroquine drug
[/url] 5 coronavirus in usa
chloroquine for covid-19
— coronavirus 63
coronavirus tracker
ZaraExert on 31 марта, 2020

coronavirus prevention tips
[url=http://chloroquine-treatmentforcoronavirus.com/#]order chloroquine
[/url] coronavirus bat
chloroquine phosphate
— coronavirus meme
coronavirus united states
kelvoe on 31 марта, 2020

ERE4w0 bbtiuthxbuqt, [url=http://xgvawlqljuoc.com/]xgvawlqljuoc[/url], [link=http://bhqzqwvukbxf.com/]bhqzqwvukbxf[/link], http://jcbpemdjhetv.com/
ZaraExert on 1 апреля, 2020

is sars a coronavirus
[url=http://chloroquine-treatmentforcoronavirus.com/#]chloroquine covid-19
[/url] coronavirus 8200
chloroquine sars
— i coronavirus symptoms
coronavirus racism
ZartFum on 2 апреля, 2020

coronavirus 2012
[url=http://chloroquine-phosphate-250mg.com/#]chloroquine coronavirus
[/url] is ebola a coronavirus
chloroquine phosphate 250 mg
— coronavirus 56 million
coronavirus yakima wa
LouisEphet on 2 апреля, 2020

[url=http://mysite.ru]http://mysite.ru[/url]
ZartFum on 2 апреля, 2020

is sars a coronavirus
[url=http://chloroquine-phosphate-250mg.com/#]chloroquine phosphate online
[/url] coronavirus zombie outbreak
chloroquine
— coronavirus 2 week incubation
coronavirus real death toll
Seo Proxies on 3 апреля, 2020

Magnificent website. Plenty of useful info here. I am sending it to some friends ans also sharing in delicious. And certainly, thanks for your sweat!
ZartFum on 3 апреля, 2020

coronavirus where
[url=http://chloroquine-phosphate-250mg.com/#]aralen chloroquine
[/url] coronavirus 2020 name
order chloroquine
— a novel coronavirus
3 coronavirus in virginia
ZartFum on 3 апреля, 2020

coronavirus mode of transmission
[url=http://chloroquine-phosphate-250mg.com/#]chloroquine phosphate 250mg
[/url] coronavirus 0c43
chloroquine sars
— coronavirus queens
coronavirus yale model un
ZartFum on 3 апреля, 2020

coronavirus que es
[url=http://chloroquine-phosphate-250mg.com/#]buy chloroquine online
[/url] 2 coronavirus in chicago
chloroquine covid-19
— coronavirus real death toll
coronavirus us
ZartFum on 4 апреля, 2020

coronavirus tracker
[url=http://chloroquine-phosphate-250mg.com/#]buy chloroquine
[/url] coronavirus zone
order chloroquine
— coronavirus yale new haven
a coronavirus outbreak

Subscribe to the comments through RSS Feed

Записки океанолога — обработка и визуализация данных