编辑-地球科学江南体育网页版堆栈交换江南电子竞技平台

你不是登录。你的编辑将被放置在一个队列,直到同行评议。

我们欢迎编辑,使读者更容易理解和更有价值。因为社区成员审查编辑,请尽量使文章更好比你怎么发现,例如,通过修正语法或添加额外的资源和超链接。

牧师

插入高斯电网常规固定网格使用双线性插值?

我想regrid GCM数据可用高斯电网常规电网,说在1度X 1度。如何做到这一点CDO或Python或MATLAB。

编辑1:数据格式是netcdf。

回答

#气候数据运营商(CDO) # # #定义网格# #我们定义一个lat-lon目标网格与1°x1°网格单元尺寸30 x30网格细胞从40°N和-10°E (= 10°W): gridtype = lonlat xsize = 30 ysize = 30 xfirst = -10 xinc = 1 yfirst = 40 yinc = 1这段文本写入一个文本文件。[1.3.2节CDO手册][1]为进一步的细节和例子。如果你已经有一个netCDF文件数据在你的目标网格,您还可以从该文件中提取网格定义:cdo griddes FILE_WITH_TARGET_GRID。数控> myGridDef myGridDef是一个文本文件。# # # #插入数据,我认为我们做了双线性插值。这是完成了“remapbil”运营商通过:cdo remapbil myGridDef INPUT_FILE。数控OUTPUT_FILE。数控的网格需要正确输入数据输入文件中定义(通常是一个问题)。# #正确定义的网格替代1 # # # # # #数据变量的空间维度下面的例子(SST)需要被命名为“朗”和“纬度”。此外,协调变量的经度和纬度的存在(时间独立!)。这些需要正确的坐标变量属性(见下面的例子)。例如:维度:经度= 30; lat = 30 ; time = UNLIMITED ; // (24 currently) variables: double lon(lon) ; lon:standard_name = "longitude" ; lon:long_name = "longitude" ; lon:units = "degrees_east" ; lon:axis = "X" ; double lat(lat) ; lat:standard_name = "latitude" ; lat:long_name = "latitude" ; lat:units = "degrees_north" ; lat:axis = "Y" ; double time(time) ; time:standard_name = "time" ; time:long_name = "time" ; time:units = "seconds since 1900-01-01 00:00:00" ; time:calendar = "standard" ; time:axis = "T" ; float SST(time, lat, lon) ; SST:long_name = "sea_surface_temperature" ; SST:units = "degree celsius" ; SST:_FillValue = NaNf ; SST:missing_value = NaNf ; SST:var_desc = "sea surface temperature" ; ### Alternative 2 ### The spatial dimensions and corresponding variables exist in the source file but are not denoted as `lon` and `lat`. Then each data variable (SST, in this example) needs to have an attribute `coordinates` that holds the names of the coordinate variables. SST:coordinates = "lon lat" ; Example: dimensions: TSTEP = UNLIMITED ; // (24 currently) COL = 112 ; ROW = 106 ; variables: " ; double lon(ROW, COL) ; lon:standard_name = "longitude" ; lon:long_name = "longitude coordinate" ; lon:units = "degrees_east" ; lon:_CoordinateAxisType = "Lon" ; double lat(ROW, COL) ; lat:standard_name = "latitude" ; lat:long_name = "latitude coordinate" ; lat:units = "degrees_north" ; lat:_CoordinateAxisType = "Lat" ; double time(TSTEP) ; time:standard_name = "time" ; time:long_name = "time" ; time:units = "seconds since 1900-01-01 00:00:00" ; time:calendar = "standard" ; float SST(TSTEP, ROW, COL) ; SST:long_name = "sea_surface_temperature" ; SST:units = "degree celsius" ; SST:coordinates = "lon lat" ; SST:var_desc = "sea surface temperature" ; ### Note: ### If the coordinate variables are missing, then they can be created via the `setgrid` operator. If we want to add a grid definition to the source file then we need the source file's grid definition in the text file (see 'define grid'). cdo setgrid,mySourceGridDef INFILE_NO_COORDS.nc INFILE_WI_COORDS.nc [1]: https://code.zmaw.de/projects/cdo/embedded/index.html#x1-150001.3.2

编辑总结

取消

\ begingroup美元有一个同样简单的方法进行设置投影坐标系的坐标?我不懂太多的文档。 \ endgroup美元
- - - - - -怪不得我
2017年11月2日上午
\ begingroup美元 @Janina你的意思是当你有一些投影,这不是项目模型网格在常规lat-lon网格?你的意思是插值或只是设置一个网格的定义?更好的你问一个新的问题(指这一个),包括您的数据集的一个例子。 \ endgroup美元
- - - - - -daniel.heydebreck
2017年11月2日,9:48
\ begingroup美元是的,没错,兰伯特正形圆锥等。从本质上讲,我想设置网格的定义,所以,例如在情节中的数据正确的地理位置。是的,最终我对regridding感兴趣。好的,我想我还是把这变成一个问题…:-) \ endgroup美元
- - - - - -怪不得我
2017年11月2日,19
\ begingroup美元 @Janina我重读我的回答。实际上解决了你的问题,如果你已经有你的坐标中可用一个netCDF文件。 \ endgroup美元
- - - - - -daniel.heydebreck
2017年11月2日,31
1

\ begingroup美元实际上,你是对的。经过反复试验进展我能够做regridding成功——谢谢你的详细的描述!但我的工作流程的有效性是有问题的。我创建并保存数据与原始坐标REPSG: 3333。我解决了运用问题通过添加投影的元数据。我转换了WGS84坐标的细胞中心坐标和拯救他们变量经度和纬度(就像你选择2)的输出griddes这个文件是一个曲线网格。 \ endgroup美元
- - - - - -怪不得我
2017年11月3日远

添加一个评论 |

正确的小错误或错误
澄清的意思,而不更改它
添加相关资源或链接
总是尊重作者的意图
不使用编辑与作者回复

创建代码用引号”或式子~篱笆
' ' '
像这样
' ' '
添加语言标识符来突出显示的代码
”“python
def函数(foo):
打印(foo)
' ' '
把段落之间的回报
在结束linebreak添加2个空间
_italic_或* *的* *
报价将>开始的
做链接(尽量使用https)
< https://example.com >
(例子)(https://example.com)
< a href = " https://example.com " > < / >示例
MathJax方程罪\ \ ^ 2θ美元

格式帮助»
回答帮助»

MathJax帮助»