R语言怎么画出固定的经纬度区域图

Python015

R语言怎么画出固定的经纬度区域图,第1张

R语言怎么画出固定的经纬度区域图：

1 如果有经纬度坐标，直接用geoda和arcgis就可以生成图了。

2 R处理数据。

3 x <- readShapePoly('bou2_4p.shp')，这个变量名x换个别的，

ggplot(china.map,aes(x=long,y=lat)) 这里x=lon

R语言做不到这个。。。

你可以用R处理经纬度数据，R本身没有（除非是某些包自带的）

你可以用百度地图的api

http://api.map.baidu.com/lbsapi/getpoint/index.html

tgα= 杆长/影长。

在地球上任何地点，只要有只表，有根竹竿，一根卷尺，就可知道当地经纬度。但表必须与该国标准时校对。

方法如下：

1、先算两分日

比如在中国某地，杆影最短时是中午13点20分，且杆长与影长之比为1，则可知该地是北纬45°（tgα＝1)，东经100°（从120°里1小时减15°，4分钟减1°）杆长与影长之比需查表求α，这里用了特殊角。

2、再算两至日经度的算法不变纬度在北半球冬至α＋23.5°，夏至α－23.5°在任意一天加减修正值即可。

3、修正值算法：就是距两分或两至日的天数差乘以94/365.比如2013年2月17日，2013年3月22日春分差33天，即太阳直射点在南纬。

33×94/365=8.5°

所以今天正午时得到的纬度是（arctgα＋8.5）°

在地球上任何地点，只要有只表，有根竹竿，一根卷尺，就可知道当地经纬度。但表必须与该国标准时校对。

1.已知地球上两点的经度、纬度：（X1,Y1), (X2,Y2)，其中X1,X2为经度，Y1,Y2为纬度。

视计算程序需要转化为弧度（＊3.1415926/180）。

地球半径为R=6371.0 km。

则两点距离d=R*arcos。

2.在地球上同一条经线上纬度相差一度，实际距离是111公里；在赤道经度相差一度，实际距离是111公里，在其他的纬线上，经度每相差一度，实际距离是111*经度数＊cos经度数，例如在北纬30度，经度相差5度，实际距离是：111*5*cos30度。