我正在尝试用Anaconda3编写一个python包，其中包含一些编译过的Fortran库。根据在线提供的示例，我需要使用numpy.distutils.core.Extension来正确构建扩展。下面是我的最小工作示例目录：
MANIFEST.in
setup.py
bin
└─── test_smat
mypkg
├─── __init__.py
├─── compiled
| ├─── __init__.py
| └─── gaussian.f90
└─── pystuff
├─── __init__.py
└─── buildma
浏览 4提问于2018-08-04得票数 1
我在Linux下找不到任何关于在swift中使用gnu科学库的资料。有可能吗？如果是这样，那么是如何实现的呢？
例如，我对用gsl_ran_gaussian生成正态分布的随机数很感兴趣。
浏览 0提问于2016-06-07得票数 0
我正在生成一个random.uniform(low=0.0, high=100.0, size=(150,150))数组。
我将其输入到生成X、x和y的函数中。
但是，如果随机测试矩阵大于100，则得到以下误差。
我试过玩弄θ值。
有人有这个问题吗？这是个虫子吗？
我正在使用python2.6和scikit-learn 0.10。我应该试试python3吗？
欢迎任何建议或意见。
谢谢。
gp.fit( XKrn, yKrn )
File "/usr/lib/python2.6/scikit_learn-0.10_git-py2.6-linux-x86_64.egg/sklearn
浏览 5提问于2012-11-13得票数 0
1回答
我编写了一个类，并开始在同一个.cpp文件中编写所有内容。然而，过了一段时间，我发现这个类越来越大了，所以我决定把它分成一个.h和一个.cpp文件。
gaussian.h文件：
class Gaussian{
private:
double mean;
double standardDeviation;
double variance;
double precision;
double precisionMean;
public:
Gaussian(double, double);
浏览 8提问于2012-01-29得票数 12
回答已采纳
我想为标量和向量实现一个一般的高斯函数。在Ada或C++中，我只是选择模板，但在Haskell中，这有点让人费解。
首先，定义一个可以应用高斯运算符的类，例如组合或计算概率：
class Gaussian g where
(!*) :: g -> g -> g
prob :: g -> a -> Float -- Here, I want a to be depending on g
data Gaussian1D = Gaussian1D Float Float
data Gaussian2D = Gaussian2D (Linear.V2 Float) L
浏览 2提问于2020-12-07得票数 3
回答已采纳
1回答
目前我的代码(在GitHub)是这样的：
module Gaussian where
import qualified Data.Complex as Complex
-- Enables toComplex.
-- Qualified to avoid clashing.
data Gaussian = (:+)
{
real :: Integer,
imaginary :: Integer
} deriving
(Show, Read, Eq)
-- The gaussian integer type.
-- A gaussian int
浏览 0提问于2019-09-04得票数 5
回答已采纳
我使用的是Scilab版本6.1.0，操作系统是Ubuntu 18.04。我正在为已经存在的库向Scilab添加一个工具箱。我已经在Scilab中分别接口了库的两个函数。我不能同时连接它们。我认为原因是"builder_gateway_cpp.sce“文件中的Linker_Flag变量，下面是该文件的内容：
mode(-1)
lines(0)
toolbox_title = "test_toolbox";
Build_64Bits = %t;
path_builder = get_absolute_file_path('builder_gateway_cpp.s
浏览 15提问于2020-06-14得票数 0
1回答
我试图使用以下算法实现图像扩展：
1. Zero pad: Put a zero pixel between two pixels
2. Convolve with a column vector of a Gaussian, and multiply by 2
3. Convolve with a row vector of a Gaussian and multiply by 2
4. Return blurred image.
下面是我使用的代码(我用3表示guassian大小)：
def expand(im, filter_size):
''' Takes
浏览 2提问于2016-12-23得票数 0
回答已采纳
1回答
我正在做一个项目，用OpenCV库将图像混合在一起。当我尝试重建金字塔的层时，我得到了一个错误，即使我认为设置大小是正确的。在代码的末尾，我得到了这个错误消息(在#Reconstructed Faces)，感谢您的帮助：
Traceback (most recent call last):
File "morph.py", line 104, in <module>
faces_reconstructed = cv2.pyrUp(faces_reconstructed, dstsize=size)
cv2.error: OpenCV(4.1.1) C:\
浏览 4提问于2019-11-29得票数 1
我正在创建一个新的库仑到一个现有的数据火花，通过搜索一个已归档的“脚本”，并返回匹配作为新列的条目。
import re as re
def sw_fix(data_str):
if re.compile(r'gaussian').search(data_str):
cleaned_str = 'gaussian'
elif re.compile(r'gromacs').search(data_str):
cleaned_str = 'gromacs&#
浏览 3提问于2017-05-22得票数 0
回答已采纳
2回答
我已经创建了一个名为FitFunction的工厂类，它在前面展示的内容之外添加了一大堆东西。标签方法pretty_string应该只返回写入的字符串。当我运行这个文件时，它会打印一个与repr一样有用的字符串。有没有人知道我会怎么去实现它？
#!/usr/bin/env python
from __future__ import print_function, absolute_import
import abc
import types
import numpy as np
class FitFunction(object):
def __init__(self, pytho
浏览 1提问于2016-09-30得票数 1
1回答
我对shell脚本很陌生，我试着理解为什么它在命令之后停止运行，该命令将一个作业提交给高斯(一个计算化学程序)。这是我的密码。
echo "hi, the script is running..."
location=$(dirname $1)
echo $location
g09 <<EEE
%chk=test2
$(cat $1)
EEE
echo $?
echo "done"
echo "need help??????" >tttt.txt
它没有显示最后三个回波，而计算程序正常结束。
我的假设是，当程序退出时，我的脚本
浏览 1提问于2016-06-21得票数 0
回答已采纳
这里是我的解析器的一部分。问题是它没有读取和打印DeltaE和强度，也许有人能帮上忙？
if (line[1:24] == "Mulliken atomic charges" or
line[1:22] == "Lowdin Atomic Charges"):
if not hasattr(self, "atomcharges"):
self.atomcharges = {}
ones = inputfile.next()
charges = []
浏览 3提问于2015-11-17得票数 1
回答已采纳
1回答
有人能给我解释一下为什么我的代码不能正常运行吗？我很久没有使用python了，忘记了很多东西，我不知道出了什么问题，请帮帮忙。下面是我的代码和结果：
def my_gaussian(x, mu=0, sigma=1):
"""
Compute the Gaussian function.
Parameters
----------
x : float or array
A single value or an array of values for which the Gaussian will be
eval
浏览 0提问于2020-10-06得票数 0
3回答
我有一个线程安全的课程：
internal class ExmplFile
{
private readonly string filename;
private int resolution;
private volatile PaddedImage gaussian;
private object lockObject=new object();
//blah, blah, blah
internal PaddedImage Gaussian()
{
if (gaussian != null)
{
浏览 5提问于2016-07-18得票数 1
回答已采纳
我使用的是Idris 2 v0.3。我不明白为什么编译器找不到Distribution的Gaussian在
interface Distribution dist where
mean : dist -> Double
data Gaussian : Type where
MkGaussian : Double -> Gaussian
Distribution Gaussian where
mean (MkGaussian mean') = mean'
KnowledgeBased : Distribution d => (d : Typ
浏览 7提问于2021-05-19得票数 1
回答已采纳
我有一个很简单的规则，我已经得到了质量乐队的问题。有人能解释一下为什么会这样吗？我读过关于的文档，但我不明白。
这是我编写的产生警告的代码：
Gaussian BornYoungPerLitterDist = new Gaussian(5.5, 0);
Gaussian BirthRateDist = new Gaussian(0.004102103451, 0.00000144);
Variable<double> BornYoungPerLitter = Variable.Random<double, Gaussian>(BornYoungPerLitterDis
浏览 1提问于2019-01-23得票数 0
回答已采纳
2回答
Python不锐化掩模
、、、、
我想在16位图像上使用unsharp mask。
图像具有640 x 480 Pixels，并保存在NumPy数组中。
我做了以下工作：
用Gaussian filter模糊图像(三种不同的方法)
然后，通过从原始图像中减去模糊图像创建一个掩码。
最后，将掩码乘以WightFaktor添加到原始图像中
但这不起作用。
下面是Python代码：
Gaussian1 = ndimage.filters.gaussian_filter(Image,sigma=10.0)
Gaussian2 = filters.gaussian_filter(Image,sigma=10.0)
G
浏览 11提问于2015-09-08得票数 15
目前，我有以下代码
def approx_binomial(n, p, size=None):
gaussian = np.random.normal(n*p, n*p*(1-p), size=size)
# Add the continuity correction to sample at the midpoint of each integral bin.
gaussian += 0.5
if size is not None:
binomial = gaussian.astype(np.int64)
else:
#
浏览 6提问于2016-02-02得票数 2
回答已采纳
1回答
我试图在BW2中填充对象，但我得到了一个完全白色的图像。为什么会发生这种情况？
clc;
clear;
grayImage = rgb2gray(imread('https://i.stack.imgur.com/afCDL.jpg'));
gaussian1 = fspecial('Gaussian', 15, 10);
gaussian2 = fspecial('Gaussian', 15, 200);
dog = gaussian1 - gaussian2;
dogFilterImage = conv2(double(grayImage
浏览 7提问于2017-01-29得票数 1
回答已采纳
我有一个模型，它在训练过程中的输入大小是可变的，以便于推广。
为了量化，我必须修正输入大小，所以我只需要用固定的输入大小重新创建模型，然后复制所有的权重和偏差，然后保存模型。然而，由于某种原因，型号的大小大致分为几个部分。注意，这是在量化或其他任何东西之前，并且参数保持不变。
以下是两个示范摘要：
型号1= 4.6MB
old_model.summary(line_length=110)
____________________________________________________________________________________________________
浏览 1提问于2019-03-13得票数 2
回答已采纳
2回答
我正在实现用于分类识别分析的“气泡算法”。在此算法中，我们必须使用高斯圆来制作掩模。圆的中心将是1(255)，并且它将在半径为0的范围内减小。我有问题，当我把圆圈在彼此之间创建蒙版时，它会在圆圈之间放一条黑线，我无法删除它。这是我的代码：
def make_gaussian(circle_center, Gaussian_base):
for t in range(circle_center[1] - radius, circle_center[1] + radius):
for tt in range(circle_center[0] - radius, circle_center[
浏览 1提问于2020-12-11得票数 0
2回答
我最近开始使用Julia来加速以前用Python编写的一些代码。我以前只有使用Python的经验，所以这是我第一次关心性能，在遍历一个结构数组时，我发现了一些奇怪的行为。我正在定义一个新的结构高斯，它表示一个2d高斯函数和一个函数强度()，它计算给定位置的函数的振幅：
struct Gaussian{T<:Float32}
x0::T
y0::T
A::T
a::T
b::T
c::T
end
function intensity(
model::Gaussian,
x::Float32,
y::Float32
浏览 5提问于2022-06-29得票数 1
回答已采纳
我有这个gaussian2 m文件：
function [g] = gaussian2(t, tau)
g = (1/sqrt(2*pi)) .*(1/tau) .* exp(- (t .^2) / (2 * tau^2));
end
我想让另一个函数m文件为RMStau：
sqrt( integral( @(t) t.^2 .* gaussian2(t,tau), -inf, inf ) / integral( @(t) gaussian2(t,tau), -inf, inf ))
如何使RMStau m文件成为gaussian2输入的函数？谢谢
浏览 2提问于2016-04-08得票数 0
有人知道scipy.ndimage.gaussian_gradient_magnitude 和scipy.ndimage.gaussian_filter function有什么区别吗？
我认为gaussian_gradient_magnitude(img, 1)和gaussian_filter(img, sigma= 1, order = 1)应该返回相同的结果，但事实并非如此。
谢谢。
浏览 14提问于2019-11-29得票数 0
1回答
我试图在Ubuntu21.10上应用gtk+3.0补丁，但失败了。但是更重要的是，即使没有应用补丁，我也无法构建软件包。
以下是几个步骤：
apt-get source "gtk+3.0"
sudo apt-get build-dep "gtk+3.0"
cd gtk+3.0-3.24.40
dpkg-buildpackage -rfakeroot -b -us -uc
但以下几个方面都失败了：
make check-local
make[6]: Entering directory '/home/gaussian/owninstalls/sway-u
浏览 0提问于2021-11-26得票数 1
2回答
我试着找到这个问题的答案，也许它很简单，这就是为什么我不能。
如果我制作了一个高斯函数，并且我想用Matplotlib.pyplot.plot绘制它，那么我如何使用浮点值来绘制它呢？即来自-20<=x<=20的值，增量为0.1
import matplotlib.pyplot as plt
import math
from math import exp
import numpy
#Parameters for the Gaussian
A=1
c=10
t=0
a=1
x=[]
p=-20.
while p<=20:
x.append(p)
p+=0.
浏览 0提问于2013-11-03得票数 0
回答已采纳
我在试着kaggle Bayesian Hyperparam Optimization of RF内核。而且我不能导入sklearn.gaussian_process.GaussianProcess。请帮帮这个可怜的新手吧。 from sklearn.gaussian_process import GaussianProcess as GP 错误： Traceback (most recent call last):
File "C:/Users/Develop/PycharmProjects/reinforcement recommandation system/BNP/bay
浏览 232提问于2019-09-03得票数 3
回答已采纳
1回答
我正在尝试编写一个方法来生成高斯整数的高斯因子序列-高斯整数要么是正常整数，要么是复数g = a + bi，其中a和b都是整数，高斯整数g的高斯因子是高斯整数d，使得g / d也是高斯整数。
我有下面的代码。
def is_gaussian_integer(c):
"""
Checks whether a given real or complex number is a Gaussian integer,
i.e. a complex number g = a + bi such that a and b are integer
浏览 23提问于2017-01-16得票数 5
回答已采纳
2回答
我试图拟合一些高斯人，对于初始参数我已经有了一个很好的想法(在这个例子中，我正在生成分布，所以我应该总是能够适应这些分布)。然而，对于两个高斯人，我似乎不知道如何强迫平均值为0。有可能吗？m.means_ = ...不起作用。
from sklearn import mixture
import numpy as np
import math
import matplotlib.pyplot as plt
from scipy import stats
a = np.random.normal(0, 0.2, 500)
b = np.random.normal(0, 2, 800)
obs
浏览 8提问于2017-12-20得票数 2
1回答
我正在努力实现的科学工具包-学习实现KDE的小输入范围。下面的代码可以工作。将除数变量增加到100，KDE挣扎：
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
sns.set_style('whitegrid')
from sklearn.neighbors import KernelDensity
# make data:
np.random.seed(0)
divisor = 1
gaussian1 = (3 * np.random.randn(1700))/divisor
浏览 1提问于2017-07-11得票数 2
回答已采纳
我正在使用ggplot2有一些带有图例的密度图，但我不能将图例添加到我的最终结果中。
m <- ggplot(chickwts, aes(x = weight))
m + geom_density(kernel = "gaussian", adjust = .3, colour="green", size=1)+
geom_density(kernel = "gaussian", adjust = 1, colour="red", size=1)+
geom_density(kernel = "gaussia
浏览 2提问于2013-10-04得票数 0
我正在为以下命令编写CMake文件：
g++ --std=c++11 -DLOG2VLEN=7 -DRUN=1 -DASSEMBLY=1 -DBITCODE=1 -I
/home/user/Qualcomm/HALIDE_Tools/2.2.06/Halide/include -fno-rtti -O3 -g
halide_gaussian_hexagon_generate.cpp /home/user/Qualcomm/HALIDE_Tools/2.2.06/Halide/tools/GenGen.cpp
-L/home/user/Qualcomm/HALIDE_Tools/2.2.0
浏览 12提问于2021-08-04得票数 0
3回答
我有下面的Matlab代码来处理两幅图像，灰度图像和RGB图像。重点是在这两个图像上应用Average、Gaussian和Laplacian过滤器。
%%Clear
clear
clc
%%Reading images
gray=imread('cameraman.tif')
[gray_table gray_map]=gray2ind(gray,256)
rgb=imread('peppers.png')
[rgb_table rgb_map]=rgb2ind(rgb,256)
%%Creating filters
average=fspecial('
浏览 5提问于2014-10-24得票数 1
回答已采纳
我试图用高斯二维卷积模糊矩阵。但在边界元素上我得到了尖锐的跃迁。
下面是我正在运行的一段代码：
// create 1D Kernel
void createGaussianKerenel_1D() {
unsigned kernelSize = 2 * kernelRad_ + 1;
gaussian1Dkernel_ = vector<double>(kernelSize);
double sigma = (double)kernelRad_;
double sum = 0.0;
for(unsigned i = 0; i < k
浏览 1提问于2015-05-22得票数 1
回答已采纳
我用的是这个密码，
import matplotlib.pyplot as plt
from numpy import exp, loadtxt, pi, sqrt
from lmfit import Model
data = loadtxt('model1d_gauss.dat')
x = data[:, 0]
y = data[:, 1] + 0.25*x - 1.0
def gaussian(x, amp, cen, wid):
"""1-d gaussian: gaussian(x, amp, cen, wid)""
浏览 2提问于2020-11-28得票数 0
1回答
我试图使用以下代码向一些图像中添加高斯噪声
import numpy as np
import cv2
import glob
mean = 0
var = 10
sigma = var ** 0.5
gaussian = np.random.normal(mean, sigma, (224, 224))
for image in glob.glob('/home/aub/myflower/flower_photos/dandelion/*.jpg'):
img = cv2.imread(image)
noisy_image = np.zeros(im
浏览 0提问于2019-03-20得票数 2
回答已采纳
1回答
我试图用双高斯曲线拟合一些数据。数据看起来几乎完全是高斯的，但是尽管我尝试，我无法得到比某个形状更合适的结果，不管我输入的最初猜测是什么。我试过使用下面列出的两个高斯方程，但两者都不太合适。总的来说，我希望它在连续体(没有“翅膀”)上更平顺，如果可能的话，它会更平滑，更接近实际的形状。
由于后续分析的性质，拟合需要是双高斯的，因为我需要拟合参数，因此不能考虑其他的拟合方法。这些数据可以在这里找到：
下面是一个代码的例子，我一直用它来拟合数据，以及输出图。
import matplotlib.pyplot as plt
import numpy as np
import scipy.stats
浏览 1提问于2020-10-02得票数 0
1回答
我是c编程和图像处理领域的新手。我试图用高斯环绕函数在空间域上用PPM图像编码单尺度的retinex。我使用的是标准的retinex公式如下：
Ri(x，y) = logIi(x，y) - logF(x，y)*Ii(x，y)
它的子指标i代表第一个光谱带，c是constant.Ri (x，y) i是Retinex输出，Ii(x，y)是I波段的输入图像分布。
我选择了3x3高斯掩码如下：
F(x，y)=e^(-2/c^2) e^(-1/c^2) e^(-2/c^2);e^(-1/c^2) 1 e^(-1/c^2);e^(-2/c^2) e^(-1/c^2) e^(-2/c^2)
浏览 2提问于2016-04-12得票数 0
1回答
我正在尝试从运行代码。此代码是为gpflow版本0.4.0编写的。我想更新gpflow 2.5.2版本的代码。我更改了下面一行
gp = gpflow.gpr.GPR(np.empty((0, 2), dtype=safe_learning.config.np_dtype),
np.empty((0, 1), dtype=safe_learning.config.np_dtype),
kernel,
mean_function=mean_function)
至
gp = gp
浏览 14提问于2022-09-13得票数 1
可能的重复：
在决定是否应该这样做时，应该使用什么标准：
<Blur Type="Gaussian", Amount=5></Blur>
或
<Blur>Gaussian, 5</Blur>
如果是#1，你会在括号里放什么？
另外，我不知道括号中的值是什么(Gaussian, 5)：inside properties
浏览 6提问于2009-09-21得票数 1
1回答
高斯拟合的返回值
、、、、
我正在尝试将高斯函数拟合到一些数据，但我想返回'center‘和'width’值并打印它们。但是，代码返回一个ValueError，因此没有绘图。从高斯函数的返回行中删除'center‘和'width’可以解决这个问题，但我不知道如何从拟合中获得这些值(以及如何让代码在拟合后将它们打印为输出)。有什么建议吗？ 代码错误： import matplotlib.pyplot as plt
import pathlib
import os
from scipy.optimize import curve_fit
#Data
data = np.loadtxt(&
浏览 35提问于2021-08-05得票数 1
我的职能如下：
tidy_multi_normal <- function(.n = 50, .mean = c(-1,0,1), .sd = c(1), .num_sims = 2){
n <- as.numeric(.n)
mu <- as.numeric(.mean)
std <- as.numeric(.sd)
num_sims <- as.numeric(.num_sims)
x <- seq(1, num_sims, 1)
ps <- seq(-n, n-1, 2
浏览 3提问于2022-01-21得票数 0
2回答
我有一个针对数据集运行的Gaussian naive bayes算法。我需要的是获取目标类的特性重要性(特性的影响性)。
这是我的密码：
from sklearn.naive_bayes import GaussianNB
from sklearn.metrics import accuracy_score
from sklearn.model_selection import train_test_split
X_train, X_test, y_train, y_test = train_test_split(inputs, target, test_size=0.2)
gaussia
浏览 1提问于2020-07-16得票数 4
1回答
我正在使用numpy的random.normal例程来创建一个具有给定均值和标准差的高斯。
array_a = an array of len(100)
gaussian = np.random.normal(loc=array_a,scale=0.1,size=len(2*array_a))
因此，我希望gaussian有一个mean=array_a和stddev=0.1，gaussian数组的大小是array_a的2倍。然而，上面的array_a返回一个与 !相同大小的数组。
如何使len(gaussian)成为2倍 len(array_a)与给定的mean和standard devia
浏览 2提问于2016-02-25得票数 0
回答已采纳
1回答
单程高斯模糊
、、
有可能实现碎片着色做一次高斯模糊？我发现了很多两通模糊(高斯模糊和方框模糊)的实现：
http://callumhay.blogspot.com/2010/09/gaussian-blur-shader-glsl.html
http://www.gamerendering.com/2008/10/11/gaussian-blur-filter-shader/
http://www.geeks3d.com/20100909/shader-library-gaussian-blur-post-processing-filter-in-glsl/
诸若此类。
我一直在考虑将高斯模糊化为卷积(实际上，
浏览 0提问于2012-04-02得票数 19
回答已采纳
我使用的是glm函数，它接受一个名为的参数，如下所示：
glm(..., family = gaussian("identity"))
我想为系列创建一个可能的选项数组，并循环它们。我最初的想法是使用引用，如下所示：
fam.array <- c( quote(gaussian("identity")), quote(gaussian("inverse"),
quote(gaussian("log")))
for(i in 1:3) glm(..., family = fam.array[i])
浏览 2提问于2015-03-05得票数 0
1回答
我正在尝试使用AutoSklearn和特定的算法列表。
Logistic Regression
Random Forest
Gaussian NB
SVC
ADA
MLP
我知道我可以用这些参数。
mdl = autosklearn.classification.AutoSklearn2Classifier(
include = {
'classifier': ["random_forest", "gaussian_nb", "libsvm_svc", "adaboost"],
浏览 11提问于2022-01-18得票数 0
回答已采纳
2回答
我试图创建一个16位灰度环模拟，由于某种原因，它只是不工作。
让我解释一下，在开始的时候，我用8位格式写的，然后我意识到我需要16位的格式。我为8位格式编写的代码工作得很好，如下所示：
from PIL import Image, ImageFilter
from scipy.ndimage import convolve, gaussian_filter
import NumPy as np
def theta(x, y, phi):
if np.angle(x - y*1j, deg=True) - phi > 180:
return 1*(np.angle
浏览 2提问于2022-01-17得票数 2
1回答
我有一个名为gaussian_array的数组，它由一系列数字组成，一旦绘制出来，就会形成一个近似的高斯数组。我需要理解这个高斯的\sigma，但我不允许使用任何形式的fit。到目前为止，我所尝试的是计算高斯的峰值，它是由阵列的第一个元素(高斯以原点为中心) gaussian_array[0]给出的，然后不知何故我认为使用\sigma和FWHM之间的众所周知的关系可能会很有用。但是，我不知道如何在python中实现这一点。我认为写一些像这样的东西可能会很有用
for i in range(len(gaussian_array)):
if gaussian[i] = FWHM:
浏览 1提问于2015-06-12得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
