Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange
3

You can do this just by using the Mosaic to New Raster geoprocessing tool. Just choose maximum for your Mosaic Operator and the new raster will use the high value anywhere that they overlap.


3

Union (geoprocessing) feature make a "union" of 2 dataset (C = A united to B). Probably you are tring to JOIN two dataset. Try this: Go to JOIN settings: (right click -> properties -> Join) Add Vector Join of the second layer (click on "+" green-plus icon) and set these parameters: Join Layer: "the second layer to join" Join field: "ID" Join Target: "ID" ...


2

It sounds like you need jiffle which is contained in the gt-process-raster module. Once you have constructed a Process you can pass in a list of source rasters and bands for it to process in the execute method. The actual operation is defined in the jiffle scripting language which is fairly powerful including loops, conditionals and good range of functions. ...


2

If the two source layers have the same fields, you could combine them together using the merge tool. Using the data you currently have, you would have to rely on the coalesce function to handle the nulls. If the field has a value, it is used. If the field is NULL, then the replacement value is used. For instance, you could replace null by 0: coalesce(...


2

calcRaster is a temporary raster, not a text string. You will need to find another way of naming your result, like outRasterName = '{}720'.format(inRaster). Python has tools for creating complete path names from path fragments (folder & file names) use os.path.join like this: CalcRasterName = os.path.join(outFolder, outRasterName).


1

There is also the NNJoin plugin, which does a nearest-neighbor join. Each point in the first layer will have the attributes added to it from the nearest point in the second layer. Similar to gianfranco di pietro's answer, but without needing a shared, unique identifier.


1

I was able to get the code working using the following approach. import arcpy,os,sys,string from arcpy import env from arcpy.sa import * arcpy.CheckOutExtension("Spatial") input_folder = r'SomeFolder' arcpy.env.workspace = input_folder output_folder = r'SomeFolder' #Set reflectance constant const = 672 # Get raster list, loop through and multiply be ...


1

If you set the extent to raster1 image, then you can use Con with IsNull functions to convert the Null data of your output raster to 0 as follows: Con(IsNull("OutputRaster.img"),0,"OutputRaster.img") The output will convert Null data to 0 of your output raster.


1

If you set your environment extent correctly to raster 1, then the extent of the output raster will be that of raster 1, but you will only get data where raster 2 has data (because Value + NoData == NoData). You can use the CellStatistics tool, with the Sum option and Ignore NoData set to "DATA" with the environment extent set to raster 1


1

I have always had a hard time trying to get rasters of differing spatial extents to output to the "Union of Inputs". You could try setting this up first before trying the steps outlined below. First try to set the processing extent to "Union of Inputs". If that does not work, try this. You can use the "Mosaic to New Raster" process. Make sure Raster 2 is ...


1

Thank you @FelixIP for the help. I managed to correct my mistake. By using the arcpy.Raster() function I, then, had to check out a Spatial Analyst license in the script. Below is my final draft of the code. It is converting Meters to Feet print "Starting Process 2: Converting Units to Feet..." #Checkout License class LicenseError(Exception): pass try: ...


1

The RemapValue is an arcpy.sa (Spatial Analyst) object. By not referencing the spatial analyst module, the script doesn't know what to do with your RemapValue statement. The easiest fix would probably be to include the following at the beginning of your script... import arcpy from arcpy.sa import * arcpy.CheckOutExtension("Spatial") Alternatively, I ...


Only top voted, non community-wiki answers of a minimum length are eligible