Spaces:
Running
Running
| import streamlit as st | |
| from streamlit_drawable_canvas import st_canvas | |
| from streamlit_image_coordinates import streamlit_image_coordinates | |
| from model.data_process.demo_data_process import process_ct_gt | |
| import numpy as np | |
| import matplotlib.pyplot as plt | |
| from PIL import Image, ImageDraw | |
| import monai.transforms as transforms | |
| from utils import show_points, make_fig, reflect_points_into_model, initial_rectangle, reflect_json_data_to_3D_box, reflect_box_into_model, run | |
| import nibabel as nib | |
| import tempfile | |
| print('script run') | |
| ############################################# | |
| # init session_state | |
| if 'option' not in st.session_state: | |
| st.session_state.option = None | |
| if 'text_prompt' not in st.session_state: | |
| st.session_state.text_prompt = None | |
| if 'reset_demo_case' not in st.session_state: | |
| st.session_state.reset_demo_case = False | |
| if 'preds_3D' not in st.session_state: | |
| st.session_state.preds_3D = None | |
| st.session_state.preds_3D_ori = None | |
| if 'data_item' not in st.session_state: | |
| st.session_state.data_item = None | |
| if 'points' not in st.session_state: | |
| st.session_state.points = [] | |
| if 'use_text_prompt' not in st.session_state: | |
| st.session_state.use_text_prompt = False | |
| if 'use_point_prompt' not in st.session_state: | |
| st.session_state.use_point_prompt = False | |
| if 'use_box_prompt' not in st.session_state: | |
| st.session_state.use_box_prompt = False | |
| if 'rectangle_3Dbox' not in st.session_state: | |
| st.session_state.rectangle_3Dbox = [0,0,0,0,0,0] | |
| if 'irregular_box' not in st.session_state: | |
| st.session_state.irregular_box = False | |
| if 'running' not in st.session_state: | |
| st.session_state.running = False | |
| if 'transparency' not in st.session_state: | |
| st.session_state.transparency = 0.25 | |
| case_list = [ | |
| 'model/asset/FLARE22_Tr_0002_0000.nii.gz', | |
| 'model/asset/FLARE22_Tr_0005_0000.nii.gz', | |
| 'model/asset/FLARE22_Tr_0034_0000.nii.gz', | |
| 'model/asset/FLARE22_Tr_0045_0000.nii.gz' | |
| ] | |
| ############################################# | |
| ############################################# | |
| # reset functions | |
| def clear_prompts(): | |
| st.session_state.points = [] | |
| st.session_state.rectangle_3Dbox = [0,0,0,0,0,0] | |
| def reset_demo_case(): | |
| st.session_state.data_item = None | |
| st.session_state.reset_demo_case = True | |
| clear_prompts() | |
| def clear_file(): | |
| st.session_state.option = None | |
| process_ct_gt.clear() | |
| reset_demo_case() | |
| clear_prompts() | |
| ############################################# | |
| st.image(Image.open('model/asset/overview back.png'), use_column_width=True) | |
| github_col, arxive_col = st.columns(2) | |
| with github_col: | |
| st.write('GitHub repo:https://github.com/BAAI-DCAI/SegVol') | |
| with arxive_col: | |
| st.write('Paper:https://arxiv.org/abs/2311.13385') | |
| # modify demo case here | |
| demo_type = st.radio( | |
| "Demo case source", | |
| ["Select", "Upload"], | |
| on_change=clear_file | |
| ) | |
| if demo_type=="Select": | |
| uploaded_file = st.selectbox( | |
| "Select a demo case", | |
| case_list, | |
| index=None, | |
| placeholder="Select a demo case...", | |
| on_change=reset_demo_case | |
| ) | |
| else: | |
| uploaded_file = st.file_uploader("Upload demo case(nii.gz)", type='nii.gz', on_change=reset_demo_case) | |
| st.session_state.option = uploaded_file | |
| if st.session_state.option is not None and \ | |
| st.session_state.reset_demo_case or (st.session_state.data_item is None and st.session_state.option is not None): | |
| st.session_state.data_item = process_ct_gt(st.session_state.option) | |
| st.session_state.reset_demo_case = False | |
| st.session_state.preds_3D = None | |
| st.session_state.preds_3D_ori = None | |
| prompt_col1, prompt_col2 = st.columns(2) | |
| with prompt_col1: | |
| st.session_state.use_text_prompt = st.toggle('Sematic prompt') | |
| text_prompt_type = st.radio( | |
| "Sematic prompt type", | |
| ["Predefined", "Custom"], | |
| disabled=(not st.session_state.use_text_prompt) | |
| ) | |
| if text_prompt_type == "Predefined": | |
| pre_text = st.selectbox( | |
| "Predefined anatomical category:", | |
| ['liver', 'right kidney', 'spleen', 'pancreas', 'aorta', 'inferior vena cava', 'right adrenal gland', 'left adrenal gland', 'gallbladder', 'esophagus', 'stomach', 'duodenum', 'left kidney'], | |
| index=None, | |
| disabled=(not st.session_state.use_text_prompt) | |
| ) | |
| else: | |
| pre_text = st.text_input('Enter an Anatomical word or phrase:', None, max_chars=20, | |
| disabled=(not st.session_state.use_text_prompt)) | |
| if pre_text is None or len(pre_text) > 0: | |
| st.session_state.text_prompt = pre_text | |
| else: | |
| st.session_state.text_prompt = None | |
| with prompt_col2: | |
| spatial_prompt_on = st.toggle('Spatial prompt', on_change=clear_prompts) | |
| spatial_prompt = st.radio( | |
| "Spatial prompt type", | |
| ["Point prompt", "Box prompt"], | |
| on_change=clear_prompts, | |
| disabled=(not spatial_prompt_on)) | |
| st.session_state.enforce_zoom = st.checkbox('Enforce zoom-out-zoom-in') | |
| if spatial_prompt == "Point prompt": | |
| st.session_state.use_point_prompt = True | |
| st.session_state.use_box_prompt = False | |
| elif spatial_prompt == "Box prompt": | |
| st.session_state.use_box_prompt = True | |
| st.session_state.use_point_prompt = False | |
| else: | |
| st.session_state.use_point_prompt = False | |
| st.session_state.use_box_prompt = False | |
| if not spatial_prompt_on: | |
| st.session_state.use_point_prompt = False | |
| st.session_state.use_box_prompt = False | |
| if not st.session_state.use_text_prompt: | |
| st.session_state.text_prompt = None | |
| if st.session_state.option is None: | |
| st.write('please select demo case first') | |
| else: | |
| image_3D = st.session_state.data_item['z_image'][0].numpy() | |
| col_control1, col_control2 = st.columns(2) | |
| with col_control1: | |
| selected_index_z = st.slider('X-Y view', 0, image_3D.shape[0] - 1, 162, key='xy', disabled=st.session_state.running) | |
| with col_control2: | |
| selected_index_y = st.slider('X-Z view', 0, image_3D.shape[1] - 1, 162, key='xz', disabled=st.session_state.running) | |
| if st.session_state.use_box_prompt: | |
| top, bottom = st.select_slider( | |
| 'Top and bottom of box', | |
| options=range(0, 325), | |
| value=(0, 324), | |
| disabled=st.session_state.running | |
| ) | |
| st.session_state.rectangle_3Dbox[0] = top | |
| st.session_state.rectangle_3Dbox[3] = bottom | |
| col_image1, col_image2 = st.columns(2) | |
| if st.session_state.preds_3D is not None: | |
| st.session_state.transparency = st.slider('Mask opacity', 0.0, 1.0, 0.25, disabled=st.session_state.running) | |
| with col_image1: | |
| image_z_array = image_3D[selected_index_z] | |
| preds_z_array = None | |
| if st.session_state.preds_3D is not None: | |
| preds_z_array = st.session_state.preds_3D[selected_index_z] | |
| image_z = make_fig(image_z_array, preds_z_array, st.session_state.points, selected_index_z, 'xy') | |
| if st.session_state.use_point_prompt: | |
| value_xy = streamlit_image_coordinates(image_z, width=325) | |
| if value_xy is not None: | |
| point_ax_xy = (selected_index_z, value_xy['y'], value_xy['x']) | |
| if len(st.session_state.points) >= 3: | |
| st.warning('Max point num is 3', icon="⚠️") | |
| elif point_ax_xy not in st.session_state.points: | |
| st.session_state.points.append(point_ax_xy) | |
| print('point_ax_xy add rerun') | |
| st.rerun() | |
| elif st.session_state.use_box_prompt: | |
| canvas_result_xy = st_canvas( | |
| fill_color="rgba(255, 165, 0, 0.3)", # Fixed fill color with some opacity | |
| stroke_width=3, | |
| stroke_color='#2909F1', | |
| background_image=image_z, | |
| update_streamlit=True, | |
| height=325, | |
| width=325, | |
| drawing_mode='transform', | |
| point_display_radius=0, | |
| key="canvas_xy", | |
| initial_drawing=initial_rectangle, | |
| display_toolbar=True | |
| ) | |
| try: | |
| print(canvas_result_xy.json_data['objects'][0]['angle']) | |
| if canvas_result_xy.json_data['objects'][0]['angle'] != 0: | |
| st.warning('Rotating is undefined behavior', icon="⚠️") | |
| st.session_state.irregular_box = True | |
| else: | |
| st.session_state.irregular_box = False | |
| reflect_json_data_to_3D_box(canvas_result_xy.json_data, view='xy') | |
| except: | |
| print('exception') | |
| pass | |
| else: | |
| st.image(image_z, use_column_width=False) | |
| with col_image2: | |
| image_y_array = image_3D[:, selected_index_y, :] | |
| preds_y_array = None | |
| if st.session_state.preds_3D is not None: | |
| preds_y_array = st.session_state.preds_3D[:, selected_index_y, :] | |
| image_y = make_fig(image_y_array, preds_y_array, st.session_state.points, selected_index_y, 'xz') | |
| if st.session_state.use_point_prompt: | |
| value_yz = streamlit_image_coordinates(image_y, width=325) | |
| if value_yz is not None: | |
| point_ax_xz = (value_yz['y'], selected_index_y, value_yz['x']) | |
| if len(st.session_state.points) >= 3: | |
| st.warning('Max point num is 3', icon="⚠️") | |
| elif point_ax_xz not in st.session_state.points: | |
| st.session_state.points.append(point_ax_xz) | |
| print('point_ax_xz add rerun') | |
| st.rerun() | |
| elif st.session_state.use_box_prompt: | |
| if st.session_state.rectangle_3Dbox[1] <= selected_index_y and selected_index_y <= st.session_state.rectangle_3Dbox[4]: | |
| draw = ImageDraw.Draw(image_y) | |
| #rectangle xz view (upper-left and lower-right) | |
| rectangle_coords = [(st.session_state.rectangle_3Dbox[2], st.session_state.rectangle_3Dbox[0]), | |
| (st.session_state.rectangle_3Dbox[5], st.session_state.rectangle_3Dbox[3])] | |
| # Draw the rectangle on the image | |
| draw.rectangle(rectangle_coords, outline='#2909F1', width=3) | |
| st.image(image_y, use_column_width=False) | |
| else: | |
| st.image(image_y, use_column_width=False) | |
| col1, col2, col3 = st.columns(3) | |
| with col1: | |
| if st.button("Clear", use_container_width=True, | |
| disabled=(st.session_state.option is None or (len(st.session_state.points)==0 and not st.session_state.use_box_prompt and st.session_state.preds_3D is None))): | |
| clear_prompts() | |
| st.session_state.preds_3D = None | |
| st.session_state.preds_3D_ori = None | |
| st.rerun() | |
| with col2: | |
| img_nii = None | |
| if st.session_state.preds_3D_ori is not None and st.session_state.data_item is not None: | |
| meta_dict = st.session_state.data_item['meta'] | |
| pred_array = st.session_state.preds_3D_ori.transpose(2, 1, 0) | |
| img_nii = nib.Nifti1Image(pred_array, affine=meta_dict['affine']) | |
| with tempfile.NamedTemporaryFile(suffix=".nii.gz") as tmpfile: | |
| nib.save(img_nii, tmpfile.name) | |
| with open(tmpfile.name, "rb") as f: | |
| bytes_data = f.read() | |
| st.download_button( | |
| label="Download result(.nii.gz)", | |
| data=bytes_data, | |
| file_name="segvol_preds.nii.gz", | |
| mime="application/octet-stream", | |
| disabled=img_nii is None | |
| ) | |
| with col3: | |
| run_button_name = 'Run'if not st.session_state.running else 'Running' | |
| if st.button(run_button_name, type="primary", use_container_width=True, | |
| disabled=( | |
| st.session_state.data_item is None or | |
| (st.session_state.text_prompt is None and len(st.session_state.points) == 0 and st.session_state.use_box_prompt is False) or | |
| st.session_state.irregular_box or | |
| st.session_state.running | |
| )): | |
| st.session_state.running = True | |
| st.rerun() | |
| # if len(st.session_state.points) > 0: | |
| # st.write(st.session_state.points) | |
| if st.session_state.running: | |
| st.session_state.running = False | |
| with st.status("Running...", expanded=False) as status: | |
| run() | |
| st.rerun() |