解决pd.read_csv出现的PermissionError

原因:路径出现错误,没有写上文件名称。

data=pd.read_csv(r"F:\anaconda",engine='python',encoding='utf-8')

修改:加上文件名称(dating.txt)后运行成功。

data=pd.read_csv(r"F:\anaconda\dating.txt",engine='python',encoding='utf-8')

 

print("样本名称:", sys.argv[1]) print("单细胞参考-表达矩阵:", sys.argv[2]) print("单细胞参考-细胞类型信息:", sys.argv[3]) print("N_cells_per_location = ", sys.argv[4], "----", sys.argv[5]) sample = sys.argv[1] ## spot的细胞核识别数量的最大值: # WMQ-605-mspl:62 # Ncellsperspot = int(sys.argv[4]) proindex= sys.argv[5] # save_path = "/data1/chengrui/workspace/zst_st/result/cell2location/" ## save_path = "/data1/zhaoshutao/projectworkspace/cell2location_mouse_atlas_20241014_ref_filter/c2lres/"+ proindex +"/" # sample_folder = "/data1/chengrui/workspace/zst_st/ourdata_main/matrix/" sample_folder = "/data1/chengrui/workspace/zst_st/all_data/" # save_file = save_path + sample save_file = save_path + sample + "/1000/" os.makedirs(save_file, exist_ok=True) countsfilepath = sys.argv[2] if countsfilepath.endswith('.csv'): data=sc.read_csv(countsfilepath) elif countsfilepath.endswith('.txt'): data=sc.read_text(countsfilepath) else: print("Unsupported file format. Please use .csv or .txt files.") # 可以选择退出程序或者抛出异常 sys.exit(1) # meta = pd.read_csv(sys.argv[3]) # meta = pd.read_table(sys.argv[3]) file_path = sys.argv[3] if file_path.endswith('.csv'): meta = pd.read_csv(file_path) elif file_path.endswith('.txt'): meta = pd.read_table(file_path) else: print("Unsupported file format. Please use .csv or .txt files.") # 可以选择退出程序或者抛出异常 sys.exit(1) ##################################################################### adata_ref=data.transpose() os.chdir(save_file) os.getcwd() if sample == "WMQ-606-mSpl": sample_path = sample_folder + "WMQ-606-mSpleen/outs" elif sample in ["WMQ-627-mSpl","WMQ-629-mLym","WMQ-636-mLung","WMQ-678-mLiver","LXZ-009","WMQ-712-mTes","WMQ-714-mTes"]: if sample == "WMQ-627-mSpl": sample_path = sample_folder + "WMQ-627-mSpleen" else: sample_path = sample_folder + sample elif sample == "WMQ-731-mOva": sample_path = sample_folder + "WMQ-731-mOVA-bu" elif sample == "WMQ-732-mOva": sample_path = sample_folder + "WMQ-732-mOVA-bu" elif sample in ["WMQ-765-mOVA","WMQ-766-mOVA","WMQ-767-mOVA","WMQ-768-mOVA"]: sample_path = "/data1/zhaoshutao/projectworkspace/rawdata_20241222_ovary_st/" + sample +"/outs" else: sample_path = sample_folder + sample +"/outs" adata_st = sc.read_visium(sample_path) adata_st if sample in ["WMQ-586-mHeart", "WMQ-584-mHeart", "WMQ-644-mHeart"]: selected_barcodes = pd.read_csv("".join(["/data1/chengrui/workspace/zst_st/eye/",sample,"_pointsinregions.csv"])) selected_barcodes = selected_barcodes["barcode"] ## 扣掉血液位置 adata_st = adata_st[~adata_st.obs.index.isin(selected_barcodes), :].copy() elif sample in ["WMQ-567-mEye", "WMQ-594-mEye", "WMQ-593-mEye"]: selected_barcodes = pd.read_csv("".join(["/data1/chengrui/workspace/zst_st/eye/",sample,"_pointsinregions.csv"])) selected_barcodes = selected_barcodes["barcode"] ## 仅保留目标位置 adata_st = adata_st[adata_st.obs.index.isin(selected_barcodes), :].copy() else: pass adata_st检查代码有没有问题,并提供注释
03-14
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值