Question

我正在尝试从月平均chl .nc卫星文件（经过处理的L2至L3产品）中提取叶绿素，该文件的地理边界坐标为-25 S至-12 N，经度-179 W至-169 E，那里1443 * 1110浮点数组元素。我必须编写一个python代码，让我输入自己的纬度和经度，然后python会打印出在该位置找到的相应chl值。

这是我的代码的起始行，错误在于整形行。

#!/usr/bin/env python

import numpy as np
import pylab

#input user's selection of latitudes and longitudes points

lat = float(input('enter your latitude here: '))
print 'lat: ',lat
lon = float(input('enter your longitude here: '))
print 'lon: ',lon 

#Call path directory, call file name, call variable 'fromfile'

path = '/home/rst/data/final_data/2002_chl_sst/Jul/Output-L3/monthly /1000m-chlor_a/mean/A20021822002212.MON.chlor_a.map.nc'
fname = open(path) 
chl = np.fromfile(fname, dtype=np.float64) 
fname.close() 

#reshape the data into 2D matrix array

nrow = 1443.0
ncol = 1110.0
chl = chl.reshape([nrow,ncol])      # this is line 21
print chl.shape

这是我代码中的第一个错误，其余代码未在此处显示。这个新数组的总数是多少？感谢任何帮助，谢谢

回溯（最近通话最近）：在第21行的文件“ prob1_extract_chl_2002JUL.py” chl = chl.reshape（[nrow，ncol]） ValueError：新数组的总大小必须保持不变

在python中使用重塑函数时如何修复valueerror

0 个答案: