[ncl-talk] NDVI from Advanced Very High-Resolution Radiometer (AVHRR)
Dennis Shea
shea at ucar.edu
Tue Jun 21 11:39:31 MDT 2022
I downloaded 5 files for a test.
%> ncl NDVI_anomaly.ncl
Copyright (C) 1995-2017 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.5.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
Variable: fils
Type: string
Total Size: 40 bytes
5 values
Number of Dimensions: 1
Dimensions and sizes: [5]
Coordinates:
(0) AVHRR-Land_v005_AVH13C1_NOAA-09_19860101_c20170612092103.nc
(1) AVHRR-Land_v005_AVH13C1_NOAA-09_19860102_c20170612093403.nc
(2) AVHRR-Land_v005_AVH13C1_NOAA-09_19861229_c20170612212717.nc
(3) AVHRR-Land_v005_AVH13C1_NOAA-09_19861230_c20170612214313.nc
(4) AVHRR-Land_v005_AVH13C1_NOAA-09_19861231_c20170612215936.nc
Variable: prc
Type: float
Total Size: 518400000 bytes
129600000 values
Number of Dimensions: 3
Dimensions and sizes: [time | 5] x [latitude | 3600] x [longitude | 7200]
Coordinates:
time: [1826..2190]
latitude: [89.975..-89.97499]
longitude: [-179.975..179.975]
Number Of Attributes: 7
long_name : NOAA Climate Data Record of Normalized Difference
Vegetation Index
units : 1
grid_mapping : crs
standard_name : normalized_difference_vegetation_index
_FillValue_original : -9999
_FillValue : -9999
========================================================================
========================================================================
========================================================================
It seems to me that you have path or file problems.
I am not sure what else I can do.
Others?
On Tue, Jun 21, 2022 at 5:49 AM Setareh Rahimi <setareh.rahimi at gmail.com>
wrote:
> Dear Dennis,
> So many thanks for your advice. I run the script on a remote server with
> at least *1 terabyte* of memory. I think this is not the case.
> I used the below commands to check the files:
> ===========================================
> fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
> print(fils)
>
> f = addfiles (fils, "r")
>
> x = short2flt(f[:]->NDVI)
>
> printVarSummary(x)
> ===========================================
> , and attached the result. Please look at the image, and kindly help me
> sort the issue out. I do not know what is wrong with the files I downloaded.
> best regards,
>
> On Sat, Jun 18, 2022 at 10:52 PM Dennis Shea <shea at ucar.edu> wrote:
>
>> Based on the provided information, I have no idea why you would get
>> "returned nothing to me"
>>
>> I know that NDVI files are 0.1 degree resolution (3600,7200). You have 36
>> years (1986-2021) of daily data.
>> So, 36*365=13140 [+ a few Leap Year days)
>>
>> prc(13140,3600,7200) => daily data
>>
>> A lot of memory is needed. Maybe your computer can not handle this?
>>
>> After the
>> xM = calculate_monthly_values (prc, "avg", 0,opt)
>> printVarSummary(xM)
>>
>> xM(432,3600,7200)
>>
>>
>>
>> On Sat, Jun 18, 2022 at 5:01 AM Setareh Rahimi via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> Dear all,
>>>
>>> I need to calculate the monthly anomaly of the Normalized Difference
>>> Vegetation Index (NDVI), from Advanced Very High-Resolution Radiometer
>>> (AVHRR) sensor. I have downloaded AVHRR data (1986-2021) from the link
>>> below:
>>>
>>> https://www.ncei.noaa.gov/metadata/geoportal/rest/metadata/item/gov.noaa.ncdc:C01558/html
>>> I first tried to prepare a script to combine all the data, and then
>>> calculate the monthly average from daily data ( the link provides daily
>>> data). After having monthly data I tried to figure out anomalies (the
>>> script has been attached). But NCL stopped and returned nothing to me. I
>>> just wonder what is wrong. Please kindly advise me on this issue.
>>> Best wishes,
>>>
>>>
>>> --
>>> S.Rahimi
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at mailman.ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>
> --
> S.Rahimi
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220621/97d58106/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NDVI_anomaly.ncl
Type: application/octet-stream
Size: 811 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20220621/97d58106/attachment.obj>
More information about the ncl-talk
mailing list