Dear Members,

        I want to make a data set with 1 for presence in point location 
(in xpoint and ypoint coordinates), 0 for absence of point location 
(inside of window W) and NA for points outside the window for each 
coordinate in my syn.ppp object of spatstat package.//I try to create a 
geotif for this//, but I have some problems that involves ppp objects 
and writeGDAL() function in rgdal package and I need members ideas for 
make this,  in my example:

require(spatstat)
require(sp)


# Two polygons 
coordinates------------------------------------------------------

#Polygon 1
x.coords1<-c(371299.9,    371266.4,    371205.6,    371111.8, 
371047.6,    371018.2,    371014.0,
               371009.3,    370983.1,    370919.7,    370853.6, 
370785.6,    370748.8,    370711.8,
               370687.8,    370696.4,    370785.9,    370885.5, 
371035.8,    371148.1,    371205.2,
               371231.7,    371236.5,    371240.3,    371285.8, 
371326.5,    371397.2,    371417.1,
               371432.9,    371445.0,    371455.7,    371466.4, 
371476.6,    371502.6,    371536.0,
               371550.0,    371546.8,    371528.3,    371470.0, 
371393.3,    371299.9,    371299.9)

y.coords1<-c(8246589,    8246560,    8246508,    8246428, 8246373,    
8246349,    8246348,
              8246352,    8246385,    8246465,    8246551, 8246638,    
8246685,    8246732,
              8246764,    8246771,    8246846,    8246932, 8247062,    
8247160,    8247209,
              8247230,    8247224,    8247221,    8247160, 8247107,    
8247016,    8246991,
              8246967,    8246939,    8246914,    8246892, 8246875,    
8246846,    8246821,
              8246809,    8246802,    8246785,    8246735, 8246669,    
8246589,    8246589)

#Polygon 2

x.coords2<-c(368382.9,    368399.4,    368394.1,    368464.7, 
368652.2,    368683.7,    368699.7,
              368714.9,    368714.9,    368714.9,    368860.9, 
368932.1,    368994.8,    368994.8,
              369015.7,    369014.5,    369014.5,    369026.8, 
369045.2,    369020.0,    368952.9,
              368951.4,    368952.9,    368951.4,    368936.8, 
368923.4,    368804.1,    368676.0,
              368178.1,    368182.4,    368202.2,    368216.2, 
368233.4,    368251.3,    368270.3,
              368281.7,    368300.6,    368323.8,    368344.8, 
368364.6,    368374.2,    368382.9)



y.coords2<-c(8249120,    8249150,    8249153,    8249244, 8249173,    
8249167,    8249154,
                 8249145,    8249145,    8249145,    8249048, 
8249026,    8249003,    8249003,
              8248993,    8248991,    8248991,    8248984, 8248970,    
8248889,    8248833,
              8248829,    8248833,    8248829,    8248790, 8248765,    
8248626,    8248526,
              8248774,    8248780,    8248815,    8248843, 8248871,    
8248903,    8248934,
              8248952,    8248978,    8249013,    8249049, 8249086,    
8249103,    8249120)

# Window creation

D1 <- unique(data.frame(x=rev(x.coords1),y=rev(y.coords1)))
W1 <- owin(poly=D1)
D2 <- unique(data.frame(x=rev(x.coords2),y=rev(y.coords2)))
W2 <- owin(poly=D2)
W   <- union.owin(W1,W2)
plot(W)

# Create plot points inside the windows
xpoint<-c(517291.6, 517290.3, 517287.3, 517286.9, 517289.5, 517290.4,    
517892.6,
             517892.5,    517892.0,    517893.7,    517893.2, 
517892.1,    517892.2)
ypoint<-c(7970141, 7970139,    7970141,    7970145, 7970147,    
7970146,    7970641,
           7970641, 7970641,    7970643,    7970643,    7970643, 7970643)
syn.ppp<-ppp(x=coordinates(xpoint),y=coordinates(ypoint),window=W)
plot(syn.ppp, main=" ")
# END

Thanks,

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)

         alexandre.santos@cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================


	[[alternative HTML version deleted]]

