抱歉我的英语不好。 我在实习,在此之前从未使用过Python。 我需要从NetCDF文件中提取数据,
我已经创建了一个创建DataFrame的循环但是当我尝试提取这个DataFrame时,我在41000上只有201个值。
import csv
import numpy as np
import pandas as pd
import netCDF4
from netCDF4 import Dataset, num2date
nc = Dataset('Q:/QGIS/2011001.nc', 'r')
chla = nc.variables['chlorophyll_a'][0]
lons = nc.variables['lon'][:]
lat = nc.variables['lat'][:]
time = nc.variables['time'][:]
nlons=len(lons)
nlat=len(lat)
第一个循环在ArcGIS python控制台中为我提供了41000个值
for i in range(0,nlat) :
dla = {'lat':lat[i],'long':lons,'chla':chla[i]}
z = pd.DataFrame(dla)
print (z)
z.to_csv('Q:/QGIS/fichier.csv', sep =',', index= True)
但是当我执行to.csv时,我只能在csv文件中获得201个值。
for y in range(0,nlat):
q[y].to_csv('Q:/QGIS/fichier.csv', sep =',', index= True)
for i in range(0,nlat):
dlo ={'lat':lat[i],'long':lons,'chla':chla[i]}
q[y] = pd.DataFrame(dlo)
print(q)
我希望你能有一个解决这个问题的答案,而且如果你有任何脚本来提取创建一个shp文件的值,我将非常感激,如果你能分享它!
祝你好运
提前谢谢
DAM