topic Re: arcpy.sa.TopoToRaster Visual Fortran error in Python Questions

arcpy.sa.TopoToRaster Visual Fortran error

RodrigoRocha — Tue, 05 Apr 2016 17:09:35 GMT

Hello everybody , I'm trying to run Topo to Raster by Arcpy, in a ArcGIS Desktop 10.2 , however, every time appears the following warning and ArcGis closes. I did an experiment by turning the Topo to Raster directly from ArcToolbox and worked perfectly. I also tried to run the tool in ModelBuilder , successfully. The error message are attached.

Re: arcpy.sa.TopoToRaster Visual Fortran error

DanPatterson_Retired — Tue, 05 Apr 2016 18:59:47 GMT

seen 3 years ago... https://community.esri.com/message/468079#comment-468079

no closed resolution

Re: arcpy.sa.TopoToRaster Visual Fortran error

SteveLynch — Tue, 05 Apr 2016 19:40:30 GMT

...we'll need info 🙂

I would guess that it is related to the output raster size and/or number of inputs. Try changing the output cellsize or limiting the number of input contours. We have done much work on trapping the "out of memory" errors so this may be non-reproducible in versions post 10.2

BTW, TopoToRaster is written in Fortran and that is most likely the reason for throwing this error.

-Steve

Re: arcpy.sa.TopoToRaster Visual Fortran error

DanPatterson_Retired — Tue, 05 Apr 2016 19:44:00 GMT

I suspect not Fortran 77 ... sigh

Re: arcpy.sa.TopoToRaster Visual Fortran error

SteveLynch — Tue, 05 Apr 2016 20:00:12 GMT

Correct, however, there are still constraints and you might want to read How large a DEM can the TopoToRaster tool create?

Re: arcpy.sa.TopoToRaster Visual Fortran error

curtvprice — Wed, 06 Apr 2016 04:27:23 GMT

If the inputs are the same, I would look at the difference in environment. The error message is low-level and makes it look like files can't be written... have you set the env workspace and scratchWorkspace the same as they are set in the ArcMap session, or are they None -- or set to a non-existent workspace (by a typo or misspelling the path by neglecting to use raw strings (r"c:\workspace")

Re: arcpy.sa.TopoToRaster Visual Fortran error

RodrigoRocha — Wed, 06 Apr 2016 19:02:23 GMT

Steve, what intrigues me is the fact to error only when running in arcpy. If I use the same inputs and cellsize in TopoToRaster directly from ArcToolbox, can generate raster without any problem !!!

The solution presented by the link: https://geonet.esri.com/thread/102133#comment-468079 would not be the cause, for my Desktop is already with such a configuration ...

This problem is occurring in ArcGIS 10.2 ...

Any help is still welcome

Thanks for help and still available.

Rodrigo

Re: arcpy.sa.TopoToRaster Visual Fortran error

RodrigoRocha — Wed, 06 Apr 2016 19:03:00 GMT

Curtis, both env.worspace as env.scratchworkspace are declared in the script, however, to get 99% of the running, the pop-up on the Visual Fortran Error appears and then closes the ArcGis!

Already I tested without declaring environment, but the problem persists!

Also revised all path and spelling, but all are correct.

Dan Patterson, read about the suggested article "How large the DEM can the tool TopoToRaster create?"(https://geonet.esri.com/external-link.jspa?url=https%3A%2F%2Fblogs.esri.com%2Fesri%2Farcgis%2F2012%2F08%2F13%2Fhow-large-a-dem-can-the-topotoraster-tool-create%2F), however, the expected DEM are within the processing capacity (Win7, Intel i5-650, 4GB RAM, ArcGis10.2, etc).

Anyway, I appreciate the ideas. The discussion is always welcome!

The following script text below:

... import arcpy

... from arcpy import env

... from arcpy.sa import *

...

... # Set environment settings

... env.workspace = r"C:\Users\rodrigotr\Documents\01_PROCESSOS\2016\indice_cartas_split_buffer"

... env.scratchWorkspace = r"C:\Users\rodrigotr\Documents\01_PROCESSOS\2016\indice_cartas_split_buffer"

...

... # Set local variables

... input = r"C:\Users\rodrigotr\Documents\01_PROCESSOS\2016\indice_cartas_split_buffer\CN_cartas_buff5mUTM23sirgas2k\CN_buff5m_26124se.shp"

... outputFolder = r"C:\Users\rodrigotr\Documents\01_PROCESSOS\2016\indice_cartas_split_buffer\outMDE"

... inContours = TopoContour([[input, 'COTA']])

... inFeatures = ([inContours])

...

... # Execute TopoToRaster

... outTTR = TopoToRaster(inFeatures)

...

... # Save the output

... outTTR.save(outputFolder + "\\" + outTTR + ".tif")

That's all folks!

Re: arcpy.sa.TopoToRaster Visual Fortran error

SteveLynch — Thu, 07 Apr 2016 15:55:14 GMT

Rodrigo

What about trying it with fewer contour lines (i.e. select about 80%)?

What about in a standalone script?

How many columns and rows does the output raster have?

-Steve

Re: arcpy.sa.TopoToRaster Visual Fortran error

PeterWilson — Sun, 10 Apr 2016 07:43:55 GMT

Hi Rodrigo

I suspect its the environment settings that are not being explicitly being set when you run the TopoToRaster from Python. When you run it from within ArcMap your current environment settings would be used. I also noted that you have not provided a Boundary Feature Class or set the output Cell Size of the DEM that you are generating. I would try to first provide a boundary Feature Class to limit the DEM extent that is being generated. Make sure they are both using same projected coordinate system.

I run TopoTo Raster on a daily basis on very large data sets that few would ever require and using the following workflow have had no errors:

Dissolve my contours by the height field (Make sure to remove the selection to create multi-parts)
Repair the geometry of the dissolved contours to remove any empty or self intersecting geometries
Provide a Boundary Feature Class to limit the extent of the processing extent
Provide a Cell Size for the output DEM based on the limitations of the input contours. As a rule of thumb I don't exceed the contour height interval by half. (.i.e. Contour vertical interval is 10m so my Cell Size would be 5m)

Hope the following helps

Re: arcpy.sa.TopoToRaster Visual Fortran error

RodrigoRocha — Sun, 10 Apr 2016 20:58:35 GMT

Hi Peter,

In previous attempts, I explored the syntax arcpy.sa.TopoToRaster the most. The vertical contour interval feature class that I am using is 10m and the cellsize I use is 5m. I also use a feature class boundary. The extent already is a parameter that generates me questions. I don't know if I declare in the script as a variable using previously arcpy.env.extent.

In another experiment, I tried to create a model in ModelBuilder, using TopoToRaster and declaring the workspace. After running the tool, captured the code with their respective parameters and variables used by the Results> Copy As Python Snippet. Later, I pasted the code in ArcMap Python Window and ran. TopoToRaster with the same parameters and variables now in arcpy, after 99% closed by error.

Some friends of nearby universities, working in the geomatics area when they heard of this problem, suggested that this is a bug in arcpy.

But I still could not reach a conclusion. I will continue performing tests. I will do another test considering Peter's suggestions.

Any news, I write. But I ask everyone to continue trying to help.

Regards, Rodrigo

Re: arcpy.sa.TopoToRaster Visual Fortran error

PeterWilson — Mon, 11 Apr 2016 14:59:36 GMT

Hi Rodrigo

If the contours file is not too large zip it and attach it to the forum postl. I'll have a look at the following to see if I can resolve the problem that you are having, Regards Peter

Re: arcpy.sa.TopoToRaster Visual Fortran error

RodrigoRocha — Tue, 04 Sep 2018 15:09:00 GMT

Hello Steve, how long did I remain without accessing this forum!

I followed your suggestion and used a region with fewer outline lines (very small).

I did not quite understand your other suggestion of using a stand-alone script.

The shapefiles I used, as well as the script are attached.

I tried to run it again in the Python terminal of Arcgis (10.2) again,

but the same error described at the beginning persists.

Re: arcpy.sa.TopoToRaster Visual Fortran error

SteveLynch — Tue, 04 Sep 2018 15:16:29 GMT

Rodrigo

Please email me (slynch@esri.com) a zipped copy of the data and the parameters that you are using.

Thanks

-Steve

Re: arcpy.sa.TopoToRaster Visual Fortran error

RodrigoRocha — Tue, 04 Sep 2018 15:22:32 GMT

Hi Peter,

After 2 years without accessing this forum, due to personal problems, the files I used in a new attempt is attached, but i'm still frustrated.

Thank you very much in advance

Re: arcpy.sa.TopoToRaster Visual Fortran error

PeterWilson — Thu, 04 Oct 2018 14:07:06 GMT

Hi Rodrigo Rocha

I'll have a look at the following and get back to you shortly. Sorry I've not been on the forums either for a while. I've changed my career direction from geospatial to data science.

Regards Peter