Why does the cloud cover percent displayed in the MTL .txt file of the force-level1-landsat program differ from the cloud cover % displayed in the log of force-level2 #352
-
I ran the following lines of codes (1 and 2) in Powershell to search and download all the LANDSAT 8 OLI sensor images in 2024 with cloud cover between 0 and 10%, clipped to the a .shp vector (I was able to make FORCE work through Windows Powershell, since it did not work in the Windows Subsystem for Linux, and I am stuck with a Windows workstation at work (I prefer Linux)): (1) OUTPUTS:
(2) OUTPUT:
Since the imagery is all between 0 and 10% cloud cover, I thought that keeping the cloud over percent as 10 (%) for the parameter MAX_CLOUD_COVER_FRAME under CLOUD DETECTION OPTIONS for the force-level2 program would not have any effect, since the paramater has been described as
However, after I run The file "LC09_L1TP_136040_20240310_20240310_02_T1_MTL.txt" in the LC09_L1TP_136040_20240310_20240310_02_T1.tar states:
However, when I check the corresponding force-level2 log files (after running
Why is the cloud cover 100% in this and skipping the processing? Is this a bug? I am attaching both the metadata and log files mentioned in this question With kind regards, |
Beta Was this translation helpful? Give feedback.
Replies: 2 comments
-
Hi @paulvpop, FORCE uses a different version of Fmask for internal cloud detection. Also, USGS metadata of cloud cover do not include cloud shadow, and also no buffer distances. FORCE cloud coverages are usually somewhat larger. Re the specific image: I cannot reproduce this, thus there still might be some issue on your end. I still end up with too much cloud cover, but this is all confusion between snow/ice and clouds, and hence less relevant for most applications. I am not sure what could cause this on your end, however. Probably the DEM (I used the Copernicus DEM)? FORCE version (I used 3.8.01)?
Best, |
Beta Was this translation helpful? Give feedback.
-
Hello David, Thanks for the clarification. Yeah, I also think that the algorithm may be confusing snow cover for clouds. I actually would need the imagery with snow cover (but less cloud cover) since I will be using it for land cover classification, including snow as a class. I checked the specific imagery in question. It indeed does have a lot of snow. I am attaching the imagery (showing the outline of the study area map touching the imagery). Based on literature survey, I found that a spline-interpolated void-filled 90 m CGIAR-updated SRTM DEM was best (in comparison to using a 30 m resolution SRTM available from any source or ASTER). So, I prepared the interpolated 30 m DEM using Google Earth Engine myself for my specific study area. I had not considered a Copernicus DEM or an ALOS World 3D which I am reading is better. So, I will try with them as well. Thanks for bringing into my attention the availability of the Copernicus DEM. I am also using v 3.8.01. With kind regards, |
Beta Was this translation helpful? Give feedback.
Hi @paulvpop,
FORCE uses a different version of Fmask for internal cloud detection. Also, USGS metadata of cloud cover do not include cloud shadow, and also no buffer distances. FORCE cloud coverages are usually somewhat larger.
Re the specific image: I cannot reproduce this, thus there still might be some issue on your end. I still end up with too much cloud cover, but this is all confusion between snow/ice and clouds, and hence less relevant for most applications. I am not sure what could cause this on your end, however. Probably the DEM (I used the Copernicus DEM)? FORCE version (I used 3.8.01)?