Subject: Marine Geospatial Ecology Tools (MGET) help
Text archives
| From: | Yao Zhao <> |
|---|---|
| To: | |
| Subject: | [mget-help] Questions with Cayula-Cornillon Fronts Detection & HYCOM Data Download |
| Date: | Mon, 23 Aug 2010 13:11:39 -0400 |
Executing: CayulaCornillonEdgeDetectionDetectEdgesInArcGISRaster C:\yzwork\test\sstraster C:\yzwork\test\edge false int(inputRaster) 3 32 16 0.65 0.25 3 0.76 0.9 0.92 1 # # # # #
Start Time: Mon Aug 23 10:34:23 2010
Running script CayulaCornillonEdgeDetectionDetectEdgesInArcGISRaster...
Finding fronts in C:\yzwork\test\sstraster using parameters: medianFilter=3, histogramWindow=32, histogramStride=16, minPropCellsWithData=0.65000000000000002, minPopProp=0.25, minPopMeanDiff=3.0, minTheta=0.76000000000000001, minSinglePopCohesion=0.90000000000000002, minGlobalPopCohesion=0.92000000000000004, threads=1.
Output rasters: fronts=C:\yzwork\test\edge, filteredImage=None, candidateCounts=None, frontCounts=None, windowStatusCodes=None, windowStatusValues=None
KeyError: '_gdal_array_1_0_3'
<type 'exceptions.KeyError'>: '_gdal_array_1_0_3'
Failed to execute (CayulaCornillonEdgeDetectionDetectEdgesInArcGISRaster).
End Time: Mon Aug 23 10:34:30 2010 (Elapsed Time: 7.00 seconds)
A similar error came out when I used the MGET tool to download HYCOM GOMl0.04 3D Variable to ArcGIS Rasters. I selected ssh variable and limited the time series from 4/1/2010 to 4/30/2010. Error message is like this:
Executing: HYCOMGOMl0043DCreateArcGISRasters ssh C:\yzwork\test\HYCOM Add %(VariableName)s;%%Y;%(VariableName)s_%%Y%%j.img # Degrees 4/1/2010 4/30/2010 1000 # true false C:\yzwork\test\HYCOM
Start Time: Mon Aug 23 11:37:43 2010
Running script HYCOMGOMl0043DCreateArcGISRasters...
RuntimeError: Failed to write a block of data of 541 columns by 385 rows at offsets x=0, y=0 to band 1 of band 1 of GDAL dataset "C:\yzwork\test\HYCOM\ssh\2010\ssh_2010091.img" with the Geospatial Data Abstraction Library (GDAL). Verify that the dataset exists, is accessible, and has the expected dimensions. Detailed error information: band.WriteArray(data, 0, 0) reported KeyError: '_gdal_array_1_0_3'.
<type 'exceptions.RuntimeError'>: Failed to write a block of data of 541 columns by 385 rows at offsets x=0, y=0 to band 1 of band 1 of GDAL dataset "C:\yzwork\test\HYCOM\ssh\2010\ssh_2010091.img" with the Geospatial Data Abstraction Library (GDAL). Verify that the dataset exists, is accessible, and has the expected dimensions. Detailed error information: band.WriteArray(data, 0, 0) reported KeyError: '_gdal_array_1_0_3'.
Failed to execute (HYCOMGOMl0043DCreateArcGISRasters).
End Time: Mon Aug 23 11:38:18 2010 (Elapsed Time: 35.00 seconds)
I also tried to repeat the example "Extracting ArcGIS rasters from the HYCOM global ocean model's 4D netCDFs". I followed the step-by-step instructions. After inputing the command "Extract_Rasters_From_HYCOM_Global_NetCDFs.py NetCDFs\archv.2003_002_00_3zt.nc temperature regional.grid.a regional.grid.b Rasters 0" in CMD shell, there is an ImportError which said "No module named netCDF3". So I just wonder if there is something that I failed to install on my computer.