Hi, I have multiple netCDF files I need to make a "netCDF raster" layer and then "project raster" and save as a ".tif" file. I have the following script but this gives me an error message. Please see the code and error message below: Thanks in advance. # Import system modules import arcpy , sys from arcpy import env from arcpy . sa import * # Input data source arcpy . env . workspace = r "C:/Users/jayaskeradl/Desktop/weekly_avg_project" arcpy . env . overwriteOutput = True # Set output folder OutputFolder = r "C:/Users/jayaskeradl/Desktop/weekly_avg_project/output" # Loop through a list of files in the workspace NCfiles = arcpy . ListFiles ( "*.nc" ) for filename in NCfiles : print ( "Processing: " + filename ) inNCfiles = arcpy . env . workspace + "//" + filename fileroot = filename [ 0 : ( len ( filename ) - 3 ) ] TempLayerFile = "SM_amount" outRaster = OutputFolder + "//" + fileroot # Process: Make NetCDF Raster Layer arcpy . MakeNetCDFRasterLayer_md ( filename , "var250" , "lon" , "lat" , TempLayerFile , "" , "" , "BY_VALUE" ) arcpy . ProjectRaster_management ( filename , TempLayerFile , "GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]" , "NEAREST" , "0.125 0.125" , "NAD_1983_To_WGS_1984_1" , "" , "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]" ) # Process: Copy Raster arcpy . CopyRaster_management ( TempLayerFile , outRaster + ".tif" , "" , "" , "" , "NONE" , "NONE" , "" ) print "***DONE!!!" print arcpy . GetMessages ( ) ExecuteError : Failed to execute . Parameters are not valid . Undefined coordinate system for input dataset . Failed to execute ( ProjectRaster ) .

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; } } }

Error in multiple netCDF file processing followed by raster projection

DuminduJayasekera

Hi,

I have multiple netCDF files I need to make a "netCDF raster" layer and then "project raster" and save as a ".tif" file. I have the following script but this gives me an error message.

Please see the code and error message below:

Thanks in advance.

<SPAN class="comment token"># Import system modules&nbsp; </SPAN>
<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</SPAN> sys 
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env&nbsp; 
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>&nbsp; 
&nbsp; 
<SPAN class="comment token"># Input data source&nbsp; </SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:/Users/jayaskeradl/Desktop/weekly_avg_project"</SPAN>&nbsp; 
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>&nbsp; 
&nbsp; 
<SPAN class="comment token"># Set output folder&nbsp; </SPAN>
OutputFolder <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:/Users/jayaskeradl/Desktop/weekly_avg_project/output"</SPAN>&nbsp; 
&nbsp; 
<SPAN class="comment token"># Loop through a list of files in the workspace&nbsp; </SPAN>
NCfiles <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListFiles<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*.nc"</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; 
&nbsp; 
<SPAN class="keyword token">for</SPAN> filename <SPAN class="keyword token">in</SPAN> NCfiles<SPAN class="punctuation token">:</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Processing: "</SPAN> <SPAN class="operator token">+</SPAN> filename<SPAN class="punctuation token">)</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; inNCfiles <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">+</SPAN> <SPAN class="string token">"//"</SPAN> <SPAN class="operator token">+</SPAN> filename&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; fileroot <SPAN class="operator token">=</SPAN> filename<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="punctuation token">(</SPAN>len<SPAN class="punctuation token">(</SPAN>filename<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; TempLayerFile <SPAN class="operator token">=</SPAN> <SPAN class="string token">"SM_amount"</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; outRaster <SPAN class="operator token">=</SPAN> OutputFolder <SPAN class="operator token">+</SPAN> <SPAN class="string token">"//"</SPAN> <SPAN class="operator token">+</SPAN> fileroot&nbsp; 
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Process: Make NetCDF Raster Layer&nbsp; </SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>MakeNetCDFRasterLayer_md<SPAN class="punctuation token">(</SPAN>filename<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"var250"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"lon"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"lat"</SPAN><SPAN class="punctuation token">,</SPAN> TempLayerFile<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"BY_VALUE"</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; 
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>ProjectRaster_management<SPAN class="punctuation token">(</SPAN>filename<SPAN class="punctuation token">,</SPAN> TempLayerFile<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"NEAREST"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"0.125 0.125"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"NAD_1983_To_WGS_1984_1"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Process: Copy Raster&nbsp; </SPAN>
&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>TempLayerFile<SPAN class="punctuation token">,</SPAN> outRaster <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"NONE"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"NONE"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; 
&nbsp; 
&nbsp; 
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"***DONE!!!"</SPAN>&nbsp; 
<SPAN class="keyword token">print</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetMessages<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN> 


‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍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>
Undefined coordinate system <SPAN class="keyword token">for</SPAN> input dataset<SPAN class="punctuation token">.</SPAN>
Failed to execute <SPAN class="punctuation token">(</SPAN>ProjectRaster<SPAN class="punctuation token">)</SPAN><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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Trouver plus de publications taguées avec

Accepted answers

IanMurray

You have the same input for Make NetCDF Layer and for Project Raster(filename or the .nc file), while you should take the output from Make NetCDF Layer (an actual raster layer, which should have a coordinate system) and use that for the input of Project Raster. The output of Project Raster would need a different name than you are currently using for the input.

All comments

DuminduJayasekera

Thanks again Ian. I have tried what you mentioned in the above comment. (Please see the code and errors above). I would like to know how I can copy the string representation of a coordinate system using the model builder as you noted in 3rd bullet.

I have opened the properties of the raster of same coordinate system but it does not provide string representation of a coordinate system. See the image below.

Again your help is highly appreciated to make the script working and moving forward.

IanMurray

Code looks much better now thanks.

Issue now is how your output coordinate system parameter is formatted, as it will not accept the string name for the coordinate system nor the WKID for the coordinate system

"The coordinate system to which the input raster will be projected. The default value is set based on the Output Coordinate System environment setting.

Valid values for this parameter are

A file with the ".prj" extension (the prj files which ship with ArcGIS can be found in "C:\Program Files\ArcGIS\Coordinate Systems").
An existing feature class, feature dataset, raster catalog (basically anything with a coordinate system).
The string representation of a coordinate system. These lengthy strings can be generated by adding a coordinate system variable to ModelBuilder, setting the variable's value as desired, then exporting the model to a Python script."

You have a few options, one is to have a .prj file somewhere with the output spatial reference you want and use that as the parameter(http://desktop.arcgis.com/en/arcmap/10.3/tools/data-management-toolbox/project-raster.htm). Similarly if you had a raster you manually projected into that coordinate system somewhere in your workspace, you could reference that as an input.

Another option I believe is to set the Output Coordinate System environment setting in your script to as it should then become the default value for the output coordinate system variable(you can do this by placing arcpy.env.outputCoordinateSystem = arcpy.SpatialReference(4269) with the rest of your environment variables(Actually I'm curious if this works, someone else had a similar error here and I wonder if you left the 3rd parameter blank the environment setting would entirely take care of things.)

DuminduJayasekera

Thanks Ian.

I appreciate your reply. I have considered your comments and modified in the script. (Please see the script and error above). My code still does not work. I get a different error now. I really appreciate if you can help me the script working.

Thanks a lot again.

IanMurray

First off, when posting code to GeoNet please use the code formatting to make your script more legible(See https://community.esri.com/people/curtvprice/blog/2014/09/25/posting-code-blocks-in-the-new-geonet?sr=search&searchId=33390d5c-862e-4489-8e39-2afae930a01c&searchIndex=2‌ for more on that).

Currently you are getting an indent error, so either something is indented where it shouldn't be, you are mixing spaces and tabs for your indenting.

I do not believe that will be the only issue since you code has some fairly serious problems. since the foreslash character is an escape character, the file paths you used in your workspace and output folder will not be read properly. This is easily fixed by placing an r character prior to the string to make it a raw string and treat the foreslash as a raw string not an escape character.

arcpy.env.workspace = r"C:/Users/Desktop/weekly_avg_project"

OutputFolder = r"C:/Users/Desktop/weekly_avg_project/output"

You have a similar problem with escape characters in the following lines:

inNCfiles = arcpy.env.workspace + "/" + filename

outRaster = OutputFolder + "/" + fileroot

they would need to be double foreslashes to end up with a single foreslash character in python(or even better, use os.path.join to concatenate filepaths instead of just string concatentation. See: https://docs.python.org/2/library/os.path.html)

inNCfiles = arcpy.env.workspace + "//" + filename

outRaster = OutputFolder + "//" + fileroot

Hopefully this gets you closer to the script working, if not repost your updated code with the code formatting and new error message and we can go from there.

Rejoindre la conversation

Les membres connectés peuvent publier, suivre les mises à jour, et plus encore. Nouveau ici ? Inscrivez-vous gratuitement.

Liens rapides

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