def testSaPaRun(self): with TempDir() as d: ratio_raster = os.path.join(d, 'test_sapa_ratio.tif') surf_raster = os.path.join(d, 'test_surface_area.tif') arcpy.env.scratchWorkspace = d surface_area_to_planar_area.main(config.bathy_raster, ratio_raster, surf_raster) self.assertTrue(os.path.exists(ratio_raster)) self.assertTrue(os.path.exists(surf_raster)) self.assertAlmostEqual(su.raster_properties(ratio_raster, "MEAN"), \ 1.0042422342676) self.assertAlmostEqual(su.raster_properties(ratio_raster, "STD"), \ 0.0058175502835692) self.assertAlmostEqual(su.raster_properties(surf_raster, "MEAN"), \ 25.119343739217) self.assertAlmostEqual(su.raster_properties(surf_raster, "STD"), \ 0.14551573347447)
def testSaPaRunCorrection(self): with TempDir() as d: ratio_raster = os.path.join(d, 'test_sapa_ratio.tif') surf_raster = os.path.join(d, 'test_surface_area.tif') arcpy.env.scratchWorkspace = d surface_area_to_planar_area.main( config.bathy_raster, ratio_raster, True, surf_raster) self.assertTrue(os.path.exists(ratio_raster)) self.assertTrue(os.path.exists(surf_raster)) self.assertAlmostEqual( su.raster_properties(ratio_raster, "MEAN"), 1.000616701751028) self.assertAlmostEqual( su.raster_properties(ratio_raster, "STD"), 0.0008445980914179108) self.assertAlmostEqual( su.raster_properties(surf_raster, "MEAN"), 25.119343091857) self.assertAlmostEqual( su.raster_properties(surf_raster, "STD"), 0.14551570953239)
def testSaPaRunCorrection(self): with TempDir() as d: ratio_raster = os.path.join(d, 'test_sapa_ratio.tif') surf_raster = os.path.join(d, 'test_surface_area.tif') arcpy.env.scratchWorkspace = d surface_area_to_planar_area.main(config.bathy_raster, ratio_raster, True, surf_raster) self.assertTrue(os.path.exists(ratio_raster)) self.assertTrue(os.path.exists(surf_raster)) self.assertAlmostEqual(su.raster_properties(ratio_raster, "MEAN"), 1.000616701751028) self.assertAlmostEqual(su.raster_properties(ratio_raster, "STD"), 0.0008445980914179108) self.assertAlmostEqual(su.raster_properties(surf_raster, "MEAN"), 25.119343091857) self.assertAlmostEqual(su.raster_properties(surf_raster, "STD"), 0.14551570953239)
def testSaPaRunNoCorrection(self): with TempDir() as d: ratio_raster = os.path.join(d, 'test_sapa_ratio.tif') surf_raster = os.path.join(d, 'test_surface_area.tif') arcpy.env.scratchWorkspace = d surface_area_to_planar_area.main(config.bathy_raster, ratio_raster, False, surf_raster) self.assertTrue(os.path.exists(ratio_raster)) self.assertTrue(os.path.exists(surf_raster)) self.assertAlmostEqual(su.raster_properties(ratio_raster, "MEAN"), 1.0042422342676) self.assertAlmostEqual(su.raster_properties(ratio_raster, "STD"), 0.0058175502835692) self.assertAlmostEqual(su.raster_properties(surf_raster, "MEAN"), 25.119343091857) self.assertAlmostEqual(su.raster_properties(surf_raster, "STD"), 0.14551573347447)
def testSaPaRunWithFgdbLocation(self): with TempDir() as d: fgdb_name = "sapa.gdb" fgdb_workspace = os.path.join(d, fgdb_name) # create a temporary File Geodatabase location arcpy.CreateFileGDB_management(d, fgdb_name) self.assertTrue(os.path.exists(fgdb_workspace)) ratio_raster = os.path.join(fgdb_workspace, 'test_sapa_ratio') surf_raster = os.path.join(fgdb_workspace, 'test_sapa_area') arcpy.env.scratchWorkspace = d surface_area_to_planar_area.main( config.bathy_raster, ratio_raster, False, surf_raster) self.assertAlmostEqual( su.raster_properties(ratio_raster, "MEAN"), 1.0042422342676) self.assertAlmostEqual( su.raster_properties(ratio_raster, "STD"), 0.0058175502835692) self.assertAlmostEqual( su.raster_properties(surf_raster, "MEAN"), 25.119343091857) self.assertAlmostEqual( su.raster_properties(surf_raster, "STD"), 0.14551573347447)
def testSaPaRunWithFgdbLocation(self): with TempDir() as d: fgdb_name = "sapa.gdb" fgdb_workspace = os.path.join(d, fgdb_name) # create a temporary File Geodatabase location arcpy.CreateFileGDB_management(d, fgdb_name) self.assertTrue(os.path.exists(fgdb_workspace)) ratio_raster = os.path.join(fgdb_workspace, 'test_sapa_ratio') surf_raster = os.path.join(fgdb_workspace, 'test_sapa_area') arcpy.env.scratchWorkspace = d surface_area_to_planar_area.main(config.bathy_raster, ratio_raster, False, surf_raster) self.assertAlmostEqual(su.raster_properties(ratio_raster, "MEAN"), 1.0042422342676) self.assertAlmostEqual(su.raster_properties(ratio_raster, "STD"), 0.0058175502835692) self.assertAlmostEqual(su.raster_properties(surf_raster, "MEAN"), 25.119343091857) self.assertAlmostEqual(su.raster_properties(surf_raster, "STD"), 0.14551573347447)