"Cannot Construct Databuffer" error despite 200Gb available RAM
Posted: Mon Nov 20, 2017 1:31 pm America/New_York
Using VIIRS data for 16 channels, I'm trying to run l2gen, followed by gpt reproject, and gpt mosaic. I'm able to run l2gen but after that, I've been getting the "Cannot Construct Databuffer" error. I allocated 200Gb of memory in the seadas/bin/seadas.vmoptions file (by specifying 200000m) but I'm still receiving the same error. This should be enough memory for this, shouldn't it? Is there an issue with my reproject xml?
I'm creating the SVM files using viirs_RDR and I'm able to open these files in seadas without any issues. I'm also able to run reproject and mosaic in seadas but would much prefer to do it from the command line.
**Workflow**
l2gen ifile=SVM01.h5 geofile=GMTCO.h5 ofile=SV_L2_file l2prod1="sst" resolution=375 ctl_pt_incr=1 ctl_pt_incr=1 proc_ocean=1 proc_sst=1 proc_land=0 atmocor=1 maskcloud=1 maskland=1 maskhilt=0 maskstlight=0
gpt.sh -e reproject.sst.xml -f netCDF4-CF -t sst.reproject.nc SV_L2_file
gpt.sh -e mosaic.sst.xml -f netCDF4-CF -t sst.mosaic.nc sst.reproject.nc
**reproject.sst.xml**
<graph id="bl2map_mosaic">
<version>1.0</version>
<node id="bandMath">
<operator>BandMaths</operator>
<sources>
<sourceProduct>${sourceProducts}</sourceProduct>
</sources>
<parameters>
<targetBands>
<targetBand>
<name>sst</name>
<expression>sst</expression>
<description>sst</description>
<type>float32</type>
<noDataValue>NaN</noDataValue>
<spectralBandIndex>0</spectralBandIndex>
</targetBand>
</targetBands>
</parameters>
</node>
<node id="reproject">
<operator>Reproject</operator>
<sources>
<source>bandMath</source>
</sources>
<parameters>
<crs>PROJCS["Mercator_1SP / World Geodetic System 1984",
GEOGCS["World Geodetic System 1984",
DATUM["World Geodetic System 1984",
SPHEROID["WGS 84", 6378135.0, 298.257223563, AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST],
AXIS["Geodetic latitude", NORTH]],
PROJECTION["Mercator_1SP"],
PARAMETER["latitude_of_origin", 0.0],
PARAMETER["central_meridian", 0.0],
PARAMETER["scale_factor", 1.0],
PARAMETER["false_easting", 0.0],
PARAMETER["false_northing", 0.0],
UNIT["m", 1.0],
AXIS["Easting", EAST],
AXIS["Northing", NORTH]]</crs>
<resampling>Nearest</resampling>
<pixelSizeX>4000</pixelSizeX>
<pixelSizeY>4000</pixelSizeY>
<orthorectify>false</orthorectify>
<noDataValue>NaN</noDataValue>
</parameters>
</node>
</graph>
Thanks!
I'm creating the SVM files using viirs_RDR and I'm able to open these files in seadas without any issues. I'm also able to run reproject and mosaic in seadas but would much prefer to do it from the command line.
**Workflow**
l2gen ifile=SVM01.h5 geofile=GMTCO.h5 ofile=SV_L2_file l2prod1="sst" resolution=375 ctl_pt_incr=1 ctl_pt_incr=1 proc_ocean=1 proc_sst=1 proc_land=0 atmocor=1 maskcloud=1 maskland=1 maskhilt=0 maskstlight=0
gpt.sh -e reproject.sst.xml -f netCDF4-CF -t sst.reproject.nc SV_L2_file
gpt.sh -e mosaic.sst.xml -f netCDF4-CF -t sst.mosaic.nc sst.reproject.nc
**reproject.sst.xml**
<graph id="bl2map_mosaic">
<version>1.0</version>
<node id="bandMath">
<operator>BandMaths</operator>
<sources>
<sourceProduct>${sourceProducts}</sourceProduct>
</sources>
<parameters>
<targetBands>
<targetBand>
<name>sst</name>
<expression>sst</expression>
<description>sst</description>
<type>float32</type>
<noDataValue>NaN</noDataValue>
<spectralBandIndex>0</spectralBandIndex>
</targetBand>
</targetBands>
</parameters>
</node>
<node id="reproject">
<operator>Reproject</operator>
<sources>
<source>bandMath</source>
</sources>
<parameters>
<crs>PROJCS["Mercator_1SP / World Geodetic System 1984",
GEOGCS["World Geodetic System 1984",
DATUM["World Geodetic System 1984",
SPHEROID["WGS 84", 6378135.0, 298.257223563, AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST],
AXIS["Geodetic latitude", NORTH]],
PROJECTION["Mercator_1SP"],
PARAMETER["latitude_of_origin", 0.0],
PARAMETER["central_meridian", 0.0],
PARAMETER["scale_factor", 1.0],
PARAMETER["false_easting", 0.0],
PARAMETER["false_northing", 0.0],
UNIT["m", 1.0],
AXIS["Easting", EAST],
AXIS["Northing", NORTH]]</crs>
<resampling>Nearest</resampling>
<pixelSizeX>4000</pixelSizeX>
<pixelSizeY>4000</pixelSizeY>
<orthorectify>false</orthorectify>
<noDataValue>NaN</noDataValue>
</parameters>
</node>
</graph>
Thanks!