Python удалить координаты со всеми специфическими переменными - netcdf - PullRequest
Новая
       12

Python удалить координаты со всеми специфическими переменными - netcdf

1 голос
/ 26 сентября 2019

У меня есть несколько файлов netcdf, и я хотел удалить данные с определенной характеристикой и сохранить их в новом файле netcdf.

Файл имеет координаты широты, долготы и времени с переменной температурой и ветром.Идея состоит в том, что если переменная ветра меньше 6 м / с, вся информация (для этой широты, долготы и времени) будет удалена для переменных ветра и температуры.Как я могу сделать это в Python?

У меня есть это пока

РЕДАКТИРОВАТЬ: файл .nc 

from netCDF4 import Dataset
from datetime import * 
import os

#-----------------------------------------------------------------------------------------------------------------------------------
path = '/home/Downloads/Arquivos_GOES_2018_OSISAF_NC/'

data_ini = datetime(2018,1,1,1,0,0) 
data_end = datetime(2018,7,1,1,0,0) 

IDval = 6 # minimo quality_level aceito

#-----------------------------------------------------------------------------------------------------------------------------------
DateTimeNow = data_ini 

while DateTimeNow<=data_end:


   namefile_net = DateTimeNow.strftime('%Y%m%d%H%M')+'00-OSISAF-L3C_GHRSST-SSTsubskin-GOES16-ssteqc_goes16_'+DateTimeNow.strftime('%Y%m%d_%H%M')+'00-v02.0-fv01.0.nc'


   if os.path.isfile(path+namefile_net) == True:


      data_net = Dataset(path+namefile_net)
      tsm = data_net.variables['sea_surface_temperature'][0,:]-273.15 
      lat = data_net.variables['lat'][:]
      lon = data_net.variables['lon'][:]
      wind = data_net.variables['wind_speed'][0,:]

      wind = np.where(qlv.mask==True,-1.0, qlv) 

      lon, lat = np.meshgrid(lon, lat) # criar uma matriz com lat e lon


      tsm_area = tsm[yid,xid] 
      lon_area = lon[yid,xid] 
      lat_area = lat[yid,xid] 
      wind_area = wind[yid,xid] 

      if np.sum(wind>=IDval) > 0:

         lon_area = np.where(qlv_area<IDval,-999.9,lon_area)
         lat_area = np.where(qlv_area<IDval,-999.9,lat_area)

         yid2, xid2 = seach_point(lon_area, lat_area, lon_point, lat_point)
         ```

Ответы [ 2 ]

1 голос
/ 27 сентября 2019

Кажется, ваши данные уже в куче.Вот идея (у вас нет вашего cdf, поэтому давайте представим, что столбец 3 - это ветер): 

>>> y = np.zeros((4,4))
>>> data = np.zeros((4,4))
>>> data[:,3] = [1,12, 4, 9]
>>> data
array([[  0.,   0.,   0.,   1.],
       [  0.,   0.,   0.,  12.],
       [  0.,   0.,   0.,   4.],
       [  0.,   0.,   0.,   9.]])
>>> data[ data[:,3] > 6]
array([[  0.,   0.,   0.,  12.],
       [  0.,   0.,   0.,   9.]])
0 голосов
/ 27 сентября 2019

Итак, я бы использовал следующую стратегию для задачи:

  • сделать копию исходного файла

  • заменить данные вновый файл

Вот пример с генератором данных в первой части: 

#!/usr/bin/env ipython
# ----------------------------------------
import numpy as np
from netCDF4 import Dataset,num2date,date2num
import os, datetime, subprocess
# ---------------------------------------
# Some input:
londata = np.arange(9.,31.,0.5);nx=np.size(londata);
latdata = np.arange(54.,66.,0.25);ny=np.size(latdata);
ntime = 4000;
timevals = np.linspace(0.,ntime*3600.0,ntime);
dataout = 0.+np.random.random((ntime,ny,nx))*25.0
# =========================================
# GENERATE SOME RANDOM FILE WITH ORIGINAL DATA:
genfile = 'data.nc'
ncout=Dataset(genfile,'w');
ncout.createDimension('lon',nx)
ncout.createDimension('lat',ny)
ncout.createDimension('time')
# ------------------------------------------
xvar = ncout.createVariable('lon','float32',('lon'));xvar[:]=londata
yvar = ncout.createVariable('lat','float32',('lat'));yvar[:]=latdata
tvar = ncout.createVariable('time','float64',('time'));tvar[:] = timevals;tvar.setncattr('units','seconds since 2019-01-01 00:00:00');
dvar = ncout.createVariable('wspeed','float32',('time','lat','lon'));
dvar[:] = dataout;
ncout.close();
# =========================================
# REMOVEDATA:
filein = genfile;
fileout = filein.replace('.nc','.modified.nc')
cstr = 'cp '+filein+' '+fileout;
subprocess.call(cstr,shell=True); # Let us make the copy of the file, so that we do not have to copy the variables on our own
# ----------------------------------------
ncout = Dataset(fileout,'a');
datain = ncout.variables['wspeed'][:];
dataout = np.array(datain,copy=True);# I am making a real copy of the data, just in case I want to compare afterwards
dataout[dataout<6.0] = np.nan; # So here I replace the values
ncout.variables['wspeed'][:] = dataout; # I replace the values in the file
ncout.close();
...