科学网—GuoKeCCC的博文

[转载]matlab求两坐标点组成直线的方位角

李祥飞 2022-7-4 11:13

参考代码：https://zhuanlan.zhihu.com/p/139760930 在其基础上进行的修改。 function =GetAzimuth_2points(point1,point2) %%功能：求Point1，相对与Point2的相对方位和距离 %point的格式： %批量操作 %point_Set= ; %point2= ...

个人分类: MATLAB|2001 次阅读|没有评论

matlab导出矢量图并插入word

李祥飞 2022-5-23 17:41

matlab导出eps矢量图 print(gcf, '-depsc', 'abc.eps'); 得到的eps文件已将四周的白边去除，可以打开并另存为pdf eps转emf 在线工具https://cloudconvert.com/emf-converter 将得到的emf直接复制粘贴到word中即可这样得到的word转为pdf后，图片不失真 ...

个人分类: MATLAB|4705 次阅读|1 个评论

CDO计算每个月的冻结融化指数(FDD/TDD)

李祥飞 2022-5-22 18:17

1. 气温单位从K转变到℃：cdo addc,-273.15 input.nc output_K2C.nc addc - 加上一个常数；input.nc为日平均温度；output_K2C.nc为输出文件；-b F32是为了转为浮点型，否则若input.nc中变量不是浮点型，可能会报numeric conversion not representable的错误。 2. 获取指示日平均温度是否大于0的mask：cdo -gtc,0 ou ...

个人分类: CDO|2456 次阅读|没有评论

根据FES2014潮汐分量计算tidal range

李祥飞 2022-5-1 17:24

【说明】并非潮汐相关专业，相关结果是在查阅相关资料后的理解，不对之处，敬请指正！ FES2014是一个全球潮汐模型，输出了34个潮汐分量（分潮，tide component），比较重要的有K1, O1, M2, S2等四个分量。首先，计算form ration F = (K1 + O1)/(M2 + S2)，对潮汐进行分类： F = 0.00 - 0.25: semidiurnal ...

个人分类: Shoreline|2986 次阅读|没有评论

CDO每隔n段时间（小时，天，月份等）统计变量（求和，平均，最大最小等）

李祥飞 2022-4-30 12:07

越来越感觉到CDO的强大了，速度快，又不占内存。我有一个大小为40G的1440x257x20x2928(lon, lat, level, time)的nc文件，想要在time维度上每隔8个计算一个平均值，命令如下： cdo timselmean,8 1980_Ta.nc 1980_Ta_daily.nc 1980_Ta.nc和1980_Ta_daily.nc分别为输入输出文件，最终输出的变量维度为1440x257x20x ...

个人分类: CDO|5015 次阅读|没有评论

CDO将曲线网格（curve linear）重采样为经纬网格（lonlat）

李祥飞 2022-4-29 16:52

一些气候模式的输出数据的格式并非规则的经纬网格，而是网格大小不一的曲线网格。比如CESM模式的海洋模块使用的网格类型（grid type）为displaced-pole grid，意为把北极点移到了格陵兰岛。此网格的代号gx1v6，水平分辨率1.125° × 0.27~0.54°。现在我们需要把数据重采样为0.25° × 0.25°的规则网格数据，使用CDO ...

个人分类: CDO|5195 次阅读|没有评论

CDO (Climate Data Operators)安装遇到的问题及解决

李祥飞 2021-11-24 21:05

CDO (Climate Data Operators)是一款处理常用气象数据格式（如netcdf）的利器，但不支持windows环境。好在win10开始内嵌了linux环境（WSL, Windows Subsystem for Linux），可与windows直接互访。其安装可参考官方网站：https://docs.microsoft.com/en-us/windows/wsl/install-manual，注意WSL1和WSL2的区别。安装好 ...

个人分类: CDO|5369 次阅读|没有评论

matlab无法打开并行

李祥飞 2021-11-12 13:10

【问题】输入 p=parpool运行后报错： Starting parallel pool (parpool) using the 'local' profile ... Error using parpool (line 103) Failed to start a parallel pool. (For information in addition to the causing error, validate the profile 'local' in the Cluster Profile Manager. ...

个人分类: MATLAB|4086 次阅读|没有评论

ArcGIS报错：The truth value of a raster is ambiguous

李祥飞 2021-8-31 17:35

使用ArcGIS的Raster Calculator时报错The truth value of a raster is ambiguous，原因是判断语句没有加括号，无法判断bool值。当涉及多个判断语句时，每个判断语句最好都加上括号。比如： Con((ESACCI-2010.tif == 160) | (ESACCI-2010.tif == 170) | (ESACCI-2010.tif == 180), 4, 1)，如果去掉括 ...

个人分类: GIS|2806 次阅读|没有评论

Matlab对矩阵元素进行分类

热度 1 李祥飞 2021-8-31 17:12

目标：对ESA CCI的土地利用类型（表1）进行重分类（表2）。代码： % 构建5*15的分类矩阵veg_class，5代表重分类后大类的总数，15表示大类中最大的子类个数，其余用NaN填充。注意：要把重分类后大类值最小（此处是0）的子类放在最后一行。 veg_class = ; % 待重分类的矩阵veg_type veg_type = ; % 找 ...

个人分类: MATLAB|3248 次阅读|2 个评论热度 1

帐号		自动登录	找回密码
密码			注册

关闭 安全验证

GuoKeCCC

关闭安全验证