Skip to Content.

mget-help - Re: [mget-help] trouble with Data Product/Nasa Ocean Color/L3

Please Wait...

Subject: Marine Geospatial Ecology Tools (MGET) help

Text archives


Re: [mget-help] trouble with Data Product/Nasa Ocean Color/L3


Chronological Thread 
  • From: Sophie Laran <>
  • To: Jason Roberts <>, "White, George" <>
  • Cc: "''" <>
  • Subject: Re: [mget-help] trouble with Data Product/Nasa Ocean Color/L3
  • Date: Wed, 29 Jan 2020 17:08:03 +0100
  • Authentication-results: mail-gw.oit.duke.edu; spf=pass ; dkim=pass header.d=univ-lr.fr header.s=default; dkim=pass header.d=univ-lr.fr header.s=default; dmarc=pass header.from=univ-lr.fr

Hi all and thanks a lot for your quick fit back and all the explanations. I'm not good to understand error message. It's more clear now.

Thanks Jason also, but I've check with my colleagues and they can wait until March for the data I have to extract, so I will follow the news on the web site. I'm taking the opportunity for tree more questions:

-I'm using ArcGIS 10.6 and most of the time the MGET toolbox desapear when I'm closing ArcCatalog

-do you have a mailing list or something to inform people of update or something like that?

- and did you plan to addapt MGET for ArcGISPro ?

Thanks a lot anycas for this wonderfull tools that is MGET


Best regards

Sophie


Le 27/01/2020 à 22:38, Jason Roberts a écrit :

Thanks, George, I can confirm that is the problem. Your post is very helpful. I was not aware of this change; I have not been keeping up on everything happening at NASA OBPG (sorry about that).

 

In order for this to work again, it is probably necessary for me to modify MGET to accept a NASA Earthdata username and password and plumb it through MGET’s HTTP stack. MGET does not use curl, at least not directly, but the Python API it relies on does allow for this sort of thing, so it shouldn’t be too complicated.

 

I was currently targeting a new version of MGET for March. If I stuck to that schedule, MGET’s ability to download NASA OceanColor data would be broken until then. This may be problematic for people. Please contact me privately if you are actively using MGET to access OceanColor data. If people are impacted, I’ll see if I can accelerate the work.

 

Best,

Jason

 

From: <> On Behalf Of White, George
Sent: Monday, January 27, 2020 1:44 PM
To:
Cc: Sophie Laran <>
Subject: RE: [mget-help] trouble with Data Product/Nasa Ocean Color/L3

 

Have you seen the recent announcement (https://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?tid=11280). EOSDIS username and password are now required to download most of the Ocean Color data products?   The error you get:

 

RuntimeError: Failed to open netCDF file c:\users\slaran\appdata\local\temp\arc5df1\GeoEco_OceanColorLevel3SMIFileSearcher_Temp_mzv5vr\A2002185.L3m_DAY_POC_poc_4km.nc (a local copy of A2002185.L3m_DAY_POC_poc_4km.nc from NASA OceanColor data archive). The file may not be in netCDF-3 or netCDF-4 format. Detailed error information: netCDF4.Dataset() reported RuntimeError: NetCDF: Unknown file format.

may result from a login failure.  This often produces a “.nc” file which contains the html EODIS login page.    You can check this by replacing “.nc” with “.html” and trying to view the page.

Many automated downloaders (including some of the Python scripts in the OCSSW Processing System) have been failing since the change.   The NASA Ocean Color Forum has several discussions, one for shell scripts using wget or curl and another for downloads of ancillary files needed for local processing, which is done using Python scripts.    For Python scripts that fail, one approach is to generate a list of the URL’s to download and then using wget or curl to download the files from the list.    New versions of the OCSSW Python scripts were released last week.   Wget and curl can store the credentials in a “.netrc” file.  If access to this file is restricted to the user, then it is not necessary to enter usernames or passwords on the command line

or hard code  them in scripts.

 

A couple changes to existing scripts that have been useful:

 

1.       For URL’s of the form https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/ using "ob" in place of "cgi" shortcuts one redirection 2.       Some data sets need the wget option “--auth-no-challenge=on” (see the wget man page) 3.       Some remote sites need to use IPv4, and fail if IPv6 is enabled (I think this is due to some intermediate site not using IPv6. DNS provides the IPv6 address but downloads fail and debugging gives an error to the effect that IPv6 is not available).  Use “wget -4” to disable IPv6.

 

All this a corollary of the rule that security is a pain in the neck (and sometimes even lower).

 

George N. White III 

 

From: <> On Behalf Of Sophie Laran
Sent: Monday, January 27, 2020 4:42 AM
To:
Subject: [mget-help] trouble with Data Product/Nasa Ocean Color/L3

 

Hello

I'm using MGET on ArcGis and I recently experiment some troubles to use

Create a Climatological raster for NASA oceanColor Lm SMI

It was runing few weeks ago and now I can't, and receive this kind of message below this e-mail.

Could you let me know if it's come from me/my computer or the tool. It's seems that none of the tools from L3 Product work.

Thanks a lot in advance for your answer and for this very nice Tools, so helpfull!!

Best regards

Sophie Laran

Executing (Create Climatological Rasters for NASA OceanColor L3 SMI Product (3)): OceanColorLevel3SMITimeSeriesCreateClimatologicalArcGISRasters Aqua Daily 4km POC_poc Mean Cumulative C:\Users\slaran\Documents\2_SIG\SPEEcov\POC Replace %(Sensor)s;%(TemporalResolution)s;%(SpatialResolution)s;%(ProductCode)s;%(ClimatologyBinType)s_Climatology;%(Sensor)s_%(ProductCode)s_%(ClimatologyBinName)s_%(Statistic)s.img # 1 1 # "-5 44.5 -0.6 47" 18/01/2019 15/02/2019 60 300 true false
Start Time: Mon Jan 27 09:23:24 2020
Running script OceanColorLevel3SMITimeSeriesCreateClimatologicalArcGISRasters...
Querying cumulative climatological mean of the aqua daily 4km POC_poc from NASA OceanColor, clipped to indices tMin = 6042, tMax = 6070, yMin = 3228, yMax = 3287, xMin = 4200, xMax = 4305.
Query complete: 0:00:00 elapsed, 1 datasets found, 0:00:00.050000 per dataset.
Importing 1 datasets into ArcGIS Folder C:\Users\slaran\Documents\2_SIG\SPEEcov\POC with mode "replace".
Import stopped before all datasets were imported: 0:00:01 elapsed, 0 datasets imported, 0:00:00 per dataset, 1 datasets not imported.
RuntimeError: Failed to open netCDF file c:\users\slaran\appdata\local\temp\arc5df1\GeoEco_OceanColorLevel3SMIFileSearcher_Temp_mzv5vr\A2002185.L3m_DAY_POC_poc_4km.nc (a local copy of A2002185.L3m_DAY_POC_poc_4km.nc from NASA OceanColor data archive). The file may not be in netCDF-3 or netCDF-4 format. Detailed error information: netCDF4.Dataset() reported RuntimeError: NetCDF: Unknown file format.
Failed script OceanColorLevel3SMITimeSeriesCreateClimatologicalArcGISRasters...

Traceback (most recent call last):
  File "C:\Program Files\GeoEco\ArcGISToolbox\Scripts\OceanColorLevel3SMITimeSeriesCreateClimatologicalArcGISRasters.py", line 5, in <module>
    ExecuteMethodFromCommandLineAsArcGISTool('GeoEco.DataProducts.NASA.OceanColor', 'OceanColorLevel3SMITimeSeries', 'CreateClimatologicalArcGISRasters')
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\ArcGISScripts.py", line 210, in ExecuteMethodFromCommandLineAsArcGISTool
    exec sourceCode in globals(), locals()
  File "<string>", line 1, in <module>
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\DataProducts\NASA\OceanColor.py", line 595, in CreateClimatologicalArcGISRasters
    workspace.ImportDatasets(collection.QueryDatasets(), mode, calculateStatistics=calculateStatistics, buildPyramids=buildPyramids)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 820, in ImportDatasets
    self._ImportDatasets(datasets, mode.lower(), reportProgress, options)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\Collections.py", line 695, in _ImportDatasets
    self._ImportDatasetsToPath(pathComponentsForPath[path], datasetsForPath[path], mode, progressReporter, options)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\ArcGIS.py", line 612, in _ImportDatasetsToPath
    self.DatasetType._ImportDatasetsToPath(os.path.join(self.Path, *pathComponents), sourceDatasets, mode, progressReporter, options)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\ArcGIS.py", line 1036, in _ImportDatasetsToPath
    elif dataset.DataType not in ArcGISRaster._GDALDataTypeForNumpyDataType.keys():
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 2673, in _GetDataType
    return self.UnscaledDataType
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 2697, in _GetUnscaledDataType
    return self.GetLazyPropertyValue('UnscaledDataType')
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 160, in GetLazyPropertyValue
    value = self._GetLazyPropertyPhysicalValue(name)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\Virtual.py", line 1970, in _GetLazyPropertyPhysicalValue
    if self._Grids[0].DataType == u'float64':
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 2671, in _GetDataType
    if self.DataIsScaled:
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 2692, in _GetDataIsScaled
    return self.GetLazyPropertyValue('ScalingFunction') is not None
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 160, in GetLazyPropertyValue
    value = self._GetLazyPropertyPhysicalValue(name)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\Virtual.py", line 563, in _GetLazyPropertyPhysicalValue
    return self._Grid.GetLazyPropertyValue(name)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 160, in GetLazyPropertyValue
    value = self._GetLazyPropertyPhysicalValue(name)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\Virtual.py", line 913, in _GetLazyPropertyPhysicalValue
    return self._Grid.GetLazyPropertyValue(name)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 160, in GetLazyPropertyValue
    value = self._GetLazyPropertyPhysicalValue(name)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\Virtual.py", line 277, in _GetLazyPropertyPhysicalValue
    return self._CachedOldestGrid.GetLazyPropertyValue(name)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\__init__.py", line 160, in GetLazyPropertyValue
    value = self._GetLazyPropertyPhysicalValue(name)
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\NetCDF.py", line 348, in _GetLazyPropertyPhysicalValue
    v = self._GetVariable()
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\NetCDF.py", line 723, in _GetVariable
    self.ParentCollection._Open()
  File "C:\Python27\ArcGIS10.6\lib\site-packages\GeoEco\Datasets\NetCDF.py", line 288, in _Open
    raise RuntimeError(_(u'Failed to open %(dn)s. The file may not be in netCDF-3 or netCDF-4 format. Detailed error information: netCDF4.Dataset() reported %(e)s: %(msg)s.') % {u'dn': self._DisplayName, u'e': e.__class__.__name__, u'msg': self._Str(e)})
RuntimeError: Failed to open netCDF file c:\users\slaran\appdata\local\temp\arc5df1\GeoEco_OceanColorLevel3SMIFileSearcher_Temp_mzv5vr\A2002185.L3m_DAY_POC_poc_4km.nc (a local copy of A2002185.L3m_DAY_POC_poc_4km.nc from NASA OceanColor data archive). The file may not be in netCDF-3 or netCDF-4 format. Detailed error information: netCDF4.Dataset() reported RuntimeError: NetCDF: Unknown file format.


Failed to execute (Create Climatological Rasters for NASA OceanColor L3 SMI Product (3)).
Failed at Mon Jan 27 09:23:37 2020 (Elapsed Time: 12,89 seconds)

 

-- 
 
Sophie Laran (PhD)
 
ADERA
 
Observatoire PELAGIS - UMS 3462
 
Université de La Rochelle - CNRS
 
Pôle Analytique - 5 allées de l'Océan
 
17 000 LA ROCHELLE - FRANCE
 
TEL : +33 (0)5 16 49 67 15 (LD)
 
TEL : +33 (0)5 46 44 99 10 (STD)
 
http://www.observatoire-pelagis.cnrs.fr
-- 
Sophie Laran (PhD)
ADERA
Observatoire PELAGIS - UMS 3462
Université de La Rochelle - CNRS
Pôle Analytique - 5 allées de l'Océan
17 000 LA ROCHELLE - FRANCE
TEL : +33 (0)5 16 49 67 15 (LD)
TEL : +33 (0)5 46 44 99 10 (STD)
http://www.observatoire-pelagis.cnrs.fr



Archive powered by MHonArc 2.6.19.

Top of Page