用下载得到的GIS数据做两张加拿大地图,第一张是用map(,'canada')画的加拿大地图,没有各省边界,第二张是某物种在加拿大的分布情况,第二张是加拿大ecozones map,用到以下语句
x= readShapePoly ('d:/cm/canada.shp')
plot(x)
问题一:与第一张相比,后两张都有些扭曲,应该是坐标系不同的关系,怎么改正让后两幅与第一幅相同?
问题二:图2中的物种分布在不同的ecozone中,如何根据所在的ecozone把seed list中的records进行分类,现在有每个record的经纬度和ecozone的GIS数据