jack_stepeny

Members
  • Content count

    101
  • Joined

  • Last visited

  • Days Won

    5

jack_stepeny last won the day on March 14

jack_stepeny had the most liked content!

Community Reputation

58 Best

1 Follower

About jack_stepeny

  • Rank
    Lecturer
  • Birthday 03/05/1984

Profile Information

  • Gender
    Male

Recent Profile Visitors

154 profile views
  1. Feel free to contact me if u encounter any problem during running the script.
  2. Here is code.. Before executing this code make sure raster library is installed in R library(raster) #setting working directory(folder path containing all rasters) setwd("G:\\new_gl_data\\4har") #listing all raster in a folder rlist=list.files(getwd(), pattern=".tif$", full.names=F) #change year here (2001:2015) for (yr in 2001:2015){ nlist = rlist[grep(yr,rlist)] stack_ras = stack(nlist) sum_ras = sum(stack_ras) sum_brk = brick(sum_ras) writeRaster(sum_brk, filename=paste("prc_",yr,sep=""), format="HFA",overwrite=TRUE) rm() gc() }
  3. If the extent of all rasters are same.. it can be done easily. 1. List all rasters in a folder 2. Loop through list to select rasters from particular year. ( It can be done by matching the year with the name of raster. As u have mentioned name include year) 3. Then sum of all the raster can be done. If u r comfortable with R I can share code that I have already written to accomplish this task.
  4. U have all the raster in same folder??? Is there any number in the raster names that can be used to identify raster which belongs to Jan which belongs to dec etc. Like MODIS follows naming convention and we can easily identify the which raster belongs to which months.. Please provide more details so that I can help you.
  5. Actually I want to use R STL (Seasonal Trend Decomposition by LOESS) function on each pixel of MODIS Time Series. Currently I am able to run STL on single pixel value which is being read from a csv file. I don't know how to loop over each pixel and get the output in three rasters say raster 1 for trend, raster 2 for seasonal component and raster 3 for remainder. Any help in this regard will be appreciated
  6. If somebody knows how to generate orthophoto from stereopair in ERDAS LPS please post step by step tutorial here. I will be highly obliged for the help.
  7. Please some help in creting Bare Earth DEM from stereopair. When I am creating DEM using stereopair it include the height of buildings and trees. I want to create DEM without trees and buildings heights.
  8. I have created a NDVI image from LANDSAT TM data in erdas. Now I want to create a Layout for the NDVI image as well as need to reclassify the NDVI image in ArcGIS using Natural Breakings. My problem is that when I opt for the option No Color for the Background it also remove the 0 values from the NDVI image. I simply want to remove 0 background values and keep the 0 values from the NDVI image. Please suggest how to do that.
  9. thanx for the information.
  10. @yousef2233 Exactly what you explained. bcoz of coordinates in arcgis parcels number which I m generating are in sequence but they starts from left bottom corner and goes to right. I just want to do it from upper left corner. right now i m using script in field calculator. Got this script from ESRI forum. Static rec As Long Dim pStart As Long Dim pInterval As Long pStart = 1 pInterval = 1 If (rec = 0) Then rec = pStart Else rec = rec + pInterval End If if you have any script please share.
  11. I have to generate unique Id sequentially for each parcel in a grid in ArcGIS. In each grid numbering will start from upper left corner to right corner then again numbering will start from left just like we type next line in microsoft word. Any help will be appreciated.