Python netcdf grid_longitude - PullRequest
Новая
       8

Python netcdf grid_longitude

0 голосов
/ 13 марта 2019

Обычно, когда я работаю с файлами, координаты указаны в lon и lat, но в этом случае они находятся в 'grid_longitude' 

import xarray as xr
dataset = xr.open_dataset(r'C:\Users\eg\Desktop\nec\sfcan20120101a20121231_rot_mask.nc')
dataset

<xarray.Dataset>
Dimensions:        (height: 1, rlat: 240, rlon: 280, time: 366)
Coordinates:
* rlon           (rlon) float64 -5.0 -4.95 -4.9 -4.85 ... 8.8 8.85 8.9 8.95
* rlat           (rlat) float64 -6.45 -6.4 -6.35 -6.3 ... 5.35 5.4 5.45 5.5
* height         (height) float64 0.0
* time           (time) datetime64[ns] 2012-01-01T06:00:00 ... 2012-12- 31T06:00:00
Data variables:
rotated_pole   |S1 ...
precipitation  (time, height, rlat, rlon) float32 ...
lon            (rlat, rlon) float32 ...
lat            (rlat, rlon) float32 ...
Attributes:
title:        AEMET High-resolution (0.05 deg) daily gridded precipitatio...
institution:  Agencia Estatal de Meteorologia (AEMET, www.aemet.es)
references:   Peral, C., Navascu�s, B., Ramos, P. Available at: http://ww...
history:      Creation year 2017
Conventions:  CF-1.7
version:      1.0

Кто-то знает, как повернуть его в градусах на восток и градусах на север ?, спасибо 

 <xarray.DataArray 'rlon' (rlon: 280)>
 array([-5.  , -4.95, -4.9 , ...,  8.85,  8.9 ,  8.95])
 Coordinates:
  * rlon     (rlon) float64 -5.0 -4.95 -4.9 -4.85 -4.8 ... 8.8 8.85 8.9 8.95
 Attributes:
long_name:      longitude in rotated pole grid
units:          degrees
standard_name:  grid_longitude
axis:           X

<xarray.DataArray 'lon' (rlat: 240, rlon: 280)>
 array([[-11.988, -11.928, -11.868, ...,   4.566,   4.625,   4.684],
   [-11.992, -11.932, -11.872, ...,   4.573,   4.633,   4.692],
   [-11.996, -11.936, -11.876, ...,   4.581,   4.64 ,   4.699],
   ...,
   [-13.136, -13.065, -12.994, ...,   6.57 ,   6.64 ,   6.71 ],
   [-13.142, -13.071, -13.   , ...,   6.58 ,   6.65 ,   6.72 ],
   [-13.148, -13.077, -13.006, ...,   6.59 ,   6.661,   6.731]],
  dtype=float32)
  Coordinates:
  * rlon     (rlon) float64 -5.0 -4.95 -4.9 -4.85 -4.8 ... 8.8 8.85 8.9 8.95
  * rlat     (rlat) float64 -6.45 -6.4 -6.35 -6.3 -6.25 ... 5.35 5.4 5.45 5.5
 Attributes:
long_name:  longitude
units:      degrees_east

То, что вы можете видеть, что они совершенно разные

...