要绘制一张赏心悦目的统计图表，坐标轴的设置至关重要。在R语言底层作图中，对坐标轴的调整主要通过调整plot函数、axis函数和title函数的一系列参数完成。

plot(x,y, ...)

axis(side,at = NULL, labels = TRUE, tick = TRUE, line = NA,

     pos= NA, outer = FALSE, font = NA, lty = "solid",

     lwd = 1, lwd.ticks = lwd, col = NULL,col.ticks = NULL,

     hadj = NA, padj = NA, ...)

title(main= NULL, sub = NULL, xlab = NULL, ylab = NULL,

      line = NA, outer = FALSE, ...)

一、plot函数的准备

　　在个性化设置坐标轴之前中，通常需调整plot函数中的ann、bty、xaxt、yaxt、xaxs和yaxs参数：

　　ann取FALSE时将不会画出标题（包括主、副标题及坐标轴标题）；

　　bty用来设置边框形式，默认值为"o"，表示四面边框都画出，其余可选值包括"l"（左下）、"7"（上右）、"c"（上下左）、"u"（左下右）、"]" （上下右）和"n"（无，即不画边框），在很多个性化绘图中，bty设为"n"，后期的边框线再使用其他函数（如axis）自行添加；

　　xaxs和yaxs 用来设置x轴和y轴的范围，默认值取“r”，表示坐标轴比给定作图范围（参数xlim和ylim给出的范围）稍微大一点儿，取”i”时表示坐标轴范围与给定作图范围完全相同，另外还可取”s”、”e”、”d”；

　　xaxt和yaxt 取”n”时，坐标轴、刻度线以及刻度值将不会画出。

 

x <- seq(-4, 4, 0.01)
y <- x^2
par(mfrow = c(2, 2), mar = c(4, 4, 1, 1))
plot(x, y)   # 未作处理
plot(x, y,  xaxs = "i", yaxs ="i")   # 绘图边框未留白
plot(x, y, bty = 'l')   # 只保留左和下两条边框
plot(x, y, ann = F, bty = "n", xaxt = "n", yaxt ="n")   # 边框、坐标轴都去掉

 



二、axis函数的用法

1、基本操作

　　side表示要操作的坐标轴，取值1、2、3、4分别代表下、左、上、右；

　　at表示刻度线及刻度值所在位置；

　　labels表示刻度值；

　　las表示坐标刻度值文字方向，las=0表示文字方向与坐标轴平行，1表示始终为水平方向，2表示与坐标轴垂直，3表示终为垂直方向。

 

x<- seq(-4, 4, 0.01)
y<- x^2
par(mfrow=c(1, 2), mar = c(4, 4, 1, 1))
#未设置文字方向
plot(x,y, ann = F, xaxt = "n", yaxt = "n")  
axis(1,-4:4, -4:4)
axis(2,seq(0, 16, 4), seq(0, 16, 4))
#文字方向水平
plot(x,y, ann = F, xaxt = "n", yaxt = "n")  
axis(1,-4:4, -4:4, las = 1)
axis(2,seq(0, 16, 4), seq(0, 16, 4), las = 1)

 



2、字体字号

　　cex.axis表示坐标轴刻度值的字号大小，

　　font.axis表示坐标轴刻度值的字体，font=1表示正体，2表示黑体，3表示斜体，4表示黑斜体。

 

x<- seq(-4, 4, 0.01)
y<- x^2
par(mfrow=c(2, 2), mar = c(4, 4, 1, 1))
#未设置字体字号
plot(x,y, ann = F, xaxt = "n", yaxt = "n")
axis(1,-4:4, -4:4)
#cex.axis = 2
plot(x,y, ann = F, xaxt = "n", yaxt = "n")
axis(1,-4:4, -4:4, cex.axis = 2)
#font.axis = 2
plot(x,y, ann = F, xaxt = "n", yaxt = "n")
axis(1,-4:4, -4:4, font.axis = 2)
#font.axis = 3
plot(x,y, ann = F, xaxt = "n", yaxt = "n")
axis(1,-4:4, -4:4, font.axis = 3)






 

 

 

 

3、颜色

　　col表示图的颜色，用在axis函数中表示坐标轴线和坐标刻度线的颜色；

　　col.axis表示坐标轴刻度值的颜色；

　　col.ticks表示坐标轴刻度线的颜色。

 

x<- seq(-4, 4, 0.01)
y<- x^2
par(mfrow=c(2, 2), mar = c(4, 4, 1, 1))
#未设置颜色
plot(x,y, ann = F, bty = "n", xaxt = "n", yaxt = "n")
axis(1,-4:4, -4:4)
#col = 2
plot(x,y, ann = F, bty = "n", xaxt = "n", yaxt = "n")
axis(1,-4:4, -4:4, col = 2)
#col.axis = 2
plot(x,y, ann = F, bty = "n", xaxt = "n", yaxt = "n")
axis(1,-4:4, -4:4, col.axis = 2)
#col.ticks = 2
plot(x,y, ann = F, bty = "n", xaxt = "n", yaxt = "n")
axis(1,-4:4, -4:4, col.ticks = 2)



 

4、位置

　　line表示坐标轴线位置与图像边框的距离，取负数时会画在图像边框以内；

　　mgp默认值为c(3, 1, 0)，三个数字分别代表坐标轴标题、刻度值和轴线与绘图边框的距离；

　　tcl默认值为-0.5，数值表示刻度线长度，负值表示刻度线朝外，正值朝里；

　　pos 表示轴线所在的位置；

　　line.outer取TRUE时，坐标轴将画在画布边缘处；

　　hadj指将刻度值沿平行坐标轴方向调整的距离；

　　padj指将刻度值沿垂直坐标轴方向调整的距离。

 

x <- seq(-4, 4, 0.01)
y <- x^2
par(mfrow= c(2, 2), mar = c(4, 4, 1, 1))
# 未设置刻度值位置
plot(x, y, ann = F, xaxt = "n", yaxt ="n")
axis(1, -4:4, -4:4)
# 使用line调整刻度值位置
plot(x, y, ann = F, xaxt = "n", yaxt ="n")
axis(1, -4:4, -4:4, line = 2)
# 使用mgp调整刻度值位置
plot(x, y, ann = F, xaxt = "n", yaxt ="n")
axis(1, -4:4, -4:4, mgp = c(3, 2, 0)) 
# 使用padj调整刻度值位置
plot(x, y, ann = F, xaxt = "n", yaxt ="n")
axis(1, -4:4, -4:4, padj = 1)



 

 

 

5、其他

　　tick取FALSE时，坐标轴线和刻度线不画出；

　　lty表示线型，用在axis函数中表示坐标轴线型；

　　lwd表示线的粗细，用在axis函数中表示坐标轴线粗细；

　　lwd.ticks表示刻度线粗细。

 

三、title函数的用法

　　main、sub、xlab和ylab分别表示主标题、副标题、x轴标题和y轴标题；

　　cex.lab表示坐标轴标题的字号大小；

　　font.lab表示坐标轴标题的字体；

　　col.lab表示坐标轴标题的颜色；

　　其余参数和axis中用法一致。

 

x<- seq(-4, 4, 0.01)
y<- x^2
par(mfrow=c(2, 2), mar = c(4, 4, 1, 1))
#无标题
plot(x,y, ann = F, xaxt = "n", yaxt = "n")  
axis(1,-4:4, -4:4)
axis(2,seq(0, 16, 4), seq(0, 16, 4))
#使用title设置标题
plot(x,y, ann = F, xaxt = "n", yaxt = "n")  
axis(1,-4:4, -4:4, las = 1)
axis(2,seq(0, 16, 4), seq(0, 16, 4), las = 1)
title(xlab= 'x', ylab = 'y')
#使用line调整位置
plot(x,y, ann = F, xaxt = "n", yaxt = "n")  
axis(1,-4:4, -4:4, las = 1)
axis(2,seq(0, 16, 4), seq(0, 16, 4), las = 1)
title(xlab= 'x', ylab = 'y', line = 2)
#调整颜色字体字号
plot(x,y, ann = F, xaxt = "n", yaxt = "n")  
axis(1,-4:4, -4:4, las = 1)
axis(2,seq(0, 16, 4), seq(0, 16, 4), las = 1)
title(xlab= 'x', ylab = 'y', col.lab = 2, font.lab = 4, cex.lab = 2)




四、刻度间隔

　　plot函数会自动给出比较合理的刻度间隔，本质上是调用了pretty函数。

　　pretty(x, n = 5, min.n = n %/% 3,  shrink.sml = 0.75,

       　　high.u.bias = 1.5, u5.bias = .5 + 1.5*high.u.bias,

       　　eps.correct = 0, ...)

　　x是一个序列，pretty函数会先对x取一个range，然后等分为大约n+1个间隔。如果在作图时对系统自动给出的间隔不满意，可以调整n的值，用pretty函数自定义合适的间隔。

 

 

 

bd09经纬度坐标。允许的值为bd09ll、bd09mc、gcj02、wgs84。bd09ll表示百度经纬度坐标，bd09mc表示百度墨卡托坐标，gcj02表示经过国测局加密的坐标，wgs84表示gps获取的坐标

 如何将下载的影像变换为西安80坐标？

                                              同步视频教程：投影转换(转西安80) 

 如何将下载的影像变换为北京54坐标？

 如何将下载的影像变换为CGCS2000坐标？

 

下载：全国路网数据、全国水系矢量 行政区划边界（省市区县）、行政地名 矢量数据  

 

                利用ARCGIS进行自定义坐标系和投影转换

 

                ARCGIS种通过三参数和其参数进行精确投影转换

     注意：投影转换成80坐标系需要下载无偏移卫星图像进行转换，有偏移的转换将导致转换后的卫星图像扭曲，坐标错误，无法配准。

     第一步：选择无偏移地图源，下载你所需要的卫星图像。

     第二步：选择BIGEMAP软件右边工具栏，选择【投影转换】，如下图所示：



      2.1   选择说明：

               1. 源文件：选择下载好的卫星图像文件（下载目录中后缀为tiff的文件）

               2. 源坐标系：打开的源文件的投影坐标系（自动读取，不需要手动填写）

               3. 输出文件：选择转换后你要保持文件的文件路径和文件名

               4. 目标坐标系：选择你要转换成的目标坐标系，如下图：



     选择上图的更多，如下图所示：



          1：选择 -Xi'an 1980  2：选择地区   3：选择分度带对应的带号（一般默认，也可以手动修改）

          选择对应的分度带或者中央子午线（请参看：如何选择分度带？），点击【确定】

      5. 重采样算法：投影转换需要将影像的像素重新排列，一次每种算法的效率不一样，一般选择【立方卷积采样】，以达到最好的效果。如下图：



     6. 指定变换参数：在不知道的情况下，可以不用填此处信息，如果√上，则如下图：



        此参数为【三参数】或者【七参数】，均为国家保密参数，需要到当地的测绘部门或者国土部门，以单位名义签保密协议进行购买，此参数各地都不一样，是严格保密的，请不要随便流通。

   

    第三步：点击【确定】，开始转换，如下图：



    第四步：完成后，打开你刚才选择的输出文件夹，里面就是转换后的卫星图像。

    第五步：如果你需要套合你手里已经有的矢量文件，请参看：【BIGEMAP无偏移影像叠加配准】

用过CASS的人都知道，野外数字测图得到的点数据（平面坐标）可以直接导入到CASS中，进一步绘制地形图。那么，带有坐标的数据能不能在ArcGIS中实现点图层的生成呢？答案是必须的！

本文以气象台站shp数据的生成为例，详细介绍ArcGIS 10.2中导入X、Y坐标（这里指的是经度、纬度），生成Shapefile点数据的流程。


相关阅读：《全站仪、RTK测量坐标数据在CASS和ArcGIS中展点的区别和联系》

目录

1、数据准备

2、添加X、Y数据

3、属性表连接（Join）

4、导出为Shapefile矢量数据

5、自定义投影转换

1、数据准备


链接：https://pan.baidu.com/s/1eaulWklowfB-Qd5S-7wnkA

提取码：4ugl




2、添加X、Y数据

（1）【文件】→【添加数据】→【添加XY数据】




（2）选择“气象台站.xlsx”中的sheet“气象台站”，X字段为经度，Y字段为纬度，Z字段为高程（此处没有），选择地理坐标系CGCS2000。

    注意：坐标系只能选择地理坐标系，而不能选择投影坐标系！！！




    点击【确定】。



   （3）结果展示



3、属性表连接（Join）

    属性表连接为可选项，可以不用连接，本文以连接各个气象台站的气温和降水数据为例。

（1）数据准备：气温降水.xlsx（点击下载）



（2）属性表连接





（3）结果展示



4、导出为Shapefile矢量数据



5、自定义投影转换

为了后面的计算面积等操作方便，现进行投影转换，本文自定义一个Albers投影，并进行投影转换。

注意：Albers投影的参数为中央经线105，标准纬线一25，标准纬线二47。



    结果展示：




刘一哥GIS：专注GIS教育，探索GIS奥秘，分享GIS价值 ！