<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
<span style="font-size: 10pt; ">Dear All</span><br><div dir="ltr"><div><br></div><div>I would really appreciate some help with a script which a colleague wrote for me, but I am having problems running (and have not been able to contact my colleague).</div><div>The script is designed to compare the area of suitable habitat in binary projections of a large number of species current and future distributions, and create an excel file detailing the total area suitable in the present and future scenario, in addition to the area which overlaps.</div><div>These maps are made from Maxent outputs, which have been reclassified to give binary projections for each species and will be used to assess how climatic change may potentially effect the total area which is suitable for each species. However despite considerable effort I cannot get the attached script to run. It creates an excel with the first species current and future areas of suitable extent detailed, but "NA" for overlap, then gives error warnings (shown below), but changes I have tried have failed to make the script run any further (though I have generated new error messages...)</div><div>Any help or advice would be very greatly appreciated, I can also send some example files of current and future distributions if anyone thinks they maybe able to help.</div><div>Thanks and best wishes</div><div><br></div><div>Alice C Hughes </div><div><br></div><div><br></div><div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">"<span style="font-size:10pt">+ png(paste(out.dir, spp,'_overlap.png', sep=''), width=nrow(spp.asc),height=ncol(spp.asc), units='px', res=300, pointsize=5, bg='white') #start the plot</span></div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">+ par(mar=c(0,0,0,0)+0.1) #remove any plot margins</div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">+ image(overlap, ann=FALSE,axes=FALSE,col=cols, zlim=c(0,3)) #plot the richness</div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">+ legend.gradient(legend.pnts,cols=cols,limits=c(0,3), title='overlap', cex=2)</div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">+ dev.off()</div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">+ }</div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)"><font color="#ff0000" style="line-height:normal">Error in nrow(spp.asc) : object 'bat.asc' not found</font></div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">In addition: There were 12 warnings (use warnings() to see them)</div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">> </div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">> write.csv(overlap.matrix, paste(out.dir, "all_species.csv",sep=''))</div><div style="line-height:17px;color:rgb(42, 42, 42);font-family:'Segoe UI', Tahoma, Verdana, Arial, sans-serif;background-color:rgb(255, 255, 255)">> "</div><br><div><div id="ecxSkyDrivePlaceholder"></div>Date: Fri, 17 Feb 2012 12:18:55 +0100<br>From: Roger.Bivand@nhh.no<br>To: michael.denslow@gmail.com<br>CC: r-sig-geo@r-project.org; skopp@esri.com<br>Subject: Re: [R-sig-Geo] R / ArcGIS<br><br><pre>On Fri, 17 Feb 2012, Michael Denslow wrote:<br> <br>> On Fri, Feb 17, 2012 at 2:15 AM, Edzer Pebesma<br>> <edzer.pebesma@uni-muenster.de> wrote:<br>>> Hi,<br>>><br>>> I might succeed meeting with some ESRI / ArcGIS folks at AAG next week;<br>>> I'm in touch with them in particular about linking R to ArcGIS.<br>>><br>>> I would like to find out if any of you has been successful (or<br>>> unsuccessful) in reading or writing ArcGIS geodatabase files directly<br>>> from R, by using the rgdal driver, or by some other means.<br>><br>> Thanks for working on this. I have had success reading .mdb files but<br>> not .gdb files. I typically use rgdal for this sort of thing. I was<br>> under the impression that GDAL support for .gdb files started with<br>> 1.9.<br> <br>This is off-topic here - the topic here is space-time data.<br> <br>Your question is answered with respect to
ef="http://article.gmane.org/gmane.comp.lang.r.geo/13571" target="_blank">http://article.gmane.org/gmane.comp.lang.r.geo/13571</a><br> <br>(link to gmane because it crosses a month boundary). You did join that <br>thread, and said there that you are using OSX. You were able to open Pgeo <br>files, I think. Are you thinking of the new file geodatabase format? If <br>you read:<br> <br><a href="http://www.gdal.org/ogr/drv_filegdb.html" target="_blank">http://www.gdal.org/ogr/drv_filegdb.html</a><br> <br>you will see that you must build GDAL linking to the proprietary FileGDB <br>API SDK, which requires an ESRI global account. Reading the OGR format <br>file would tell you that, no?<br> <br>Roger<br> <br>><br>> I am very interested in this issue since I am a regular user of ArcGIS<br>> (teaching) and R (research).<br>><br>> Michael<br>><br>><br>>><br>>> Another issue is reading/writing spatio-temporal features (netcdf does<br>>> spatio-temporal grids/arrays pretty well it seems); rumours go that<br>>> ArcGIS might, in the near future, support reading/writing SpatiaLite.<br>>><br>>> With best regards,<br>>> --<br>>> Edzer Pebesma<br>>> Institute for Geoinformatics (ifgi), University of Münster<br>>> Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251<br>>> 8333081, Fax: +49 251 8339763 <a href="http://ifgi.uni-muenster.de" target="_blank">http://ifgi.uni-muenster.de</a><br>>> <a href="http://www.52north.org/geostatistics" target="_blank">http://www.52north.org/geostatistics</a> e.pebesma@wwu.de<br>>><br>>> _______________________________________________<br>>> R-sig-Geo mailing list<br>>> R-sig-Geo@r-project.org<br>>> <a href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-sig-geo</a><br>><br>><br>><br>><br> <br>-- <br>Roger Bivand<br>Department of Economics, NHH Norwegian School of Economics,<br>Helleveien 30, N-5045 Bergen, Norway.<br>voice: +47 55 95 93 55; fax +47 55 95 95 43<br>e-mail: Roger.Bivand@nhh.no<br></pre><br>_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo</div></div> </div> </div></body>
</html>