Mercurial > repos > goeckslab > rename_tiff_channels
view scale_cell_coordinates.py @ 6:4763362914f4 draft default tip
planemo upload for repository https://github.com/goeckslab/tools-mti/tree/main/tools/mti-utils commit 5c324d06ad54dca76aa2b4c8118df5654f49da66
author | goeckslab |
---|---|
date | Mon, 15 Jul 2024 21:37:39 +0000 |
parents | 61fd94c7ce42 |
children |
line wrap: on
line source
import argparse import json import warnings import anndata as ad def main(inputs, output): """ inputs : str File path to galaxy tool JSON inputs config file output: str File path to save the output h5ad file """ warnings.simplefilter('ignore') # read inputs JSON with open(inputs, 'r') as param_handler: params = json.load(param_handler) # read input anndata file adata = ad.read_h5ad(params['anndata']) # scale coords unit = params['unit'] new_col_names = [] for c in [params['x_coord'], params['y_coord']]: scaled_col_name = f'{c}_{unit}' adata.obs[scaled_col_name] = adata.obs[c] * params['resolution'] new_col_names.append(scaled_col_name) # overwrite adata.obsm['spatial'] with scaled coordinates adata.obsm['spatial'] = adata.obs[scaled_col_name].values # write out anndata to h5ad file adata.write_h5ad(output) if __name__ == '__main__': aparser = argparse.ArgumentParser() aparser.add_argument( "-i", "--inputs", dest="inputs", required=True) aparser.add_argument( "-o", "--output", dest="output", required=False) args = aparser.parse_args() main(args.inputs, args.output)