DEM error when using ARCSI for processing Worldview data

[Gwawr Jones]

Hi,

I get the following error when trying to use the --dem function in ARCSI to process Worldview-2 and Worldview-3 data:

(arcsienv2) gwawr.jones@gwawrj-ubuntu:~/Documents/057190571010_01_P001_MUL$ arcsi.py -s wv2 -f KEA --stats -p RAD TOA SREF STDSREF METADATA --aeropro Maritime --atmospro Tropical --aot 0.05 --surfacealtitude 0.122 --dem srtm_24_09.tif -o ./OutputsAOT005 -t ./tmp -i 17SEP20142742-M2AS-057190571010_01_P001.XML --imagefile 17SEP20142742-M2AS-057190571010_01_P001_Ortho_UTM20N.tif ARCSI 2.2.24 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:38:45, Thu Nov 11 2017.

 

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: WV2_20170920_lat0lon6749

 

Checking Input Images are valid

 

Converting to Radiance

Opening: /home/gwawr.jones/Documents/057190571010_01_P001_MUL/17SEP20142742-M2AS-057190571010_01_P001_Ortho_UTM20N.tif

New image width = 6717 height = 10445 bands = 8

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Setting Band Names...

Calculating Statistics...

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Calculating Image Pyramids.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

 

Converting to TOA

Open ./OutputsAOT005/WV2_20170920_lat0lon6749_rad.kea

New image width = 6717 height = 10445 bands = 8

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Setting Band Names...

Calculating Statistics...

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Calculating Image Pyramids.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

Started .0..10..20..30..40..50..60..70..80..90..100. Complete.

 

Started.0..10..20..30..40..50..60..70..80..90. Complete.

Subset and reproject DEM...

Output DEM:  ./OutputsAOT005/WV2_20170920_lat0lon6749_dem.kea

Input Header: '17SEP20142742-M2AS-057190571010_01_P001.XML'

Output Basename: 'WV2_20170920_lat0lon6749'

Error: maskImage() argument 2 must be str, not None

Error: The following products were not generated:

SREF STDSREF TOPOSHADOW METADATA

Input Header: '17SEP20142742-M2AS-057190571010_01_P001.XML'

Output Basename: 'WV2_20170920_lat0lon6749'

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

I've used 3 different DEM's that I've got available and I get the same error each time. Any ideas?

Thanks,

Gwawr

[Daniel Clewley]

Hi Gwawr,

Can you try running the same command again but add the `--debug` flag and post the output. The debug option will print more info about where the error is occurring. 

Thanks,

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

[Gwawr Jones]

Hi Dan,

Here's the error with the debug option:

Started.0..10..20..30..40..50..60..70..80..90. Complete.

Subset and reproject DEM...

Output DEM:  ./OutputsAOT005/WV2_20170920_lat0lon7349_dem.kea

Input Header: '17SEP20142742-M2AS-057190571010_01_P001.XML'

Output Basename: 'WV2_20170920_lat0lon7349'

Error: maskImage() argument 2 must be str, not None

Error: The following products were not generated:

TOPOSHADOW METADATA SREF STDSREF

Input Header: '17SEP20142742-M2AS-057190571010_01_P001.XML'

Output Basename: 'WV2_20170920_lat0lon7349'

Traceback (most recent call last):

  File "/home/gwawr.jones/miniconda3/envs/arcsienv2/bin/arcsi.py", line 605, in <module>

    arcsilib.arcsirun.runARCSI(args.inputheader, args.imagefile, args.cloudmask, args.sensor, args.inwkt, args.format, args.outpath, args.outbasename, args.outwkt, args.outproj4, args.projabbv, args.ximgres, args.yimgres, args.prods, args.stats, args.aeropro, args.atmospro, args.aeroimg, args.atmosimg, args.grdrefl, args.surfacealtitude, args.atmosozone, args.atmoswater, atmosOZoneWaterSpecified, args.aerowater, args.aerodust, args.aerooceanic, args.aerosoot, aeroComponentsSpecified, args.aot, args.vis, args.tmpath, args.minaot, args.maxaot, args.lowaot, args.upaot, args.dem, args.demnodata, args.aotfile, (not args.localdos), args.dosout, args.simpledos, args.debug, args.scalefac, args.interp, args.interpresamp, args.cs_initdist, args.cs_initminsize, args.cs_finaldist, args.cs_morphop, args.fullimgouts, args.checkouts, args.classmlclouds, args.cloudtrainclouds, args.cloudtrainother, args.resample2lowres, args.keepfileends)

  File "/home/gwawr.jones/miniconda3/envs/arcsienv2/lib/python3.5/site-packages/arcsilib/arcsirun.py", line 1307, in runARCSI

    prepareDEM(paramsObj)

  File "/home/gwawr.jones/miniconda3/envs/arcsienv2/lib/python3.5/site-packages/arcsilib/arcsirun.py", line 991, in prepareDEM

    rsgislib.imageutils.maskImage(outDEMNameTmp, paramsObj.validMaskImage, paramsObj.outDEMName, paramsObj.outFormat, rsgislib.TYPE_32FLOAT, paramsObj.demNoDataVal, 0)

TypeError: maskImage() argument 2 must be str, not None

Thanks,

Gwawr

[Daniel Colson]

Hi,

I've been getting the same issue when trying to process WorldView-2 imagery. I was wondering if this was resolved, or if anyone has found a solution?

Many thanks

Daniel

[Daniel Colson]

Apologies, I meant to say the version I was running. I have tried this with both ARCSI version 2.2.24 and version 3.1.6, with a couple of DEMs (local higher-res and SRTM) in both .tif and .kea formats. Any help is much appreciated.

Many thanks

Daniel