Hi Helen,
I can help you get this working. It looks like NASA finally reprocessed the CZCS data to be consistent with the rest of the ocean color products. Previously it
was stored in a different way that required us to special case some things about it in the tool. Now that is no longer necessary, but the tool is doing it anyway, which prevents it from working with the new format.
I fix our code and send you an updated file in a little while. When you try it, you should use CHL_chlor_a as the product code, although the documentation for
the tool will not reflect this (we can only update the documentation in a new release of MGET). Also, CZCS is only available at 9km resolution, so you should choose that.
More soon,
Jason
From: helen perivier [mailto:]
Sent: Thursday, February 05, 2015 6:28 AM
To:
Subject: Fwd:
(please ignore previous message. I hit the send button by mistake)
I am new to MGET and making use of its great tools to plug into an Arc-GIS ocean project. Thanks for providing this to the public.
I am writing now to call attention to a possible glitch in accessing the CZCS data from the NASA Ocean Color L3 SMI site with the Data Products tool set, either to interpolate at points or to create rasters.
The tool doesn't find the data and I think the problem could be that the source may have changed its data organization and product codes. CHLO comes up empty. I also tried alternative codes CHLOR and CHL_chlor_a because some of the CZCS files have those titles.
My sciprt got fartheest with the CHL_chlor_a product code.
Any suggestions would be great. I am working with data from the late 1970s and 1980s and it seems the CZCS chlorophyll data is the only set going back that far. I'm trying to download and apply the data into
Arc directly but so far without a lot of success.
Here is one example of the script I get back... (n this case using CHLO as the product code)
Executing: OceanColorLevel3SMITimeSeriesInterpolateAtArcGISPoints CZCS 8day 4km CHLO Extract_Stn_Lpt2 Chl_CZCS Min_of_Date Nearest # # # 60 300 # 256 16 16 3
Start Time: Thu Feb 05 10:55:47 2015
Running script OceanColorLevel3SMITimeSeriesInterpolateAtArcGISPoints...
CollectionIsEmptyError: The NASA OceanColor data archive contains no datasets matching the _expression_ Sensor = 'czcs' AND TemporalResolution = '8day' AND SpatialResolution = '4km' AND ProductCode = 'chlo'.
Failed script OceanColorLevel3SMITimeSeriesInterpolateAtArcGISPoints...
Traceback (most recent call last):
File "C:\Program Files\GeoEco\ArcGISToolbox\Scripts\OceanColorLevel3SMITimeSeriesInterpolateAtArcGISPoints.py", line 5, in <module>
ExecuteMethodFromCommandLineAsArcGISTool('GeoEco.DataProducts.NASA.OceanColor', 'OceanColorLevel3SMITimeSeries', 'InterpolateAtArcGISPoints')
File "C:\Python27\ArcGIS10.2\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.2\lib\site-packages\GeoEco\DataProducts\NASA\OceanColor.py", line 615, in InterpolateAtArcGISPoints
grid = cls(sensor, temporalResolution, spatialResolution, product, timeout, maxRetryTime, cacheDirectory)
File "C:\Python27\ArcGIS10.2\lib\site-packages\GeoEco\DataProducts\NASA\OceanColor.py", line 527, in __init__
super(OceanColorLevel3SMITimeSeries, self).__init__(collection, _expression_=_expression_, reportProgress=False)
File "C:\Python27\ArcGIS10.2\lib\site-packages\GeoEco\Datasets\Virtual.py", line 65, in __init__
raise CollectionIsEmptyError(collection.DisplayName, _expression_)
CollectionIsEmptyError: The NASA OceanColor data archive contains no datasets matching the _expression_ Sensor = 'czcs' AND TemporalResolution = '8day' AND SpatialResolution = '4km' AND ProductCode = 'chlo'.
Failed to execute (OceanColorLevel3SMITimeSeriesInterpolateAtArcGISPoints).
Failed at Thu Feb 05 10:55:50 2015 (Elapsed Time: 2.96 seconds)
|