Python script to read and visualize HDF CALIPSO files
Python script to read and visualize HDF CALIPSO files
Hello everyone,
I'm trying to open CALIPSO hdf files and analyze the Vertical Feature Mask of them, and to do this I want to plot (via Python) a vertical profile of the data (with the Feature Type flags). I found the Compreensive Examples site:
https://www.hdfeos.org/zoo/index_openLaRC_Examples.php
But with the python scripts from there I can only plot the vertical profile of a single Feature Type, and it's a little bit confuse (I am studying to understand beter the data and how to handle with them) to adapt to what I need.
How can I make a full vertical Feature Type plot? I'm trying to do it, but it has been difficult. Thank you very much!
(OBS.: Here is an example of hdf file I'm using: CAL_LID_L2_VFM-Standard-V4-20.2006-06-13T04-54-47ZN.hdf)
I'm trying to open CALIPSO hdf files and analyze the Vertical Feature Mask of them, and to do this I want to plot (via Python) a vertical profile of the data (with the Feature Type flags). I found the Compreensive Examples site:
https://www.hdfeos.org/zoo/index_openLaRC_Examples.php
But with the python scripts from there I can only plot the vertical profile of a single Feature Type, and it's a little bit confuse (I am studying to understand beter the data and how to handle with them) to adapt to what I need.
How can I make a full vertical Feature Type plot? I'm trying to do it, but it has been difficult. Thank you very much!
(OBS.: Here is an example of hdf file I'm using: CAL_LID_L2_VFM-Standard-V4-20.2006-06-13T04-54-47ZN.hdf)
Tags:
-
- Posts: 17
- Joined: Tue Nov 19, 2019 11:41 am America/New_York
Re: Python script to read and visualize HDF CALIPSO files
Dear Ben-hur Martins,
Thank you for your interest in HDF CALIPSO files. I have reached out to our CALIPSO SME. We should have an answer for you soon.
Thank you for your interest in HDF CALIPSO files. I have reached out to our CALIPSO SME. We should have an answer for you soon.
-
- Posts: 17
- Joined: Tue Nov 19, 2019 11:41 am America/New_York
Re: Python script to read and visualize HDF CALIPSO files
Dear Ben-Hur Martins,
I have received the following information from our subject matter expert:
The sample code at hdfeos.org is tailored to only plot a single feature type; type=2(cloud).
# Extract Feature Type only (1-3 bits) through bitmask.
data = data & 7
……
# Focus on cloud (=2) data only.
data[data > 2] = 0;
data[data < 2] = 0;
data[data == 2] = 1;
Hopefully this information will help you. If you need further assistance please let us know.
I have received the following information from our subject matter expert:
The sample code at hdfeos.org is tailored to only plot a single feature type; type=2(cloud).
# Extract Feature Type only (1-3 bits) through bitmask.
data = data & 7
……
# Focus on cloud (=2) data only.
data[data > 2] = 0;
data[data < 2] = 0;
data[data == 2] = 1;
Hopefully this information will help you. If you need further assistance please let us know.
-
- User Services
- Posts: 33
- Joined: Mon Nov 23, 2020 3:57 pm America/New_York
Re: Python script to read and visualize HDF CALIPSO files
Hi @benhurmp,
Regarding your question on how to make a full Vertical Feature plot, are you looking for an image of all types plotting as a function of altitude (see example attached)?
If you'd like to connect with us virtually and walk us through what you are hoping to plot we'd be glad to take a look and offer up further assistance. Just let us know if you'd prefer this approach and we can set something up.
In the meantime, here are some additional resources that may assist you in plotting full vertical Feature Types:
VOCAL tool http://nasa-develop.github.io/VOCAL/dev ... index.html
How to Access and Visualize CALIPSO data http://meteothink.org/examples/meteoinf ... lipso.html
The NASA CALIPSO website also offers browse imagery where you can explore multiple images of CALIPSO profile observations https://www-calipso.larc.nasa.gov/produ ... 2020-07-01. This also includes the Vertical Feature Mask (VFM) browse image.
Hope this is helpful.
-Joe
Regarding your question on how to make a full Vertical Feature plot, are you looking for an image of all types plotting as a function of altitude (see example attached)?
If you'd like to connect with us virtually and walk us through what you are hoping to plot we'd be glad to take a look and offer up further assistance. Just let us know if you'd prefer this approach and we can set something up.
In the meantime, here are some additional resources that may assist you in plotting full vertical Feature Types:
VOCAL tool http://nasa-develop.github.io/VOCAL/dev ... index.html
How to Access and Visualize CALIPSO data http://meteothink.org/examples/meteoinf ... lipso.html
The NASA CALIPSO website also offers browse imagery where you can explore multiple images of CALIPSO profile observations https://www-calipso.larc.nasa.gov/produ ... 2020-07-01. This also includes the Vertical Feature Mask (VFM) browse image.
Hope this is helpful.
-Joe
- Attachments
-
- 2020-01-03_01-12-02_V4.10_3_6.png (85.27 KiB) Not viewed yet
Re: Python script to read and visualize HDF CALIPSO files
Hello @joseph.f.koch,
Your materials are great, and I'm trying to do exactly the kind of plot you attached. In the Meteothink website you sent me there is a script to plot all the Feature Types (the 3rd code), but when I try to run the script, the compiler shows a message saying that there is an inconsistency on line 32 (a ']' missing, perhaps). Can you help me, please?
Thank you! :)
Your materials are great, and I'm trying to do exactly the kind of plot you attached. In the Meteothink website you sent me there is a script to plot all the Feature Types (the 3rd code), but when I try to run the script, the compiler shows a message saying that there is an inconsistency on line 32 (a ']' missing, perhaps). Can you help me, please?
Thank you! :)
-
- User Services
- Posts: 33
- Joined: Mon Nov 23, 2020 3:57 pm America/New_York
Re: Python script to read and visualize HDF CALIPSO files
Hi @benhurmp,
I have reviewed line 32 and see no inconsistencies or syntax errors, so not sure why the script you are running is not able to execute, especially considering it is the same code that was used to generate the image listed on the Meteothink webpage.
You could try adding a "]" to where it may be missing, but I would recommend reaching out to Yaqiang Wang at yaqiang.wang@gmail.com to see if they can further assist since they wrote and maintain this code and we (ASDC) do not.
Hope this helps.
Thanks,
Joe
I have reviewed line 32 and see no inconsistencies or syntax errors, so not sure why the script you are running is not able to execute, especially considering it is the same code that was used to generate the image listed on the Meteothink webpage.
You could try adding a "]" to where it may be missing, but I would recommend reaching out to Yaqiang Wang at yaqiang.wang@gmail.com to see if they can further assist since they wrote and maintain this code and we (ASDC) do not.
Hope this helps.
Thanks,
Joe