Dear R-Sig-Geo Members,
?? I create a virtual machine (VM) in Google Cloud with Ubuntu 18.04
with 8 CPU and 30 RAM memory and R 3.6.0 version, but I try to improve
my spatial analysis without success or same a more faster process. If I
use packages snow and doMC with all the 8 CPU's in a
Dear R-Sig-Geo Members,
?? Many years ago I used enfa() function of adehabitat package with my
*asc files with environmental variables and my script works very well,
but now adehabitat is deprecated and there are many versions of
adehabitat. I install the adehabitatHS but for doesn't work
Thanks for your help Marcelino e for the careful explanation Rolf Turner
and so sorry about my post script configuration, I expected that I
solved that in my new post.
First my variable area is a marked point (attribute or auxiliary
information about my point process - page 7 and not a spatial
Dr. Rolf Turner,
Thanks again and your explanations as a personal spatial course for
me :) You are very clear, but there are some questions yet. First, these
data set is not artificial, I have real coordinates and area of nests of
leaf-cutting ants in my example (in my complete data set I
Dear members,
I've like to create a raster using my event coordinates and information
about the size in this coordinates as a rule for my raster creation.
First, I make:
library(raster)
library(spatstat)
#Points coordinates in UTM
xp<-c(371278.588,371250.722,371272.618,371328.421,371349.974,
3
Dear R-Sig-Geo Members,
I have the hypothetical point process situation:
library(spatstat)
set.seed(2019)
A <- rpoispp(100) ## First event
B <- rpoispp(50) ## Second event
C <- rpoispp(50) ## Third event
plot(A, pch=16)
plot(B, col="red", add=T)
plot(C, col="blue", add=T)
I've like to know an a
; boxplot(permout$ABd - permout$ACd)
> points(1, mean(ABd) - mean(ACd), col="red")
>
> table(abs(mean(ABd) - mean(ACd)) >= abs(permout$ABd - permout$ACd))
> # FALSE TRUE
> # 573 426
>
> sum(abs(mean(ABd) - mean(ACd)) >= abs(permout$ABd - permout$ACd)) / npe
Dear r-sig-geo members,
I've to know if there is a function for sample for extract coordinates
and marks of a marked planar point pattern object? In my example:
library(spatstat)
## Using a longleaf data set
longleaf
#Marked planar point pattern: 584 points
#marks are numeric, of storage type
Dear R-Sig-Geo members,
I've like to fitting an Inhomogeneous Poisson model and I have one
question about better approaches and/or steps order for this.
### Simulating Inhomogeneous Poisson process
library(spatstat)
cI <- rpoispp(function(x, y) {200 * x^2}, lmax = 1000)
cII <- rpoispp(function(x
Dear r-sig-geo members,
??? I?ve like to know the possibility of make some modifications in
resample() function in raster package. First, in "bilinear" method by
default is assigns a weighted average of the four nearest cells, and
I?ll like to change for five and six nearest cells, too, is pos
https://www.eecbio.ufg.br/)
> Federal University of Goiás | Brazil
> RG: https://www.researchgate.net/profile/Frederico_Faleiro
>
>
> Em seg., 23 de mar. de 2020 às 10:41, ASANTOS via R-sig-Geo
> mailto:r-sig-geo@r-project.org>> escreveu:
>
> Dear r-sig-geo members,
>
>
Dear r-sig-geo Members,
??? I'll like to read several shapefiles, count occurrences in the same
coordinate and create a final shapefile with a threshold number of
occurrences. I try to convert the shapefiles in ppp object (because I
have some part of my data set in shapefile and another in ppp
.N), by=.(x,y)]
> res.xy1 = res1[target_sub1, on=c("x","y")]
>
> all.equal(res.xy, res.xy1, check.attributes=FALSE) # should return TRUE
>
> If you're using SQL then you just join the raw table with the grouped table
> and you should get the table coordinates
Dear R-Sig_Geo Members,
I'd like to find a more simple way to filtering a set of points in a
"ppp" object by minimum or maximum distance.
In my example using a ants ("ppp" object) in spatstat package:
#Packages
library(spatstat)
library(sp)
#Point process example
data(ants)
str(ants)
#- attr(*
Dear R-Sig-Geo Members,
I'd like to find any way to filtering a set of points in a "ppp" object by
minimum distance just only between different marks. In my example:
#Package
library(spatstat)
#Point process example - ants
data(ants)
ants.ppp<-ppp(x=ants$x,y=ants$y,marks=rep("antNests",length(a
/researcher/3085587/alexandre-dos-santos/
--
Em 16/09/2020 03:18, Marcelino de la Cruz Rot escreveu:
Hi Alexandre,
may be this?
ddd <- nndist(insects.ppp, by=factor(insects.ppp$marks))
subset(insects.ppp, marks=="termiNests" & ddd[,"termiNests"] >20)
Cheers,
Ma
publons.com/researcher/3085587/alexandre-dos-santos/
--
Em 16/09/2020 03:18, Marcelino de la Cruz Rot escreveu:
Hi Alexandre,
may be this?
ddd <- nndist(insects.ppp, by=factor(insects.ppp$marks))
subset(insects.ppp, marks=="termiNests" & ddd[,"termiNests"] >20)
Dear R-sig-geo Members,
I'd like to use de extract() function using the raster package for
calculate a proportion of pixel with "1"s values inside a buffer given
some coordinates in a raster. I try to create a function for this
without success, in my hypothetical example:
#Package
library(rast
Dear r-sig-geo Members,
I'd like to find any way to create 1 (total 9 pixels) and 2 pixels
(total 25 pixels) surrounding the neighborhood of each pixel (ant) in my
plot image (antscount).
In my example:
#Packages
library(spatstat)
library(raster)
#Selection of ants data set
data(ants)
geo.f
t;antscount" as the same
variable. Maybe? It's not clear to me.
Finally, when you call raster::adjacent() you haven't defined the value for
the cells argument when computing e1 and e2. What do you want them to be?
I wish I could be more helpful.
Cheers,
Ben
On Mon, Nov 9, 2020 at 1
oking "antscounts" and "antscount" as the same
variable. Maybe? It's not clear to me.
Finally, when you call raster::adjacent() you haven't defined the value for
the cells argument when computing e1 and e2. What do you want them to be?
I wish I could be more hel
Dear r-sig-geo Members,
I'd like to calculate the centroid of each cluster (considering points
with < 10m distance to the same cluster) of spatial points using mean
operation (spdplyr package) for coordinates and another operation (sum)
for the attribute(area) without success.
In my example:
Dear members,
I will try to calculate each polygon percentage inside a circles
given an arbitrary radius in a shapefile object with the code below and
my output needs to be (Two first columns with center os circle
coordinates and values of each polygon percentage):
"pts$x" "pts$y" "I
Dear Rolf Turner,
It's much better a clean code with a minimum packages, thank
you very much for your answer. But "pct" object give me a total polygon
percentage around each point and I need too an identification (in
columns) of individual contribution of each polygon. In my simulat
Dear Rolf Turner,
Amazing! I will never put R in doubt.
My problem was solved with success, thanks very much again!!!
Best wishes,
Alexandre
--
==
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto
Dear R-sig Member
I'd like to extract mean pixels values using x,y coordinates by
defining a buffer or area surrounding individual point coordinate. My
circular sample represents the radius of a circular region around each
point. I try unsuccessfully this using extract function with rast
Dear Members,
I created a georeferenced PDF using a shapefile and OK no problem, but
I'd like to know if is possible to insert a title inside a output pdf
(in top position "My geo PDF"), points of interest (pts.sampling) and
values too.
In my example:
#Packages
|require(rgdal)require(maptoo
Dear R Sig Geo Members,
I've like to know if there are any function in any package for
estimation density in a marked point process (e.g. geographic position
and size of ants nests in square meters). My goal will be represents the
density of ants nest estimated, but use nests sizes as cova
via R-sig-Geo wrote:
Dear R Sig Geo Members,
I've like to know if there are any function in any package for
estimation density in a marked point process (e.g. geographic
position and size of ants nests in square meters). My goal will be
represents the density of ants nest estimated
Dear R-Sig-Geo Members,
I've like to extract xy coordinates from raster of 24 neigborhood
cells that surround the 4 given points (pts_s) pixels and for this I used:
#Packages
require(raster)
require(sp)
## Create a raster
r <- raster(nc=30, nr=30)
r <- setValues(r, round(runif(ncell(r))*
- data.frame(xd,yd)
> pts_s<- SpatialPoints(pts)
> projection(pts_s) <- projection(r)
> points(pts_s, col="red", pch=16)
>
> ## Find pixels center of each point
> N_cells <- cellFromXY(r, pts_s)
>
> e <- adjacent(r, N_cells , directions=
Dear R-sig-geo Members,
I've like to identify hotspots points (centroid/geometric
center of distances(r) when CSR is not satisfied), in my study case,
centroids with points around 0.75 radius. This thinking in the map
representation, for this objective I make:
#Package
library(spatst
Dear R-sig-geo Members,
I've like to create Kernel image results as *tif using an object of
density() function output in spatstat package. But in my example,
doesn't work when I try:
#Packages
library(spatstat)
library(raster)
library(rgdal)
#Swedishpines's data set in spatstat package
d
,
instead of converting the image to a SpatialPixelDataFrame converting
to a raster object might be an alternative.
r_pines<-raster(d_pines)
writeRaster(r_pines, "Pines.tif")
Regards,
Florian
Am 26.09.2018 um 22:25 schrieb ASANTOS via R-sig-Geo:
Dear R-sig-geo Members,
I'
Dear R-sig-geo Members,
I've like to use R and raster package (extract() function) for
pixel calculations (SD, Skewness and Kurtosis in 1 unit radius of a
buffer) in a geoTiff raster, but my original geoTiff image has
6244(nrow), 8721(ncol) and 54453924 (ncell) dimensions. This image size
35 matches
Mail list logo