ACSI error

[fabiano.c...@gmail.com]

Dear All,

I am new here and I just started with ARCSI.

I am trying tu run the following command line:

arcsi.py -s sen2 -p RAD  --stats -f KEA -o .\Fbcc\Desktop\arcsi\Outputs -i .\Desktop\arcsi\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738.SAFE\MTD_MSIL1C.xml

And I get the following  error message:

Start Time: 15:06:16, Wed Sep 09 2017.
Reading header file
Found Sentinel-2 image file product 'S2MSI1C'

Input imagery as listed in header file is present.

Get aero profile from image...
Converting Image Projection
Input Header: 'C:\Users\Fbcc\Desktop\arcsi\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738.SAFE\MTD_MSIL1C.xml'
Error: Coordinate System Transformation failed (Min).
ARCSI took in less than a second to process the input image. Thank you for using ARCSI.

What could I do to solve this?

Best regards

Fabiano

[Pete Bunting [pfb]]

Hi Fabiano,

Most likely the GDAL_DATA environment variable needs to be defined so GDAL can find the files it needs to convert coordinate systems.

I you have installed arcsi/gdal just into the base miniconda location then the path will be miniconda3/share/gdal/ if you have created an environment then it’ll be miniconda3/envs/<environment name>/share/gdal/

Best wishes,

Pete

****************************************************
* Dr Pete Bunting
* Reader in Remote Sensing
* Earth Observation and Ecosystem Dynamics Group
* Department of Geography and Earth Sciences
* Aberystwyth University
* Aberystwyth
* Ceredigion
* SY23 3DB
* UK
*
* Ph: +44 (0) 1970 622615
* Mob: +44 (0) 7917 842743
* Email: p...@aber.ac.uk
* ORCID: http://orcid.org/0000-0002-7435-0148
****************************************************
"Please consider the environment before printing this email or any documents attached”

> --
> 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.
> Visit this group at https://groups.google.com/group/rsgislib-support.
> For more options, visit https://groups.google.com/d/optout.



--------------------------------------------------------------------
Un o’r 5 prifysgol uchaf yn y DU a’r orau yng Nghymru am fodlonrwydd myfyrwyr.
(Arolwg Cenedlaethol y Myfyrwyr 2017)
www.aber.ac.uk

Top 5 UK university and best in Wales for student satisfaction
(National Student Survey 2017)
www.aber.ac.uk

[Fabiano Costantini]

Dear Pete,

Thank you very much it helped a lot. now I managed to run the process and i got this error message :.

Start Time: 15:19:55, Thu Sep 09 2017.

Reading header file
Found Sentinel-2 image file product 'S2MSI1C'

Input imagery as listed in header file is present.

Get aero profile from image...

Aerosol Profile =  Maritime

Get atmos profile from image...
Northern Hemisphere
Atmosphere Profile =  Tropical

Image Base Name: SEN2_20170224_lat15lon897_T15PZS

Checking Input Images are valid
Resampling image bands to match one another.
Started.0..10..20..30..40..50..60..70..80..90..100. Complete.

Input Header: 'C:\Users\Fbcc\Desktop\arcsi\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738.SAFE\MTD_MSIL1C.xml'

Output Basename: 'SEN2_20170224_lat15lon897_T15PZS'
Error: The HDF5ImageDataset driver does not support update access  to existing datasets.

Error: The following products were not generated:
RAD

Input Header: 'C:\Users\Fbcc\Desktop\arcsi\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738.SAFE\MTD_MSIL1C.xml'

Output Basename: 'SEN2_20170224_lat15lon897_T15PZS'



ARCSI took 00:00:01 to process the input image. Thank you for using ARCSI.

Kind Regards

Fabiano

> 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.

[Pete Bunting [pfb]]

Dear Fabiano,

It looks like the KEA file format isn't properly setup within your installation.

Can you run:

gdalinfo --formats (it might be gdalinfo.exe on windows) and check whether the KEA format is there?

Secondly, can you check which gdal version you have install from conda? Do you just have the one gdal installation on your system (I.e., from conda-forge)?

Best wishes

Pete

Get Outlook for iOS

---

From: Fabiano Costantini <fabiano.c...@gmail.com>
Sent: Thursday, September 14, 2017 3:33:29 PM
To: Pete Bunting [pfb]
Cc: RSGISLib Support
Subject: Re: [rsgislib-support] ACSI error

 

[Fabiano Costantini]

Dear Pete,

Thank you for your prompt reply.

You are right KEA format is not among the ones listed with gdalinfo --formats.

My gdal version is: GDAL 2.1.4, released 2017/06/23

Which one should I use instead?

Best ragards

Fabiano

[Pete Bunting [pfb]]

Dear Fabiano, 

I am using gdal 2.2.1 but in terms of your setup, I assume you have installed GDAL through conda?

Can you look at your package list and check the version but also check that it is from the conda-forge channel. If you are using the version of gdal from conda-forge then is should have KEA support build-in. 

Is conda the only GDAL installation you have on your system?

Best wishes, 

Pete

****************************************************
* Dr Pete Bunting
* Reader in Remote Sensing
* Earth Observation and Ecosystem Dynamics Group
* Department of Geography and Earth Sciences
* Aberystwyth University
* Aberystwyth
* Ceredigion
* SY23 3DB
* UK
* 
* Ph: +44 (0) 1970 622615
* Mob: +44 (0) 7917 842743
* Email: p...@aber.ac.uk
* ORCID: http://orcid.org/0000-0002-7435-0148
****************************************************

"Please consider the environment before printing this email or any documents attached”

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.

Visit this group at https://groups.google.com/group/rsgislib-support.
For more options, visit https://groups.google.com/d/optout.

[Daniel Clewley]

Hi Fabiano,

Assuming you installed using conda, did you install ARCSI into a separate environment?

Most of the environmental variables needed by GDAL are set when you run the activate script, which isn’t called if you run in the root environment.

Thanks,

Dan

[Fabiano Costantini]

Dear Pete and Daniel,

thank you very much. I installed gdal through conda forge and now KEA format is present when I type gdalinfo --formats.

Now when i process my command line I have the following error message:

arcsi.py -s sen2 -p RAD  --stats -f KEA -o C:\Users\Fbcc\Desktop\arcsi\Outputs -i C:\Users\Fbcc\Desktop\arcsi\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738.SAFE\MTD_MSIL1C.xml
ARCSI 1.4.2 Copyright (C) 2017 Pete Bunting, Dan Clewley
This program comes with ABSOLUTELY NO WARRANTY.
This is free software, and you are welcome to redistribute it
under certain conditions; See website (http://www.rsgislib.org/arcsi).
Bugs are to be reported to rsgislib...@googlegroups.com.

Start Time: 10:24:42, Fri Sep 09 2017.

Reading header file
Found Sentinel-2 image file product 'S2MSI1C'

Input imagery as listed in header file is present.

Get aero profile from image...
Aerosol Profile =  Maritime

Get atmos profile from image...
Northern Hemisphere
Atmosphere Profile =  Tropical

Image Base Name: SEN2_20170224_lat15lon897_T15PZS

Checking Input Images are valid
Resampling image bands to match one another.
Started.0..10..20..30..40..50..60..70..80..90..100. Complete.
Input Header: 'C:\Users\Fbcc\Desktop\arcsi\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738.SAFE\MTD_MSIL1C.xml'
Output Basename: 'SEN2_20170224_lat15lon897_T15PZS'

Error: Object given is not a Python function

Error: The following products were not generated:
RAD
Input Header: 'C:\Users\Fbcc\Desktop\arcsi\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738.SAFE\MTD_MSIL1C.xml'
Output Basename: 'SEN2_20170224_lat15lon897_T15PZS'



ARCSI took 00:00:01 to process the input image. Thank you for using ARCSI.

Thank you for the patience and I hope to get there soon !!

Fabiano

[Fabiano Costantini]

 Dear All,

I tried again many times and it seems I cannot find the solution of my problem (after making sure that arcsi,gdal and the other dependecies are under the same conda-forge enviroment).

Do you have and suggestion for that?

Best regards

Fabiano

2017-09-15 10:27 GMT+01:00 Fabiano Costantini <fabiano.c...@gmail.com>:

Dear Pete and Daniel,

thank you very much. I installed gdal through conda forge and now KEA format is present when I type gdalinfo --formats.

Now when i process my command line I have the following error message:

arcsi.py -s sen2 -p RAD  --stats -f KEA -o C:\Users\Fbcc\Desktop\arcsi\Outputs -i C:\Users\Fbcc\Desktop\arcsi\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738.SAFE\MTD_MSIL1C.xml
ARCSI 1.4.2 Copyright (C) 2017 Pete Bunting, Dan Clewley
This program comes with ABSOLUTELY NO WARRANTY.
This is free software, and you are welcome to redistribute it
under certain conditions; See website (http://www.rsgislib.org/arcsi).

Bugs are to be reported to rsgislib-support@googlegroups.com.

Fabiano

[Pete Bunting [pfb]]

Dear Fabiano, 

We’re shortly going to release a new version of ARCSI so hopefully that will fix the issue but the new version debug option so hopefully we’ll be able to get some more info on the error if it still doesn’t work. 

This is the first version we’ve released for windows so it sounds like we’ve still got some issues to work out. In the mean time if you use linux or MacOS then it is all working as far as I am aware. 

Best wishes, 

Pete

****************************************************
* Dr Pete Bunting
* Reader in Remote Sensing
* Earth Observation and Ecosystem Dynamics Group
* Department of Geography and Earth Sciences
* Aberystwyth University
* Aberystwyth
* Ceredigion
* SY23 3DB
* UK
* 
* Ph: +44 (0) 1970 622615
* Mob: +44 (0) 7917 842743
* Email: p...@aber.ac.uk
* ORCID: http://orcid.org/0000-0002-7435-0148
****************************************************

"Please consider the environment before printing this email or any documents attached”

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.

Visit this group at https://groups.google.com/group/rsgislib-support.
For more options, visit https://groups.google.com/d/optout.

[Fabiano Costantini]

Dear Pete,

Thanks for the suggestion ! I will move to Ubuntu and I will try again when the new Window version will be realesed.

Have a nice

Fabiano

[PI TI]

Hi everyone,

I do have the same problem for sentinel2. --> Error: Object given is not a Python function

My installation is on windows and I use ARCSI version 3.2.4.

I have dig into code and I think I have found where the problems comes from.

Basially this is a GDAL error.

Within miniconda3\envs\arcsienv\Lib\site-packages\rsgislib\imageutils\__init__.py file and specifically on line 507 there is a gdal call.

This is

gdal.Warp(outFile, inFile, options=wrpOpts)

Now, gettting into wrpOpts I found that  a callback function is passed among other info. But for some reason the callback is not called properly and this is causing the error.


More info about the case may be found here --> https://gis.stackexchange.com/questions/237479/using-callback-with-python-gdal-rasterizelayer

It is not the same but the problem is equal to the one we face here.

However, removing the callback form  the gdal call solves the problem and arcsi contiue to work withoun problems.

So in line 507 of the file mentioned I do change the mentioned line to this:

gdal.Warp(outFile, inFile)

For my self it seems to work as expected but I would be more than grateful if someone could verify that this change will not create any kind of faulty data.


I ll keep experimenting with that so if I come up with something interesting I ll let you know.








Thanks


Pavlos  

On 18 Sep 2017, at 16:30, Fabiano Costantini <fabiano....@gmail.com> wrote:

 Dear All,

I tried again many times and it seems I cannot find the solution of my problem (after making sure that arcsi,gdal and the other dependecies are under the same conda-forge enviroment).

Do you have and suggestion for that?

Best regards

Fabiano

2017-09-15 10:27 GMT+01:00 Fabiano Costantini <fabiano....@gmail.com>:

Dear Pete and Daniel,

thank you very much. I installed gdal through conda forge and now KEA format is present when I type gdalinfo --formats.

Now when i process my command line I have the following error message:

arcsi.py -s sen2 -p RAD  --stats -f KEA -o C:\Users\Fbcc\Desktop\arcsi\Outputs -i C:\Users\Fbcc\Desktop\arcsi\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738\S2A_MSIL1C_20170224T162331_N0204_R040_T15PZS_20170224T163738.SAFE\MTD_MSIL1C.xml
ARCSI 1.4.2 Copyright (C) 2017 Pete Bunting, Dan Clewley
This program comes with ABSOLUTELY NO WARRANTY.
This is free software, and you are welcome to redistribute it
under certain conditions; See website (http://www.rsgislib.org/arcsi).

Bugs are to be reported to rsgisli...@googlegroups.com.

Fabiano

From: Fabiano Costantini <fabiano....@gmail.com>

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

> Visit this group at https://groups.google.com/group/rsgislib-support.
> For more options, visit https://groups.google.com/d/optout.



--------------------------------------------------------------------
Un o’r 5 prifysgol uchaf yn y DU a’r orau yng Nghymru am fodlonrwydd myfyrwyr.
(Arolwg Cenedlaethol y Myfyrwyr 2017)
www.aber.ac.uk

Top 5 UK university and best in Wales for student satisfaction
(National Student Survey 2017)
www.aber.ac.uk

--------------------------------------------------------------------
Un o’r 5 prifysgol uchaf yn y DU a’r orau yng Nghymru am fodlonrwydd myfyrwyr.
(Arolwg Cenedlaethol y Myfyrwyr 2017)
www.aber.ac.uk

Top 5 UK university and best in Wales for student satisfaction
(National Student Survey 2017)
www.aber.ac.uk

--
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 rsgisli...@googlegroups.com.
To post to this group, send email to rsgisli...@googlegroups.com.

Visit this group at https://groups.google.com/group/rsgislib-support.
For more options, visit https://groups.google.com/d/optout.

--------------------------------------------------------------------
Un o’r 5 prifysgol uchaf yn y DU a’r orau yng Nghymru am fodlonrwydd myfyrwyr.
(Arolwg Cenedlaethol y Myfyrwyr 2017)
www.aber.ac.uk

Top 5 UK university and best in Wales for student satisfaction
(National Student Survey 2017)
www.aber.ac.uk

--
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 rsgisli...@googlegroups.com.
To post to this group, send email to rsgisli...@googlegroups.com.

Visit this group at https://groups.google.com/group/rsgislib-support.
For more options, visit https://groups.google.com/d/optout.

--
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...@googlegroups.com.
To post to this group, send email to rsgisli...@googlegroups.com.

Visit this group at https://groups.google.com/group/rsgislib-support.
For more options, visit https://groups.google.com/d/optout.