me and my student are using Landsat-7 ETM imagery to create high-resolution ice-water maps.
For this we create a product based on bands 2 to 4 which we read from the *_MTL.txt files
We had problems saving the product and/or a subset of the created bands / products and were told that we need to re-sample the bands first before we do our computations. So even though one does computations based purely on the bands with 30 m pixel size SNAP is mixing up information from all bands - also the panchromatic one, which has 15 m pixel size. Fine, we did the re-sampling and were able to safe the product properly for the Arctic. By properly we mean: the product, latitude, and longitude for each pixel.
Now, the same procedure seems not to work for images from Antarctica. Neither are we able to create a “subset” in the “export” mode of SNAP nor are we able to properly safe the latitude and longitude coordinates along with the product - no matter whether we do or not do the re-sampling. Either the coordinates are not written at all (which happens when we re-sample the bands before we do our computations) or the coordinates are zero everywhere.
What we tried:
Import data -> re-sampling -> our computations -> export into netCDF
Import data -> re-sampling -> our computations -> export as BEAM-DIMAP or geoTIFF -> re-import -> export into netCDF
Import data -> our computations -> deletion of all bands we don’t need -> export into netCDF
Import data -> our computations -> deletion of all bands we don’t need -> re-sampling -> export into netCDF
All these did not work and we don’t understand while it is fine for the Arctic but not fine for the Antarctic.
Perhaps somebody has an idea into what / which information we could look else to solve the riddle? Thank you so much!
p.s.: As a side note: Even though Landsat-5 TM also has bands with different pixel resolutions we did not need to do any re-sampling and working with the data of bands 2 to 4 with SNAP was simply a pleasure.