標簽[healpy] - 堆棧內存溢出

[英]In healpy how to split mask/map into 2 parts?

我有一個適合文件的面具/地圖，可在此處下載，我可以使用以下代碼使用 healpy plot 這個面具：它產生這個數字：我想做的是將面具分成 2 個新面具 - 一個應該只包含 A 部分，另一個應該只包含 B 部分。我不知道該怎么做，我想對所有點做一些條件表達式which longitude + ...

使用不同於 l_max 的 alm 大小運行 healpy.sphtfunc.alm2map 時出錯

[英]Error running healpy.sphtfunc.alm2map with alm size different from l_max required

我運行 healpy.sphtfunc.alm2map，將 alm 作為 l_max = 256 的數組傳遞，並詢問 output 和 l_max = 191 的 map，但似乎 function 沒有正確接受新的 l_max。這是我正在做的示例代碼，其中包含人工生成的施舍：我收到此錯誤：誰能 ...

plot healpix binning我應該怎么處理

[英]How should I plot healpix binning

我正在繪制等面積分箱的天球。我的第一個意圖是根據我從大學一年級微積分 class 中學到的基本微積分知識進行分箱。以下是我所做的。有了這個我會得到這樣的東西在繪制我制作的 plot 后，我了解到有一種名為 healpix 的分箱方法，它可以提供不同形狀的等面積分箱。我還了解到有一個名為 he ...

Healpy 中的正交投影

[英]Orthographic Projection in Healpy

我正在嘗試使用 healpy 地圖，但遇到了正交投影問題。我想在 healpy mollview 投影圖中分配“N”個像素，假設值為 100。但是，當我對地圖進行正交投影時，並非此處所有“N”個像素的值為 100。可以有人幫我嗎？ healpy中正交投影的像素面積是多少？找到下面的偽代碼。根 ...

為什么 python function 即使不應該並行工作？

[英]Why does a python function work in parallel even if it should not?

我正在使用healpy package 運行這段代碼。我沒有使用多處理，我需要它在單核上運行。它工作了一段時間，但是，當我現在運行它時，function healpy.projector.GnomonicProj.projmap占用了所有可用的內核。這是有罪的代碼塊：有人知道為什么這個 fu ...

從 3 個 healpy 貼圖創建 RGB 圖像

[英]Creating a RGB image out of 3 healpy maps

我有三個 healpy 映射（對應於不同頻率的三個 Stokes 參數），我試圖將它們 plot 放在同一個 map 中，以突出差異。我想做一些類似於第一個 plot 中的大 map 中所做的事情： https://sci.esa.int/web/planck/-/56323-planck-leg ...

將非中心高斯與 healpy map 卷積的快速方法？

[英]Fast way to convolve a non-centered gaussian with a healpy map?

我有一個healpix map（nside = 256），我試圖將它與一個不居中的高斯卷積，這是模擬“望遠鏡指向”中的偏移。有沒有一種簡單的方法可以做到這一點？我知道 `healpy.smoothing' 適用於通常的高斯函數，但在這種情況下找不到任何東西。 ...

將徑向亮度分布轉換為二維強度分布

[英]Converting radial brightness profile to 2D intensity distribution

我有 2 列表，將光源的光強度分布表示為其 angular 半徑的 function。這些表是將本圖中第一個plot和第三個plot數字化得到的。 x-aix 是兩個源的 angular 半徑。 y 軸是源的強度。現在我想將一維分布轉換為二維分布，如經度和緯度的 function。另外，我 ...

HealPix 中不同的 Nside 值對應不同的 a_lm 系數值

[英]Different values of a_lm coefficients for different values of Nside in HealPix

我正在使用 HealPix 來計算粒子通過具有一定半徑的球體的流出率，並且我正在嘗試確定球諧函數的系數（基本上用球諧函數擬合一些數據）。基本步驟是：但是，如果我為相同的 lmax 值更改 map 的分辨率（更改 nside），我會得到不同的系數值（不同的 a_lm 值）。我嘗試從系數中復制初始 ...

如何從 Healpy 中的完整天空地圖獲取方形地圖（較小的補丁）

[英]How to obtain square maps (smaller patches) from a full sky map in Healpy

我有一個完整的天空收斂圖。我想從這張地圖上剪下更小的方形補丁。我想知道這怎么可能。如果我閱讀地圖並繪制它，它是這樣的：此后我不知道如何切割較小的方形補丁。 ...

健康中的 FITS 數據處理

[英]FITS data handling in healpy

我正在嘗試在 python 中映射一個 FITS 數據文件，並開始遵循他們網站上提供的健康文檔。但是我無法讀取我必須為其創建地圖的 FITS 文件，因為它給了我一些我無法理解的錯誤。我將提供下面的代碼 - 我提供的錯誤是 - 有人可以告訴我我做錯了什么嗎？ ...

Healpy：在 hp.ma() 中將掩碼與 map 一起旋轉與單獨旋轉產生不同的結果？

[英]Healpy: Rotate a mask together with the map in hp.ma() vs. separately produce different results?

我想更好地理解 rotate_map_pixel() 在這種特殊情況下的行為。下面我嘗試提供一個最小的玩具示例：這會產生略有不同的 map.mask 和 mask。 “不同”是指和要么都產生不同的空切。我猜這是由於插值程序？如果我想使用面膜作選擇用途，我應該怎么做？ ...

***檢測到堆棧粉碎***：使用 HEALPix C 子例程終止中止（核心轉儲）

[英]*** stack smashing detected ***: terminated Aborted (core dumped) using HEALPix C subroutines

我正在嘗試將 python 代碼轉換為 C 代碼。所以我試圖翻譯那行：import healpy as hp OldTobySensitivMap = hp.read_map('file.fits', dtype=numpy.float64) 這將使用 HEALPIx 投影打開一個 FIT ...

Healpix - 將經度、緯度轉換為銀河坐標

[英]Healpix - Converting longitude, latitude to galactic coordinates

我有一個 Healpix 像素協調文件https://wwwmpa.mpa-garching.mpg.de/~ensslin/research/data/faraday2020.html 這包含我想在天空圖中實現的緯度和經度值。然而，這些值需要首先轉換為銀河坐標 (l,b)。我的代碼是：有沒有 ...

Healpix / Healpy：我如何使用 healpy 對天空進行 SHT 處理？

[英]Healpix / Healpy: How can I use healpy to do SHT for a patch on the sky?

Healpix 是一個非常有用的軟件，可以對球體進行球面分析。例如，我們可以使用map2alm分析 Healpix RING 有序 map 並返回球面諧波。這里map的參數應該是一個具有$Npix = 12*Nside^2$的數組。如果我只想對一個補丁進行分析，而不是對一個完整的天空進行分 ...

復值 map 的球諧變換？

[英]Spherical harmonic transform of complex-valued map?

healpy 可以計算復值 map 的球諧變換嗎？當我使用healpy.sphtfunc.map2alm嘗試此操作時，沒有警告，但 function 僅針對 m>0 給出 a_{l,m}。這對於實值映射是有意義的，其中 a_{l,-m} = (-1)^m * a_{l,m}^*。但是 ...

將healpy（或healpix）map 旋轉到任意方向

[英]Rotating a healpy(or healpix) map to an arbitrary direction

我有一個治愈的 map，其中包含一個 object，中心位於 Dec=0 度和 RA=0 度的位置。現在我想將 object 的 position 移動到任意點，例如，Dec=20 度，RA=10 度。是否有一種簡單而精確的旋轉方式？我紅了治療文檔，但我仍然不清楚如何旋轉 map。假設我有一 ...

healpy anafast 導致大的 Cl 的

[英]healpy anafast leads to large Cl's

我正在使用 HEALPy 的 anafast 從普朗克地圖（數據 map [1] 或模擬地圖）中提取 Cl。似乎即使在我應用了他們的 Intensity common mask [2] 之后，我的功率譜也比他們的發布 [3] 大五倍。另一方面，如果我自己使用 synfast 合成一個 map，然 ...

Healpy 可以計算雙譜來研究 CMB map 的非高斯性嗎？

[英]Can Healpy compute a bispectrum to study non-gaussianity of the CMB map?

大家見鬼 Healpy 可以計算 CMB map 的雙譜嗎？ Healpy 庫中似乎沒有這樣的內置 function。謝謝！一切順利， ...

HEALpy - 從像素、l 和 m 的 function 獲取球面諧波

[英]HEALpy - Getting spherical harmonics from as function of pixel, l and m

我需要計算球面諧波的積分及其在球體的一小部分上的不同階數（l，m）的乘積（為了簡單起見，我們假設在半個球體上）。 Healpix/Healpy package 使得計算球面諧波變換變得非常容易和快速，將 function 乘以球面諧波 (SH) 的積分替換為預定義像素中的函數 * SH 值的總和（ ...