Error module 'acolite.aerlut' has no attribute 'interp2d' for WorldView

[hnthang]

Hi colleges,

I applied the acolite_mr for a scene of WorldView but got the error:

Code: Select all

AttributeError: module 'acolite.aerlut' has no attribute 'interp2d'

The details are as follow:

Code: Select all

File "/.../acolite_mr/acolite/acolite/acolite_mr_ac.py", line 96, in acolite_mr_ac
    ac.worldview.worldview_ac(bundle, output=output, limit=limit,
  File "/.../acolite_mr/acolite/worldview/worldview_ac.py", line 103, in worldview_ac
    ttg = ac.ac.gas_transmittance(sza, vza, uoz=uoz, uwv=uwv)
  File "/.../acolite_mr/acolite/ac/gas_transmittance.py", line 20, in gas_transmittance
    wv_o2_hs, tt_o2_hs = ac.ac.o2lut_interp(sza, vza)
  File "/.../acolite_mr/acolite/ac/o2lut_interp.py", line 19, in o2lut_interp
    iw = ac.aerlut.interp2d(o2lut[:,:,par_id,:], ths_id, thv_id)
AttributeError: module 'acolite.aerlut' has no attribute 'interp2d'

I checked the code at line 19 of o2lut_interp.py:

Code: Select all

iw = ac.aerlut.interp2d(o2lut[:,:,par_id,:], ths_id, thv_id)

but not sure this error regarding the lack of any package/import or the module itself.

Many thanks,
Thang

[quinten]

Dear Thang

The interp2d function should be there (see: https://github.com/acolite/acolite_mr/b ... nterp2d.py)

Do you have another "acolite" package in your Python path that gets imported before the one from acolite_mr?

Quinten

[hnthang]

Dear Quinten,

I checked the file interp2d and yes, it is located in the folder /.../acolite_mr/acolite/aerlut/interp2d.py. It seems there is not any acolite package imported before. I used acolite for Sentinel 2 processing a few days ago, but I created another environment (i.e acolite_mr) to run acolite_mr independently. However, the error is still there:

Code: Select all

    iw = ac.aerlut.interp2d(o2lut[:,:,par_id,:], ths_id, thv_id)
AttributeError: module 'acolite.aerlut' has no attribute 'interp2d'

I clone the latest version of acolite_mr, cleaning all the command in the bash but it did not work. Do you have any ideas for further check? Should I test on another computer or check any path in the system? I acknowledge not a fatal error but the acolite_mr just stop.

Kind reagrds,
Thang

[quinten]

Hi Thang

Could you run:

Code: Select all

import sys
print(sys.path)

Quinten

[hnthang]

Hi Quinten,

Here is the path for acolite_mr env:

Code: Select all

>>> import sys
>>> print(sys.path)
['', '/home/metal/anaconda3/envs/acolite_mr/lib/python38.zip', '/home/metal/anaconda3/envs/acolite_mr/lib/python3.8', '/home/metal/anaconda3/envs/acolite_mr/lib/python3.8/lib-dynload', '/home/metal/anaconda3/envs/acolite_mr/lib/python3.8/site-packages']

the path of the system:

Code: Select all

>>> import sys
>>> print(sys.path)
['', '/home/metal/anaconda3/lib/python37.zip', '/home/metal/anaconda3/lib/python3.7', '/home/metal/anaconda3/lib/python3.7/lib-dynload', '/home/metal/anaconda3/lib/python3.7/site-packages']

and the rest for another env:

Code: Select all

>>> import sys
>>> print(sys.path)
['', '/home/metal/anaconda3/envs/thang/lib/python38.zip', '/home/metal/anaconda3/envs/thang/lib/python3.8', '/home/metal/anaconda3/envs/thang/lib/python3.8/lib-dynload', '/home/metal/anaconda3/envs/thang/lib/python3.8/site-packages']

Many thanks,
Thang

[quinten]

Hi Thang

My apologies for not checking before, it seems there is a line missing in acolite/aerlut/__init__.py

Could you add:

Code: Select all

 from .interp2d import interp2d

and try again?

Quinten

[hnthang]

Hi Quinten,

Thank you very much, it worked by adding a new line of code. It's now, another problem with the RHOS values with only NaN values whilst there's no problem with the RHOT.

This is the output code:

Code: Select all

Launching ACOLITE MR processing.
Processing worldview file
2020-05-06T23:54:32 - Started processing /home/metal/Downloads/Compressed/010217881010_01_P002_MUL
2020-05-06T23:54:32 - Imported LUTs
PONDER-LUT-201704-MOD1 0.001 0.048881449723840634 0 4
PONDER-LUT-201704-MOD2 0.001 0.04889167508863663 0 4
2020-05-06T23:54:46 - Fitted model C, band pair COASTAL:RED, taua 550: 0.001
2020-05-06T23:54:46 - Writing lat/lon
Processing tile R1C1 [0, 0]
2020-05-06T23:55:22 - Finished writing /home/metal/Downloads/Compressed/010217881010_01_P002_MUL/ac_wv/WorldView2_2012_07_11_15_34_18_P002_L2R.nc

I tested several time but (1) acolite_mr only work with the format .TIF, not .tif and (2) the the RHOS output only has NaN values.

My folder contains the following files:

Code: Select all

12JUL11153418-M2AS_R1C1-010217881010_01_P002.TIF
12JUL11153418-M2AS-010217881010_01_P002.XML

From your guidance:

Code: Select all

The directory containing the subtile .tif files and .xml metadata should be passed as $in.

but if I run with the extension .tif, error out:

Code: Select all

Launching ACOLITE MR processing.
Processing worldview file
2020-05-07T00:26:40 - Started processing /home/metal/Downloads/Compressed/010217881010_01_P002_MUL
2020-05-07T00:26:48 - Imported LUTs
Traceback (most recent call last):
  File "acolite_mr.py", line 106, in <module>
    run_acolite_mr()
  File "acolite_mr.py", line 91, in run_acolite_mr
    ac.acolite.acolite_mr_ac(args.input, output=args.output, limit=args.limit,
  File "/home/metal/Documents/phd/atmospheric/acolite_mr/acolite/acolite/acolite_mr_ac.py", line 96, in acolite_mr_ac
    ac.worldview.worldview_ac(bundle, output=output, limit=limit,
  File "/home/metal/Documents/phd/atmospheric/acolite_mr/acolite/worldview/worldview_ac.py", line 164, in worldview_ac
    rhod[band_names[b]]['rhod'] = np.sort(rhod[band_names[b]]['rhod'])
KeyError: 'COASTAL'

And will only 1 file .nc is created? It's pretty different from the processing of PlanetScope (wth .nc and RGB of RHOS, RHOT).

Thang

[quinten]

Dear Thang

Just passing the directory for processing is the way to go. As to why you have only NaNs, this could be due to processing set up or an input file issue. The aerosol optical depth is computed as 0.001 (the minimum in the LUTs), indicating that the dark spectrum (or the whole image) may be lower than Rayleigh reflectance. A couple of questions to help the troubleshooting:

• Is there a dark spectrum plot generated?
• Is your target at high elevation?
• Do you have a TOA image, or has an atmospheric correction already been applied?
• Could you perhaps share the image or the metadata file?

Quinten

[hnthang]

Dear Quinten,

Regrading your questions:

1. Is there a dark spectrum plot generated?
No, only a .NC file generated (WorldView2_2012_07_11_15_34_18_P002_L2R.nc)

2. Is your target at high elevation?
Yes, it is in the mountainous areas.

3. Do you have a TOA image, or has an atmospheric correction already been applied?
No, the input image is radiometric corrected and not applied atmospheric correction before.

4. Could you perhaps share the image or the metadata file?
We can not share the image (due to private project) but please have a look at the attached XML file (I paste here due to a invalid extension problem when attaching). Is there a change to raise any relevant errors in this case in the python processing windows?

Code: Select all

<PAMDataset>
  <SRS>PROJCS["WGS 84 / UTM zone 19S",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-69],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","32719"]]</SRS>
  <GeoTransform> 2.1406799997402000e+005, 2.0000000000000000e+000, 0.0000000000000000e+000, 8.5563420000015702e+006, 0.0000000000000000e+000,-2.0000000000000000e+000</GeoTransform>
  <Metadata domain="IMD">
    <MDI key="bandId">"Multi"</MDI>
    <MDI key="BAND_B.absCalFactor">1.260825e-02</MDI>
    <MDI key="BAND_B.effectiveBandwidth">5.430000e-02</MDI>
    <MDI key="BAND_B.LLHAE">2332.70</MDI>
    <MDI key="BAND_B.LLLat">-13.22847154</MDI>
    <MDI key="BAND_B.LLLon">-71.63832808</MDI>
    <MDI key="BAND_B.LRHAE">2332.70</MDI>
    <MDI key="BAND_B.LRLat">-13.22909668</MDI>
    <MDI key="BAND_B.LRLon">-71.57710406</MDI>
    <MDI key="BAND_B.TDILevel">14</MDI>
    <MDI key="BAND_B.ULHAE">2332.70</MDI>
    <MDI key="BAND_B.ULLat">-13.04555359</MDI>
    <MDI key="BAND_B.ULLon">-71.63637407</MDI>
    <MDI key="BAND_B.URHAE">2332.70</MDI>
    <MDI key="BAND_B.URLat">-13.04616980</MDI>
    <MDI key="BAND_B.URLon">-71.57519527</MDI>
    <MDI key="BAND_C.absCalFactor">9.295654e-03</MDI>
    <MDI key="BAND_C.effectiveBandwidth">4.730000e-02</MDI>
    <MDI key="BAND_C.LLHAE">2332.70</MDI>
    <MDI key="BAND_C.LLLat">-13.22847154</MDI>
    <MDI key="BAND_C.LLLon">-71.63832808</MDI>
    <MDI key="BAND_C.LRHAE">2332.70</MDI>
    <MDI key="BAND_C.LRLat">-13.22909668</MDI>
    <MDI key="BAND_C.LRLon">-71.57710406</MDI>
    <MDI key="BAND_C.TDILevel">24</MDI>
    <MDI key="BAND_C.ULHAE">2332.70</MDI>
    <MDI key="BAND_C.ULLat">-13.04555359</MDI>
    <MDI key="BAND_C.ULLon">-71.63637407</MDI>
    <MDI key="BAND_C.URHAE">2332.70</MDI>
    <MDI key="BAND_C.URLat">-13.04616980</MDI>
    <MDI key="BAND_C.URLon">-71.57519527</MDI>
    <MDI key="BAND_G.absCalFactor">9.713071e-03</MDI>
    <MDI key="BAND_G.effectiveBandwidth">6.300000e-02</MDI>
    <MDI key="BAND_G.LLHAE">2332.70</MDI>
    <MDI key="BAND_G.LLLat">-13.22847154</MDI>
    <MDI key="BAND_G.LLLon">-71.63832808</MDI>
    <MDI key="BAND_G.LRHAE">2332.70</MDI>
    <MDI key="BAND_G.LRLat">-13.22909668</MDI>
    <MDI key="BAND_G.LRLon">-71.57710406</MDI>
    <MDI key="BAND_G.TDILevel">14</MDI>
    <MDI key="BAND_G.ULHAE">2332.70</MDI>
    <MDI key="BAND_G.ULLat">-13.04555359</MDI>
    <MDI key="BAND_G.ULLon">-71.63637407</MDI>
    <MDI key="BAND_G.URHAE">2332.70</MDI>
    <MDI key="BAND_G.URLat">-13.04616980</MDI>
    <MDI key="BAND_G.URLon">-71.57519527</MDI>
    <MDI key="BAND_N.absCalFactor">1.224380e-02</MDI>
    <MDI key="BAND_N.effectiveBandwidth">9.890000e-02</MDI>
    <MDI key="BAND_N.LLHAE">2332.70</MDI>
    <MDI key="BAND_N.LLLat">-13.22847154</MDI>
    <MDI key="BAND_N.LLLon">-71.63832808</MDI>
    <MDI key="BAND_N.LRHAE">2332.70</MDI>
    <MDI key="BAND_N.LRLat">-13.22909668</MDI>
    <MDI key="BAND_N.LRLon">-71.57710406</MDI>
    <MDI key="BAND_N.TDILevel">10</MDI>
    <MDI key="BAND_N.ULHAE">2332.70</MDI>
    <MDI key="BAND_N.ULLat">-13.04555359</MDI>
    <MDI key="BAND_N.ULLon">-71.63637407</MDI>
    <MDI key="BAND_N.URHAE">2332.70</MDI>
    <MDI key="BAND_N.URLat">-13.04616980</MDI>
    <MDI key="BAND_N.URLon">-71.57519527</MDI>
    <MDI key="BAND_N2.absCalFactor">9.042234e-03</MDI>
    <MDI key="BAND_N2.effectiveBandwidth">9.960000e-02</MDI>
    <MDI key="BAND_N2.LLHAE">2332.70</MDI>
    <MDI key="BAND_N2.LLLat">-13.22847154</MDI>
    <MDI key="BAND_N2.LLLon">-71.63832808</MDI>
    <MDI key="BAND_N2.LRHAE">2332.70</MDI>
    <MDI key="BAND_N2.LRLat">-13.22909668</MDI>
    <MDI key="BAND_N2.LRLon">-71.57710406</MDI>
    <MDI key="BAND_N2.TDILevel">24</MDI>
    <MDI key="BAND_N2.ULHAE">2332.70</MDI>
    <MDI key="BAND_N2.ULLat">-13.04555359</MDI>
    <MDI key="BAND_N2.ULLon">-71.63637407</MDI>
    <MDI key="BAND_N2.URHAE">2332.70</MDI>
    <MDI key="BAND_N2.URLat">-13.04616980</MDI>
    <MDI key="BAND_N2.URLon">-71.57519527</MDI>
    <MDI key="BAND_R.absCalFactor">1.103623e-02</MDI>
    <MDI key="BAND_R.effectiveBandwidth">5.740000e-02</MDI>
    <MDI key="BAND_R.LLHAE">2332.70</MDI>
    <MDI key="BAND_R.LLLat">-13.22847154</MDI>
    <MDI key="BAND_R.LLLon">-71.63832808</MDI>
    <MDI key="BAND_R.LRHAE">2332.70</MDI>
    <MDI key="BAND_R.LRLat">-13.22909668</MDI>
    <MDI key="BAND_R.LRLon">-71.57710406</MDI>
    <MDI key="BAND_R.TDILevel">10</MDI>
    <MDI key="BAND_R.ULHAE">2332.70</MDI>
    <MDI key="BAND_R.ULLat">-13.04555359</MDI>
    <MDI key="BAND_R.ULLon">-71.63637407</MDI>
    <MDI key="BAND_R.URHAE">2332.70</MDI>
    <MDI key="BAND_R.URLat">-13.04616980</MDI>
    <MDI key="BAND_R.URLon">-71.57519527</MDI>
    <MDI key="BAND_RE.absCalFactor">4.539619e-03</MDI>
    <MDI key="BAND_RE.effectiveBandwidth">3.930000e-02</MDI>
    <MDI key="BAND_RE.LLHAE">2332.70</MDI>
    <MDI key="BAND_RE.LLLat">-13.22847154</MDI>
    <MDI key="BAND_RE.LLLon">-71.63832808</MDI>
    <MDI key="BAND_RE.LRHAE">2332.70</MDI>
    <MDI key="BAND_RE.LRLat">-13.22909668</MDI>
    <MDI key="BAND_RE.LRLon">-71.57710406</MDI>
    <MDI key="BAND_RE.TDILevel">24</MDI>
    <MDI key="BAND_RE.ULHAE">2332.70</MDI>
    <MDI key="BAND_RE.ULLat">-13.04555359</MDI>
    <MDI key="BAND_RE.ULLon">-71.63637407</MDI>
    <MDI key="BAND_RE.URHAE">2332.70</MDI>
    <MDI key="BAND_RE.URLat">-13.04616980</MDI>
    <MDI key="BAND_RE.URLon">-71.57519527</MDI>
    <MDI key="BAND_Y.absCalFactor">5.101088e-03</MDI>
    <MDI key="BAND_Y.effectiveBandwidth">3.740000e-02</MDI>
    <MDI key="BAND_Y.LLHAE">2332.70</MDI>
    <MDI key="BAND_Y.LLLat">-13.22847154</MDI>
    <MDI key="BAND_Y.LLLon">-71.63832808</MDI>
    <MDI key="BAND_Y.LRHAE">2332.70</MDI>
    <MDI key="BAND_Y.LRLat">-13.22909668</MDI>
    <MDI key="BAND_Y.LRLon">-71.57710406</MDI>
    <MDI key="BAND_Y.TDILevel">24</MDI>
    <MDI key="BAND_Y.ULHAE">2332.70</MDI>
    <MDI key="BAND_Y.ULLat">-13.04555359</MDI>
    <MDI key="BAND_Y.ULLon">-71.63637407</MDI>
    <MDI key="BAND_Y.URHAE">2332.70</MDI>
    <MDI key="BAND_Y.URLat">-13.04616980</MDI>
    <MDI key="BAND_Y.URLon">-71.57519527</MDI>
    <MDI key="bitsPerPixel">16</MDI>
    <MDI key="compressionType">"None"</MDI>
    <MDI key="generationTime">2019-06-05T09:20:20.000000Z</MDI>
    <MDI key="imageDescriptor">"ORStandard2A"</MDI>
    <MDI key="IMAGE_1.attitudeKnowledgeSrc">"R"</MDI>
    <MDI key="IMAGE_1.avgLineRate">5000.00</MDI>
    <MDI key="IMAGE_1.CatId">"103001001A50C800"</MDI>
    <MDI key="IMAGE_1.cloudCover">0.000</MDI>
    <MDI key="IMAGE_1.colUncertainty">522.53</MDI>
    <MDI key="IMAGE_1.exposureDuration">0.00020000</MDI>
    <MDI key="IMAGE_1.firstLineTime">2012-07-11T15:34:18.154272Z</MDI>
    <MDI key="IMAGE_1.maxCollectedColGSD">1.957</MDI>
    <MDI key="IMAGE_1.maxCollectedRowGSD">1.939</MDI>
    <MDI key="IMAGE_1.maxCrossTrackViewAngle">-9.0</MDI>
    <MDI key="IMAGE_1.maxInTrackViewAngle">-9.3</MDI>
    <MDI key="IMAGE_1.maxOffNadirViewAngle">13.3</MDI>
    <MDI key="IMAGE_1.maxSatAz">233.8</MDI>
    <MDI key="IMAGE_1.maxSatEl">75.7</MDI>
    <MDI key="IMAGE_1.maxSunAz">28.8</MDI>
    <MDI key="IMAGE_1.maxSunEl">50.2</MDI>
    <MDI key="IMAGE_1.meanCollectedColGSD">1.957</MDI>
    <MDI key="IMAGE_1.meanCollectedGSD">1.947</MDI>
    <MDI key="IMAGE_1.meanCollectedRowGSD">1.937</MDI>
    <MDI key="IMAGE_1.meanCrossTrackViewAngle">-9.2</MDI>
    <MDI key="IMAGE_1.meanInTrackViewAngle">-9.5</MDI>
    <MDI key="IMAGE_1.meanOffNadirViewAngle">13.2</MDI>
    <MDI key="IMAGE_1.meanSatAz">232.2</MDI>
    <MDI key="IMAGE_1.meanSatEl">75.2</MDI>
    <MDI key="IMAGE_1.meanSunAz">28.7</MDI>
    <MDI key="IMAGE_1.meanSunEl">50.1</MDI>
    <MDI key="IMAGE_1.minCollectedColGSD">1.956</MDI>
    <MDI key="IMAGE_1.minCollectedRowGSD">1.934</MDI>
    <MDI key="IMAGE_1.minCrossTrackViewAngle">-9.3</MDI>
    <MDI key="IMAGE_1.minInTrackViewAngle">-9.7</MDI>
    <MDI key="IMAGE_1.minOffNadirViewAngle">13.1</MDI>
    <MDI key="IMAGE_1.minSatAz">230.9</MDI>
    <MDI key="IMAGE_1.minSatEl">74.8</MDI>
    <MDI key="IMAGE_1.minSunAz">28.7</MDI>
    <MDI key="IMAGE_1.minSunEl">50.0</MDI>
    <MDI key="IMAGE_1.mode">"FullSwath"</MDI>
    <MDI key="IMAGE_1.PNIIRS">3.0</MDI>
    <MDI key="IMAGE_1.positionKnowledgeSrc">"R"</MDI>
    <MDI key="IMAGE_1.resamplingKernel">"MTF"</MDI>
    <MDI key="IMAGE_1.revNumber">14463</MDI>
    <MDI key="IMAGE_1.rowUncertainty">406.46</MDI>
    <MDI key="IMAGE_1.satId">"WV02"</MDI>
    <MDI key="IMAGE_1.scanDirection">"Reverse"</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.colSpacing">2.00</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.datumName">"WE"</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.datumOffset">(0.000,0.000,0.000)</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.DEMCorrection">"Base Elevation"</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.earliestAcqTime">2012-07-11T15:34:18.841772Z</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.inverseFlattening">298.257223563</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.latestAcqTime">2012-07-11T15:34:18.841772Z</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.LLH">2332.70</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.LLX">214068.99997304</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.LLY">8536091.00000158</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.LRH">2332.70</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.LRX">220708.99997728</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.LRY">8536091.00000133</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.mapHemi">"S"</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.mapProjCode">1</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.mapProjName">"UTM"</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.mapProjParam">(0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000,0.000000000)</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.mapZone">19</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.numGCP">0</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.orientationAngle">0.0</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.originX">214068.99997402</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.originY">8556341.00000157</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.productGSD">2.00</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.productUnits">"M"</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.rowSpacing">2.00</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.semiMajorAxis">6378137.0000</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.terrainHae">2332.70</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.ULH">2332.70</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.ULX">214068.99997402</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.ULY">8556341.00000157</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.URH">2332.70</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.URX">220708.99997812</MDI>
    <MDI key="MAP_PROJECTED_PRODUCT.URY">8556341.00000131</MDI>
    <MDI key="numberOfLooks">1</MDI>
    <MDI key="numColumns">3321</MDI>
    <MDI key="numRows">10126</MDI>
    <MDI key="outputFormat">"GeoTIFF"</MDI>
    <MDI key="panSharpenAlgorithm">"None"</MDI>
    <MDI key="productCatalogId">"None"</MDI>
    <MDI key="productLevel">"LV2A"</MDI>
    <MDI key="productOrderId">"010217881010_01_P002"</MDI>
    <MDI key="productType">"Standard"</MDI>
    <MDI key="radiometricEnhancement">"Off"</MDI>
    <MDI key="radiometricLevel">"Corrected"</MDI>
    <MDI key="version">"28.4"</MDI>
  </Metadata>
  <Metadata domain="RPC">
    <MDI key="HEIGHT_OFF">2899</MDI>
    <MDI key="HEIGHT_SCALE">1110</MDI>
    <MDI key="LAT_OFF">-13.1387</MDI>
    <MDI key="LAT_SCALE">0.0961</MDI>
    <MDI key="LINE_DEN_COEFF">+1.000000E+00 +1.340142E-05 +2.280973E-04 -2.590785E-04 +1.004087E-08 +4.063676E-07 -2.241707E-07 -1.301510E-07 +4.136808E-08 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +4.875542E-08 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00</MDI>
    <MDI key="LINE_NUM_COEFF">+4.580781E-04 -3.815467E-03 -1.014704E+00 -1.738947E-02 +1.184011E-05 -2.643526E-05 -2.630275E-04 +2.588546E-05 -2.288126E-04 -2.970462E-06 +4.987988E-07 +0.000000E+00 +1.943760E-08 -3.950414E-08 -4.044091E-08 -5.194577E-08 +0.000000E+00 -6.440836E-08 -5.186995E-07 +0.000000E+00</MDI>
    <MDI key="LINE_OFF">5090</MDI>
    <MDI key="LINE_SCALE">5243</MDI>
    <MDI key="LONG_OFF">-71.6073</MDI>
    <MDI key="LONG_SCALE">0.0355</MDI>
    <MDI key="SAMP_DEN_COEFF">+1.000000E+00 -3.979801E-04 -1.747896E-04 -8.680570E-04 +6.207230E-08 -1.472495E-08 +2.592397E-07 -7.160463E-08 +5.117194E-07 -7.848415E-07 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00</MDI>
    <MDI key="SAMP_NUM_COEFF">+7.701323E-04 +1.034443E+00 -2.990413E-02 +6.173273E-02 +2.348054E-04 +7.720960E-04 -6.760084E-04 -4.246901E-04 +1.116831E-04 +2.606434E-05 -3.500904E-08 +0.000000E+00 -1.035875E-06 +1.702233E-07 -1.193321E-07 -1.829311E-08 -7.570366E-07 -3.600665E-07 +4.359219E-07 -2.843736E-08</MDI>
    <MDI key="SAMP_OFF">1688</MDI>
    <MDI key="SAMP_SCALE">1863</MDI>
  </Metadata>
  <Metadata domain="IMAGERY">
    <MDI key="ACQUISITIONDATETIME">2012-07-11 15:34:18</MDI>
    <MDI key="CLOUDCOVER">0</MDI>
    <MDI key="SATELLITEID">WV02</MDI>
  </Metadata>
  <Metadata>
    <MDI key="METADATATYPE">DG</MDI>
  </Metadata>
  <Metadata domain="Esri">
    <MDI key="PyramidResamplingType">NEAREST</MDI>
  </Metadata>
</PAMDataset>

[quinten]

Hi Thang

If the lake is at high altitude, this needs to be taken into account in the atmospheric correction. This is currently not supported for WV, but I can add that in the near future. I'll inform you here if I have made the edits to the code base.

Quinten