Dear All
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).
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.
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...)
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.
Thanks and best wishes
Alice C Hughes
"+ 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
+ par(mar=c(0,0,0,0)+0.1) #remove any plot margins
+ image(overlap, ann=FALSE,axes=FALSE,col=cols, zlim=c(0,3)) #plot the richness
+ legend.gradient(legend.pnts,cols=cols,limits=c(0,3), title='overlap',
cex=2)
+ dev.off()
+ }
Error in nrow(spp.asc) : object 'bat.asc' not found
In addition: There were 12 warnings (use warnings() to see
them)>
> write.csv(overlap.matrix, paste(out.dir,
"all_species.csv",sep=''))
> "
Date: Fri, 17 Feb 2012 12:18:55 +0100
From: Roger.Bivand at nhh.no
To: michael.denslow at gmail.com
CC: r-sig-geo at r-project.org; skopp at esri.com
Subject: Re: [R-sig-Geo] R / ArcGIS
On Fri, 17 Feb 2012, Michael Denslow wrote:
> On Fri, Feb 17, 2012 at 2:15 AM, Edzer Pebesma
> <edzer.pebesma at uni-muenster.de> wrote:
>> Hi,
>>
>> I might succeed meeting with some ESRI / ArcGIS folks at AAG next week;
>> I'm in touch with them in particular about linking R to ArcGIS.
>>
>> I would like to find out if any of you has been successful (or
>> unsuccessful) in reading or writing ArcGIS geodatabase files directly
>> from R, by using the rgdal driver, or by some other means.
>
> Thanks for working on this. I have had success reading .mdb files but
> not .gdb files. I typically use rgdal for this sort of thing. I was
> under the impression that GDAL support for .gdb files started with
> 1.9.
This is off-topic here - the topic here is space-time data.
Your question is answered with respect to Personal geodatabases on the
Windows platform in thread:
http://article.gmane.org/gmane.comp.lang.r.geo/13571
(link to gmane because it crosses a month boundary). You did join that
thread, and said there that you are using OSX. You were able to open Pgeo
files, I think. Are you thinking of the new file geodatabase format? If
you read:
http://www.gdal.org/ogr/drv_filegdb.html
you will see that you must build GDAL linking to the proprietary FileGDB
API SDK, which requires an ESRI global account. Reading the OGR format
file would tell you that, no?
Roger
>
> I am very interested in this issue since I am a regular user of ArcGIS
> (teaching) and R (research).
>
> Michael
>
>
>>
>> Another issue is reading/writing spatio-temporal features (netcdf does
>> spatio-temporal grids/arrays pretty well it seems); rumours go that
>> ArcGIS might, in the near future, support reading/writing SpatiaLite.
>>
>> With best regards,
>> --
>> Edzer Pebesma
>> Institute for Geoinformatics (ifgi), University of M?nster
>> Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251
>> 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de
>> http://www.52north.org/geostatistics e.pebesma at wwu.de
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>
>
>
--
Roger Bivand
Department of Economics, NHH Norwegian School of Economics,
Helleveien 30, N-5045 Bergen, Norway.
voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no
_______________________________________________ R-sig-Geo mailing list R-sig-Geo
at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
