sentinel2 images lost related geographical information after correction

44 views

Skip to first unread message

Roberta Fagandini

unread,

Jul 6, 2017, 4:14:27 AM7/6/17

to RSGISLib Support

Hi all!

I used arcsi to perform atmospheric correction of Sentinel2 images, I tried both the example described in the arcsi.py file [1] and the one shown here [2]. Everything works fine in both cases but the output images lost their geographical information. I'm using a Sentinel2 L1C image of Italy (RS wgs84-utm32 epsg:32632).

Following you can find one of the command line that I launched:

arcsi.py --sensor sen2 -i S2A_MSIL1C_20170527T102031_N0205_R065_T32TMQ_20170527T102301.SAFE/MTD_MSIL1C.xml -o S2A_MSIL1C_20170527T102031_N0205_R065_T32TMQ_20170527T102301.SAFE/OutputsAOTInv --tmpath ./tmp -f KEA --stats -p RAD DOSAOTSGL SREF --aeroimg Documents/arcsi-1.4.2/data/WorldAerosolParams.kea --atmosimg Documents/arcsi-1.4.2/data/WorldAtmosphereParams.kea --dem S2A_MSIL1C_20170527T102031_N0205_R065_T32TMQ_20170527T102301.SAFE/dem_10m_utm32 --minaot 0.05 --maxaot 0.6 --simpledos

Is it possible to avoid this problem? How?

Thank you in advance!

Roberta

[1] Example: arcsi.py -s sen2 -p RAD TOA DOS --stats -f KEA -o ./Outputs \

                  -i ./S2A_MSIL1C_20170105T112441_N0204_R037_T30UVD_20170105T112439.SAFE/MTD_MSIL1C.xml

[2] https://spectraldifferences.wordpress.com/2015/02/27/inverting-for-aerosol-optial-thickness-aot-using-arcsi/

Daniel Clewley

unread,

Jul 6, 2017, 8:38:48 AM7/6/17

to Roberta Fagandini, RSGISLib Support

The images should retain all their geographic information but they will be in KEA format. Are you trying to open them in a program which doesn't have support for the KEA format, if so it will open as a HDF5 file and not read the spatial information.

You can convert to another format (e.g., GeoTiff) after processing using gdal_translate and this should copy across the spatial information.

Dan

--
You received this message because you are subscribed to the Google Groups "RSGISLib Support" group.
To unsubscribe from this group and stop receiving emails from it, send an email to rsgislib-support+unsubscribe@googlegroups.com.
To post to this group, send email to rsgislib-support@googlegroups.com.
Visit this group at https://groups.google.com/group/rsgislib-support.
For more options, visit https://groups.google.com/d/optout.

Roberta Fagandini

unread,

Jul 6, 2017, 9:40:53 AM7/6/17

to RSGISLib Support, robifa...@gmail.com

Thank you so much for your answer Dan!

Now it works..the problem was that I used gdal_translate outside the conda environment.

Thanks again!

Roberta

Il giorno giovedì 6 luglio 2017 14:38:48 UTC+2, Dan Clewley ha scritto:

The images should retain all their geographic information but they will be in KEA format. Are you trying to open them in a program which doesn't have support for the KEA format, if so it will open as a HDF5 file and not read the spatial information.

You can convert to another format (e.g., GeoTiff) after processing using gdal_translate and this should copy across the spatial information.

Dan

On 6 July 2017 at 09:14, Roberta Fagandini <robifa...@gmail.com> wrote:

Hi all!
I used arcsi to perform atmospheric correction of Sentinel2 images, I tried both the example described in the arcsi.py file [1] and the one shown here [2]. Everything works fine in both cases but the output images lost their geographical information. I'm using a Sentinel2 L1C image of Italy (RS wgs84-utm32 epsg:32632).

Following you can find one of the command line that I launched:

arcsi.py --sensor sen2 -i S2A_MSIL1C_20170527T102031_N0205_R065_T32TMQ_20170527T102301.SAFE/MTD_MSIL1C.xml -o S2A_MSIL1C_20170527T102031_N0205_R065_T32TMQ_20170527T102301.SAFE/OutputsAOTInv --tmpath ./tmp -f KEA --stats -p RAD DOSAOTSGL SREF --aeroimg Documents/arcsi-1.4.2/data/WorldAerosolParams.kea --atmosimg Documents/arcsi-1.4.2/data/WorldAtmosphereParams.kea --dem S2A_MSIL1C_20170527T102031_N0205_R065_T32TMQ_20170527T102301.SAFE/dem_10m_utm32 --minaot 0.05 --maxaot 0.6 --simpledos

Is it possible to avoid this problem? How?

Thank you in advance!

Roberta

[1] Example: arcsi.py -s sen2 -p RAD TOA DOS --stats -f KEA -o ./Outputs \
-i ./S2A_MSIL1C_20170105T112441_N0204_R037_T30UVD_20170105T112439.SAFE/MTD_MSIL1C.xml

[2] https://spectraldifferences.wordpress.com/2015/02/27/inverting-for-aerosol-optial-thickness-aot-using-arcsi/

--
You received this message because you are subscribed to the Google Groups "RSGISLib Support" group.

To unsubscribe from this group and stop receiving emails from it, send an email to rsgislib-suppo...@googlegroups.com.
To post to this group, send email to rsgislib...@googlegroups.com.

Reply all

Reply to author

Forward

0 new messages