Re: [R-sig-Geo] Anisotropic point processes are inhomogeneous?

2017-11-28 Thread Domenico Giusti 

Thank you Ege for your reply,


In general it seems like a difficult problem to separate inhomogeneity and 
anisotropy from each other


this point is indeed not very clear to me.

If the homogeneous Poisson point process is defined by the properties of 
i) homogeneity, ii) independence and iii) Poisson distribution; the 
inhomogeneous Poisson process is a modification of the former with i) 
intensity function, ii) independence and iii) Poisson distribution.


An anisotropic point process, defined as directional dependent process, 
could be seen as inhomogeneous process having an oriented intensity 
function, or violates the assumption of independence in both homogeneous 
and inhomogeneous processes?


Similarly, is anisotropy violating the assumptions of independent 
labels/components in inhomogeneous multitype point processes?


Thanks,

On 11/21/2017 10:41 PM, Ege Rubak wrote:
stationarity of your process I would think that the results from 
something like Ksector can be misleading.


___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] Anisotropic point processes are inhomogeneous?

2017-11-21 Thread Domenico Giusti 

Dear all,

I have an anisotropic point pattern, analyzed by means of directional 
correlograms, wavelets and Ksector functions.


I wonder if anisotropy prevents the use of inhomogeneous K-functions.

Thank you in advance,

--
Domenico Giusti <dnc...@gnewarchaeology.it>
GPG keyID: 2048R/A3AB7054F6E5D778

___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Re: [R-sig-Geo] Wavelet analysis of anisotropy in point patterns

2017-07-19 Thread Domenico Giusti 

Thank you Marcelino,

I already had a look at those functions and, although they do the job, I 
was wondering how to analyze anisotropy in point pattern using wavelet 
analysis.


I tried to use the WaveletComp package, but with no success. It looks to 
me that packages for wavelet analysis work fine with temporal data, but 
I couldn't find the way to use it with spatial point patterns.


How could I get in R the same analysis provided by the PASSaGE software?

Thank you!

On 07/18/2017 06:17 PM, Marcelino de la Cruz Rot wrote:

Hi Domenico,

In spatstat you have several functions to detect anysotropy in point 
patterns:

- the sector K-function, Ksector();
- the pair orientation distribution, pairorient();
- the anisotropic pair correlation function, Kmeasure().

See the help pages of these functios or, even better, consult pages 
236:242 in  Adrian Baddeley, Ege Rubak, Rolf Turner (2015). Spatial 
Point  Patterns: Methodology and Applications with R. London: Chapman 
and   Hall/CRC Press, 2015.
http://www.crcpress.com/Spatial-Point-Patterns-Methodology-and-Applications-with-R/Baddeley-Rubak-Turner/9781482210200/ 



Cheers,

Marcelino

El 18/07/2017 a las 16:30, Domenico Giusti escribió:

Hi all,

I'm running wavelet analysis of anisotropy in point patterns using 
the PASSaGE software by Rosenberg et al.


Rosenberg, M.S., and C.D. Anderson (2011) PASSaGE: Pattern Analysis, 
Spatial Statistics and Geographic Exegesis. Version 2. Methods in 
Ecology & Evolution 2(3):229-232.

http://www.passagesoftware.net

Rosenberg, M.S. 2004. Wavelet analysis for detecting anisotropy in 
point patterns. Journal of Vegetation Science 15:277-284.

DOI: 10./j.1654-1103.2004.tb02262.x

I wonder how I can do the same analysis using R.

Thanks,

Dome





___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] Wavelet analysis of anisotropy in point patterns

2017-07-18 Thread Domenico Giusti 

Hi all,

I'm running wavelet analysis of anisotropy in point patterns using the 
PASSaGE software by Rosenberg et al.


Rosenberg, M.S., and C.D. Anderson (2011) PASSaGE: Pattern Analysis, 
Spatial Statistics and Geographic Exegesis. Version 2. Methods in 
Ecology & Evolution 2(3):229-232.

http://www.passagesoftware.net

Rosenberg, M.S. 2004. Wavelet analysis for detecting anisotropy in point 
patterns. Journal of Vegetation Science 15:277-284.

DOI: 10./j.1654-1103.2004.tb02262.x

I wonder how I can do the same analysis using R.

Thanks,

Dome

--
Domenico Giusti <dnc...@gnewarchaeology.it>
GPG keyID: 2048R/A3AB7054F6E5D778

___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Re: [R-sig-Geo] R 3.3.0 - spatstat 1.45-2 - envelope -- Error in nrank%%1 : non-numeric argument to binary operator

2016-05-18 Thread Domenico Giusti 
Thank you!

P.S. I use to type the options, to remind me what I'm doing :)

On 05/17/2016 11:08 PM, Rolf Turner wrote:
> On 17/05/16 23:44, Domenico Giusti wrote:
>> Dear all,
>>
>> I get "Error in nrank%%1 : non-numeric argument to binary operator"
>> running
>>
>>> data("amacrine")
>>> E <- envelope(amacrine, Kcross, "on", "off", nsim=19, global=FALSE)
>>
>> Envelope works fine running
>>
>>> E <- alltypes(amacrine, Kcross, nsim=19, envelope=TRUE, global=FALSE)
>>
>> Thanks,
> 
> Try:
> 
> E <- envelope(amacrine, Kcross, funargs=list(i="on",j="off"),
>   nsim=19,global=FALSE)
> 
> cheers,
> 
> Rolf Turner
> 
> P. S.  Note that you *do not need* to set "global=FALSE"; this is the
> default.
> 
> R. T.
> 
> 

-- 
Domenico Giusti <dnc...@gnewarchaeology.it>
GPG keyID: 2048R/A3AB7054F6E5D778

___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] R 3.3.0 - spatstat 1.45-2 - envelope -- Error in nrank%%1 : non-numeric argument to binary operator

2016-05-17 Thread Domenico Giusti 
Dear all,

I get "Error in nrank%%1 : non-numeric argument to binary operator" running

> data("amacrine")
> E <- envelope(amacrine, Kcross, "on", "off", nsim=19, global=FALSE)

Envelope works fine running

> E <- alltypes(amacrine, Kcross, nsim=19, envelope=TRUE, global=FALSE)

Thanks,

-- 
Domenico Giusti <dnc...@gnewarchaeology.it>
GPG keyID: 2048R/A3AB7054F6E5D778

___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Re: [R-sig-Geo] How to get the borders of scattered spacial points

2015-11-19 Thread Domenico Giusti 
Hi,

convexhull is probably what you need.

Best,

On 11/18/2015 11:56 PM, Zhong-Yuan Zhang wrote:
> Dear All:
> 
>As a freshman, I am now analyzing some spacial data.
> 
> I want to get the border of scattered spacial points. Are there
> 
> any libraries that I can use?  Also is there any library that can
> 
> use Baidu map API?  I highly appreciate your help and
> 
> suggestions.
> 
>    Best Regards Always.
> 

-- 
Domenico Giusti <dnc...@gnewarchaeology.it>
GPG keyID: 2048R/A3AB7054F6E5D778

___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] R 3.1.3 - GNU/Linux - readOGR -- Layer not found

2015-03-30 Thread Domenico Giusti 
Dear list,

I cannot read some of the layers (table and view) from a
Postgres/PostGIS DB.

Error in ogrInfo.

 ogrListLayers(PG:dbname=dbname)

doesn't list all the tables and views. e.g. I have two views with both a
3D point geom, but only one of them is in the list.

What could be the problem?

Thanks,

 sessionInfo()
R version 3.1.3 (2015-03-09)
Platform: i486-pc-linux-gnu (32-bit)
Running under: Debian GNU/Linux 7 (wheezy)

locale:
 [1] LC_CTYPE=en_US.UTF-8   LC_NUMERIC=C
LC_TIME=en_US.UTF-8
 [4] LC_COLLATE=en_US.UTF-8 LC_MONETARY=en_US.UTF-8
LC_MESSAGES=en_US.UTF-8
 [7] LC_PAPER=en_US.UTF-8   LC_NAME=C  LC_ADDRESS=C

[10] LC_TELEPHONE=C LC_MEASUREMENT=en_US.UTF-8
LC_IDENTIFICATION=C

attached base packages:
[1] stats graphics  grDevices utils datasets  methods   base

other attached packages:
[1] rgdal_0.9-2 sp_1.0-17

loaded via a namespace (and not attached):
[1] grid_3.1.3  lattice_0.20-30 tools_3.1.3

-- 
Domenico Giusti dnc...@gnewarchaeology.it
GPG keyID: 2048R/A3AB7054F6E5D778

___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo