求助；把等密度面插值出来？

呆妹小霸王

怎么才能从三维(lon-lat-depth)密度场中，把想要的等密度面的坐标找出来呢？（感觉就跟画等值线图一样，把某一根线找出来）
可以给我个大致思路或者某一个函数，让我知道怎么去做。(最好让我看下程序截图)

zqzq

使用find函数
k = find(X) 返回一个包含数组 X 中每个非零元素的线性索引的向量。
如果 X 为向量，则 find 返回方向与 X 相同的向量。
如果 X 为多维数组，则 find 返回由结果的线性索引组成的列向量。
如果 X 包含非零元素或为空，则 find 返回一个空数组。
[row,col] = find(___) 使用前面语法中的任何输入参数返回数组 X 中每个非零元素的行和列下标。

find不只可以用于查找数值为0的点的坐标
例如find(isnan(sst))等
也可以查找某一特定值的位置

呆妹小霸王

zqzq 发表于 2019-12-13 15:37
使用find函数
k = find(X) 返回一个包含数组 X 中每个非零元素的线性索引的向量。
如果 X 为向量，则 fin ...

谢谢你的回复。终于有人回我了，我太难了。find只能找数组中确实存在的一个数。
我现在是用interp1(density，depth，[26.5])；这样去插值处某一密度所处的深度，需要先把密度场中nan的值剔除出来。
我把密度随深度变化的线plot出来，看上去没有什么大问题。
还是感谢你的回复。

kaiofchina

我是新手来学习的

萌新

新手，求楼主请教

呆妹小霸王

萌新 发表于 2021-9-26 15:57
新手，求楼主请教

啊，什么问题

萌新

楼主，有代码分享吗

萌新

呆妹小霸王 发表于 2019-12-13 19:40
谢谢你的回复。终于有人回我了，我太难了。find只能找数组中确实存在的一个数。
我现在是用interp1(dens ...

楼主，您用的是python还是什么软件

呆妹小霸王

萌新 发表于 2022-2-16 21:07
楼主，您用的是python还是什么软件

我用的是matlab

代码就是我前面回复的这个interp1(density，depth，[26.5])