topic Re: Invalid floating point operation. arcpy in Python Questions

Invalid floating point operation. arcpy (ArcSDE only)

MarcGoetzke — Thu, 26 Feb 2015 13:50:36 GMT

I always get this Error when i do simple raster calculations with arcpy. i use PyScripter 2.7.5 and ArcGIS 10.2.1 :

Invalid floating point operation

Exception class: EInvalidOp

Exception message: Invalid floating point operation.

Exception address: 00007FF93BB646BE

------------------------------------------------------------------------------

PyScripter version : 2.5.3.0 x64

Python DLL : python27.dll

Python Engine : peInternal

I think its an ArcSDE only memory/cache problem. I tried things like arcpy.ClearWorkspaceCache_management() but it didnt help.

Any suggestions?

Thx

My Program:

#simple raster calculation of a list of rasters

outRas = Int(Raster('test1') + Raster('test2'))

outRas.save(destination + '/' + 'test3')

Re: Invalid floating point operation. arcpy

DanPatterson_Retired — Thu, 26 Feb 2015 13:58:48 GMT

Where is destination defined? something is incomplete, does your first posted line work? have you tried

outRas.save() without specifying a location?

Re: Invalid floating point operation. arcpy

VandanaRaghunathan — Fri, 27 Feb 2015 00:30:47 GMT

Hello Marc,

Can you try adding a .tif extension to your output raster?

outRas = Int(Raster('test1') + Raster('test2'))

outRas.save(destination + '\\' + 'test3.tif')

-or-

import os

outRas.save(os.path.join(destination, 'test3.tif')

Thanks,

Vandana

Re: Invalid floating point operation. arcpy

DanPatterson_Retired — Fri, 27 Feb 2015 01:34:45 GMT

without the extension, it will produce an esri grid, which is preferable to tifs

Re: Invalid floating point operation. arcpy

VandanaRaghunathan — Fri, 27 Feb 2015 01:46:36 GMT

Hello Dan,

That is correct, without extension it will be an esri grid. However, esri grid has the following limitations and if they are not met, it will error out:

Storage limitation

The name of a grid is limited as follows:

It cannot be stored using spaces.
It cannot start with a number.
It cannot be longer than 13 characters (a multiband grid is allowed up to 9 characters).

There is a limit to the number of files that can be stored in an INFO directory for both coverages and grids. This total is approximately 10,000. Therefore, this limits the number of grids you can store in a workspace. For example, the following lists the theoretical maximum number of grid datasets that can be stored in a single workspace directory:

Fewer than 5,000 floating point grids, or
Fewer than 3,333 integer grids, with VATs (fewer than 5,000 if no VATs), or
Fewer than 10,000 grid stacks

In addition, there is a limit on the number of characters allowed for a GRID path in windows. The path cannot be longer than 120 characters (including the GRID name)

Esri Grid format—Help | ArcGIS for Desktop

35317 - Do GRID paths have a limit on the number of characters allowed on Windows platforms?

Thanks,

Vandana

Re: Invalid floating point operation. arcpy

XanderBakker — Fri, 27 Feb 2015 01:46:42 GMT

Apart from the output destination being unclear, the inputs are also "incomplete". Running your code from PyScripter, will require the import of arcpy and defining where the inputs can be found. If not specifying the entire path to the input data you should defined the arcpy.env.workspace where the data resides. May need to checkout the spatial analyst extension too.

So please post the entire code you are running. If possible attach (part of) the inputs or provide more information on the type of raster, size, projection, etc. These could all influence.

Did you try to run the code from the python window inside ArcMap? If so, is the same error occurring?

Re: Invalid floating point operation. arcpy

DanPatterson_Retired — Fri, 27 Feb 2015 02:19:29 GMT

and none of these are an issue if you know what you are doing... sorry ... you are barking at the wrong old dog

Re: Invalid floating point operation. arcpy

MarcGoetzke — Fri, 27 Feb 2015 11:59:25 GMT

Sorry this is the whole program:

import arcpy, os

from arcpy import env

from arcpy.sa import *

arcpy.CheckOutExtension("spatial")

arcpy.env.overwriteOutput = True

def main():

OD_path = 'S:/Project/Raster1'

DD_path = 'S:/Project/Raster2'

destination = 'S:/Project/Output'

liste = ('n','nh4','no3')

for chem in liste:

liste2 = ('te','tu','ta')

for item in liste2:

OD = OD_path + '/' + '%s_s%s09'%(chem, item)

DD = DD_path + '/' + '%s_t%s09'%(chem, item)

Out = '%s_f%s09'%(chem, item)

outRas = Raster(OD) + Raster(DD)

outRas.save(destination + '/' + Out)

arcpy.ClearWorkspaceCache_management()

if __name__ == '__main__':

main()

Re: Invalid floating point operation. arcpy

MarcGoetzke — Fri, 27 Feb 2015 11:59:51 GMT

thx for the replies

Re: Invalid floating point operation. arcpy

MarcGoetzke — Fri, 27 Feb 2015 12:11:19 GMT

My program works for "n" but not for "nh4" and "no3", i always get the Invalid floating point operation Error. When i do the calculation in ArcGIS manually it works fine. I dont have many files but i get this error from time to time in different programs and sometimes this error appears and sometimes not in the same program

Re: Invalid floating point operation. arcpy

XanderBakker — Sat, 11 Dec 2021 17:56:15 GMT

The problem might be data related, but could you try and replace the + with using arcpy.sa.Plus?

import arcpy, os
from arcpy.sa import *

def main():
    OD_path = r'S:\Project\Raster1'
    DD_path = r'S:\Project\Raster2'
    destination = r'S:\Project\Output'

    liste = ('n','nh4','no3')
    liste2 = ('te','tu','ta')

    arcpy.env.overwriteOutput = True
    arcpy.CheckOutExtension("SPATIAL")

    for chem in liste:
        for item in liste2:
            OD = os.path.join(OD_path, "{0}_s{1}09".format(chem, item))
            DD = os.path.join(DD_path, "{0}_t{1}09".format(chem, item))

            Out = "{0}_f{1}09".format(chem, item)

            # outRas = Raster(OD) + Raster(DD)
            outRas = Plus(OD, DD)
            outRas.save(os.path.join(destination, Out))

    # arcpy.ClearWorkspaceCache_management()
    arcpy.CheckInExtension("SPATIAL")

if __name__ == '__main__':
    main()

BTW: you can increase readability of the code by using syntax highlighting for your code (see: Posting Code blocks in the new GeoNet )

Re: Invalid floating point operation. arcpy

MarcGoetzke — Fri, 27 Feb 2015 14:49:16 GMT

@Xander: Thx for the help. i tried your optimized program but i get the same error. i forgot to mention that two errors alternate. The floating point error and this one:

Access violation at address 00007FF959AF8AB3 in module 'ntdll.dll'. Read of address 0000000000000029.

Someone gave me the advice that it could really be a memory/ram problem. Hmmm

Re: Invalid floating point operation. arcpy

MarcGoetzke — Fri, 13 Mar 2015 12:40:02 GMT

Okay i think i solved the problem. I simply put "arcpy.ClearWorkspaceCache_management()" at the beginning of the program so the arcpy work cache cant overflow. If you have many calculation steps put it in between. The access violation problem occures because arcpy has a limited work cache space assigned.

Thx for the help everyone

Re: Invalid floating point operation. arcpy

DanPatterson_Retired — Sun, 15 Mar 2015 10:38:00 GMT

Please indicate in your original post that this error applies ArcSDE only

Clear workspace cache