Final working version of SDSS mockup

Author: Ezequiel92

Manifest.toml

@@ -2,7 +2,7 @@
 
 julia_version = "1.12.4"
 manifest_format = "2.0"
-project_hash = "202ccab0baded9db29714373162e784d21a62543"
+project_hash = "1560840b2eb579fcf2e15ac9956ef29a6618bfd4"
 
 [[deps.AbstractFFTs]]
 deps = ["LinearAlgebra"]
@@ -99,6 +99,17 @@ version = "7.22.0"
 uuid = "56f22d72-fd6d-98f1-02f0-08ddc0907c33"
 version = "1.11.0"
 
+[[deps.AstroAngles]]
+git-tree-sha1 = "bc188d9a6507511e7360444d54ed57d0a9d6cf91"
+uuid = "5c4adb95-c1fc-4c53-b4ea-2a94080c53d2"
+version = "0.2.0"
+
+[[deps.AstroImages]]
+deps = ["AbstractFFTs", "AstroAngles", "ColorSchemes", "DimensionalData", "FITSIO", "FileIO", "ImageAxes", "ImageBase", "ImageIO", "ImageShow", "MappedArrays", "PlotUtils", "PrecompileTools", "Printf", "RecipesBase", "Statistics", "Tables", "UUIDs", "WCS"]
+git-tree-sha1 = "b036ab0541311fee6f278e1469b7839cc3af8f19"
+uuid = "fe3fc30c-9b16-11e9-1c73-17dabf39f4ad"
+version = "0.5.1"
+
 [[deps.Automa]]
 deps = ["PrecompileTools", "SIMD", "TranscodingStreams"]
 git-tree-sha1 = "a8f503e8e1a5f583fbef15a8440c8c7e32185df2"
@@ -143,6 +154,18 @@ git-tree-sha1 = "389ad5c84de1ae7cf0e28e381131c98ea87d54fc"
 uuid = "fa961155-64e5-5f13-b03f-caf6b980ea82"
 version = "0.5.0"
 
+[[deps.CFITSIO]]
+deps = ["CFITSIO_jll"]
+git-tree-sha1 = "8c6b984c3928736d455eb53a6adf881457825269"
+uuid = "3b1b4be9-1499-4b22-8d78-7db3344d1961"
+version = "1.7.2"
+
+[[deps.CFITSIO_jll]]
+deps = ["Artifacts", "Bzip2_jll", "JLLWrappers", "LibCURL_jll", "Libdl", "Zlib_jll"]
+git-tree-sha1 = "15e80be798d7711411f4ac4273144cdb2a89eb2f"
+uuid = "b3e40c51-02ae-5482-8a39-3ace5868dcf4"
+version = "4.6.2+0"
+
 [[deps.CPUSummary]]
 deps = ["CpuId", "IfElse", "PrecompileTools", "Preferences", "Static"]
 git-tree-sha1 = "f3a21d7fc84ba618a779d1ed2fcca2e682865bab"
@@ -386,6 +409,36 @@ git-tree-sha1 = "c55f5a9fd67bdbc8e089b5a3111fe4292986a8e8"
 uuid = "927a84f5-c5f4-47a5-9785-b46e178433df"
 version = "1.6.6"
 
+[[deps.DimensionalData]]
+deps = ["Adapt", "ArrayInterface", "ConstructionBase", "DataAPI", "Dates", "Extents", "Interfaces", "IntervalSets", "InvertedIndices", "IteratorInterfaceExtensions", "LinearAlgebra", "OrderedCollections", "PrecompileTools", "Random", "RecipesBase", "Statistics", "TableTraits", "Tables"]
+git-tree-sha1 = "147961441e5cb35da0af404aa4684d2e74ec68eb"
+uuid = "0703355e-b756-11e9-17c0-8b28908087d0"
+version = "0.29.25"
+
+    [deps.DimensionalData.extensions]
+    DimensionalDataAbstractFFTsExt = "AbstractFFTs"
+    DimensionalDataAlgebraOfGraphicsExt = "AlgebraOfGraphics"
+    DimensionalDataCategoricalArraysExt = "CategoricalArrays"
+    DimensionalDataChainRulesCoreExt = "ChainRulesCore"
+    DimensionalDataDiskArraysExt = "DiskArrays"
+    DimensionalDataMakieExt = "Makie"
+    DimensionalDataNearestNeighborsExt = "NearestNeighbors"
+    DimensionalDataPythonCallExt = "PythonCall"
+    DimensionalDataSparseArraysExt = "SparseArrays"
+    DimensionalDataStatsBaseExt = "StatsBase"
+
+    [deps.DimensionalData.weakdeps]
+    AbstractFFTs = "621f4979-c628-5d54-868e-fcf4e3e8185c"
+    AlgebraOfGraphics = "cbdf2221-f076-402e-a563-3d30da359d67"
+    CategoricalArrays = "324d7699-5711-5eae-9e2f-1d82baa6b597"
+    ChainRulesCore = "d360d2e6-b24c-11e9-a2a3-2a2ae2dbcce4"
+    DiskArrays = "3c3547ce-8d99-4f5e-a174-61eb10b00ae3"
+    Makie = "ee78f7c6-11fb-53f2-987a-cfe4a2b5a57a"
+    NearestNeighbors = "b8a86587-4115-5ab1-83bc-aa920d37bbce"
+    PythonCall = "6099a3de-0909-46bc-b1f4-468b9a2dfc0d"
+    SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
+    StatsBase = "2913bbd2-ae8a-5f71-8c99-4fb6c76f3a91"
+
 [[deps.Distances]]
 deps = ["LinearAlgebra", "Statistics", "StatsAPI"]
 git-tree-sha1 = "c7e3a542b999843086e2f29dac96a618c105be1d"
@@ -480,6 +533,12 @@ git-tree-sha1 = "6d6219a004b8cf1e0b4dbe27a2860b8e04eba0be"
 uuid = "f5851436-0d7a-5f13-b9de-f02708fd171a"
 version = "3.3.11+0"
 
+[[deps.FITSIO]]
+deps = ["CFITSIO", "Printf", "Reexport", "Tables"]
+git-tree-sha1 = "f57de3f533590c785210893030736dc11c4a4afb"
+uuid = "525bcba6-941b-5504-bd06-fd0dc1a4d2eb"
+version = "0.17.5"
+
 [[deps.FileIO]]
 deps = ["Pkg", "Requires", "UUIDs"]
 git-tree-sha1 = "d60eb76f37d7e5a40cc2e7c36974d864b82dc802"
@@ -882,6 +941,11 @@ deps = ["Markdown"]
 uuid = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
 version = "1.11.0"
 
+[[deps.Interfaces]]
+git-tree-sha1 = "331ff37738aea1a3cf841ddf085442f31b84324f"
+uuid = "85a1e053-f937-4924-92a5-1367d23b7b87"
+version = "0.3.2"
+
 [[deps.Interpolations]]
 deps = ["Adapt", "AxisAlgorithms", "ChainRulesCore", "LinearAlgebra", "OffsetArrays", "Random", "Ratios", "SharedArrays", "SparseArrays", "StaticArrays", "WoodburyMatrices"]
 git-tree-sha1 = "65d505fa4c0d7072990d659ef3fc086eb6da8208"
@@ -2045,6 +2109,18 @@ git-tree-sha1 = "d1d9a935a26c475ebffd54e9c7ad11627c43ea85"
 uuid = "3d5dd08c-fd9d-11e8-17fa-ed2836048c2f"
 version = "0.21.72"
 
+[[deps.WCS]]
+deps = ["ConstructionBase", "WCS_jll"]
+git-tree-sha1 = "c12065744b66adfed32d24c2a13a3053cc235ea7"
+uuid = "15f3aee2-9e10-537f-b834-a6fb8bdb944d"
+version = "0.6.3"
+
+[[deps.WCS_jll]]
+deps = ["Artifacts", "JLLWrappers", "Libdl", "Pkg"]
+git-tree-sha1 = "947bfa11fcd65dac9e9b2e963504fba6b4971d31"
+uuid = "550c8279-ae0e-5d1b-948f-937f2608a23e"
+version = "7.7.0+0"
+
 [[deps.WeakRefStrings]]
 deps = ["DataAPI", "InlineStrings", "Parsers"]
 git-tree-sha1 = "b1be2855ed9ed8eac54e5caff2afcdb442d52c23"

Project.toml

@@ -3,6 +3,7 @@ uuid = "990b1489-e520-4e93-b9dc-df8612ec7e1d"
 authors = ["Ezequiel Lozano <elozano@df.uba.ar>"]
 
 [deps]
+AstroImages = "fe3fc30c-9b16-11e9-1c73-17dabf39f4ad"
 CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
 CairoMakie = "13f3f980-e62b-5c42-98c6-ff1f3baf88f0"
 ColorSchemes = "35d6a980-a343-548e-a6ea-1d62b119f2f4"
@@ -33,6 +34,7 @@ UnitfulAstro = "6112ee07-acf9-5e0f-b108-d242c714bf9f"
 WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
 
 [compat]
+AstroImages = "0.5.1"
 Distances = "0.10.12"
 Interpolations = "0.16.2"
 StaticArrays = "1.9.15"

src/GalaxyInspector.jl

@@ -14,7 +14,8 @@
 
 module GalaxyInspector
 
-using CSV,
+using AstroImages,
+    CSV,
     CairoMakie,
     ColorSchemes,
     Colors,
@@ -131,6 +132,7 @@ export stellarDensityMaps
 export gasDensityMaps
 export gasFractionsEvolution
 export evolutionVideo
+export SDSSMockup
 export snapshotReport
 export simulationReport
 export quantityReport

src/analysis/compute_quantities/masses.jl

@@ -1375,7 +1375,7 @@ function computeNumberDensity(data_dict::Dict, component::Symbol)::Vector{<:Numb
 end
 
 """
-    computeNumber(data_dict::Dict, component::Symbol)::Vector{Int64}
+    computeNumber(data_dict::Dict, component::Symbol)::Vector{Float64}
 
 Compute the number of a given `component`.
 
@@ -1414,14 +1414,14 @@ Compute the number of a given `component`.
 
 L. Blitz et al. (2006). *The Role of Pressure in GMC Formation II: The H2-Pressure Relation*. The Astrophysical Journal, **650(2)**, 933. [doi:10.1086/505417](https://doi.org/10.1086/505417)
 """
-function computeNumber(data_dict::Dict, component::Symbol)::Vector{Int64}
+function computeNumber(data_dict::Dict, component::Symbol)::Vector{Float64}
 
     if component ∉ COMPONENTS
         throw(ArgumentError("computeNumber: `component` can only be one of the elements of \
         `COMPONENTS` (see `./src/constants/globals.jl`), but I got :$(component)"))
     end
 
-    if component ∈ [:stellar, :dark_matter, :black_hole. :gas]
+    if component ∈ [:stellar, :dark_matter, :black_hole, :gas]
 
         ############################################################################################
         # Number as the amount of particles
@@ -1445,7 +1445,7 @@ function computeNumber(data_dict::Dict, component::Symbol)::Vector{Int64}
 
         M = computeMass(data_dict, component)
 
-        N = M / (2.0 * Unitful.mp)
+        N = ustrip.(Unitful.NoUnits, M / (2.0 * Unitful.mp))
 
     else
 
@@ -1455,7 +1455,7 @@ function computeNumber(data_dict::Dict, component::Symbol)::Vector{Int64}
 
         M = computeMass(data_dict, component)
 
-        N = M / Unitful.mp
+        N = ustrip.(Unitful.NoUnits, M / Unitful.mp)
 
     end
 
@@ -1467,7 +1467,7 @@ function computeNumber(data_dict::Dict, component::Symbol)::Vector{Int64}
 
     end
 
-    return ustrip.(Unitful.NoUnits, N)
+    return N
 
 end
 
@@ -2086,7 +2086,7 @@ function quantity3DProjection(
 
         end
 
-        voxel_values = histogram3D(positions, qty_values, grid; empty_nan, tall=true)
+        voxel_values = histogram3D(positions, qty_values, grid; empty_nan)
 
     end