Hi All,

I read a lot about loop formula but this morning I tried to applied the loop
for the following formulas   

I know it's easy for an expertise, but sorry I am working with R in this
mounts and I change completely my GIS philosophy.

 

Thanks for help

 

give a coordinates attribution

 

coordinates(centre_plot1) =~X+Y 

coordinates(centre_plot2) =~X+Y

coordinates(centre_plot3) =~X+Y

coordinates(centre_plot4) =~X+Y

......................

.....................

(I HAVE 500 plot)

 

proj4string(centre_plot1) <- CRS("+init=epsg:26911")

proj4string(centre_plot2) <- CRS("+init=epsg:26911")

proj4string(centre_plot3) <- CRS("+init=epsg:26911")

proj4string(centre_plot4) <- CRS("+init=epsg:26911")

 

writeOGR(centre_plot1 ["Z"], " centre_plot1.shp", " centre_plot1",
driver="ESRI Shapefile")

writeOGR(centre_plot2 ["Z"], " centre_plot2.shp", " centre_plot2",
driver="ESRI Shapefile")

writeOGR(centre_plot3 ["Z"], " centre_plot3.shp", " centre_plot3",
driver="ESRI Shapefile")

writeOGR(centre_plot4 ["Z"], " centre_plot4.shp", " centre_plot4",
driver="ESRI Shapefile")

 

Ale


	[[alternative HTML version deleted]]


