почему gdal дал мне разные максимальные и минимальные значения переменной из R?

Эти три строки откроют и прочитают одну переменную из файла netcdf:

f=open.ncdf("C:\\BR_Ji1.nc")
A = get.var.ncdf(nc=f,varid="date",verbose=TRUE)
B = get.var.ncdf(nc=f,varid="GPP",verbose=TRUE)
[1] "get.var.ncdf: entering. Here is varid:"
[1] "GPP"                 
[1] "checking to see if passed varid is actually a dimvar"
[1] "entering vobjtodimname with varid= GPP"
[1] "vobjtodimname: is a character type varid.  This file has 6 dims"
[1] "vobjtodimname: no cases found, returning FALSE"
[1] "get.var.ncdf: isdimvar: FALSE"
[1] "vobjtovarid: entering with varid=GPP"
[1] "Variable named GPP found in file with varid= 17"
[1] "vobjtovarid: returning with varid deduced from name; varid= 17"
[1] "get.var.ncdf: ending up using varid= 17"
[1] "ndims: 2"            
[1] "get.var.ncdf: varsize:" 
[1]     1 17520          
[1] "get.var.ncdf: start:" 
[1] 1 1                     
[1] "get.var.ncdf: count:"
[1]     1 17520            
[1] "get.var.ncdf: totvarsize: 17520"
[1] "Getting var of type 3  (1=short, 2=int, 3=float, 4=double, 5=char, 6=byte)"
[1] "get.var.ncdf: C call returned 0"
[1] "count.nodegen: 17520    Length of data: 17520"
[1] "get.var.ncdf: final dims of returned array:"
[1] 17520       
[1] "varid: 17" 
[1] "nc$varid2Rindex: 1"  "nc$varid2Rindex: 2"  "nc$varid2Rindex: 0"  "nc$varid2Rindex: 0"  "nc$varid2Rindex: 3" 
[6] "nc$varid2Rindex: 0"  "nc$varid2Rindex: 4"  "nc$varid2Rindex: 5"  "nc$varid2Rindex: 6"  "nc$varid2Rindex: 7" 
[11] "nc$varid2Rindex: 8"  "nc$varid2Rindex: 9"  "nc$varid2Rindex: 10" "nc$varid2Rindex: 11" "nc$varid2Rindex: 12" 
[16] "nc$varid2Rindex: 13" "nc$varid2Rindex: 14" "nc$varid2Rindex: 15" "nc$varid2Rindex: 16" "nc$varid2Rindex: 17"
[21] "nc$varid2Rindex: 18" "nc$varid2Rindex: 19" "nc$varid2Rindex: 20"
[1] "nc$varid2Rindex[varid]: 14"
[1] "get.var.ncdf: setting missing values to NA"
[1] "missval: -9999   tol: 0.09999"
[1] "get.var.ncdf: implementing add_offset ( FALSE ) and scale_factor ( FALSE )"
[1] "var has NEITHER add_offset nor scale_factor"
**K = get.var.ncdf(nc=f,varid="Qle",verbose=TRUE)
write.table(t(rbind(A,B,K)),"C:\\Ji1-gpp-lat.txt")**

Когда я посмотрел на извлеченные данные, я обнаружил, что у меня очень маленькие и отрицательные значения. Затем я использовал gdalinfo и обнаружил, что максимальное значение равно 400, а минимальное равно 0 (что нормально, и именно так должны выглядеть значения). Любые идеи, почему я получаю номер неправильно?


r netcdf rgdal gdal
person Sami Yemein    schedule 16.10.2012    source источник
comment
Вы уверены, что это не проблема дисплея? Например, если gdalinfo форматирует числа как целые числа, отрицательные числа ничтожной величины будут отображаться как 0.   -  person Carl Witthoft    schedule 16.10.2012
comment
Повторите gdalinfo с -stats, чтобы принудительно восстановить статистику. Что сообщил R? Еще один хороший инструмент, который вы можете использовать, если знаете Python/NumPy, — это netcdf4-python.   -  person Mike T    schedule 23.10.2012

Ответы (1)


arrow_upward
0
arrow_downward

См. этот вопрос: Учитывает ли ComputeBandStats отсутствие данных?. Статистика, содержащаяся в метаданных файла, может быть устаревшей или точно рассчитанной. Если вам нужны точные статистические данные, вы должны рассчитать их самостоятельно или убедиться, что флаг расчета с использованием приближений установлен на False.

person Benjamin    schedule 08.02.2017