CN112102430B - 单要素变化曲线至地理地图的投影计算方法 - Google Patents
单要素变化曲线至地理地图的投影计算方法 Download PDFInfo
- Publication number
- CN112102430B CN112102430B CN202010831650.6A CN202010831650A CN112102430B CN 112102430 B CN112102430 B CN 112102430B CN 202010831650 A CN202010831650 A CN 202010831650A CN 112102430 B CN112102430 B CN 112102430B
- Authority
- CN
- China
- Prior art keywords
- scale
- axis
- size
- coordinates
- ratio
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 71
- 230000008859 change Effects 0.000 title claims abstract description 44
- 238000002372 labelling Methods 0.000 claims description 23
- 238000000034 method Methods 0.000 claims description 17
- 240000004282 Grewia occidentalis Species 0.000 claims description 7
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明公开了一种单要素变化曲线至地理地图的投影计算方法,包括以下步骤:获取待投影要素曲线的初始计算参数;获取绘图区域的基本信息率;根据初始计算参数和基本信息计算待投影要素曲线坐标轴的左下角坐标;计算X轴和Y轴上的各刻度间距的比率;计算待投影要素中各要素值的坐标;计算绘图边框的四角坐标;计算要素变化曲线投影后的X轴刻度坐标和Y轴刻度坐标;计算X轴和Y轴的各刻度值标注坐标和单位标注坐标;输出计算结果。本发明步骤清晰、计算简便、易于实现,输出结果随用户比例尺变化而自动变化。
Description
技术领域
本发明涉及投影变换技术领域,更具体的说是涉及一种单要素变化曲线至地理地图的投影计算方法。
背景技术
传统图集能够充分展示某个要素的时空变化状况。基于经纬度方区统计的要素变化曲线,展示的是要素随时间的变化情况,是图集的一种重要类型。
目前,在传统的图集制作过程中,该要素变化曲线是采用手工叠加的方式进行的,费时、费力。GIS软件所进行的专题图绘制,虽然能够实现曲线的自动绘制,但其曲线与地理地图坐标系是脱离的,在地图比例尺切换时,图形的位置会发生变化,在进行图形操作和发布时存在诸多障碍。
因此,如何提供一种能够在地图比例尺切换时不存在位置偏移的单要素变化曲线至地理地图的投影计算方法,是本领域技术人员亟需解决的技术问题。
发明内容
有鉴于此,本发明提供了一种单要素变化曲线至地理地图的投影计算方法及系统,能够解决单要素变化曲线至地理地图的投影计算问题,将基于经纬度方区统计的要素变化曲线投影至地理地图,以解决地图比例尺切换存在的问题和地图后期发布问题。
为了实现上述目的,本发明采用如下技术方案:
一种单要素变化曲线至地理地图的投影计算方法,包括以下步骤:
步骤S1、获取待投影要素曲线的初始计算参数;
步骤S2、获取绘图区域的基本信息,所述绘图区域位于待投影垂向剖面所处的经纬度方区内,且相对于所述经纬度方区的边缘预留有四个方向的空白区;所述基本信息包括四个方向的空白区相对于方区大小的比率,X轴刻度相对于方区大小的比率,Y轴刻度相对于方区大小的比率,X轴标注相对于方区大小的比率,以及Y轴标注相对于方区大小的比率;
步骤S3、根据步骤S1中的初始计算参数和步骤S2中的基本信息计算待投影要素曲线坐标轴的左下角坐标;
步骤S4、计算X轴和Y轴上的各刻度间距的比率;
步骤S5、计算待投影要素中各要素值的坐标;
步骤S6、计算绘图边框的四角坐标;
步骤S7、计算要素变化曲线投影后的X轴刻度坐标和Y轴刻度坐标;
步骤S8、计算X轴和Y轴的各刻度值标注坐标和单位标注坐标;
步骤S9、对步骤S5-S8的计算结果进行分类和整理,输出计算结果。
经由上述的技术方案可知,与现有技术相比,本发明对单要素随时间的变化情况,不再需要手工绘制,仅需输入待投影的单要素曲线的计算参数以及要素所处的经纬度方区的经纬度信息,便能直接绘制到地理地图的经纬网格中,能够直观的显示该要素随时间的变化信息,能够将基于经纬度方区统计的要素变化值,投影至地理地图,进一步直接绘制到地理地图的统计网格中,可广泛应用于GIS专题图的制作。通过对单要素变化曲线投影至地理地图进行计算,将基于经纬度方区统计的要素变化曲线投影至地理地图,以解决地图比例尺切换存在的问题和地图后期发布问题。
优选的,在上述一种单要素变化曲线至地理地图的投影计算方法中,步骤S1中的初始计算参数包括经纬度方区的左下角经度lon和纬度lat、经纬度方区大小gsize、要素的数据值序列l,要素的时间序列z,X轴标注Xlab,X轴单位xStr,Y轴单位yStr,其中z值的个数设为nz,l的个数设为nv。
优选的,在上述一种单要素变化曲线至地理地图的投影计算方法中,步骤S2中,四个方向的空白区分别为左空白区、右空白区、上空白区和下空白区;左空白区相对于经纬度方区大小的比率lScale为0.10;右空白区对于经纬度方区大小的比率rScale为0.10;下空白区相对于经纬度方区大小的比率bScale为0.10;上空白区相对于经纬度方区大小的比率tScale为0.13;X轴刻度相对于方区大小的比率xtickScale和Y轴刻度相对于方区大小的比率ytickScale均为0.02;X轴标注相对于方区大小的比率xlabScale和Y轴标注相对于方区大小的比率ylabScale均为0.05。本发明采用上述数值,使在经纬度方区内绘图时更加合理和美观。
优选的,在上述一种单要素变化曲线至地理地图的投影计算方法中,步骤S3的计算公式如下:
x0=lon+(lScale+xtickScale+xlabScale)*gsize;
y0=lat+(bScale+ytickScale+ylabScale)*gsize;
(x0,y0)表示待投影要素曲线坐标轴的左下角坐标。
优选的,在上述一种单要素变化曲线至地理地图的投影计算方法中,步骤S4中X轴刻度间距比率的计算公式如下;
monScale=(1-lScale-rScale-ytickScale-xlabScale)*gsize/(nz-1);
Y轴刻度间距比率的计算公式如下;
zScale=(1-tScale-bScale-xtickScale-ylabScale)*gsize/(nv-1)。
优选的,在上述一种单要素变化曲线至地理地图的投影计算方法中,步骤S5中各要素值坐标的计算公式如下:
xn(k)=x0+(k-1)*monScale;
yn(k)=y0+(l(k)-l(1))*zScale;
其中,k∈(1,nz);(xn(k),yn(k))表示待投影要素第k个要素值坐标;l(k)表示Y轴的第k个刻度值;l(1)表示轴的第一个刻度。
优选的,在上述一种单要素变化曲线至地理地图的投影计算方法中,步骤S6中绘图边框坐标的计算公式如下:
xmin=xn(1);
ymin=y0;
xmax=xn(nz);
ymax=y0+(l(end)-l(1))*zScale;
上式中,通过上式得出绘图边框的左上角坐标(xmin,ymax),左下角坐标(xmin,ymin),右下角坐标(xmax,ymin)和右上角坐标(xmax,ymax)。
优选的,在上述一种单要素变化曲线至地理地图的投影计算方法中,步骤S7中X轴刻度的计算公式如下:
xt=x0+(n-1)*monScale;
yt1=y0;
yt2=y0-xtickScale*gsize;
第n个X轴刻度短线的起点坐标为(xt,yt1),终点坐标为(xt,yt2);
其中,n∈(1,nv);
Y轴刻度的计算公式如下:
xt1=x0;
xt2=x0-ytickScale*gsize;
yt=y0+(l(k)-l(1))*zScale;
第k个Y轴刻度短线的起点坐标为(xt1,yt),终点坐标为(xt2,yt);
其中,k∈(1,nz)。
优选的,在上述一种单要素变化曲线至地理地图的投影计算方法中,步骤S8中X轴刻度值标注坐标的计算公式如下:
xb=x0+(k-1)*monScale;
yb=y0-xtickScale*gsize;
X轴刻度值标注坐标为(xb,yb),标注值为Xlab(k);
其中,k∈(1,nz);
X轴单位标注坐标的计算公式如下:
xbn=x0+(nz-1)*monScale;
ybn=y0-5*xtickScale*gsize;
X轴单位标注的坐标为(xbn,ybn),标注值为xStr;
Y轴的刻度值标注坐标的计算公式如下:
xc=x0-ytickScale*gsize;
yc=y0+(l(n)-l(1))*zScale;
Y轴的刻度值标注的坐标为(xc,yc),标注值为l(n);
其中,n∈(1,nv);
Y轴单位标注坐标的计算公式如下:
xcn=x0(n);
ycn=y0(n)+(l(end)-l(1))*zScale+4*ytickScale*gsize;
Y轴单位标注的坐标为(xcn,ycn),标注值为yStr。
优选的,在上述一种单要素变化曲线至地理地图的投影计算方法中,步骤S7的输出结果包括要素曲线、绘图边框线、刻度线和标注点;其中要素曲线为依次连接步骤S5中各要素值坐标形成的线文件;绘图边框线为连接步骤S6得出的绘图边框四角坐标形成的线文件;刻度线为循环连接步骤S7计算得出的X轴和Y轴刻度的起点和终端形成的线文件;标注点为步骤S8计算得出的X轴和Y轴刻度标注和单位标注的点文件。本发明的输出成果清晰,包括绘图边框线、要素曲线、X轴和Y轴刻度线、X轴标注点和Y轴标注点共4个文件,易于成果在GIS软件中的符号化、渲染和使用;其要素曲线、边框、X轴和Y轴刻度线随用户比例尺变化而自动变化,无需进行再次编辑与调整。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1附图为本发明提供的垂向剖面的要素值至地理地图的投影计算方法的流程图;
图2附图为本发明提供的计算结果输出效果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明实施例公开了一种单要素变化曲线至地理地图的投影计算方法,包括以下步骤:
步骤S1、获取待投影要素曲线的初始计算参数;
步骤S2、获取绘图区域的基本信息,所述绘图区域位于待投影垂向剖面所处的经纬度方区内,且相对于所述经纬度方区的边缘预留有四个方向的空白区;所述基本信息包括四个方向的空白区相对于方区大小的比率,X轴刻度相对于方区大小的比率,Y轴刻度相对于方区大小的比率,X轴标注相对于方区大小的比率,以及Y轴标注相对于方区大小的比率;
步骤S3、根据步骤S1中的初始计算参数和步骤S2中的基本信息计算待投影要素曲线坐标轴的左下角坐标;
步骤S4、计算X轴和Y轴上的各刻度间距的比率;
步骤S5、计算待投影要素中各要素值的坐标;
步骤S6、计算绘图边框的四角坐标;
步骤S7、计算要素变化曲线投影后的X轴刻度坐标和Y轴刻度坐标;
步骤S8、计算X轴和Y轴的各刻度值标注坐标和单位标注坐标;
步骤S9、对步骤S5-S8的计算结果进行分类和整理,输出计算结果。
本发明对单要素随时间的变化情况,不再需要手工绘制,仅需输入待投影的单要素曲线的计算参数以及要素所处的经纬度方区的经纬度信息,便能直接绘制到地理地图的经纬网格中,能够直观的显示该要素随时间的变化信息,能够将基于经纬度方区统计的要素变化值,投影至地理地图,进一步直接绘制到地理地图的统计网格中,可广泛应用于GIS专题图的制作。通过对单要素变化曲线投影至地理地图进行计算,将基于经纬度方区统计的要素变化曲线投影至地理地图,以解决地图比例尺切换存在的问题和地图后期发布问题。
具体的,步骤S1、获取计算参数:
获取步骤S1中的初始计算参数包括经纬度方区的左下角经度lon和纬度lat、经纬度方区大小gsize、要素的数据值序列l,要素的时间序列z,X轴标注Xlab,X轴单位xStr,Y轴单位yStr,其中z值的个数设为nz,l的个数设为nv。
步骤S2、获取绘图区域的基本信息:
设定四个方向的空白区分别为左空白区、右空白区、上空白区和下空白区;分别设定四个空白区相对于经纬度方区大小的比率lScale、rScale、bScale、tScale;设定X轴刻度和Y轴刻度相对于方区大小的比率相对于方区大小的比率xScale、yScale;设定X轴标注和Y轴标注相对于方区大小的比率xlabScale、ylabScale。经反复计算,为保证经纬度方区绘图的合理性和美观性,将各参数的数值设定如下:
lScale=0.10;
rScale=0.10;
tScale=0.13;
bScale=0.10;
xScale=0.02;
yScale=0.02;
xlabScale=0.05;
ylabScale=0.05。
步骤S3、计算待投影要素曲线坐标轴的左下角坐标:
x0=lon+(lScale+xtickScale+xlabScale)*gsize;
y0=lat+(bScale+ytickScale+ylabScale)*gsize;
(x0,y0)表示待投影要素曲线坐标轴的左下角坐标。
步骤S4、计算X轴和Y轴上的各刻度间距的比率:
X轴刻度间距比率的计算公式如下;
monScale=(1-lScale-rScale-ytickScale-xlabScale)*gsize/(nz-1);
Y轴刻度间距比率的计算公式如下;
zScale=(1-tScale-bScale-xtickScale-ylabScale)*gsize/(nv-1)。
步骤S5、计算待投影要素中各要素值的坐标;
xn(k)=x0+(k-1)*monScale;
yn(k)=y0+(l(k)-l(1))*zScale;
其中,k∈(1,nz);(xn(k),yn(k))表示待投影要素第k个要素值坐标;l(k)表示Y轴的第k个刻度值;l(1)表示轴的第一个刻度。
步骤S6、计算绘图边框的四角坐标;
xmin=xn(1);
ymin=y0;
xmax=xn(nz);
ymax=y0+(l(end)-l(1))*zScale;
上式中,通过上式得出绘图边框的左上角坐标(xmin,ymax),左下角坐标(xmin,ymin),右下角坐标(xmax,ymin)和右上角坐标(xmax,ymax)。
步骤S7、计算要素变化曲线投影后的X轴刻度坐标和Y轴刻度坐标;
X轴刻度的计算公式如下:
xt=x0+(n-1)*monScale;
yt1=y0;
yt2=y0-xtickScale*gsize;
第n个X轴刻度短线的起点坐标为(xt,yt1),终点坐标为(xt,yt2);
其中,n∈(1,nv);
Y轴刻度的计算公式如下:
xt1=x0;
xt2=x0-ytickScale*gsize;
yt=y0+(l(k)-l(1))*zScale;
第k个Y轴刻度短线的起点坐标为(xt1,yt),终点坐标为(xt2,yt);
其中,k∈(1,nz)。
步骤S8、计算X轴和Y轴的各刻度值标注坐标和单位标注坐标;
X轴刻度值标注坐标的计算公式如下:
xb=x0+(k-1)*monScale;
yb=y0-xtickScale*gsize;
X轴刻度值标注坐标为(xb,yb),标注值为Xlab(k);
其中,k∈(1,nz);
X轴单位标注坐标的计算公式如下:
xbn=x0+(nz-1)*monScale;
ybn=y0-5*xtickScale*gsize;
X轴单位标注的坐标为(xbn,ybn),标注值为xStr;
Y轴的刻度值标注坐标的计算公式如下:
xc=x0-ytickScale*gsize;
yc=y0+(l(n)-l(1))*zScale;
Y轴的刻度值标注的坐标为(xc,yc),标注值为l(n);
其中,n∈(1,nv);
Y轴单位标注坐标的计算公式如下:
xcn=x0(n);
ycn=y0(n)+(l(end)-l(1))*zScale+4*ytickScale*gsize;
Y轴单位标注的坐标为(xcn,ycn),标注值为yStr。
步骤S9、对步骤S5-S8的计算结果进行分类和整理,输出计算结果。
步骤S7的输出结果包括要素曲线、绘图边框线、刻度线和标注点;其中要素曲线为依次连接步骤S5中各要素值坐标形成的线文件;绘图边框线为连接步骤S6得出的绘图边框四角坐标形成的线文件;刻度线为循环连接步骤S7计算得出的X轴和Y轴刻度的起点和终端形成的线文件;标注点为步骤S8计算得出的X轴和Y轴刻度标注和单位标注的点文件。
步骤S9实现了对步骤S5-S8的计算结果的分类与整理,使计算结果更加有序,易于进行下一步的图形绘制。本发明的输出成果清晰,包括绘图边框线、要素曲线、X轴和Y轴刻度线、X轴标注点和Y轴标注点共4个文件,易于成果在GIS软件中的符号化、渲染和使用;其要素曲线、边框、X轴和Y轴刻度线随用户比例尺变化而自动变化,无需进行再次编辑与调整
如图2所示,本发明输出成果清晰,包括绘图边框线文件、要素曲线、刻度线文件和标注点文件四个文件,易于成果在GIS软件中的符号化、渲染和使用。
图2的四角为经纬度坐标,框内为温度随时间变化的曲线,X轴刻度为时间,单位月份,Y轴刻度为温度,单位℃。通过该图可以直观判断该经纬度方区内温度随时间的变化情况。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (2)
1.一种单要素变化曲线至地理地图的投影计算方法,其特征在于,包括以下步骤:
步骤S1、获取待投影要素曲线的初始计算参数;初始计算参数包括经纬度方区的左下角经度lon和纬度lat、经纬度方区大小gsize、要素的数据值序列l,要素的时间序列z,X轴标注Xlab,X轴单位xStr,Y轴单位yStr,其中z值的个数设为nz,l的个数设为nv;
步骤S2、获取绘图区域的基本信息,所述绘图区域位于待投影垂向剖面所处的经纬度方区内,且相对于所述经纬度方区的边缘预留有四个方向的空白区,分别为左空白区、右空白区、上空白区和下空白区;所述基本信息包括四个方向的空白区相对于方区大小的比率lScale、rScale、bScale、tScale,X轴刻度相对于方区大小的比率xtickScale,Y轴刻度相对于方区大小的比率ytickScale,X轴标注相对于方区大小的比率xlabScale,以及Y轴标注相对于方区大小的比率ylabScale;
步骤S3、根据步骤S1中的初始计算参数和步骤S2中的基本信息计算待投影要素曲线坐标轴的左下角坐标;计算公式为:
x0=lon+(lScale+xtickScale+xlabScale)*gsize;
y0=lat+(bScale+ytickScale+ylabScale)*gsize;
(x0,y0)表示待投影要素曲线坐标轴的左下角坐标;
步骤S4、计算X轴和Y轴上的各刻度间距的比率,X轴刻度间距比率的计算公式如下:
monScale=(1-lScale-rScale-ytickScale-xlabScale)*gsize/(nz-1);
Y轴刻度间距比率的计算公式如下:
zScale=(1-tScale-bScale-xtickScale-ylabScale)*gsize/(nv-1);
步骤S5、计算待投影要素中各要素值的坐标;各要素值坐标的计算公式如下:
xn(k)=x0+(k-1)*monScale;
yn(k)=y0+(l(k)-l(1))*zScale;
其中,k∈(1,nz);(xn(k),yn(k))表示待投影要素第k个要素值坐标;l(k)表示Y轴的第k个刻度值;l(1)表示轴的第一个刻度;
步骤S6、计算绘图边框的四角坐标;绘图边框坐标的计算公式如下:
xmin=xn(1);
ymin=y0;
xmax=xn(nz);
ymax=y0+(l(end)-l(1))*zScale;
上式中,通过上式得出绘图边框的左上角坐标(xmin,ymax),左下角坐标(xmin,ymin),右下角坐标(xmax,ymin)和右上角坐标(xmax,ymax);
步骤S7、计算要素变化曲线投影后的X轴刻度坐标和Y轴刻度坐标;
X轴刻度的计算公式如下:
xt=x0+(n-1)*monScale;
yt1=y0;
yt2=y0-xtickScale*gsize;
第n个X轴刻度短线的起点坐标为(xt,yt1),终点坐标为(xt,yt2);
其中,n∈(1,nv);
Y轴刻度的计算公式如下:
xt1=x0;
xt2=x0-ytickScale*gsize;
yt=y0+(l(k)-l(1))*zScale;
第k个Y轴刻度短线的起点坐标为(xt1,yt),终点坐标为(xt2,yt);
其中,k∈(1,nz);
步骤S8、计算X轴和Y轴的各刻度值标注坐标和单位标注坐标;X轴刻度值标注坐标的计算公式如下:
xb=x0+(k-1)*monScale;
yb=y0-xtickScale*gsize;
X轴刻度值标注坐标为(xb,yb),标注值为Xlab(k);
其中,k∈(1,nz);
X轴单位标注坐标的计算公式如下:
xbn=x0+(nz-1)*monScale;
ybn=y0-5*xtickScale*gsize;
X轴单位标注的坐标为(xbn,ybn),标注值为xStr;
Y轴的刻度值标注坐标的计算公式如下:
xc=x0-ytickScale*gsize;
yc=y0+(l(n)-l(1))*zScale;
Y轴的刻度值标注的坐标为(xc,yc),标注值为l(n);
其中,n∈(1,nv);
Y轴单位标注坐标的计算公式如下:
xcn=x0(n);
ycn=y0(n)+(l(end)-l(1))*zScale+4*ytickScale*gsize;
Y轴单位标注的坐标为(xcn,ycn),标注值为yStr;
步骤S9、对步骤S5-S8的计算结果进行分类和整理,输出计算结果;输出结果包括要素曲线、绘图边框线、刻度线和标注点;其中要素曲线为依次连接步骤S5中各要素值坐标形成的线文件;绘图边框线为连接步骤S6得出的绘图边框四角坐标形成的线文件;刻度线为循环连接步骤S7计算得出的X轴和Y轴刻度的起点和终端形成的线文件;标注点为步骤S8计算得出的X轴和Y轴刻度标注和单位标注的点文件。
2.根据权利要求1所述的一种单要素变化曲线至地理地图的投影计算方法,其特征在于,步骤S2中,左空白区相对于经纬度方区大小的比率lScale为0.10;右空白区对于经纬度方区大小的比率rScale为0.10;下空白区相对于经纬度方区大小的比率bScale为0.10;上空白区相对于经纬度方区大小的比率tScale为0.13;X轴刻度相对于方区大小的比率xtickScale和Y轴刻度相对于方区大小的比率ytickScale均为0.02;X轴标注相对于方区大小的比率xlabScale和Y轴标注相对于方区大小的比率ylabScale均为0.05。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010831650.6A CN112102430B (zh) | 2020-08-18 | 2020-08-18 | 单要素变化曲线至地理地图的投影计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010831650.6A CN112102430B (zh) | 2020-08-18 | 2020-08-18 | 单要素变化曲线至地理地图的投影计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112102430A CN112102430A (zh) | 2020-12-18 |
CN112102430B true CN112102430B (zh) | 2023-11-21 |
Family
ID=73754539
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010831650.6A Active CN112102430B (zh) | 2020-08-18 | 2020-08-18 | 单要素变化曲线至地理地图的投影计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112102430B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB454654A (en) * | 1935-01-03 | 1936-10-05 | Arthur Thomas Keeble | A great circle graphic calculator for use in solving problems in navigation and astronomy |
JP2002366027A (ja) * | 2001-06-05 | 2002-12-20 | Sokken:Kk | 経度間隔の違う地図の、経度目盛りの設定方法と既知経緯度座標を地図上に特定する方法及び、この方法を使った装置。 |
CN201060594Y (zh) * | 2007-04-03 | 2008-05-14 | 姬成强 | 陶瓷天球仪 |
CN104537606A (zh) * | 2014-12-09 | 2015-04-22 | 中国农业大学 | 地理坐标投影变化方法 |
CN106204684A (zh) * | 2016-07-13 | 2016-12-07 | 国家海洋信息中心 | 一种基于任务流的海洋专题图件自动化编制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7038681B2 (en) * | 2000-03-29 | 2006-05-02 | Sourceprose Corporation | System and method for georeferencing maps |
US9336240B2 (en) * | 2011-07-15 | 2016-05-10 | Apple Inc. | Geo-tagging digital images |
-
2020
- 2020-08-18 CN CN202010831650.6A patent/CN112102430B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB454654A (en) * | 1935-01-03 | 1936-10-05 | Arthur Thomas Keeble | A great circle graphic calculator for use in solving problems in navigation and astronomy |
JP2002366027A (ja) * | 2001-06-05 | 2002-12-20 | Sokken:Kk | 経度間隔の違う地図の、経度目盛りの設定方法と既知経緯度座標を地図上に特定する方法及び、この方法を使った装置。 |
CN201060594Y (zh) * | 2007-04-03 | 2008-05-14 | 姬成强 | 陶瓷天球仪 |
CN104537606A (zh) * | 2014-12-09 | 2015-04-22 | 中国农业大学 | 地理坐标投影变化方法 |
CN106204684A (zh) * | 2016-07-13 | 2016-12-07 | 国家海洋信息中心 | 一种基于任务流的海洋专题图件自动化编制方法 |
Non-Patent Citations (2)
Title |
---|
Assisting video surveillance in micro-spatial environments with a GIS approach;Xiang Li 等;《15th International Conference on Geoinformatics》;1-9 * |
基于GIS海洋环境统计分析系统设计与实现;张增健 等;《地理空间信息》;第21卷(第4期);130-133 * |
Also Published As
Publication number | Publication date |
---|---|
CN112102430A (zh) | 2020-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107192393B (zh) | 一种道路地图数据的维护系统及方法 | |
CN106409129B (zh) | 路况绘制方法及装置 | |
CN107065504B (zh) | 包括具有交互式时间转换的视觉图形指示符的时钟对象 | |
CN107577750B (zh) | 绘制导航数据矢量路口的方法和系统 | |
CN102789361B (zh) | 一种矢量电子地图信息的显示风格设置方法 | |
CN101852619A (zh) | 导航显示方法及装置 | |
CN112102430B (zh) | 单要素变化曲线至地理地图的投影计算方法 | |
CN115130168A (zh) | 基于Blender+BIM+无人机的施工进度模拟与管控方法 | |
CN105809616B (zh) | 基于ArcGIS的1:1万国家标准地形图的制作方法 | |
CN112102425B (zh) | 垂向剖面的要素值至地理地图的投影计算方法及系统 | |
CN104700349B (zh) | 一种挑花图案数字化精确采集方法 | |
CN103345565B (zh) | 一种用于数据展示基于Flash技术的仪表盘构建方法 | |
CN105608127A (zh) | 一种农地权属公示图自动生成方法 | |
CN103345769A (zh) | 全路径遥感影像专题制图法 | |
CN108133502B (zh) | 一种业务流程图的展示方法 | |
CN104142154A (zh) | 一种导航方法、装置和终端 | |
Stuckey | A Comparison of ArcGIS and QGIS for Animation | |
CN113051648A (zh) | 基于知户型的彩色平面方案生成系统及其生成方法 | |
CN113850993A (zh) | 一种交通道路路口的信号渠化渲染方法及系统 | |
JPH07129630A (ja) | 類似部品のcad図面作画方法 | |
CN112214559A (zh) | 一种历史地图动态显示方法 | |
JPH04157560A (ja) | 表作成方法 | |
JPH05120273A (ja) | 表生成方式 | |
CN102194028A (zh) | 机械铭表标牌制作计算机绘图方法 | |
JP3701267B2 (ja) | 図形編集方法および図形編集装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |