guilhem_lavaux/vide_public
0
0
Fork 0
mirror of https://bitbucket.org/cosmicvoids/vide_public.git synced 2025-07-04 15:21:11 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

Allow mask of any resolution

Browse source
This commit is contained in:
Guilhem Lavaux 2021-02-16 12:56:00 +02:00
parent affb56ff48
commit 85e0f36c99
1 changed files with 9 additions and 13 deletions
Download patch file Download diff file Expand all files Collapse all files

22
c_tools/mock/generateFromCatalog.cpp
View file

@ -377,10 +377,10 @@ void generateSurfaceMask(generateFromCatalog_info& args ,
output_data.uniqueID.push_back(-1); output_data.uniqueID.push_back(-1);
insertion++; insertion++;
fprintf(fp, "%e %e %e\n", // fprintf(fp, "%e %e %e\n",
(p.xyz[0]), // (p.xyz[0]),
(p.xyz[1]), // (p.xyz[1]),
(p.xyz[2])); // (p.xyz[2]));
} }
} }
} }
@ -445,10 +445,10 @@ void generateSurfaceMask(generateFromCatalog_info& args ,
output_data.redshift.push_back(-1); output_data.redshift.push_back(-1);
output_data.uniqueID.push_back(-1); output_data.uniqueID.push_back(-1);
insertion++; insertion++;
fprintf(fp, "%e %e %e\n", // fprintf(fp, "%e %e %e\n",
(p.xyz[0]), // (p.xyz[0]),
(p.xyz[1]), // (p.xyz[1]),
(p.xyz[2])); // (p.xyz[2]));
} }
fclose(fp); fclose(fp);
// END PMS // END PMS
@ -606,7 +606,6 @@ int main(int argc, char **argv)
cout << "Loading data " << args_info.catalog_arg << "..." << endl; cout << "Loading data " << args_info.catalog_arg << "..." << endl;
vector<NYU_Data> data; vector<NYU_Data> data;
Healpix_Map<float> o_mask;
vector<int> pixel_list; vector<int> pixel_list;
vector<int> full_mask_list; vector<int> full_mask_list;
ParticleData output_data; ParticleData output_data;
@ -615,12 +614,9 @@ int main(int argc, char **argv)
cout << "Loading mask..." << endl; cout << "Loading mask..." << endl;
read_Healpix_map_from_fits(args_info.mask_arg, o_mask);
Healpix_Map<float> mask; Healpix_Map<float> mask;
mask.SetNside(128, RING); read_Healpix_map_from_fits(args_info.mask_arg, mask);
mask.Import(o_mask);
computeContourPixels(mask,pixel_list); computeContourPixels(mask,pixel_list);
computeMaskPixels(mask,full_mask_list); computeMaskPixels(mask,full_mask_list);