Loading [MathJax]/extensions/tex2jax.js
LSST Applications g04a91732dc+cc8870d3f5,g07dc498a13+5aa0b8792f,g0fba68d861+80045be308,g1409bbee79+5aa0b8792f,g1a7e361dbc+5aa0b8792f,g1fd858c14a+f64bc332a9,g208c678f98+1ae86710ed,g35bb328faa+fcb1d3bbc8,g4d2262a081+47ad8a29a8,g4d39ba7253+9633a327c1,g4e0f332c67+5d362be553,g53246c7159+fcb1d3bbc8,g60b5630c4e+9633a327c1,g668ecb457e+25d63fd678,g78460c75b0+2f9a1b4bcd,g786e29fd12+cf7ec2a62a,g7b71ed6315+fcb1d3bbc8,g8852436030+8b64ca622a,g89139ef638+5aa0b8792f,g89e1512fd8+04325574d3,g8d6b6b353c+9633a327c1,g9125e01d80+fcb1d3bbc8,g989de1cb63+5aa0b8792f,g9f33ca652e+b196626af7,ga9baa6287d+9633a327c1,gaaedd4e678+5aa0b8792f,gabe3b4be73+1e0a283bba,gb1101e3267+71e32094df,gb58c049af0+f03b321e39,gb90eeb9370+2807b1ad02,gcf25f946ba+8b64ca622a,gd315a588df+a39986a76f,gd6cbbdb0b4+c8606af20c,gd9a9a58781+fcb1d3bbc8,gde0f65d7ad+4e42d81ab7,ge278dab8ac+932305ba37,ge82c20c137+76d20ab76d,gfe73954cf8+a1301e4c20,w.2025.11
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
test_gaussian.py
Go to the documentation of this file.
1# This file is part of gauss2d.
2#
3# Developed for the LSST Data Management System.
4# This product includes software developed by the LSST Project
5# (https://www.lsst.org).
6# See the COPYRIGHT file at the top-level directory of this distribution
7# for details of code ownership.
8#
9# This program is free software: you can redistribute it and/or modify
10# it under the terms of the GNU General Public License as published by
11# the Free Software Foundation, either version 3 of the License, or
12# (at your option) any later version.
13#
14# This program is distributed in the hope that it will be useful,
15# but WITHOUT ANY WARRANTY; without even the implied warranty of
16# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17# GNU General Public License for more details.
18#
19# You should have received a copy of the GNU General Public License
20# along with this program. If not, see <https://www.gnu.org/licenses/>.
21
22import pytest
23
24import lsst.gauss2d as g2d
25
26prefix_namespace = "lsst.gauss2d."
27
28
30 gauss0 = g2d.Gaussian()
31 assert gauss0 == g2d.Gaussian()
32 centroid = g2d.Centroid()
33 ellipse = g2d.Ellipse()
34 gauss1 = g2d.Gaussian(centroid, ellipse)
35 gauss2 = g2d.Gaussian(centroid, ellipse, g2d.GaussianIntegralValue(2))
36 gauss2.integral_value = gauss1.integral_value
37 assert gauss1 == gauss2
38
39 centroid.x = 1
40 assert gauss1.centroid is centroid
41 kwargs = {
42 "str_values_cen": "CentroidValues(x=1.000000e+00, y=0.000000e+00)",
43 "str_values_ell": "EllipseValues(sigma_x=0.000000e+00, sigma_y=0.000000e+00, rho=0.000000e+00)",
44 "str_integral": "GaussianIntegralValue(value=1.000000e+00)",
45 }
46 str_gauss_format = (
47 "{prefix_namespace}Gaussian(centroid={prefix_namespace}Centroid(data={prefix_namespace}"
48 "{str_values_cen}), ellipse={prefix_namespace}Ellipse(data={prefix_namespace}{str_values_ell}), "
49 "integral={prefix_namespace}{str_integral})"
50 )
51 assert str(gauss1) == str_gauss_format.format(prefix_namespace="", **kwargs)
52 assert repr(gauss1) == str_gauss_format.format(prefix_namespace=prefix_namespace, **kwargs)
53
54
56 gauss1 = g2d.Gaussian()
57 gs = g2d.Gaussians([gauss1, gauss1])
58 assert gs.at(0) == gs.at(1)
59 with pytest.raises(IndexError):
60 gs.at(2)
61 for idx in range(len(gs)):
62 assert gs.at(idx) == gauss1
63
64 str_g1 = str(gauss1)
65 str_gs = f"Gaussians(data=[{str_g1}, {str_g1}])"
66 assert str(gs) == str_gs
67
68
69@pytest.fixture(scope="module")
71 gauss1 = g2d.Gaussian(centroid=None, ellipse=g2d.Ellipse(3, 6, 0))
72 gauss2 = g2d.Gaussian(centroid=None, ellipse=g2d.Ellipse(4, 8, 0))
73 return gauss1, gauss2
74
75
76@pytest.fixture(scope="module")
77def convolved_gaussian(gaussians):
78 gauss_conv = g2d.ConvolvedGaussian(gaussians[0], gaussians[1])
79 return gauss_conv
80
81
82def test_ConvolvedGaussian(convolved_gaussian, gaussians):
83 assert convolved_gaussian.source is gaussians[0]
84 assert convolved_gaussian.kernel is gaussians[1]
85
86
87def test_ConvolvedGaussians(convolved_gaussian, gaussians):
88 conv_gaussian2 = g2d.ConvolvedGaussian(gaussians[1], gaussians[0])
89 conv_list = [convolved_gaussian, conv_gaussian2]
90 n_conv = len(conv_list)
91 conv_gaussians = g2d.ConvolvedGaussians(conv_list)
92 assert len(conv_gaussians) == n_conv
93 assert conv_gaussians.size == n_conv
94 assert conv_gaussians[0] is conv_list[0]
95 assert conv_gaussians.at(1) is conv_list[1]
96 with pytest.raises(IndexError):
97 conv_gaussians[n_conv]
A 2D coordinate representing the center of a plane figure.
Definition centroid.h:114
A convolution of a Gaussian source and kernel.
Definition gaussian.h:220
A collection of ConvolvedGaussian objects.
Definition gaussian.h:249
An Ellipse with sigma_x, sigma_y, and rho values.
Definition ellipse.h:283
A 2D Gaussian with a Centroid, Ellipse, and integral.
Definition gaussian.h:99
A GaussianIntegral storing a float value.
Definition gaussian.h:72
An array of Gaussian objects.
Definition gaussian.h:175
convolved_gaussian(gaussians)
test_ConvolvedGaussian(convolved_gaussian, gaussians)
test_ConvolvedGaussians(convolved_gaussian, gaussians)