Latest Contributions by ShouvikJha

Reply from Dan Patterson‌ If the rasters aren't huge convert to points, then Ordinary Least Squares (OLS)—Help | ArcGIS for Desktop I think the report allows for a graph ... View more

Thank you. I am not getting your points. i am trying plot two variable of raster like this ...  ... View more

Reply from Dan Patterson‌ before you go off into a side tangent of graphing and trying to 'correlate', 'predict' or extrapolate see...    Geographically Weighted Regression (GWR)—Help | ArcGIS for Desktop  Only if you are working with interval/ratio and there is some real predictive behaviour between the two variables.  Looking for correlations, relationships between two variables via a magical statistical test is overused.  But if you are looking for an association then the above and other substitutes will save a flurry of activity in the wrong direction   For nominal data...  Combine—Help | ArcGIS for Desktop   gives you the basis for tests of Association such as the chi-squared test.  And at worst it gives you a tabular output for discussion   Make your intent clear and specific... otherwise people with keep throwing suggestions out that will meet vague goals ... not necessary specific ones. ... View more

Dan Patterson‌  Thank you. I want to do  on  cell by cell basis. I have two different parameters of raster data, so one parameter i want to plot  in X axis and another one Y axis to make it relation between two parameters.  How do i convert it in row and column  format. can i convert it in  ASCII format?  Data is tiff format same as which we got from our above scrip suggested by xander_bakker‌ ... View more

xander_bakker‌, Thank you for the suggestion. I am extremely sorry. Next time surely i will follow your suggestion. Sorry for the inconvenience.  ... View more

Thank you. I am not getting your points. i am trying plot two variable of raster like this ...    ... View more

xander_bakker‌. Thank you very much. Above Code running perfectly. Regarding  error issue which you got different error massage, because i corrected that line before running script.   ... View more

Thank you very much  all of you for your kind cooperation. Xander Bakker  I ran the above code , but i am getting error. I think not having problem with code. Might be we need more explanation what i am trying to do, please follow the below  Message     File Name     Line     Position     Traceback                        <module>     <module1>     20              CellStatistics     C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py     2951              swapper     C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Utils.py     53              Wrapper     C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py     2947          RuntimeError: ERROR 999998: Unexpected Error.                    ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ I am giving here  what I am trying to do, I have three sub-folder (Name : NPP_2008, NPP_2009, NPP_2010), NPP_2008      NPP_2009       NPP_2010       MEAN_NPP (output Folder) r001_NPP.tif  r001_npp.tif   r001_NPP.tif   r001_MEAN_NPP.tif (mean of all r001_NPP) r002_NPP.tif r002_npp.tif    r002_NPP.tif   r002_MEAN_NPP.tif (mean of all r002_NPP) So on .............. Every sub-folder hold monthly raster data for different year with the name of e.g r001_NPP, r002_NPP……..r012_NPP.  So I want calculate mean of those raster which is r001_NPP series in every sub-folder, similarly program will loop through till r012_NPP series of all raster to calculate mean and save it another subfolder (Name : Mean_NPP) with output name 001_MEAN_NPP(which is mean of all r001_NPP), 002_MEAN NPP…..012_MEAN_NPP. Here with i have attached reference raster which i am using for my raster calculation ... View more

Dan , thanks . I have solved the indentation error, and i updated my code but still i could not solve the problem, updated code below  error masage ile "<module1>", line 18     sum.save(os.path.join(out_workspace, 'r{0}_MEAN_NPP.TIF'.format(ras_num)                                                                            ^ SyntaxError: invalid syntax‍‍‍‍ Updated Code  import arcpy arcpy.CheckOutExtension("spatial") from arcpy.sa import * # generate list of rasters workspace = r"D:\TEST" out_workplace = "D:\TEST\MEAN_NPP" for ras_num in [str(i).zfill(3) for i in range(1, 12)]: rasters = [] walk = arcpy.da.Walk(workspace, datatype="RasterDataset", type="TIF") for dirpath, dirnames, filenames in walk:     for filename in filenames:      rasters.append(os.path.join(dirpath, filename)) # calculate mean sum = CellStatistics(rasters, "MEAN") sum.save(os.path.join(out_workspace, 'r{0}_MEAN_NPP.TIF'.format(ras_num)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

xander_bakker‌, We need your cooperation to solve above problem. Kindly cooperate  with us.  ... View more

Dan . thanks , you are right. I didn't post the error massage. Line no 13 i am getting the Indentation Error massage. I tried to fix it but not solved yet. error massage below  File "<module1>", line 13     rasters.append(os.path.join(dirpath, filename))           ^ IndentationError: expected an indented block ‍‍‍‍ ... View more

Thank you. I have updated my suggested by you, problem still persist, i am not understanding how to iterate through each folder, take all series of 001_NPP from all folder and calculate average oh them and save it as 001_NPP_MEAN, below is my updated code and along with code i have attached reference raster also which i am using for my raster calculation,  import arcpy arcpy.CheckOutExtension("spatial") from arcpy.sa import * # generate list of rasters workspace = r"D:\TEST" out_workplace = "D:\TEST\MEAN_NPP" rasters = [] walk = arcpy.da.Walk(workspace, datatype="RasterDataset", type="TIF") for dirpath, dirnames, filenames in walk:     for filename in filenames:     rasters.append(os.path.join(dirpath, filename)) # calculate mean sum = CellStatistics(rasters, "MEAN") sum.save(os.path.join(out_workspace, "npp_sum.tif")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more