GDAL中的地理坐标系、投影坐标系及其相互转换

最新推荐文章于 2024-05-15 07:38:09 发布
最新推荐文章于 2024-05-15 07:38:09 发布
阅读量563 收藏 4
点赞数 5
分类专栏： 遥感开发学习 文章标签： c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_35662333/article/details/137022136
版权

目录

地理坐标系

国内常用地理坐标系

投影坐标系

国内常用投影坐标系（不推荐使用）

坐标转换

地理坐标转为投影坐标

投影坐标转为地理坐标

地理坐标系

原理参考这篇文章：
地理坐标系与投影坐标系区别与联系
https://yunxingluoyun.blog.csdn.net/article/details/123970678

国内常用地理坐标系

#include <cstdio>
#include "gdal_priv.h"
#include <iostream>

int main()
{
	OGRSpatialReference oSRS1;
	oSRS1.SetGeogCS("自定义地理坐标系",// 定义地理坐标系名称
		"WGS_1984",// 基准面
		"WGS84 椭球",// 椭球体名称
		SRS_WGS84_SEMIMAJOR, SRS_WGS84_INVFLATTENING,// WGS84椭球体的长半轴，WGS84椭球体扁率的倒数
		"Greenwich", 0.0,// 格林尼治子午线
		"degree", 0.0174532925199433);// 角度度量单位
	char* WGS84_WTK = NULL;
	oSRS1.exportToPrettyWkt(&WGS84_WTK);// 以规整的WTK格式输出地理坐标系的信息
	std::cout << WGS84_WTK << std::endl;
	// CGCS2000
	OGRSpatialReference oSRS2;
	oSRS2.SetWellKnownGeogCS("EPSG:4490");
	char* CGCS2000_WTK = NULL;
	oSRS2.exportToPrettyWkt(&CGCS2000_WTK);
	std::cout << CGCS2000_WTK << std::endl;
	// Beijing_1954
	OGRSpatialReference oSRS3;
	oSRS3.SetWellKnownGeogCS("EPSG:4214");
	char* Beijing_1954_WTK = NULL;
	oSRS3.exportToPrettyWkt(&Beijing_1954_WTK);
	std::cout << Beijing_1954_WTK << std::endl;
	// Xian_1980
	OGRSpatialReference oSRS4;
	oSRS4.SetWellKnownGeogCS("EPSG:4610");
	char* Xian_1980_WTK = NULL;
	oSRS4.exportToPrettyWkt(&Xian_1980_WTK);
	std::cout << Xian_1980_WTK << std::endl;

}

结果：

GEOGCS["自定义地理坐标系",
    DATUM["WGS_1984",
        SPHEROID["WGS84 椭球",6378137,298.257223563]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST]]
    
GEOGCS["China Geodetic Coordinate System 2000",
    DATUM["China_2000",
        SPHEROID["CGCS2000",6378137,298.257222101,
            AUTHORITY["EPSG","1024"]],
        AUTHORITY["EPSG","1043"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST],
    AUTHORITY["EPSG","4490"]]
    
GEOGCS["Beijing 1954",
    DATUM["Beijing_1954",
        SPHEROID["Krassowsky 1940",6378245,298.3,
            AUTHORITY["EPSG","7024"]],
        AUTHORITY["EPSG","6214"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST],
    AUTHORITY["EPSG","4214"]]
   
GEOGCS["Xian 1980",
    DATUM["Xian_1980",
        SPHEROID["IAG 1975",6378140,298.257,
            AUTHORITY["EPSG","7049"]],
        AUTHORITY["EPSG","6610"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST],
    AUTHORITY["EPSG","4610"]]

投影坐标系

设置横轴墨卡托投影的函数SetTM()有六个参数：

国内常用投影坐标系（不推荐使用）

推荐使用例3中的方式使用坐标系，EPSG请查看网址：

EPSG.io: Coordinate Systems Worldwide
注意：这种方式定义的坐标轴与例3相反，如果使用这样方法，后面的poCT->Transform(1, &x, &y)影改为poCT->Transform(1, &y, &x)。

#include <cstdio>
#include "gdal_priv.h"
#include <iostream>

int main()
{
	OGRSpatialReference oSRS1;
	oSRS1.SetProjCS("UTM 17(WGS84) in northern hemisphere.");// 设置投影坐标系名称
	oSRS1.SetWellKnownGeogCS("EPSG:4326");// 设置地理坐标系
	oSRS1.SetUTM(17, TRUE);//设置UTM投影的投影参数
	char* UTM17_WTK = NULL;
	oSRS1.exportToPrettyWkt(&UTM17_WTK);
	std::cout << UTM17_WTK << std::endl;

	OGRSpatialReference oSRS2;
	oSRS2.SetProjCS("CGCS2000 / Gauss-Kruger CM 117E");
	oSRS2.SetWellKnownGeogCS("EPSG:4490");
	oSRS2.SetTM(0,117,1,500000,0);//起始纬度，中央经线，比例因子，东偏移量，北偏移量
	char* CGCS2000_GK_117E_WTK = NULL;
	oSRS2.exportToPrettyWkt(&CGCS2000_GK_117E_WTK);
	std::cout << "CGCS2000 / Gauss-Kruger CM 117E" << std::endl;
	std::cout <<"是否为投影坐标系：" << oSRS2.IsProjected() << std::endl;
	std::cout << CGCS2000_GK_117E_WTK << std::endl;

	OGRSpatialReference oSRS3;
	oSRS3.SetProjCS("CGCS2000 3 Degree GK Zone 39");
	oSRS3.SetWellKnownGeogCS("EPSG:4490");
	oSRS3.SetTM(0, 117, 1, 39500000, 0);//起始纬度，中央经线，比例因子，东偏移量，北偏移量
	char* CGCS2000_GK_Zone39_WTK = NULL;
	oSRS3.exportToPrettyWkt(&CGCS2000_GK_Zone39_WTK);
	std::cout << "CGCS2000 3 Degree GK Zone 39" << std::endl;
	std::cout << "是否为投影坐标系：" << oSRS3.IsProjected() << std::endl;
	std::cout << CGCS2000_GK_Zone39_WTK << std::endl;

	OGRSpatialReference oSRS4;
	oSRS4.SetProjCS("Beijing 1954 3 Degree GK CM 117E");
	oSRS4.SetWellKnownGeogCS("EPSG:4214");
	oSRS4.SetTM(0, 117, 1, 500000, 0);//起始纬度，中央经线，比例因子，东偏移量，北偏移量
	char* Beijing1954_GK_117E_WTK = NULL;
	oSRS4.exportToPrettyWkt(&Beijing1954_GK_117E_WTK);
	std::cout << "Beijing 1954 3 Degree GK CM 117E" << std::endl;
	std::cout << "是否为投影坐标系：" << oSRS4.IsProjected() << std::endl;
	std::cout << Beijing1954_GK_117E_WTK << std::endl;

	OGRSpatialReference oSRS5;
	oSRS5.SetProjCS("Xian 1980 3 Degree GK Zone 39");
	oSRS5.SetWellKnownGeogCS("EPSG:4610");
	oSRS5.SetTM(0, 117, 1, 39500000, 0);//起始纬度，中央经线，比例因子，东偏移量，北偏移量
	char* Xian_1980_GK_Zone39_WTK = NULL;
	oSRS5.exportToPrettyWkt(&Xian_1980_GK_Zone39_WTK);
	std::cout << "Xian 1980 3 Degree GK Zone 39" << std::endl;
	std::cout << "是否为投影坐标系：" << oSRS5.IsProjected() << std::endl;
	std::cout << Xian_1980_GK_Zone39_WTK << std::endl;
}

 结果：

PROJCS["UTM 17(WGS84) in northern hemisphere.",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-81],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH]]
    
CGCS2000 / Gauss-Kruger CM 117E
是否为投影坐标系：1
PROJCS["CGCS2000 / Gauss-Kruger CM 117E",
    GEOGCS["China Geodetic Coordinate System 2000",
        DATUM["China_2000",
            SPHEROID["CGCS2000",6378137,298.257222101,
                AUTHORITY["EPSG","1024"]],
            AUTHORITY["EPSG","1043"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4490"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",117],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH]]
    
CGCS2000 3 Degree GK Zone 39
是否为投影坐标系：1
PROJCS["CGCS2000 3 Degree GK Zone 39",
    GEOGCS["China Geodetic Coordinate System 2000",
        DATUM["China_2000",
            SPHEROID["CGCS2000",6378137,298.257222101,
                AUTHORITY["EPSG","1024"]],
            AUTHORITY["EPSG","1043"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4490"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",117],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",39500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH]]
    
Beijing 1954 3 Degree GK CM 117E
是否为投影坐标系：1
PROJCS["Beijing 1954 3 Degree GK CM 117E",
    GEOGCS["Beijing 1954",
        DATUM["Beijing_1954",
            SPHEROID["Krassowsky 1940",6378245,298.3,
                AUTHORITY["EPSG","7024"]],
            AUTHORITY["EPSG","6214"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4214"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",117],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH]]
    
Xian 1980 3 Degree GK Zone 39
是否为投影坐标系：1
PROJCS["Xian 1980 3 Degree GK Zone 39",
    GEOGCS["Xian 1980",
        DATUM["Xian_1980",
            SPHEROID["IAG 1975",6378140,298.257,
                AUTHORITY["EPSG","7049"]],
            AUTHORITY["EPSG","6610"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4610"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",117],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",39500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH]]

坐标转换

地理坐标转为投影坐标

#include <cstdio>
#include "gdal_priv.h"

int main()
{
    OGRSpatialReference oSourceSRS, oTargetSRS;
    OGRCoordinateTransformation* poCT;
    oSourceSRS.importFromEPSG(4610); // 地理坐标系：Xian 1980
    oTargetSRS.importFromEPSG(2383); // 投影坐标系：Xian 1980 / 3-degree Gauss-Kruger CM 114E
    
    char* Xian_1980 = NULL;
    oSourceSRS.exportToPrettyWkt(&Xian_1980);
    char* Xian_1980_GK_114E_WTK = NULL;
    oTargetSRS.exportToPrettyWkt(&Xian_1980_GK_114E_WTK);
    printf("%s,%s\n", "oSourceSRS:\n", Xian_1980);
    printf("%s,%s\n", "oTargetSRS:\n", Xian_1980_GK_114E_WTK);


    poCT = OGRCreateCoordinateTransformation(&oSourceSRS, &oTargetSRS);
    double x, y;
    x = 38.8; //纬度
    y = 113.6; //经度
    printf("经纬度坐标：%.9lf	%.9lf", x, y);
    if (poCT == NULL || !poCT->Transform(1, &x, &y))
    {
        printf("Transformation failed.\n");
    }
    else
    {
        printf("\n平面坐标：%.9lf	%.9lf", x, y);
    }
    
}

 结果：

oSourceSRS:
,GEOGCS["Xian 1980",
    DATUM["Xian_1980",
        SPHEROID["IAG 1975",6378140,298.257,
            AUTHORITY["EPSG","7049"]],
        AUTHORITY["EPSG","6610"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST],
    AUTHORITY["EPSG","4610"]]
oTargetSRS:
,PROJCS["Xian 1980 / 3-degree Gauss-Kruger CM 114E",
    GEOGCS["Xian 1980",
        DATUM["Xian_1980",
            SPHEROID["IAG 1975",6378140,298.257,
                AUTHORITY["EPSG","7049"]],
            AUTHORITY["EPSG","6610"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4610"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",114],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Northing",NORTH],
    AXIS["Easting",EAST],
    AUTHORITY["EPSG","2383"]]
经纬度坐标：38.800000000        113.600000000
平面坐标：4296379.277209882     465252.023887851

投影坐标转为地理坐标

#include <cstdio>
#include "gdal_priv.h"

int main()
{
    OGRSpatialReference oSourceSRS, oTargetSRS;
    OGRCoordinateTransformation* poCT;
    oSourceSRS.importFromEPSG(2383); // Xian 1980 / 3-degree Gauss-Kruger CM 114E
    oTargetSRS.importFromEPSG(4610); // Xian 1980 

    char* Xian_1980_GK_114E_WTK = NULL;
    oSourceSRS.exportToPrettyWkt(&Xian_1980_GK_114E_WTK);
    char* Xian_1980 = NULL;
    oTargetSRS.exportToPrettyWkt(&Xian_1980);

    printf("%s,%s\n", "oSourceSRS:\n", Xian_1980);
    printf("%s,%s\n", "oTargetSRS:\n", Xian_1980_GK_114E_WTK);


    poCT = OGRCreateCoordinateTransformation(&oSourceSRS, &oTargetSRS);
    double x, y;  
    x = 4296379.277209882; //纬度
    y = 465252.023887851; //经度
    printf("\n平面坐标：%.9lf	%.9lf", x, y);
    
    if (poCT == NULL || !poCT->Transform(1, &x, &y))
    {
        printf("Transformation failed.\n");
    }
    else
    {
        printf("\n经纬度坐标：%.9lf  %.9lf", x, y);
    }
    
}

结果：

oSourceSRS:
,GEOGCS["Xian 1980",
    DATUM["Xian_1980",
        SPHEROID["IAG 1975",6378140,298.257,
            AUTHORITY["EPSG","7049"]],
        AUTHORITY["EPSG","6610"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AXIS["Latitude",NORTH],
    AXIS["Longitude",EAST],
    AUTHORITY["EPSG","4610"]]
oTargetSRS:
,PROJCS["Xian 1980 / 3-degree Gauss-Kruger CM 114E",
    GEOGCS["Xian 1980",
        DATUM["Xian_1980",
            SPHEROID["IAG 1975",6378140,298.257,
                AUTHORITY["EPSG","7049"]],
            AUTHORITY["EPSG","6610"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4610"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",114],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Northing",NORTH],
    AXIS["Easting",EAST],
    AUTHORITY["EPSG","2383"]]

平面坐标：4296379.277209882     465252.023887851
经纬度坐标：38.800000000  113.600000000

参考资料：

《GDAL源码剖析与开发指南》
ogr_srs_api.h:空间参考系C api — GDAL 文档

去往火星
关注
  • 5
    点赞
  • 踩
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
解决python gdal投影坐标系转换的问题
01-20
要将xian80地理坐标系转换成投影坐标系： xian1980 = GEOGCS[GCS_Xian_1980, DATUM[Xian_1980, SPHEROID[Xian_1980,6378140.0,298.257]], PRIMEM[Greenwich,0.0], UNIT[Degree,0.0174532925199433]] ...
基于GDAL的地理坐标与投影坐标互转
HeyLoong的博客
03-05 3295
在地理数据处理过程中，往往需要进行坐标之间的转换，将地理坐标转换为投影坐标是最常见的需求。然而由于在指定投影参数下将地理坐标转换为投影坐标时涉及过多参数，如何通过C++快速实现坐标转换成为亟需解决的难题。 GDAL是处理地理数据的高效开源库，支持C++、Python等多种语言，可实现绝大部分数据处理需求，因此本人采用C++语言基于GDAL库进行地理坐标转换为投影坐标处理。 操作如下： void projection(double &dLong, double &dLat) { //投影坐标
地理坐标系与投影坐标系区别与联系
陨星落云的博客
04-05 4293
地理坐标系与投影坐标系区别 空间参考中主要包含大地水准面、地球椭球体、投影坐标系等几部分内容。地图投影就是把地球表面的任意点，利用一定数学法则，转换到地图平面上的理论和方法，一般有两种坐标系来进行表示，分别是地理坐标系和投影坐标系。如下图所示，描述了地理坐标系与投影坐标系之间的关系。 参考椭球体 表1 常见的参考椭球体及其参数 椭球名称 长半轴（米） 短半轴（米） 扁率的倒数 Clarke1866 6378206.4 6356583.79999898 294.9786982 Besse
GDAL开发：地理坐标与投影坐标的转换
Tjmies的博客
06-30 3136
最近运用到QGis的开发，通过gdal库（操作各种栅格地理数据格式的库）解析shp文件得到想要的地理数据，开发过程中，想得到某段道路的折线长度，网上找很多资料都没有类似的代码示例（主要是讲解地理坐标和投影坐标的定义，很清晰、易理解，但没有找到自己想要的示例代码）。 开发环境：visual studio2013 （windows）gdal-2.3.1 例：需要得到图中这段道路的长度（gdal开发） 如图所示，gdal代码需要获取道路的长度921.130米 OGRLineString *lin...
GDAL库进行坐标转换实例代码
欧特GO
05-19 1274
常用到的地理坐标系有WGS84，beijing54，xian80，CGCS2000这四种。GDAL可以像proj4那样自定义坐标系，也可以仅通过字符串定义一些常用的坐标系，但本文认为最方便的还是通过EPSG数据库定义的代码来定义一个地理坐标系统；
3857投影坐标和4326地理坐标互转的python代码
01-29
相当于wgs84和墨卡托的相互转换 点击文件，输入的文件名改成自己的，直接运行集合输出 程序实现地理坐标的相互转换，可以学习到GIS领域gdal库的相关运用，同时学习基础的python语法
Qt+GDAL库制作经纬度坐标转换工具，支持CGCS2000/WGS84/XIAN80/BEIJING54四种常用坐标系
04-14
5、支持经纬度转投影坐标以及投影坐标转经纬度。 详情查看博客：...
proj4j-0.1.1.jar GIS开发空间坐标转换、投影可以使用
11-01
maven依赖 <groupId>io.jeo <artifactId>proj4j <version>0.1.1 </dependency>
C++ 容器（二）——容器操作
鲸落博客
05-09 517
1、动态数组Vector没有提供头部增删（push_front()、pop_front()）功能（和数组一样，删除除尾部索引的需要修改后面元素的索引值，很繁琐），链表list有。代码示例：通过上述函数，实现动态数组vector的增删改查。
C++初探_虚基类的声明
qq_43634157的博客
05-09 137
存在的原因：希望在多重继承时（即：派生类有多个直接基类，而这些直接基类是从同一个祖先基类派生而来），派生类只继承一个共同的祖先基类对象。分析：禁止信息通过中间类自动传递给虚基类。因此，在派生类中，必须通过显式地调用虚基类构造函数进行初始化。
初识C语言——第十四天
硬件菜鸟的成长之路
05-07 742
【代码】初识C语言——第十四天。
C语言如何在链表中的指定位置插⼊结点？
最新发布
PAP
05-15 136
C语言如何在链表中的指定位置插⼊结点？
【C++ list所有函数举例如何使用】
m0_45463480的博客
05-09 354
类中一些常用函数的示例用法。您可以使用这些示例来操作列表，例如添加、删除、访问元素等。请记住，在实际应用中，应该根据具体情况选择合适的函数来操作。，以确保数据结构的正确性和高效性。希望这些示例对您有所帮助。如果您有任何疑问，请随时提出。是一个双向链表，提供了在列表中添加、删除、访问元素等操作的方法。
C++高精度算法-加法
weixin_60673089的博客
05-12 439
当一个数很大无法被普通变量所表示（不包括特殊关键字，只谈论STL范围内的）那么，就要用一点歪招了，那就是用数组存储每个数位上的内容。如果强制表示，就会溢出，这里的溢出大家可以自行百度，反正就是会出一些-5665434之类的数。这的确是一个令人非常头疼的问题，不过只要按照，向下一位进位的法则，也就没有那么难了。在C++的运算中，难免会出现很大很大的数，下面是各个关键字的表示范围。1.因为最后一位也有进位，那么可能长度会+1,所以要写一个特判。这里需要注意的是，字符的下标是从0开始的，所以不要写成1了。
【C++】多态 — 多态的细节补充（下篇）
weixin_42976474的博客
05-09 291
我们学习了多态的形式和如何使用多态，这一章我们将来讲一讲多态的原理…
C++中内部类详解
qq_54933168的博客
05-13 392
举个通俗易懂的例子，A相当于一个大的农场，B相当于果园，当陌生人问我做什么的时候，我说我有一个大的农场，农村里的果子很好吃，我只需要让A外想吃果子的人可以吃到果子，置于我农场里的果园B怎么种果子是由我决定的，种果子的方法是商业机密，不能外传，这个时候我们就可以用内部类。内部类是一个独立的类，他不属于外部类，更不能通过外部类的对象去访问内部类的成员。但是内部类可以通过外部类的对象参数来访问外部类中的所有成员。2.注意内部类可以直接访问外部类的static成员，不需要外部类的对象或类名。
C++：STL-vector
2302_79054145的博客
05-09 573
vector是表示可变大小数组的序列容器。就像数组一样，vector也采用的连续存储空间来存储元素。也就是意味着可以采用下标对vector的元素进行访问，和数组一样高效。但是又不像数组，它的大小是可以动态改变的，而且它的大小会被容器自动处理。本质讲，vector使用动态分配数组来存储它的元素。当新元素插入时候，这个数组需要被重新分配大小为了增加存储空间。其做法是，分配一个新的数组，然后将全部元素移到这个数组。
python gdal 影像 地理坐标转投影坐标
06-02
在 GDAL 中，可以通过 `osr` 模块中的 `SpatialReference` 和 `CoordinateTransformation` 类来实现地理坐标与投影坐标的转换。下面是一个简单的例子： ```python from osgeo import gdal, osr # 定义源影像和目标影像的文件路径 src_file = 'input.tif' dst_file = 'output.tif' # 定义目标投影 dst_proj = 'EPSG:32650' # 例如，设置目标坐标系为 WGS84 UTM Zone 50N # 打开源影像文件 src_ds = gdal.Open(src_file) # 获取源影像的地理信息和投影信息 src_geo = src_ds.GetGeoTransform() src_proj = osr.SpatialReference() src_proj.ImportFromWkt(src_ds.GetProjection()) # 定义目标投影 dst_proj = osr.SpatialReference() dst_proj.SetFromUserInput(dst_proj) # 创建坐标转换对象 coord_trans = osr.CoordinateTransformation(src_proj, dst_proj) # 计算转换后的坐标 x, y, _ = coord_trans.TransformPoint(src_geo[0], src_geo[3]) # 输出转换后的坐标 print('X: {}'.format(x)) print('Y: {}'.format(y)) # 关闭影像文件 src_ds = None ``` 在上面的例子中，我们打开了一个名为 `input.tif` 的影像文件，并获取了其地理信息和投影信息。然后，我们定义了一个名为 `dst_proj` 的目标投影，创建了一个名为 `coord_trans` 的坐标转换对象，并使用 `coord_trans.TransformPoint` 方法计算了转换后的坐标。最后，我们输出了转换后的坐标，并关闭了影像文件。 需要注意的是，在上面的例子中，我们使用了 `osr.SpatialReference.SetFromUserInput` 方法来设置目标投影，该方法可以根据用户输入的字符串自动识别投影坐标系。如果你已经知道目标投影的 EPSG 代码或 WKT 字符串，也可以使用 `osr.SpatialReference.ImportFromEPSG` 或 `osr.SpatialReference.ImportFromWkt` 方法来设置目标投影。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值