Hi, I am trying to extract a band (byte 3) from the Quality_Assurance SDS of MODIS Cloud Mask swath product (MOD35-MYD35). Overall, I can not extract a band or array using gdalwarp since the dimensions of the array seems to be [Along, Across, Bytes] [2040x1354x10] and not [Bytes, Along, Across]. Therefore, I am trying to use gdalmdimtranslate to first transpose the array (transpose=[2,0,1]) and then do the extraction as GTiff.
Despite this, I am not totally sure about code writing or logic for this; maybe some of you can help me. My code is the following for the enclosed file below gdalmdimtranslate HDF4_EOS:EOS_SWATH:MOD35_L2.A2017161.1525.061.2017315035809.hdf:mod35:Quality_Assurance example.tif -scaleaxes 'QA_Dimension:mod35(3)' I am running this from the cmd as a test using gdal 3.1.1. I always get an ERROR 5: GDALDataset::GetRasterBand(-1) - Illegal band I also try to call the dimension in different ways out results. If you have a different approach to extract the Byte 3 of this SDS that is also welcome. I also enclosed a .txt as a results of gdalinfo MOD35_L2.A2017161.1525.061.2017315035809.hdf Best regards, https://drive.google.com/file/d/1fMZrfMLl18e_D6CjJbscwhLvoFFGORJp/view?usp=sharing https://drive.google.com/file/d/1f7oda-N9YT_cyf7YRw80BNQZr2O0RUMe/view?usp=sharing -- _______________________ J. Antonio Guzmán Q., *M.Sc* Ph.D Candidate and Vanier Scholar Earth and Atmospheric Sciences Department
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
