Latest Contributions by JeremyProville

‎08-06-2013

Thanks Luke - that and a combination of the other workarounds seems to have done the trick. I ended up making a blank 'extent' raster with -9999 values in every cell, and using that to set the extent. I also had to toss in a constant (-9998) in the final raster calculator formula, for HighestPosition. I put this in last place so that I could filter out these NoData (now -9998) areas. Appreciate everyone's help!

‎08-05-2013

Are you saying the output extent is wrong (this is what the environment controls) or just the particular cells are only getting a value if all the inputs from the stack are data cells? If the latter, Luke's suggestion, Con(IsNull(raster), -999, raster) should work. It is the latter - I've tried the code you mentioned and variations of it, but even by forcing NoData values into all input datasets, it still only gives me the intersection of all inputs, and does not consider all areas that were previously NoData. It is very strange. The closest I've gotten to understanding this is by going in ArcCatalog and forcing in NoData values to -9999. But it still is erratic, I'll illustrate with an example: Rasters A and B contain data. I force NoData values in these rasters to be -9999, either in Catalog or by using raster calculator with the code above. If I now use HighestPosition on these, it still only returns intersecting data in areas where values were anything but -9999... i.e. it seems to know where those NoData areas were, and omits them. However, if I put in Raster A, and a constant dummy value of e.g. -1 as the inputs to HighestPosition, it works and includes all '-9999' areas in the output. Putting in Rasters A, B and the dummy constant makes it fail again, and the output ignores -9999 areas in it's calculation. Any insight much appreciated.

‎08-05-2013

Thanks to both for your responses. TKONG, indeed lpinner's first answer was right, I was looking for the HighestPosition tool and not the max tool. lpinner, appreciate your code suggestion - I tried a great deal of variations on it, but I continue to be baffled by the fact that I only keep getting the intersection of the inputs rather than the union. No matter what I seem to try (e.g. tinkering with the GP extent, con() isnull() and setnull() in raster calculator), I get the same result. I've also tried forcing -999 NoData values through Catalog. In any case, thanks for all your help and I'll post back if I figure out the problem. Let me know if you have any further ideas.

‎07-31-2013

Use the HighestPosition function/tool. Used to be called UPos in earlier versions. Thanks lpinner, this was spot-on! In theory it is doing what I need, however in practice it's not working out exactly as planned. I have a follow up question I'm hoping I can get some help with: Running the tool with three or more input datasets seems to include only results that are an intersection of all the inputs. All other non-intersecting areas are NoData. What I would like is to have a dataset represented in the results, even if it is the only one with a value in a specific cell. I tried get around this by adding in a constant value as an input, or a blank raster with constant values everywhere, such that there would be a 'dummy' for comparison for every cell. This didn't seem to work and I'd love to hear about any ideas or further insight into how this tool works. Much appreciated!

‎07-30-2013

Hi, I have a pile of 23 different rasters, on which I'm performing various calculations. One of them is a simple Maximization, using Cell Statistics, and one thing I'd like to see is a raster showing me which of the 23 layers contained the maximum value that was chosen. I'm finding this is pretty hard to do. In Arc 9.3, I solved this problem using a long string of Con() and Max() statements in the raster calculator, to assign an integer index value to each raster that was selected in the maximization. This is no longer possible in Arc 10 since Max() has been removed. I tried an approach in Python using the code below, but it does not really work. What I'm looking for is a way to return a raster with either index values from 1-23, or ideally the names of each input raster. Each name would be assigned to the cell if it was chosen as a maximum value. Hopefully this makes sense, I can elaborate if not. Any ideas here would be hugely helpful. Here's some sample code I've been playing with, but there has to be a simpler, more elegant solution - perhaps an obvious tool I'm missing? Thanks! opt2 = Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "blank",-9999, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "CottonLB" ,1, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "FodderLB" ,2, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "FruitsLB" ,3, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "GroundnutLB" ,4, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "JatrophaLB" ,5, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "MaizeLB" ,6, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "MilletLB" ,7, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "MiscanthusLB" ,8, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "NarcostimulantsLB" ,9, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "OilpalmLB" ,10, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "PulsesLB" ,11, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "RapeseedLB" ,12, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "RiceLB" ,13, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "RootsLB" ,14, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "SorghumLB" ,15, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "SoybeanLB" ,16, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "SugarbeetLB" ,17, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "SugarcaneLB" ,18, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "SunflowerLB" ,19, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "SwitchgrassLB" ,20, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "TubersLB" ,21, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "VegetablesLB" ,22, Con(CellStatistics(["CottonLB","FodderLB","FruitsLB","GroundnutLB","JatrophaLB","MaizeLB","MilletLB","MiscanthusLB","NarcostimulantsLB","OilpalmLB","PulsesLB","RapeseedLB","RiceLB","RootsLB","SorghumLB","SoybeanLB","SugarbeetLB","SugarcaneLB","SunflowerLB","SwitchgrassLB","TubersLB","VegetablesLB","WheatLB","blank"],"MAXIMUM","DATA") == "WheatLB" ,23,-9999))))))))))))))))))))))))

‎07-30-2013

Hi, I've seen some of the other threads about this, and I'm aware that the Max() statement has been stripped from the Raster Calculator as of ArcGIS 10. This was extremely handy for running quick Maximization routines and the reason why functionality like this would be stripped rather than left in eludes me.... Leaving that aside, I have some RC code I was using in Arc 9.3 with Max() statements embedded within Con()s. I'd love to hear any creative ideas of how I might be able to get this to work in Arc 10, because I can't see how using Cell Statistics or Raster Properties GUI tools can be a solution here. The code basically looks through a pile of ~23 rasters, and assigns an index value ranging from 1 to 23 depending on which of these rasters contained the maximum value. My plan is to try and replace specific terms to make the syntax work with Python but any other ideas would be hugely helpful. Thanks! Sample code below: Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) < -9997,-9999, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [CottonLB] ,1, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [FodderLB] ,2, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [FruitsLB] ,3, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [GroundnutLB] ,4, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [JatrophaLB] ,5, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [MaizeLB] ,6, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [MilletLB] ,7, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [MiscanthusLB] ,8, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [NarcostimulantsLB] ,9, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [OilpalmLB] ,1LB, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [PulsesLB] ,11, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [RapeseedLB] ,12, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [RiceLB] ,13, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [RootsLB] ,14, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [SorghumLB] ,15, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [SoybeanLB] ,16, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [SugarbeetLB] ,17, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [SugarcaneLB] ,18, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [SunflowerLB] ,19, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [SwitchgrassLB] ,2LB, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [TubersLB] ,21, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [VegetablesLB] ,22, Con(max([CottonLB],[FodderLB],[FruitsLB],[GroundnutLB],[JatrophaLB],[MaizeLB],[MilletLB],[MiscanthusLB],[NarcostimulantsLB],[OilpalmLB],[PulsesLB],[RapeseedLB],[RiceLB],[RootsLB],[SorghumLB],[SoybeanLB],[SugarbeetLB],[SugarcaneLB],[SunflowerLB],[SwitchgrassLB],[TubersLB],[VegetablesLB],[WheatLB]) == [WheatLB] ,23))))))))))))))))))))))))

‎11-01-2011

Can you just run your regression on a subset of your data? Take a systematic sample of your rasters and run your anayses on those. Hi passiveoptical, Thanks, your suggestion is a good idea and I'm sure would have been a good solution in other circumstances. I've been running into problems when trying that approach: since the dataset is global in scale and only one of many regressions I'll have to do, the sampling would be very time intensive. I suppose I could script it but there's another issue in that after getting the expected values, I have to reconvert to raster. The regression outputs need to be combined with other datasets subsequently (using the raster calculator). I might just have to try the batch scripting method if there's no other ideas...I do find it odd that the regression tools don't accept rasters of identical dimensions though, it seems like a pretty straightforward and necessary operation. Thanks for the advice!

‎10-28-2011

Hello, I've got two raster layers. What I want to do is use one layer, that covers the whole globe, to compute the expected values of another raster in places where there is no data. What are the tools or processes that I should look into in order to do this? Thank you!

‎10-28-2011

Hi, I've got two raster layers to do with a specific crop: Yields Suitability Index for growing What I'd like to do is use the suitability layer, which covers the entire globe, to impute the potential (expected) yield in places where the crop is not grown. I would do this using the actual yields layer. I can't use rasters as inputs for the OLS or GWR tool. I tried converting them to vector format, doing a spatial join, and then running a regression but it won't work as the datasets are simply too big when I convert them to polygons. How would I run a regression on rasters, to impute the expected values of one based on another? Thanks!

‎10-27-2011

Hi, I've got two raster layers to do with a specific crop: Yields Suitability Index for growing What I'd like to do is use the suitability layer, which covers the entire globe, to impute the potential yield in places where the crop is not grown. I would do this using the actual yields layer. I'm thinking that what I need to do is a geographically weighted regression, and fill in the expected values using the output of that. However, I can't use rasters as inputs for the GWR tool and I don't want to convert them to polylines as I would have to make them integer type first and I'd lose all resolution in my data. Any ideas on how I might proceed? Thanks!

Online Status	Offline
Date Last Visited	‎11-11-2020 02:24 AM

My Ideas

Latest Contributions by JeremyProville

Re: Returning names of rasters with maximum value

Re: Returning names of rasters with maximum value

Re: Returning names of rasters with maximum value

Re: Returning names of rasters with maximum value

Returning names of rasters with maximum value

Max() in ArcGIS 10

Re: OLS or GWR regression on rasters

Building missing values in a raster with another

OLS or GWR regression on rasters

Imputing raster values based on another raster