How to Regrid the netcdf data in order to achieve smaller/greater spatial resolution?

Question

BN 2019-10-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483799-how-to-regrid-the-netcdf-data-in-order-to-achieve-smaller-greater-spatial-resolution

评论： Suman 2023-11-14

I have a netcdf file with 0.5*0.5 degree (Lat*Lon) grided resolution for precipitation from 1982-2015 in all around the globe. now I want to interpolate (regrid) this file to 0.25*0.25 degree without changing the time. In fact I want to interpolate the precipitation data in 0.25*0.25 lat and lon.

please help me with this issue. If it is possible a sample code can really helpful for me.

Best Regards,

Behzad

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Umi Pollmann 2021-5-11

Has this issue been resolved? I face the same issue at the moment. Thank you.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Chad Greene 2019-10-7

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483799-how-to-regrid-the-netcdf-data-in-order-to-achieve-smaller-greater-spatial-resolution#answer_395212

在 MATLAB Online 中打开

Hey Behzad,

This is pretty straightforward: Just use interp2.

When reading the data from the NetCDF file you might need to rotate and flip the grid to get the orientation right, but then let longitude act as the x variable and let latitude act as the y variable.

You can create the quarter-degree grid in one step with cdtgrid like this:

[lati,loni] = cdtgrid(1/4);

or you can create it manually by

[loni,lati] = meshgrid(min(lon):1/4max(lon),max(lat):-1/4:min(lat));

Then use interp2 like

zi = interp2(lon,lat,z,loni,lati);

Hope that helps.

3 个评论
显示 1更早的评论隐藏 1更早的评论

BN 2019-10-7

在 MATLAB Online 中打开

Dear Dr. Chad Greene

Thank you for your answer. I do wath you say and I faced this error. do you have any idea?

my script is:

clc
clear
format compact
close all
filename='precip.mon.total.v2018.nc';
ncdisp(filename)       
lat = ncread(filename,'lat');     
lon = ncread(filename,'lon');  
precip = ncread(filename,'precip');
t = ncread(filename,'time');
lat=double(lat);
lon=double(lon);
precip=double(precip);
t=double(t);
index=find(precip==-9.969209968386869e+36);
precip(index)=nan;
[lati,loni] = cdtgrid(1/4); 
zi = interp2(lon,lat,precip,loni,lati); 

command window show:

Source:
           C:\Users\Behzad\Desktop\precip.mon.total.v2018.nc
Format:
           netcdf4_classic
Global Attributes:
           Original_Source = 'http://www.dwd.de/en/FundE/Klima/KLIS/int/GPCC/GPCC.htm
                              is the webpage and the data is at ftp://ftp.dwd.de/pub/data/gpcc/download.html'
           Reference       = 'Users of the data sets are kindly requested to give feed back and to refer to GPCC publications on this webpage: http://www.dwd.de/bvbw/appmanager/bvbw/dwdwwwDesktop/?_nfpb=true&_pageLabel=_dwdwww_klima_umwelt_datenzentren_wzn&T12404518261141645246564gsbDocumentPath=Content%2FOeffentlichkeit%2FKU%2FKU4%2FKU42%2Fteaser__product__access.html&_state=maximized&_windowLabel=T12404518261141645246564&lastPageLabel=_dwdwww_klima_umwelt_datenzentren_wzn'
           original_source = 'ftp://ftp-anon.dwd.de/pub/data/gpcc/html/fulldata_download.html'
           Conventions     = 'CF 1.0'
           dataset_title   = 'Global Precipitation Climatology Centre (GPCC)'
           References      = 'https://www.esrl.noaa.gov/psd/data/gridded/data.gpcc.html'
           title           = 'GPCC Full Data Reanalysis Version 2018 0.5x0.5 Monthly Total'
           history         = 'Created 09/2018  based on V2018 data obtained via ftp'
           data_modified   = '2019-03-12'
           _NCProperties   = 'version=2,netcdf=4.6.3,hdf5=1.10.5'
Dimensions:
           lat   = 360
           lon   = 720
           nbnds = 2
           time  = 1512  (UNLIMITED)
Variables:
    lat   
           Size:       360x1
           Dimensions: lat
           Datatype:   single
           Attributes:
                       long_name          = 'Latitude'
                       units              = 'degrees_north'
                       standard_name      = 'latitude'
                       axis               = 'Y'
                       coordinate_defines = 'point'
                       actual_range       = [89.75       -89.75]
    lon   
           Size:       720x1
           Dimensions: lon
           Datatype:   single
           Attributes:
                       long_name          = 'Longitude'
                       units              = 'degrees_east'
                       standard_name      = 'longitude'
                       actual_range       = [0.25        359.75]
                       axis               = 'X'
                       coordinate_defines = 'point'
    time  
           Size:       1512x1
           Dimensions: time
           Datatype:   double
           Attributes:
                       long_name          = 'Time'
                       units              = 'days since 1800-1-1 00:00:00'
                       delta_t            = '0000-01-00 00:00:00'
                       avg_period         = '0000-01-00 00:00:00'
                       standard_name      = 'time'
                       axis               = 'T'
                       coordinate_defines = 'start'
                       actual_range       = [33237  79227]
    precip
           Size:       720x360x1512
           Dimensions: lon,lat,time
           Datatype:   single
           Attributes:
                       missing_value = -9.969209968386869e+36
                       units         = 'mm'
                       var_desc      = 'Precipitation'
                       level_desc    = 'Surface'
                       parent_stat   = 'Observations'
                       long_name     = 'GPCC Monthly total of precipitation'
                       valid_range   = [0  8000]
                       statistic     = 'Total'
                       level         = 'Surface'
                       actual_range  = [0      4552.4302]
                       dataset       = 'GPCC Precipitation 0.5degree V2018 Full Reanalysis'
Error using  .' 
Transpose on ND array is not defined. Use PERMUTE instead.
Error in interp2 (line 122)
        V = V.';
Error in test100 (line 23)
zi = interp2(lon,lat,precip,loni,lati); 

Thank you and I looking forward to your answer.

Best Regards

O.Hubert 2022-3-16

在 MATLAB Online 中打开

You probably figured it out by now, but the error you have arises because you try to transpose a 3-dimensional array (longitude, latitude, time) using a shorthand for 2-dimensional arrays. If you want to transpose longitude and latitude without touching the time dimension, use:

new_array = permute(old_array,[2 1 3]);

Suman 2023-11-14

Where to insert this code?

请先登录，再进行评论。

Answer 2

abbas zoulikha 2020-11-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483799-how-to-regrid-the-netcdf-data-in-order-to-achieve-smaller-greater-spatial-resolution#answer_541453

i have the same error in my script how did you resolved this errors ?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

How to Regrid the netcdf data in order to achieve smaller/greater spatial resolution?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

回答（2 个）

3 个评论
显示 1更早的评论隐藏 1更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

Community Treasure Hunt

How to Regrid the netcdf data in order to achieve smaller/greater spatial resolution?

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

回答（2 个）

3 个评论 显示 1更早的评论隐藏 1更早的评论

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

Community Treasure Hunt

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论