Geographic projection missing in output

[nmatton]

Sorry to hear about that. But it seems to be a gdal issue.
Be sure to have the JP2000 driver installed (look for JP2OpenJPEG (most common) in output of gdalinfo --formats )

[unnic]

thanks for the hint. Indeed the JP2000 driver was missing on the machine where sen2cor 2.2.1 was installed.
I removed sen2cor and installed GDAL individually. JPEG2000 driver is there and .jp2 can be read (gdalinfo). Unfortunately installing sen2cor 2.2.1 overwrites the former gdal installation thus removing the driver. I tried to remove the gdal installation from the setup.py file but wasn’t successfull.

Any hints are welcome

[AEA]

Dear unnic,
I really appreciate your reply and the time invested in making it work. However, if you don’t mine I’ll wait for you and the rest of pros to solve it before trying…
Thank you!

[Jan]

@unnic

Hi,

I’ve raised it with my colleagues, and it appears that it is a known issue. As to when it is sorted; that I cannot say with any certainty. Sorry.

I’ve asked my colleagues to look at this thread.

Cheers

Jan

[Rea]

Hello @nmatton,

I tried to georefrence L2A images using your script on windows x64 system but script proccesing isn’t getting any progress at all. All i see is 0.0% at any product i tried to process.

Any help?

[Rea]

Hello,

i am trying to use nmatton’s script @nmatton but the progress always remains at 0.0%, no error appears and nothing happens. I had exactly the same result using the script on windows and linux system. @unnic the script worked in your case? What could be the cause? Maybe i should use a specific version of python and gdal ? Well, i was thinking about doing the projection using ground control points in ArcGis but i think nmatton’s code could help me a lot.

[unnic]

@ Rea

I suppose it’s the same problem as described here: Geographic projection missing in output . You can try by checking if the jp2 driver (JP2OpenJPEG) is missing in the output of gdalinfo --formats

@nmatton s approach isn’t much different from gdalcopyproj.py (using gdal.Open() ) so it’s also not working in this case because the driver is missing.

I suggest uninstalling everything (Python, Anaconda …) and reinstalling sen2cor 2.0.6 which does not remove the needed driver.

[Rea]

@unnic
Thanks for your answer.

I installed Anaconda on Ubuntu 16.04 and then manually i installed osgeo gdal package 1.11.2. I checked the installations and all packages respond correctly. Also jp2 driver is supported by GDAL as you can see in the image bellow. But still script is not running it stucks on progress 0.0%…

I don’t know what else to do

[unnic]

sry, I can’t speak for nmattons script, it did not work for me but gdalcopyproj.py does.

[unnic]

@nmatton
So I managed to reinstall GDAL with a driver for jp2 files.
running gdalinfo filename.jp2 from the command line returns the expected output but in python

import gdal
data=gdal.Open(filename)

does not return anything.

I assume python is not directing to the right GDAL installation Is there a way to set the installation path of gdal within python?

[Val]

Hi,

I was in the belief that the most recent version of sen2cor (2.3.0) fixed the problem with missing geo reference in the L2A output. Anyhow, for me the problem still occurs … the only thing new is the output of *geo.jp2 files, which I can’t seem to open or use.

Am I missing something obvious or can someone shed some light on this for me?

Much appreciated,

Val

[unnic]

Hey Val,

i just processed a recent acquisition with sen2cor 2.3.0. The resulting .jp2 files are projected correctly.
Check if L2A_process -h indicates that the correct version is used; e.g.: Sentinel-2 Level 2A Prototype Processor (Sen2Cor). Version: 2.3.0,

[Val]

Hi unnic,

thanks for the answer.

I just checked again, and I’m running definitely version 2.3.0.

Sentinel-2 Level 2A Prototype Processor (Sen2Cor). Version: 2.3.0, created:
2016.11.18, supporting Level-1C product version: 14.

That’s what’s confusing me … also the *geo.jp2 files. Are these supposed to be there or maybe a sign of an unfinished process?

[unnic]

So your sen2cor output .jp2 files do open but are simply unprojected (like it was with sen2cor <2.3) ? Any error messages?

[Val]

Precisely … with the only difference of having additionally files like

S2A_USER_MSI_L2A_TL_MTI__20160709T114209_A005466_T36PWT_B02_10m_geo.jp2

as well.

No error messages during the processing, the only thing I’ve noticed is that sen2cor ocassionally “pauses” until I hit the return button … which I tried to sort out here:

Sen2cor stuck?

Hi, I'm running the most recent version of sen2cor on a Win7 x64 machine in cmd. Everything seems to be working fine, but occasionally the process seems to be "stuck" , meaning there is no error or malfunction, it's just not progressing any more. It's simply resolved by hitting the return key once, which is a bit of a problem for automating the process. I'm wondering if anyone experienced similar issues or if you have any other tips for me. Thanks! Val

[unnic]

I can check if I get the same result for your granule.
Did you try reinstalling sen2cor?

[Val]

Doesn’t seem to be a Granule issue … i processed somewhat 30 Gs and the issue is the same in all of them.

I guess re-installing is the next logical step. I just wanted to make sure if I’m missing something here … Thanks for your efforts!

Much appreciated

[Val]

Well well well … reinstalling sen2cor actually did the trick. I feel like I’m back in the 90s, when the only solution on hand was giving the ol’ TV a proper smack.

Anyways, thanks for the support @unnic

[m.anelli]

Hi,
I read your discussion and I found many similarities with my issue. After the latest update of SNAP, I can’t obtain georeferenced products if I run sen2cor. I tried reinstalling the plugin but I still have this problem.

Could you suggest me a solution?
Thanks in advance!

[Val]

Sorry, I can’t say more than make sure you’re running the most recent version of sen2cor (2.3.0) and if so, uninstall and re-install it following the instructions.

This did the trick for me.

Best