|
POST
|
This question has been branched from this question: https://community.esri.com/thread/198969-trend-analysis-through-time-series-of-raster-data to give it more visibility Dear Xander Bakker and Dan Patterson thank you for the awesome discussion. I have different case. Let say I have raster timeseries for monthly rainfall anomaly data from 1981-2017 for 1 province. And then I have sea surface temperature anomaly data for the same periods and from 1 location only (text data). I want to see the correlation between both data (rainfall and temperature anomaly) for each pixel in raster data. Is it possible to do this using above script? I want to have sets of raster (slope, intercept, r value, p value and std error) as an output. From this output I hope I can see if the temperature increase in the Pacific, which area will experience more/less rainfall. Thanks, Benny
... View more
11-01-2018
08:12 PM
|
0
|
6
|
3683
|
|
POST
|
Awesome! Yes, I confirmed it takes less than 30 seconds. I tested using Indonesia data (whole country), and modified line 67 and 74. I guess its the lower left coordinates. Thank you very much Dan Patterson you made my day. FYI, I will use this script to see the relation between max rain+date, flood frequency and extreme rainfall more than threshold (generated using recurrence interval). This question is still related with my previous question about raster percentile from timeseries data
... View more
07-30-2018
07:02 PM
|
0
|
1
|
1015
|
|
POST
|
Perfectly match with manual process using CellStatistics MAX. Nearest mm its acceptable. Thank you Dan!
... View more
07-30-2018
03:08 AM
|
0
|
1
|
3199
|
|
POST
|
Done! I try to make it small, its only daily rainfall for Bali island 1981-1982
... View more
07-30-2018
01:23 AM
|
0
|
0
|
3199
|
|
POST
|
Thank you Dan, I understand the logic but still hard for me to think the next step. Let me try to do it manually for 1 year using Excel, and will see how long it takes to do it.
... View more
07-29-2018
08:26 PM
|
0
|
2
|
3199
|
|
POST
|
Hi Dan, thank you for reply I am not familiar with numpy and friend, I can only use script, without knowledge write a code. You and Xander Bakker always helping me solve my previous problem My latest question to you and Xander about "percentile from raster timeseries". From your suggestion, which tool should I use to compare? Is it doable using Model Builder raster iteration? Benny
... View more
07-29-2018
08:11 AM
|
0
|
4
|
3199
|
|
POST
|
Hi, I have daily rainfall data in a year in geotiff format with file naming convention "rainfall_yyyymmdd.tif" and I would like to get maximum rainfall value per year per pixel and information about which date the max value came from. Getting the max value is doable using Cell Statistics, but how to get information about date of the max value? I am planning to do the calculation for 36 years data (1981 - 2017). Thank you, Benny raster timeseries arcpy python cellstatistics precipitation rainfall
... View more
07-29-2018
06:36 AM
|
0
|
12
|
4895
|
|
POST
|
Thank you Dan, finally I managed to add the spatial ref. I added info on line 23-27 to line 30 out = arcpy.NumPyArrayToRaster(n_90, arcpy.Point(114.4317366,-8.8492618), x_cell_size=0.049242234706753, y_cell_size=0.049242234706753, value_to_nodata=0)
... View more
01-17-2018
08:02 AM
|
0
|
1
|
818
|
|
POST
|
Hi Dan, Pardon of my ignorance. After reading carefully, I can run the script. For line 22-30, I added arcpy.NumPyArrayToRaster(n_90,
arcpy.Point(114.4317366,-8.8492618),
0.049242234706753,
0.049242234706753,
0)
out = arcpy.NumPyArrayToRaster(n_90, value_to_nodata=0) But the geotiff file still have not spatial reference. Which part I missed it? Thank you.
... View more
01-17-2018
03:47 AM
|
0
|
3
|
6128
|
|
POST
|
One more thing, the TIF output don't have any spatial reference. How to make it available automatically on the output during the calculation, so it's ready to overlay with other data.
... View more
01-17-2018
02:23 AM
|
0
|
5
|
6128
|
|
POST
|
Wow excellent! The result almost same with manual version, slightly different on the minimum value (your script is 109.391 and manual version is 105.947, for the max is same. See pictures attached (map and statistics info from gdal). The raster extent is not exactly same with yours, because its from excel - convert to point - convert to raster and clip with boundary. Here's the file: n_90_manual.zip And yes you are correct, its pretty coarse grid with 0.05 deg spatial resolution. The data came from CHIRPS global precipitation data CHG - Data and this data is the best available global precipitation data with temporal resolution daily, 5days, 10days, monthly, 2monthly, 3monthly and annual from 1981 until now and free for public. So far this data are very useful for me when getting weather station data from govt is very difficult. Thank you Dan Patterson
... View more
01-16-2018
10:48 PM
|
0
|
7
|
6128
|
|
POST
|
Thank you Dan Patterson and Xander Bakker for explanation. Unfortunately it's difficult for me to modified above python script with numpy. So far I have try using gdal, I found a script from StackExchange "gdal_calc.py -A stack.vrt allBands=A --calc='nanpercentile(A.astype(int16),85,axis=0)' --outfile out.tif" and arcpy script mentioned in this discussion Pool of raster values to calculate percentile I am getting an output using both script, however looks entirely wrong if I compared using manual calculation using Excel (I convert all raster data to csv, and calculate the percentile for each point using formula below) The idea is to get the rainfall threshold and rainfall-triggered flood with return period, let say for 10 and 25 year. If my data is in excel format, I used =PERCENTILE(DATARANGE,0.9) for calculating the rainfall threshold for 10 year return period of rainfall-triggered flood. I have upload example data of rainfall below. Rainfall Timeseries 1981 - 2016
... View more
01-16-2018
05:52 PM
|
0
|
0
|
6128
|
|
POST
|
Branched from: https://community.esri.com/thread/95403 Hi Xander Bakker thank you for above script. I have tested and it works using my data. But I have question. Let say I have timeseries raster data from 1981 to 2016, its the max of daily rainfall in a year, so its just 1 raster data for 1 year. For each cell in x,y location from 1981 to 2016, I would like to calculate the percentile of i.e 90%. How to calculate the percentile for each cell from timeseries raster using your script?
... View more
01-16-2018
01:26 AM
|
0
|
12
|
9804
|
|
POST
|
Hi, I have calculate the SPI using monthly CHIRPS CHG - Data (1981 - present). Here's the step: 1. Calculate rainfall average and standard deviation by month (I used 1981 - 2014 as a base period), you can used Cell Statistics tool in the Spatial Analyst toolbox, and choose MEAN and STD to calculate the average and standard deviation for JAN, FEB, ..., DEC 1981 - 2014. 2. Calculate the SPI (example Jul 2015): SPI = (Rainfall Jul 2015 - Rainfall AVG Jul 1981-2014) / Standard Deviation of Rainfall in Jul 1981-2014
... View more
08-11-2015
02:03 AM
|
0
|
1
|
4137
|
|
POST
|
Hi, I use Windows 7 Pro 64-bit running in Parallels Desktop 10 for Mac and using Microsoft Office 2013 64bit. Now just installed 64-bit version of ESRI Maps for Office 3.0 and I can't log in to the ESRI maps in Microsoft Excel 2013 - The sign in window always loading. (See attached picture). But when I try to log in Microsoft Powerpoint 2013, it's worked. Any thoughts on what's going on here? Best, Benny
... View more
06-19-2015
12:12 AM
|
0
|
1
|
3040
|
| Online Status |
Offline
|
| Date Last Visited |
04-07-2022
04:30 AM
|