2

I have been trying to mask clouds with the quality control band with landsat 8 imagery. The problem I continue to run into is that the quality control bands mistake tops of buildings and concrete as clouds.

Have you guys had better success with Fmask? I tried using it earlier, but couldn't find out where it put the file. The website said it would be called XXXFmask. Would I even be able to open Fmask files in Arcmap?

These are the steps I've taken. 1) Built raster attribute table for QC band (see below). 2) Selected only the pixels whose values said there were clouds present. 3) Used Raster to Polygon tool to convert selected pixels to a shapefile. 4) Selected all the cloud polygons and used mask tool on image analysis window to remove data from the cloud cells. 5) Now I plan to run my image classification tools. I've included an image of Greenville, TN before and after the raster is masked. NoData cells are displayed in pink for visibility.

QC layer's attribute tableComparison of the city of Greenville, TN with and without masking.  Masked cells with NoData displayed in pink.

  • 1
    what procedure are you applying to use the QC band? – nickves Oct 31 '14 at 6:32
  • Ok I just edited my post to include more detail. Just let me know what you think. Thanks! – Hgotwald Nov 1 '14 at 19:02
  • 1
    Please take a look at Landsat QBinfo and an example on how to create bitmasks with GRASS. Unfortunally i don't have arcgis to provide a step-by-step procedure. Probably you get a pure 'cloud mask' using the expression con("QC_BAND" & 49152 == 49152,1,0) with map calculator. Take into account that there is the possibility to 'false' positives where you must manually decide if they are clouds or not. BTW can you provide the scene name? – nickves Nov 1 '14 at 19:32
  • By scene name do you mean the Landsat 8 tile I downloaded? It's downloaded .tar file's name was LC80180352014263LGN00. The date is from 20-Sept-2014, and it's path 18 and row 35. – Hgotwald Nov 1 '14 at 20:14
1

What I ended up doing:

1) Download and decompress your landsat tile. Do not open and look at the images in ArcMap because ArcMap will create .proj on the fly and this will confuse poor fmask and prevent it from running.

2) Once your landsat tile of choice's file is decompressed/unzipped/etc. Follow the directions on fmask's website for masking clouds. I used their stand alone program so I didn't have to download Matlab.

3) After fmask runs there will be a file titled xxxfmask in with all the bands for the landsat imagry.

4) This file can only be opened in ENVI NOT in ArcGIS

5) I was lucky enough to have access to ENVI so I loaded the files onto one of my school's computers and opened them in ENVI and converted them to .shp files

6) Do whatever you want with the .shp cloud masks in Arcmap

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.