______________________________________________________________ Je rappelle que ce qu'on prevoit maintenant suite a diverses discussions c est d'utiliser HEALPIX pour le codage et le stockage des images d'un survey donné à différentes résolutions. Une fois converti en HEALPIX on n'a plus de problemes pour passer d'une coin du ciel au voisin puisqu'on a une projection valable partout, pas centrée en fait, et qui nous donne des pixels quadrangulaires de surface egale. reprojection des survey en mode Healpix ------------------------------------- Pour un survey donné et une résolution donnée on parcourt les pixels healpix et on calcule l'alpha delta central de ce pixel... On peut alors aller chercher ce pixel dans les ficheirs d'origine du survey (en utilisant eventuellement le serveur classique) et ecrire dans le (les) fichier(s) HEALPIX la valeur trouvée. Affichage: ----------- En affichage on se choisit une grille ecran dans une projection donnée. La cartesienne a l'avantage de ne pas necessiter de calcul trigonometrique pour passer du pixel a la coordonée alpha delta. elle n'est pas adaptée aux tres garnds champs (superieurs a 20/30 deg) et aux tres grandes ballades sauf a recalculer autour d'un nouveau centre. On peut bien sur aussi choisr une autre projection en sortie avec un cout calcul cependant. Pour chaqye pixel une fois qu'on a son alpha delta on peut calculer le pixel HEALPIX a la resolution consideree et recuperer la valeur qu'il contient pour l'affecter a ce pixel ecran.... Stockage -------- Il y a un probleme d'organisation: Si on met les fichiers a toutes les resolutions en FITS/HEALPIX simple on ne gagne rien puisqu'on aura des fichiers gigantesques pour les petites resolutions. le mode "ring" n est pas interessant non plus puisque il ne permet pas de faire des zones dans le plan mais seulement de parcourir des anneaux. En fait le modes nested en partant d'un découpage en 12 losanges permet de faire une organisation impriquée puisaue chaque losange a une resolution contient 4 losanges a la resolution inferieure. la numerotation des pixels est telle que si on est face au pixel losange "N" a une resolution Nside donnée, les 4 pixels inclus dans ce losange a la resolution Nside + 1 sont numerotés de 4*(N-1)+1 à 4*N. Si on considere deux resolutions eloignées de plus d'un cran (Nside2-Nside1) diff de 1, alors les numeros de pixels pour la resolution Nside2 varient de (4**(Nside2-Nside1))*(N-1)+1 a (4**(Nside2-Nside1))*N Pour voir comment organiser la base en terme de fichiers a chaque resolution on va faire quelques suppositions. - On veut descendre en résolution par des facteurs 4 sur le coté du pixel (16 en surface, tests d'Anais de l'ete dernier) - On souhaite avoir une image à l'ecran du type 1024x1024. Pour la faire a partir d'un nombre de morceaux raisonnables on peut donc faire un decoupage en 1024x1024 (4 morceaux en moyenne) ou 512x512 ( 9 morceaux en moyenne) pour realiser cela le "losange fichier" doit avoir un Nside inférieur de 9 (1024x1024) ou de 8(512x512) au Nside du "losange pixel". Voici un petit tableau des tailles de losange (coté) en fonction du Nside Nside side Nside side 0 full sky 12 1.71 arcmin 1 58.65 deg 13 0.85 arcmin 2 29.82 14 25.76 arcsec 3 14.66 15 12.88 arcsec 4 7.32 16 6.44 5 3.66 17 3.22 6 1.83 18 1.61 7 0.92 19 0.81 8 0.46 20 0.40 9 13.74 arcmin 21 0.20 10 6.87 22 0.10 11 3.44 Avec une organisation zoom4/1024x1024 çà nous donnerait la succession suivante "coté" du fichier "coté" du pixel full sky 13.74arcmin 29.82 deg 3.44arcmin 7.33 deg 0.85 arcmin 1.83 deg 12.88 arcsec 0.46 deg 3.22 arcsec 6.86 arcmin 0.81 arcsec 1.71 arcmin 0.20 arcsec Avec une organisation zoom4/512x512 çà nous donnerait la succession suivante "coté" du fichier "coté" du pixel full sky 0.46 deg 29.82 deg 6.87arcmin 7.33 deg 1.71 arcmin 1.83 deg 25.76 arcsec 0.46 deg 6.44 arcsec 6.86 arcmin 1.61 arcsec 1.71 arcmin 0.40 arcsec les resolutions 0.20/0.40 arcsec sont adaprtées pour l'ACS 0.81 pour Mama et DSS2, 2MASS, DENIS 1.71 arcmin pour IRAS, WENSS 6.87 arcmin pour Planck High 0.46 deg pour Planck Low