question about how to use python visualize cloudsat product.

[GuangYou-Zhou]

Recently i read about your moonbook project "Get the latitude and longitude of a 2B-CLDCLASS-lidar granule with python". I am a undergraduate student study about it. But i'm not familiar with python so i was wonder if you can share the code about how to visualize that 2B-CLDCLASS-LIDAR cloud phase product. The picture was on the bottom of that moonbook project. Thanks a lot!

[benjamin-hg-marchant]

I created a jupyter notebook showing how to read and plot 2B-CLDCLASS LIDAR Product. It can be found here: https://en.moonbooks.org/Jupyter/Read_and_Plot_CloudSat_2B_CLDCLASS_LIDAR_Product/ I use the following file as an example: https://drive.google.com/file/d/1tcdqc1VjRjhA3i392Gq3B7CVWRjzaeBy/view?usp=sharing Note: To run the code, you can use python anaconda distribution ( https://www.anaconda.com/products/distribution) and install the following modules: conda install -c anaconda jupyter conda install -c anaconda ipykernel python -m ipykernel install --user --name=worklab conda install -c anaconda numpy conda install -c conda-forge pyhdf conda install -c anaconda netcdf4 conda install -c conda-forge netcdf4

…

On Wed, Mar 22, 2023 at 8:42 AM GuangYou-Zhou ***@***.***> wrote: Recently i read about your moonbook project "Get the latitude and longitude of a 2B-CLDCLASS-lidar granule with python". I am a undergraduate student study about it. But i'm not familiar with python so i was wonder if you can share the code about how to visualize that 2B-CLDCLASS-LIDAR cloud phase product. The picture was on the bottom of that moonbook project. Thanks a lot! — Reply to this email directly, view it on GitHub <#1>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AA7BRNSKFH4XDKWFHCP5E53W5LXUBANCNFSM6AAAAAAWDYNBYQ> . You are receiving this because you are subscribed to this thread.Message ID: ***@***.*** com>

[GuangYou-Zhou]

thanks a lot! That really helps. Actually I'm trying to plot cloud phase product. But your code seems can only plot cloud mask. if you have code about how to plot it i'll be appreciate a lot. But thanks for your code again. I may try to do it myself.

[benjamin-hg-marchant]

GuangYou, see below the code to plot the cloud phase: f = plt.figure() ax = f.add_subplot(111) for i in range(cldclass_lidar_start_idx,cldclass_lidar_end_idx): nb_cloud_layer = np.where(CloudLayerBase[i,:] < 0 )[0][0] for j in range(nb_cloud_layer): if CloudLayerBase[i,j] > 0 and CloudLayerTop[i,j] > 0.0: bar_xi = i bar_width = 1.0 bar_base = CloudLayerBase[i,j] bar_height = CloudLayerTop[i,j] - bar_base bar_color = '1.0' bar_edgecolor = '1.0' if CloudPhase[i,j] == 3: bar_color = 'blue' bar_edgecolor = 'blue' if CloudPhase[i,j] == 2: bar_color = 'salmon' bar_edgecolor = 'salmon' plt.bar(bar_xi, height=bar_height, width=bar_width, bottom=bar_base, color=bar_color, edgecolor=bar_edgecolor) ax.set_facecolor('xkcd:lightblue') plt.bar(cldclass_lidar_start_idx, height=0, width=0, bottom=0, color='1.0', edgecolor='1.0', label='ice clouds') plt.bar(cldclass_lidar_start_idx, height=0, width=0, bottom=0, color='salmon', edgecolor='salmon', label='mixed clouds') plt.bar(cldclass_lidar_start_idx, height=0, width=0, bottom=0, color='blue', edgecolor='blue', label='liquid clouds') plt.legend() plt.xlim(cldclass_lidar_start_idx,cldclass_lidar_end_idx) plt.ylim(0,20) xticks_idx = [int(i) for i in np.linspace(cldclass_lidar_start_idx,cldclass_lidar_end_idx,10)] xticks_latitude = [ "{0:.2f}".format(latitude_table[i][0]) for i in xticks_idx] plt.xticks(xticks_idx, xticks_latitude,rotation=45) plt.xlabel('Latitude') plt.ylabel('Altitude (in km)') plt.title('Cloud thermodynamic phase') #plt.savefig( 'plot_'+tag+'_cldclass_lidar_cloud_layers_and_thermodynamic_phase_method_01.png', dpi=100, bbox_inches='tight' ) plt.show() plt.close()

…

On Fri, Mar 31, 2023 at 2:15 AM GuangYou-Zhou ***@***.***> wrote: thanks a lot! That really helps. Actually I'm trying to plot cloud phase product. But your code seems can only plot cloud mask. if you have code about how to plot it i'll be appreciate a lot. But thanks for your code again. I may try to do it myself. — Reply to this email directly, view it on GitHub <#1 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AA7BRNSR6YSU2ZFI5NKJPKLW6ZZBFANCNFSM6AAAAAAWDYNBYQ> . You are receiving this because you commented.Message ID: <benjamin-hg-marchant/ml-dl-modis-multilayer-clouds/issues/1/1491355638@ github.com>

[GuangYou-Zhou]

It works! Thanks a lot! Next i will try to write Remarks and learn how it works. I notice you are working on the cloud machine learning project. i'd like to read them maybe Cite them! Appreciate it again!

[benjamin-hg-marchant]

Good!, I am glad you were able to make it work GuangYou.

…

On Thu, Apr 6, 2023 at 11:26 PM GuangYou-Zhou ***@***.***> wrote: It works! Thanks a lot! Next i will try to write Remarks and learn how it works. I notice you are working on the cloud machine learning project. i'd like to read them maybe Cite them! Appreciate it again! — Reply to this email directly, view it on GitHub <#1 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AA7BRNSUDQS2PPVKLL45BWDW76CM7ANCNFSM6AAAAAAWDYNBYQ> . You are receiving this because you commented.Message ID: <benjamin-hg-marchant/ml-dl-modis-multilayer-clouds/issues/1/1499892818@ github.com>