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 ) Reference_NetCDF_RSTR.zip

I took the code provided by Dan Patterson and changed it a bit for the file naming part. However it threw me an error. Doing a manual conversion threw the error 000243: 000243: The units value is either invalid or not supported. Description The specified unit is either invalid or not supported by the UDUNITS library You can find the valid UDUNITS here: http://www.unidata.ucar.edu/software/udunits/udunits.txt When I clear the units variable (line 10) it runs without problems. Are you sure the units are OK (I'm running an old 10.3.1 version of ArcMap). # Import system modules import arcpy import os # Set local variables ws_in = r "C:\GeoNet\NPP_2010" ws_out = r "C:\GeoNet\NPP_2010\results" variable = "NPP" units = "" # gcsqm seems to be unvalid (ERROR: 000243) XDimension = "x" YDimension = "y" bandDimension = "" arcpy . env . workspace = ws_in rasters = arcpy . ListRasters ( "*" , "TIF" ) for rasname in rasters : raspath = os . path . join ( ws_in , rasname ) name , ext = os . path . splitext ( rasname ) netcdf = os . path . join ( ws_out , name ) # Process: RasterToNetCDF print "Exporting {0} to {1}" . format ( rasname , netcdf ) arcpy . RasterToNetCDF_md ( raspath , netcdf , variable , units , XDimension , YDimension , bandDimension ) Prints: Exporting r001_NPP.TIF to C:\GeoNet\NPP_2010\results\r001_NPP Exporting r002_NPP.TIF to C:\GeoNet\NPP_2010\results\r002_NPP Exporting r003_NPP.TIF to C:\GeoNet\NPP_2010\results\r003_NPP Exporting r004_NPP.TIF to C:\GeoNet\NPP_2010\results\r004_NPP Exporting r005_NPP.TIF to C:\GeoNet\NPP_2010\results\r005_NPP Content of output folder:

nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

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<SPAN class="punctuation token">:</SPAN> Failed to execute<SPAN class="punctuation token">.</SPAN> Parameters are <SPAN class="operator token">not</SPAN> valid<SPAN class="punctuation token">.</SPAN>
ERROR <SPAN class="number token">000840</SPAN><SPAN class="punctuation token">:</SPAN> The value <SPAN class="keyword token">is</SPAN> <SPAN class="operator token">not</SPAN> a Raster Layer<SPAN class="punctuation token">.</SPAN>
ERROR <SPAN class="number token">000840</SPAN><SPAN class="punctuation token">:</SPAN> The value <SPAN class="keyword token">is</SPAN> <SPAN class="operator token">not</SPAN> a Raster Catalog<SPAN class="punctuation token">.</SPAN>
Failed to execute <SPAN class="punctuation token">(</SPAN>RasterToNetCDF<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

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

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env

<SPAN class="comment token"># Set environment settings</SPAN>
env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"D:\2012A"</SPAN>

<SPAN class="comment token"># Set local variables</SPAN>
inRaster <SPAN class="operator token">=</SPAN> <SPAN class="string token">"D:\2012A"</SPAN>
outNetCDFFile <SPAN class="operator token">=</SPAN> <SPAN class="string token">"D:\2012A/nppnetcdf.nc"</SPAN>
variable <SPAN class="operator token">=</SPAN> <SPAN class="string token">"NPP"</SPAN>
units <SPAN class="operator token">=</SPAN> <SPAN class="string token">"gCSqm"</SPAN>
XDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"x"</SPAN>
YDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"y"</SPAN>
bandDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">""</SPAN>

<SPAN class="comment token"># Process: RasterToNetCDF</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>RasterToNetCDF_md<SPAN class="punctuation token">(</SPAN>inRaster<SPAN class="punctuation token">,</SPAN> outNetCDFFile<SPAN class="punctuation token">,</SPAN> variable<SPAN class="punctuation token">,</SPAN> units<SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; XDimension<SPAN class="punctuation token">,</SPAN> YDimension<SPAN class="punctuation token">,</SPAN> bandDimension<SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Reference_NetCDF_RSTR.zip

Encontre mais publicações marcadas com

Accepted answers

XanderBakker

I took the code provided by Dan Patterson and changed it a bit for the file naming part. However it threw me an error. Doing a manual conversion threw the error 000243:

000243: The units value is either invalid or not supported.
Description
The specified unit is either invalid or not supported by the UDUNITS library

You can find the valid UDUNITS here: http://www.unidata.ucar.edu/software/udunits/udunits.txt

When I clear the units variable (line 10) it runs without problems. Are you sure the units are OK (I'm running an old 10.3.1 version of ArcMap).

<SPAN class="comment token"># Import system modules</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> os

<SPAN class="comment token"># Set local variables</SPAN>
ws_in <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\GeoNet\NPP_2010"</SPAN>
ws_out <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\GeoNet\NPP_2010\results"</SPAN>

variable <SPAN class="operator token">=</SPAN> <SPAN class="string token">"NPP"</SPAN>
units <SPAN class="operator token">=</SPAN> <SPAN class="string token">""</SPAN> <SPAN class="comment token"># gcsqm seems to be unvalid (ERROR: 000243)</SPAN>
XDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"x"</SPAN>
YDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"y"</SPAN>
bandDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">""</SPAN>

arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws_in
rasters <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">for</SPAN> rasname <SPAN class="keyword token">in</SPAN> rasters<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; raspath <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_in<SPAN class="punctuation token">,</SPAN> rasname<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; name<SPAN class="punctuation token">,</SPAN> ext <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>splitext<SPAN class="punctuation token">(</SPAN>rasname<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; netcdf <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_out<SPAN class="punctuation token">,</SPAN> name<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Process: RasterToNetCDF</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Exporting {0} to {1}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>rasname<SPAN class="punctuation token">,</SPAN> netcdf<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>RasterToNetCDF_md<SPAN class="punctuation token">(</SPAN>raspath<SPAN class="punctuation token">,</SPAN> netcdf<SPAN class="punctuation token">,</SPAN> variable<SPAN class="punctuation token">,</SPAN> units<SPAN class="punctuation token">,</SPAN> XDimension<SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; YDimension<SPAN class="punctuation token">,</SPAN> bandDimension<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Prints:

Exporting r001_NPP.TIF to C:\GeoNet\NPP_2010\results\r001_NPP
Exporting r002_NPP.TIF to C:\GeoNet\NPP_2010\results\r002_NPP
Exporting r003_NPP.TIF to C:\GeoNet\NPP_2010\results\r003_NPP
Exporting r004_NPP.TIF to C:\GeoNet\NPP_2010\results\r004_NPP
Exporting r005_NPP.TIF to C:\GeoNet\NPP_2010\results\r005_NPP<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Content of output folder:

All comments

XanderBakker

I took the code provided by Dan Patterson and changed it a bit for the file naming part. However it threw me an error. Doing a manual conversion threw the error 000243:

000243: The units value is either invalid or not supported.
Description
The specified unit is either invalid or not supported by the UDUNITS library

You can find the valid UDUNITS here: http://www.unidata.ucar.edu/software/udunits/udunits.txt

When I clear the units variable (line 10) it runs without problems. Are you sure the units are OK (I'm running an old 10.3.1 version of ArcMap).

<SPAN class="comment token"># Import system modules</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> os

<SPAN class="comment token"># Set local variables</SPAN>
ws_in <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\GeoNet\NPP_2010"</SPAN>
ws_out <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\GeoNet\NPP_2010\results"</SPAN>

variable <SPAN class="operator token">=</SPAN> <SPAN class="string token">"NPP"</SPAN>
units <SPAN class="operator token">=</SPAN> <SPAN class="string token">""</SPAN> <SPAN class="comment token"># gcsqm seems to be unvalid (ERROR: 000243)</SPAN>
XDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"x"</SPAN>
YDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"y"</SPAN>
bandDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">""</SPAN>

arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> ws_in
rasters <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">for</SPAN> rasname <SPAN class="keyword token">in</SPAN> rasters<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; raspath <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_in<SPAN class="punctuation token">,</SPAN> rasname<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; name<SPAN class="punctuation token">,</SPAN> ext <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>splitext<SPAN class="punctuation token">(</SPAN>rasname<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; netcdf <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>ws_out<SPAN class="punctuation token">,</SPAN> name<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Process: RasterToNetCDF</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Exporting {0} to {1}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>rasname<SPAN class="punctuation token">,</SPAN> netcdf<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>RasterToNetCDF_md<SPAN class="punctuation token">(</SPAN>raspath<SPAN class="punctuation token">,</SPAN> netcdf<SPAN class="punctuation token">,</SPAN> variable<SPAN class="punctuation token">,</SPAN> units<SPAN class="punctuation token">,</SPAN> XDimension<SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; YDimension<SPAN class="punctuation token">,</SPAN> bandDimension<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Prints:

Exporting r001_NPP.TIF to C:\GeoNet\NPP_2010\results\r001_NPP
Exporting r002_NPP.TIF to C:\GeoNet\NPP_2010\results\r002_NPP
Exporting r003_NPP.TIF to C:\GeoNet\NPP_2010\results\r003_NPP
Exporting r004_NPP.TIF to C:\GeoNet\NPP_2010\results\r004_NPP
Exporting r005_NPP.TIF to C:\GeoNet\NPP_2010\results\r005_NPP<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Content of output folder:

DanPatterson_Retired

fixed... we need to pool our resources and get ArcThingy for our iPads

RebeccaStrauch__GISP

I'm coming late to the party , but doesn't the ListRaster command need/like 2args to use it this way? I think both are optional, but maybe try

rasters = arcpy.ListRasters("*", "TIF")

and then print or make sure it is returning something to rasters

ListRasters—Help | ArcGIS Desktop

I'm on an iPad so can't test anything, but help sample...

import arcpy
# Set the current workspace
arcpy.env.workspace = "c:/data/DEMS"
# Get and print a list of GRIDs from the workspace
rasters = arcpy.ListRasters("*", "GRID")
for raster in rasters:
print(raster)

Change "GRID" to "TIF"

DanPatterson_Retired

<SPAN class="comment token"># Import system modules</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env

<SPAN class="comment token"># Set environment settings</SPAN>
env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA"</SPAN>

<SPAN class="comment token"># Set local variables</SPAN>
inRaster <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NPP_2012A"</SPAN>
<SPAN class="comment token">#outNetCDFFile = r"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NC_FILE"</SPAN>
rasters <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">"TIF"</SPAN><SPAN class="punctuation token">)</SPAN>

c <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
<SPAN class="keyword token">for</SPAN> filename <SPAN class="keyword token">in</SPAN> rasters<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; variable <SPAN class="operator token">=</SPAN> <SPAN class="string token">"NPP"</SPAN>
&nbsp;&nbsp;&nbsp; units <SPAN class="operator token">=</SPAN> <SPAN class="string token">"gcsqm"</SPAN>
&nbsp;&nbsp;&nbsp; XDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"x"</SPAN>
&nbsp;&nbsp;&nbsp; YDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"y"</SPAN>
&nbsp;&nbsp;&nbsp; bandDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">""</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Process: RasterToNetCDF</SPAN>
&nbsp;&nbsp;&nbsp; out <SPAN class="operator token">=</SPAN> <SPAN class="string token">"{}{}.nc"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>r<SPAN class="string token">"c:\test\file_0"</SPAN><SPAN class="punctuation token">,</SPAN> c<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>RasterToNetCDF_md<SPAN class="punctuation token">(</SPAN>filename<SPAN class="punctuation token">,</SPAN> out<SPAN class="punctuation token">,</SPAN> variable<SPAN class="punctuation token">,</SPAN> units<SPAN class="punctuation token">,</SPAN>XDimension<SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; YDimension<SPAN class="punctuation token">,</SPAN> bandDimension<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; c <SPAN class="operator token">+=</SPAN> <SPAN class="number token">1</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Despite the path being named... r"c:\test\file_0", c) ... this is totally untested

ShouvikJha

Dan_Patterson‌ xander_bakker‌, can it be make a time series NetCDF (using my 12 months raster) using arcPy???

ShouvikJha

Dan_Patterson‌ xander_bakker‌ @Lan Murry , i have tried in my best to do it but i won't understanding the which way i have to solve the issue to make the time series NetCDF

ShouvikJha

I have posted my updated code reply with @Lan Murry comment. and that Stake overflow code is not posted by me , I mostly using geonet. My one of research junior also trying to solve the same issue. Thank you.

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

<SPAN class="comment token"># Import system modules</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env

<SPAN class="comment token"># Set environment settings</SPAN>
env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA"</SPAN>

<SPAN class="comment token"># Set local variables</SPAN>
inRaster <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NPP_2012A"</SPAN>
outNetCDFFile <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NC_FILE"</SPAN>
rasters <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*.tif"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">for</SPAN> filename <SPAN class="keyword token">in</SPAN> rasters<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; variable <SPAN class="operator token">=</SPAN> <SPAN class="string token">"NPP"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; units <SPAN class="operator token">=</SPAN> <SPAN class="string token">"gcsqm"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; XDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"x"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; YDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"y"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; bandDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">""</SPAN>

<SPAN class="comment token"># Process: RasterToNetCDF</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>RasterToNetCDF_md<SPAN class="punctuation token">(</SPAN>filename<SPAN class="punctuation token">,</SPAN> outNetCDFFile<SPAN class="punctuation token">,</SPAN> variable<SPAN class="punctuation token">,</SPAN> units<SPAN class="punctuation token">,</SPAN>XDimension<SPAN class="punctuation token">,</SPAN> YDimension<SPAN class="punctuation token">,</SPAN> bandDimension<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

DanPatterson_Retired

Hmmm it seems you haven't responded to iamurray‌ comments or implemented or commented on previous suggestions. Perhaps you could show what you have tried in more detail then since the repetition of the error message really isn't providing any new information to help solve your problem. Also, maybe you will get some suggestions over on Stack Overflow, that may help people here or even prevent them from reinventing a suggestion.

ShouvikJha

Dan_Patterson‌, I am trying to find the solution, following by above comments but its not solved yet.

DanPatterson_Retired

moved on to try here python - How to convert multiple Raster to NetCDF using ArcPy - Stack Overflow

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.

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

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.

<SPAN class="comment token"># Import system modules</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env

<SPAN class="comment token"># Set environment settings</SPAN>
env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA"</SPAN>

<SPAN class="comment token"># Set local variables</SPAN>
inRaster <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NPP_2012A"</SPAN>
outNetCDFFile <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"D:\NPP_GUJARAT_NEW\MODIS_BKU_DATA\NC_FILE\nppnetcdf.nc"</SPAN>
variable <SPAN class="operator token">=</SPAN> <SPAN class="string token">"NPP"</SPAN>
units <SPAN class="operator token">=</SPAN> <SPAN class="string token">"gcsqm"</SPAN>
XDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"x"</SPAN>
YDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">"y"</SPAN>
bandDimension <SPAN class="operator token">=</SPAN> <SPAN class="string token">""</SPAN>

<SPAN class="comment token"># Process: RasterToNetCDF</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>RasterToNetCDF_md<SPAN class="punctuation token">(</SPAN>inRaster<SPAN class="punctuation token">,</SPAN> outNetCDFFile<SPAN class="punctuation token">,</SPAN> variable<SPAN class="punctuation token">,</SPAN> units<SPAN class="punctuation token">,</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; XDimension<SPAN class="punctuation token">,</SPAN> YDimension<SPAN class="punctuation token">,</SPAN> bandDimension<SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

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.

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

ShouvikJha

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

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&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;File Name&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Line&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Position&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Traceback&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&nbsp;&nbsp;&nbsp; <SPAN class="operator token"><</SPAN>module<SPAN class="operator token">></SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="operator token"><</SPAN>module1<SPAN class="operator token">></SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="number token">19</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&nbsp;&nbsp;&nbsp; RasterToNetCDF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;C<SPAN class="punctuation token">:</SPAN>\Program Files\ArcGIS\Desktop10<SPAN class="number token">.3</SPAN>\ArcPy\arcpy\md<SPAN class="punctuation token">.</SPAN>py&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="number token">253</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
ExecuteError<SPAN class="punctuation token">:</SPAN> Failed to execute<SPAN class="punctuation token">.</SPAN> Parameters are <SPAN class="operator token">not</SPAN> valid<SPAN class="punctuation token">.</SPAN>
ERROR <SPAN class="number token">000840</SPAN><SPAN class="punctuation token">:</SPAN> The value <SPAN class="keyword token">is</SPAN> <SPAN class="operator token">not</SPAN> a Raster Layer<SPAN class="punctuation token">.</SPAN>
ERROR <SPAN class="number token">000840</SPAN><SPAN class="punctuation token">:</SPAN> The value <SPAN class="keyword token">is</SPAN> <SPAN class="operator token">not</SPAN> a Raster Catalog<SPAN class="punctuation token">.</SPAN>
Failed to execute <SPAN class="punctuation token">(</SPAN>RasterToNetCDF<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

DanPatterson_Retired

you have set your environment workspace (ie a folder or a geodatabase) as \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

Participe da Conversa

Membros conectados podem postar, seguir atualizações e mais. Novo aqui? Registre uma conta gratuita.

Links Rápidos

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.