基于Arcgis的投影坐标与经纬度坐标的相互转换代码
2012-02-07 15:46
435 查看
转自:http://www.ggiiss.com/gis/9e/338.html
现在好多人在使用ArcObject的时候都可能需要作经纬度和平面坐标的相互转换。由于经纬度是球面坐标,平面坐标是X-Y的笛卡尔坐标系统,所以这是一个看起来比较难的问题。
好多人一上来就搬出地图学、地图投影学或者测绘学中的投影公式(如高斯投影 或 墨卡托投影),又是基准坐标又是角度的搞的人头都大了。
实际上要想实现这个功能非常easy。 ARC Engine中的IPoint就可以进行投影和反投影运算了。
投影过程(C#):
/// flatref 投影的坐标系统,这里的54013是世界投影 ,世界投影所有经纬度都可以转换为平面坐标,但是由于投影面积大失真也会比较大(相当于把整个地球劈成一片片的,然后拉伸最后贴到平面上,失真当然大了)。当然也可以选择精度更高的平面如:esriSRProjCS_Beijing1954GK_23N 对应数值21483 仅仅把北京附近的地球平面拉伸铺在平面上,由于投影面积变小,所以投影经度提高。但是由于面积变小,所以有些经度纬度不能转换,比如所美国的精度纬度用北京投影就投不了。
flatref = pfactory.CreateProjectedCoordinateSystem(54013);
//没什么说的,标准大地经纬度,可以将X-Y逆投影为经度和纬度
earthref = pfactory.CreateGeographicCoordinateSystem((int)esriSRGeoCSType.esriSRGeoCS_NAD1983);
/// 将经纬度点转换为平面坐标。
private IPoint GetProject(double x, double y)
{
IPoint pt = new PointClass();
pt.PutCoords(x, y);
IGeometry geo = (IGeometry)pt;
geo.SpatialReference = earthref;
geo.Project(flatref);
return pt;
}
/// 将平面坐标转换为经纬度。
private IPoint GetGeo(double x, double y)
{
IPoint pt = new PointClass();
pt.PutCoords(x, y);
IGeometry geo = (IGeometry)pt;
geo.SpatialReference = flatref;
geo.Project(earthref);
double xx = pt.X;
return pt;
}
实际上IPoint的投影和任何地图都没什么大关系,完全可以不用地图,直接调用IPoint进行投影的转换
现在好多人在使用ArcObject的时候都可能需要作经纬度和平面坐标的相互转换。由于经纬度是球面坐标,平面坐标是X-Y的笛卡尔坐标系统,所以这是一个看起来比较难的问题。
好多人一上来就搬出地图学、地图投影学或者测绘学中的投影公式(如高斯投影 或 墨卡托投影),又是基准坐标又是角度的搞的人头都大了。
实际上要想实现这个功能非常easy。 ARC Engine中的IPoint就可以进行投影和反投影运算了。
投影过程(C#):
/// flatref 投影的坐标系统,这里的54013是世界投影 ,世界投影所有经纬度都可以转换为平面坐标,但是由于投影面积大失真也会比较大(相当于把整个地球劈成一片片的,然后拉伸最后贴到平面上,失真当然大了)。当然也可以选择精度更高的平面如:esriSRProjCS_Beijing1954GK_23N 对应数值21483 仅仅把北京附近的地球平面拉伸铺在平面上,由于投影面积变小,所以投影经度提高。但是由于面积变小,所以有些经度纬度不能转换,比如所美国的精度纬度用北京投影就投不了。
flatref = pfactory.CreateProjectedCoordinateSystem(54013);
//没什么说的,标准大地经纬度,可以将X-Y逆投影为经度和纬度
earthref = pfactory.CreateGeographicCoordinateSystem((int)esriSRGeoCSType.esriSRGeoCS_NAD1983);
/// 将经纬度点转换为平面坐标。
private IPoint GetProject(double x, double y)
{
IPoint pt = new PointClass();
pt.PutCoords(x, y);
IGeometry geo = (IGeometry)pt;
geo.SpatialReference = earthref;
geo.Project(flatref);
return pt;
}
/// 将平面坐标转换为经纬度。
private IPoint GetGeo(double x, double y)
{
IPoint pt = new PointClass();
pt.PutCoords(x, y);
IGeometry geo = (IGeometry)pt;
geo.SpatialReference = flatref;
geo.Project(earthref);
double xx = pt.X;
return pt;
}
实际上IPoint的投影和任何地图都没什么大关系,完全可以不用地图,直接调用IPoint进行投影的转换
相关文章推荐
- 【代码】ArcGIS For Flex 由地理坐标转换成自定义投影坐标系统参考
- ArcEngine下投影坐标和经纬度坐标的相互转换
- ArcGIS Engine 下投影坐标和经纬度坐标的相互转换
- ArcEngine下投影坐标和经纬度坐标的相互转换
- ArcGIS for Android 如何将经纬度坐标转换成投影坐标
- ArcEngine下投影坐标和经纬度坐标的相互转换
- WGS84经纬度坐标与BJ54坐标之间的转换(Arcgis)(转)
- 经纬度到平面坐标的相互转换
- 经纬度点与投影坐标系下坐标点的互相转换
- 【ArcGIS for Android】经纬度坐标、地图投影坐标、屏幕坐标互相转换
- 经纬度到平面坐标的相互转换
- 使用ArcGIS实现WGS84经纬度坐标到北京54高斯投影坐标的转换
- 【转】ArcGIS投影转换与坐标转换
- 经纬度坐标与高斯坐标的转换代码(转)
- 火星坐标、百度坐标、WGS-84坐标相互转换及墨卡托投影坐标转经纬度JavaScript版
- 火星坐标、百度坐标、WGS-84坐标相互转换及墨卡托投影坐标转经纬度JavaScript版 (网络搜集)
- 经纬度坐标与高斯坐标的转换代码
- 经纬度点与投影坐标系下坐标点的互相转换
- 【ArcGIS for Android】经纬度坐标、地图投影坐标、屏幕坐标互相转换
- 利用米勒投影法实现经纬度和平面坐标的相互转换