Jump to content

Photo

Spatial analysis, regression and correlation between 2 rasters


  • Please log in to reply
12 replies to this topic

#1 Arhanghelul

Arhanghelul

    Sergeant

  • Members
  • 195 posts
  • LocationKingdom of Heaven (RO, east Europe)

Posted 06 May 2013 - 10:40 AM

Hello !

I would like to make some spatial analysis regarding regressions and correlation coeficients between 2 rasters.

I have 2 rasters:
  • Land use/ Land cover extracted from aerial images
  • Land Surface Temperature (LST) derived from LANDSAT thermal band.
Both raster images are co-registered, in the same projection, having a RMSE error <=0.5 m.

Now I would like to see if there are some spatial correlations between these 2 images.
What software can I use for this task ?
I know that there are some specialised softwares like R and SPSS. But how to import and use rasters ?

Where can I find some tutorials, lab exercises or web-links about this ?

Thank you, my friends.


Arhanghelul

#2 3dbu

3dbu

    Captain

  • Members
  • 1,314 posts
  • Locationfrom hell

Posted 06 May 2013 - 11:24 AM

did you check all toolbox of spatial analyst ...regards

#3 SIGIS

SIGIS

    Private

  • Members
  • 46 posts

Posted 06 May 2013 - 12:14 PM

ENVI is the software for raster tasks as long as those two images have the same dimensions (row/line).
But it need some practice to handling it.

You can also check IDRISI a good tool

#4 yousef2233

yousef2233

    Sergeant

  • Uploader
  • 328 posts

Posted 06 May 2013 - 02:21 PM

here is the the code for spatial cross-correlation (Matlab)

there are 2 surfaces we are going to compare or find correlations
maybe you should convert them to raster first

A = imread ('...'/firstraster.jpg');
B = imread ('.../secondraster.jpg');

Result = normxcorr2(B,A);
size if B should be equal or lesser than size of A
the result is an arrey of correlation as kernel ( B moves across the Image A )

max(max(Result)) % the max xross corr found by kernel
imshow (Result) % xross surface

using ArcGIS can be very easy and helpful :)
ArcToolbox > Spatial Analyst > Multivariate > Band Collection Statistics ;)
  • HTC likes this

#5 Arhanghelul

Arhanghelul

    Sergeant

  • Members
  • 195 posts
  • LocationKingdom of Heaven (RO, east Europe)

Posted 07 May 2013 - 08:33 AM

here is the the code for spatial cross-correlation (Matlab)

there are 2 surfaces we are going to compare or find correlations
maybe you should convert them to raster first

A = imread ('...'/firstraster.jpg');
B = imread ('.../secondraster.jpg');

Result = normxcorr2(B,A);
size if B should be equal or lesser than size of A
the result is an arrey of correlation as kernel ( B moves across the Image A )

max(max(Result)) % the max xross corr found by kernel
imshow (Result) % xross surface

using ArcGIS can be very easy and helpful :)
ArcToolbox > Spatial Analyst > Multivariate > Band Collection Statistics ;)



Thank you !

#6 thelostking

thelostking

    Private

  • Members
  • 11 posts

Posted 08 June 2013 - 11:46 PM

Hi guys .. it might be late to answer this but I will explain the best way (as I can see)
first, u should use (extract by mask) which is one of the tool in ArcGIS where u can mask the first layer (i.e.LULC) using (LST) or vise versa
After finishing the masking, u can check in the properties that both layers have the same number of row and columns.
Second, use (raster to ASCII) tool which also available in ArcGIS to convert both layers into text format then open both of them using notepad or any other alternative program (I prefer Editpad Lite) and delete the header of them
Note: the header file can also show u the number of rows and columns of the resulted image.
Third, from file menu in SPSS, click on read text file and just click next until the part where u should define the number of columns where u should choose 1
do this process to the both layers then apply LR or any other algorithm u like to use.
  • HTC likes this

#7 Arhanghelul

Arhanghelul

    Sergeant

  • Members
  • 195 posts
  • LocationKingdom of Heaven (RO, east Europe)

Posted 08 June 2013 - 11:52 PM

Hi guys .. it might be late to answer this but I will explain the best way (as I can see)
first, u should use (extract by mask) which is one of the tool in ArcGIS where u can mask the first layer (i.e.LULC) using (LST) or vise versa
After finishing the masking, u can check in the properties that both layers have the same number of row and columns.
Second, use (raster to ASCII) tool which also available in ArcGIS to convert both layers into text format then open both of them using notepad or any other alternative program (I prefer Editpad Lite) and delete the header of them
Note: the header file can also show u the number of rows and columns of the resulted image.
Third, from file menu in SPSS, click on read text file and just click next until the part where u should define the number of columns where u should choose 1
do this process to the both layers then apply LR or any other algorithm u like to use.


Thank you for your reply !

Where can I find some tutorials, lab exercises for using SPSS ? (about importing and creating new projects, how to create regressions and graphics ? )

#8 HTC

HTC

    Sergeant

  • Inactive Member
  • 214 posts

Posted 09 June 2013 - 01:43 AM

Regressions in SPSS is Normal. We can choose: Linear, Logistics, Curbic (Non - Linear) for simulating correlation between 2 raster.
step1: You make sample (than 30 sample) from Land use and Land Surface Temperature.
step2: You extract value of point from sample.
Step3: You use: SPSS or Statistic software (R software, statgraphics, spss....) to make correlation.
Step4: You use Gis software for compute raster (Multiply, split, add, subtract follow equation in step3).
Good luck!
  • yousef2233 likes this

#9 yousef2233

yousef2233

    Sergeant

  • Uploader
  • 328 posts

Posted 09 June 2013 - 03:01 AM

SPSS Statistics Essential Training

http://s6.p30download.com/users/606/tutorial/office/spss-statistics/Lynda.SPSS.Statistics.Essential.Training.Full_p30download.com.part1.rar
http://s6.p30download.com/users/606/tutorial/office/spss-statistics/Lynda.SPSS.Statistics.Essential.Training.Full_p30download.com.part2.rar
pass = www.p30download.com
disable your antivirus, there are nothing wrong with the files.
  • HTC likes this

#10 Arhanghelul

Arhanghelul

    Sergeant

  • Members
  • 195 posts
  • LocationKingdom of Heaven (RO, east Europe)

Posted 09 June 2013 - 01:25 PM

SPSS Statistics Essential Training
 

http://s6.p30download.com/users/606/tutorial/office/spss-statistics/Lynda.SPSS.Statistics.Essential.Training.Full_p30download.com.part1.rar
http://s6.p30download.com/users/606/tutorial/office/spss-statistics/Lynda.SPSS.Statistics.Essential.Training.Full_p30download.com.part2.rar
pass = www.p30download.com
disable your antivirus, there are nothing wrong with the files.

 

The links are not working :(

 

I dont understand arabic language. :blink:



#11 yousef2233

yousef2233

    Sergeant

  • Uploader
  • 328 posts

Posted 09 June 2013 - 09:19 PM

works for me

http://irs6.p30download.com/users/606/tutorial/office/spss-statistics/Lynda.SPSS.Statistics.Essential.Training.Full_p30download.com.part1.rar
http://irs6.p30download.com/users/606/tutorial/office/spss-statistics/Lynda.SPSS.Statistics.Essential.Training.Full_p30download.com.part2.rar

  • HTC likes this

#12 HTC

HTC

    Sergeant

  • Inactive Member
  • 214 posts

Posted 10 June 2013 - 12:47 AM

Yes, the file is live. We can disable virus before downloading it!



#13 rahmansunbeam

rahmansunbeam

    Sergeant

  • GIS Analyst
  • 478 posts
  • LocationDhaka

Posted 10 June 2013 - 10:03 AM

http://thepiratebay.sx/torrent/6688193/LYNDA.COM_-_SPSS_STATISTICS_ESSENTIAL_TRAINING_[thethingy]

torrent is power !


  • Arhanghelul likes this




0 user(s) are reading this topic

0 members, 0 guests, 0 anonymous users