Conversion multiple Raster to NetCDF using Arcpy

[ShouvikJha]

Dear all,

I have multiples raster files, i am trying to convert it in NetCDF file using Arcpy script . below script i am using for conversion. While  i ran the below script , i am getting the error, 

Below i have attached my raster file also, which i trying to convert in NetCDF format. 

Error Code 

ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000840: The value is not a Raster Layer.
ERROR 000840: The value is not a Raster Catalog.
Failed to execute (RasterToNetCDF).‍‍‍‍‍‍‍‍

Below script i am using for conversion of raster to NetCDF. 

import arcpy
from arcpy import env

# Set environment settings
env.workspace = "D:\2012A"

# Set local variables
inRaster = "D:\2012A"
outNetCDFFile = "D:\2012A/nppnetcdf.nc"
variable = "NPP"
units = "gCSqm"
XDimension = "x"
YDimension = "y"
bandDimension = ""

# Process: RasterToNetCDF
arcpy.RasterToNetCDF_md(inRaster, outNetCDFFile, variable, units,
                        XDimension, YDimension, bandDimension)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

[DanPatterson_Retired]

you have set your environment workspace (ie a folder or a geodatabase) as D:\2012A  then you reassigned it as a raster.  This won't work.  And you have to start using raw formatting in any situation ... with the 'r' prefix to paths

r"c:\yourfolder\yoursubfolder\yourgdb.gdb"  for example

[ShouvikJha]

Dan Patterson‌, Thank you. I have updated my code by the add of r comment. and my raster file in .tif format. but unfortunately i am getting the error massage, below is error code, 

can you please explain, why i have to add .gdb  and is it possible to run without  gdb file ?

Message     File Name     Line     Position     
Traceback                    
    <module>     <module1>     19          
    RasterToNetCDF     C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\md.py     253          
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000840: The value is not a Raster Layer.
ERROR 000840: The value is not a Raster Catalog.
Failed to execute (RasterToNetCDF).
                    
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

[ShouvikJha]

Xander Bakker‌, I have posted one of my problem, Please cooperate with us to solve the issues.  

[DanPatterson_Retired]

line 5 and 8 are wrong... read my post, putting an r in front of those won't work. Post your new code

[XanderBakker]

As Dan Patterson mentions, revise the definition of the input data. If the "raster" is indeed a raster (Esri grid format with the name "2012A" in the root of drive D) then you should rename the raster (since you should not start the name of an Esri grid with a number).

Another general advise; run the tool manually with the data, check the result and if correct, copy the Python snippet to verify the syntax of the python code that you should use. 

[ShouvikJha]

Dan Patterson‌ Xander Bakker‌, Thank you very much all of you. I have updated my code, still problem persist, and and my raster name , r001_NPP.tif, r002_NPP.tif ........... r012_NPP.tif , and i checked one raster to converted in  NetCDF using Arc GIS tool, i got perfect output from there, but i want to convert by the script , i have many of raster to convert. 

# Import system modules
import arcpy
from arcpy import env

# Set environment settings
env.workspace = r"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA"

# Set local variables
inRaster = r"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NPP_2012A"
outNetCDFFile = r"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NC_FILE\nppnetcdf.nc"
variable = "NPP"
units = "gcsqm"
XDimension = "x"
YDimension = "y"
bandDimension = ""

# Process: RasterToNetCDF
arcpy.RasterToNetCDF_md(inRaster, outNetCDFFile, variable, units,
                        XDimension, YDimension, bandDimension)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

[IanMurray]

your inRaster variable is currently a folder I believe and not the actual raster file.  The input for the tool has to be a raster file(in your case it looks like .tif).  If you have multiple rasters in that directory, you will need to make a list of them all and iterate through them, using each individual raster as input for the Raster to NetCDF Tool(See the link about the ListRasters).   Please also note you will need unique output names for each output NetCDF file, or you will just be overwriting the prior file.

http://pro.arcgis.com/en/pro-app/arcpy/functions/listrasters.htm

[IanMurray]

Just as a note, I would make sure to run the ListRasters and print all the raster names as the example states, to make sure that the tool is using the correct files for input into the Raster to NetCDF.  Also if the rasters are in the NPP_2012A folder as I think they are from your code, you will need to make that your workspace when you use ListRasters.

[ShouvikJha]

Dan Patterson, @Lan Murray, Xander Bakker  i have tried using below command , but wont understand how to complete this code, and i want to take all 12 rasters in time series NetCDF file. Please guide on same 

# Import system modules
import arcpy
from arcpy import env

# Set environment settings
env.workspace = r"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA"

# Set local variables
inRaster = r"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NPP_2012A"
outNetCDFFile = r"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NC_FILE"
rasters = arcpy.ListRasters("*.tif")

for filename in rasters:
     variable = "NPP"
     units = "gcsqm"
     XDimension = "x"
     YDimension = "y"
     bandDimension = ""

# Process: RasterToNetCDF
arcpy.RasterToNetCDF_md(filename, outNetCDFFile, variable, units,XDimension, YDimension, bandDimension)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍