我试图将低于阈值的数组成员设置为nan.这是QA/QC过程的一部分,并且输入数据可能已经具有nan的时隙.
因此,作为示例,我的阈值可能是-1000,因此我想在以下数组中将-3000设置为nan
x = np.array([np.nan,1.,2.,-3000.,np.nan,5.])
以下内容:
x[x < -1000.] = np.nan
产生正确的行为,但也是RuntimeWarning,但是禁用警告的开销
warnings.filterwarnings("ignore")
...
warnints.resetwarnings()
有点沉重,可能有点不安全.
尝试使用花式索引进行两次索引,如下所示不会产生任何影响:
nonan = np.where(~np.isnan(x))[0]
x[nonan][x[nonan] < -1000.] = np.nan
我认为这是因为整数索引或使用索引两次进行复制.
有人有一个相对简单的解决方案吗？在这个过程中使用一个蒙面数组会很好,但最终的产品必须是一个ndarray,我不能引入新的依赖.谢谢.
我在Sphinx中记录了一个python模块.我有一个完整的模块使用示例的源代码文件.我想引用这个文件.内联作为连续代码太长了.有没有办法创建一个完整源文件的链接,格式化为代码友好的方式(即文字或行号)？
谢谢.
我正在尝试使用contourf在matplotlib中生成一个填充的等高线图.在图底部附近的锯齿状图案中缺少数据.轮廓图不仅在原始数据被遮盖的地方变成空白,而且在口袋中也是空白的,其中轮廓算法不能干净地插入,因为没有足够的良好数据邻域.
我知道如何扩展数据集以在这些口袋中产生合理的轮廓.但是,如果我绘制扩展数据,我会在任何地方获得轮廓填充.我想掩盖原始数据缺失的黑色或白色区域.
在之前的一个帖子中,我通过绘制第一张图像然后用掩盖坏区域的另一张图像覆盖它来学习如何为图像做到这一点.模拟将是下面的代码片段,但它不适用于轮廓......我无法通过bad_data imshow来掩盖扩展的contourf图.可能吗？
谢谢,Eli
import matplotlib.pyplot as plt
lev = [0.0,0.1,0.2,0.5,1.0,2.0,4.0,8.0,16.0,32.0]
norml = colors.BoundaryNorm(lev, 256)
# this is the contour plot, using extended_data so that the contours are plausibly extended
cs = plt.contourf(x,z,extended_data,levels = lev, cmap = cm.RdBu_r,norm = norml)
# now the attempt to cover it up -- but imshow will not cover up the original plot as it will with another image
bad_data = np.ma.masked_where(~data.mask, data.mask, copy=True)
plt.imshow(bad_data, interpolation='nearest', aspect = 'auto', cmap=cm.gray)
plt.show()
我希望我在sphinx中的toctree包含main(index.rst)页面.文档建议包含'self'会这样做,但正如其他人在stackoverflow上所说,这只包括主页面的标题,因此如果深度> 1,主页面的行为与链接页面不同.
我也看过为这个问题建议的内容指令,对我来说,这只是创建一个单独的本地内容列表,并没有帮助toctree.我真的在寻找能够在toctree中产生一致行为的解决方案.有吗？
我与一个开发基于MPI的C++数字应用程序的团队合作.该小组使用cxxtest来构建单个单元测试或小套件,但是1)在cxxtest的常用功能中存在一些并发目录,并且2)有一些集成测试通过启动mpirun从外部"更容易"实现一个python线程.
我们想使用py.test作为将它固定在一起的粘合剂,因为它宣称自己能够运行非python测试(我可以说服它跳到鼻子上).
谁能让我开始做这个的最佳实践？再次,因为它似乎是py.test的广告功能之一,我喜欢以最初设想的方式去做.
谢谢,
礼
将我的计算机从 1.6 升级到 1.8.3 (SVN 1.8.4) 的 Tortoise SVN 版本后，我遇到了可怕的 403 禁止错误。
我们的服务器是 1.8.1 版本，当时我们也从 Apache 2.2 切换到与 SVN 服务器捆绑的 CollabNet 版本。我确信存储库已升级到 1.8 版。我在存储库中被列为具有 rw 权限。
就防火墙或身份验证而言，1.8 有什么不同之处吗？Apache 服务器使用 html 并使用网络密码进行身份验证。什么是下一步？
I have hourly data records that were recorded in local daylight time (for me this is US/Pacific). These will be read in through csv. A gap exists at the beginning of DST at 02:00 when we spring forward. In fall, I believe that the data collected at 01:00 PDT is labeled 01:00 and the next hour is labeled 02:00 (and assumes PST).
I would like to translate the timestamps so they play well with other data stored in PST. Below …
我想在matplotlib(Win7,1.1.0)中创建一个填充的等高线图.我想突出显示某些值,并且级别更接近于log而不是线性.
有许多适合我的色彩图,但我忽略了我对cmap的选择.
我是否需要创建自定义"规范化"？如果是这样,每个轮廓根据其边缘值着色,然后用相同的颜色填充到下一个较低的值？为什么这个症状会忽略我的颜色图...这是一个例子,在构造期间被捕获并且我的请求被默默地忽略了？
我的原始数据缺少值.我一直在玩这个大小的......在每种情况下我都试过掩盖它们而不是屏蔽"外部"值.我也使用默认级别和规范尝试了所有排列.
lev = [0.1,0.2,0.5,1.0,2.0,4.0,8.0,16.0,32.0]
norml = colors.normalize(0,32)
cs = plt.contourf(x,z,data,cmap=cm.gray, levels=lev, norm = norml)
我希望这个片段足以至少开始对话.
谢谢,Eli