【小沐学GIS】基于C++绘制三维数字地球Earth（OpenGL、glfw、glut）第一期

已于 2024-02-15 18:53:58 修改
阅读量8.5k 收藏 103
点赞数 41
分类专栏： GIS OpenGL C/C++ 文章标签： 3d地球 opengl globe earth 瓦片地图 opengl地球 数字地球
于 2022-03-03 23:23:12 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/hhy321/article/details/123264181
版权
C/C++ 同时被 3 个专栏收录
66 篇文章 45 订阅
订阅专栏
GIS
25 篇文章 46 订阅
订阅专栏
OpenGL
21 篇文章 3 订阅
订阅专栏
🍺三维数字地球系列相关文章如下🍺：
1【小沐学GIS】基于C++绘制三维数字地球Earth（456:OpenGL、glfw、glut）第一期
2【小沐学GIS】基于C++绘制三维数字地球Earth（456:OpenGL、glfw、glut）第二期
3【小沐学GIS】基于C++OpenSceneGraph(OSG)绘制三维数字地球Earth（7:OpenGL）
4【小沐学GIS】基于C++QT绘制三维数字地球Earth（8:OpenGL）
5【小沐学GIS】基于C++绘制太阳系SolarSystem（9:OpenGL、glfw、glut）
6【小沐学GIS】基于C#绘制三维数字地球Earth（10:OpenGL）
7【小沐学GIS】基于Python绘制三维数字地球Earth（11:OpenGL）
8【小沐学GIS】基于Android绘制三维数字地球Earth（12:OpenGL）
9【小沐学GIS】基于WebGL绘制三维数字地球Earth（13:OpenGL）

文章目录

在这里插入图片描述

1、计算公式

球的三维坐标表示为：
x 2 + y 2 + z 2 = R 2 x^2 + y^2 + z^2 = R^2 x2+y2+z2=R2

引入球的参数坐标方程进行离散化。以(u,v)表示球面上某一点的坐标，且u,v的取值范围为[0 , 1]，定义(u,v)到(x,y,z)的转换如下：

{ x = R × s i n ( π × v ) c o s ( 2 π × u ) y = R × s i n ( π × v ) s i n ( 2 π × u ) z = R × c o s ( π × v ) \begin{cases} x=R×sin(\pi×v)cos(2\pi×u) \\ y=R×sin(\pi×v)sin(2\pi×u) \\ z=R×cos(\pi×v) \end{cases} ⎩ ⎨ ⎧​x=R×sin(π×v)cos(2π×u)y=R×sin(π×v)sin(2π×u)z=R×cos(π×v)​

2、绘图接口

在OpenGL中，所有图形都是通过分解成三角形的方式进行绘制。
glDrawArrays 和 glDrawElements 的作用都是从一个数据数组中提取数据渲染基本图元。

在这里插入图片描述

  • GL_POINTS：把每一个顶点作为一个点进行处理，顶点n即定义了点n，共绘制N个点。
  • GL_LINES：连接每两个顶点作为一个独立的线段，顶点2n－1和2n之间共定义了n条线段，总共绘制N/2条线段。
  • GL_LINE_STRIP：绘制从第一个顶点到最后一个顶点依次相连的一组线段，第n和n+1个顶点定义了线段n，总共绘制n－1条线段。
  • GL_LINE_LOOP：绘制从第一个顶点到最后一个顶点依次相连的一组线段，然后最后一个顶点和第一个顶点相连，第n和n+1个顶点定义了线段n，总共绘制n条线段。
  • GL_TRIANGLES：把每三个顶点作为一个独立的三角形，顶点3n－2、3n－1和3n定义了第n个三角形，总共绘制N/3个三角形。
  • GL_TRIANGLE_STRIP：绘制一组相连的三角形，对于奇数n，顶点n、n+1和n+2定义了第n个三角形；对于偶数n，顶点n+1、n和n+2定义了第n个三角形，总共绘制N-2个三角形。
  • GL_TRIANGLE_FAN：绘制一组相连的三角形，三角形是由第一个顶点及其后给定的顶点确定，顶点1、n+1和n+2定义了第n个三角形，总共绘制N-2个三角形。

其中：
GL_TRIANGLES：V0V1V2, V3V4V5, V6V7V8……
GL_TRIANGLE_FAN：V0V1V2, V0V2V3, V0V3V4……
GL_TRIANGLE_STRIP：V0V1V2, V1V2V3, V2V3V4……

2.1 glDrawArrays

The glDrawArrays function specifies multiple primitives to render.

void WINAPI glDrawArrays(
   GLenum  mode,
   GLint   first,
   GLsizei count
);

mode:
GL_POINTS, GL_LINE_STRIP, GL_LINE_LOOP, GL_LINES, GL_TRIANGLE_STRIP, GL_TRIANGLE_FAN, GL_TRIANGLES, GL_QUAD_STRIP, GL_QUADS, and GL_POLYGON.

first:
The starting index in the enabled arrays.

count:
The number of indexes to render.

2.2 glDrawElements

The glDrawElements function renders primitives from array data.
The glDrawElements function is only available in OpenGL version 1.1 or later.

void WINAPI glDrawElements(
         GLenum  mode,
         GLsizei count,
         GLenum  type,
   const GLvoid  *indices
);

mode:
The kind of primitives to render. It can assume one of the following symbolic values: GL_POINTS, GL_LINE_STRIP, GL_LINE_LOOP, GL_LINES, GL_TRIANGLE_STRIP, GL_TRIANGLE_FAN, GL_TRIANGLES, GL_QUAD_STRIP, GL_QUADS, and GL_POLYGON.

count:
The number of elements to be rendered.

type:
The type of the values in indices. Must be one of GL_UNSIGNED_BYTE, GL_UNSIGNED_SHORT, or GL_UNSIGNED_INT.

indices:
A pointer to the location where the indices are stored.

绑定VAO时也会自动绑定EBO。
不要在解绑VAO之前解绑EBO（GL_ELEMENT_ARRAY_BUFFER）。

3、代码实现

#pragma once
#include "Point3d.h"
#include <vector>

class EarthGrid3d
{
public:
	EarthGrid3d();

	void generalSphereGrid();
	void generalSphereLines();

	int _uStepsNum;
	int _vStepNum;

	float *_verticesArr;
	int _verticesNum;

	int *_indicesArr;
	int _indicesNum;

	float *_linesArr;
	int _linesNum;
};

#include "EarthGrid3d.h"

EarthGrid3d::EarthGrid3d()
{
	_uStepsNum = 120;
	_vStepNum = 120;
}

void EarthGrid3d::generalSphereGrid()
{
	double ustep = 1 / (double)_uStepsNum, vstep = 1 / (double)_vStepNum;
	int np = 0;
	int nf = 0;

	_verticesNum = (_vStepNum*_uStepsNum + 1)*(3 + 2);
	_verticesArr = new float[_verticesNum];
	
	_indicesNum = (_uStepsNum * 2 + _uStepsNum * 2 * (_vStepNum - 2)) * 3;
	_indicesArr = new int[_indicesNum];
	
	//
	// 北极的一个点
	{
		Point3d pt0 = getSpherePoint(0, 0);
		_verticesArr[5 * np] = pt0.x;
		_verticesArr[5 * np + 1] = pt0.y;
		_verticesArr[5 * np + 2] = pt0.z;
		_verticesArr[5 * np + 3] = 0.5;
		_verticesArr[5 * np + 4] = 0;
		np++;
	}

	// 中间的点
	for (int j = 1; j < _vStepNum; j++) //v-1个点（共v+1个点）
	{
		for (int i = 0; i <= _uStepsNum; i++) //u+1个点 (共u+1个点)
		{
			Point3d pt = getSpherePoint(ustep*i, vstep*j);
			_verticesArr[5 * np] = pt.x;
			_verticesArr[5 * np + 1] = pt.y;
			_verticesArr[5 * np + 2] = pt.z;
			_verticesArr[5 * np + 3] = 1 - ustep * i;
			_verticesArr[5 * np + 4] = vstep * j;
			np++;
		}
	}

	// 南极的一个点
	{
		Point3d pt1 = getSpherePoint(1, 1);
		_verticesArr[5 * np] = pt1.x;
		_verticesArr[5 * np + 1] = pt1.y;
		_verticesArr[5 * np + 2] = pt1.z;
		_verticesArr[5 * np + 3] = 0.5;
		_verticesArr[5 * np + 4] = 1;
		//np++;
	}

	//
	// 上下2行的三角形组
	for (int i = 0; i < _uStepsNum; i++) { 
		//第一层u个三角形
		_indicesArr[nf++] = 0;
		_indicesArr[nf++] = 1 + i;
		_indicesArr[nf++] = 2 + i;
		
		//最后一层u个三角形
		_indicesArr[nf++] = np - 1;
		_indicesArr[nf++] = np - 2 - i;
		_indicesArr[nf++] = np - 3 - i;
	}

	// 中间的v-2行的三角形组
	for (int j = 0; j < _vStepNum - 2; j++) {  //共v-2行
		for (int i = 0; i < _uStepsNum; i++) { //共u列

			/*
			*       |\
			*       | \
			*       |__\
			*/
			_indicesArr[nf++] = (_uStepsNum + 1)*j + 1 + i;
			_indicesArr[nf++] = (_uStepsNum + 1)*(j + 1) + 1 + i;
			_indicesArr[nf++] = (_uStepsNum + 1)*(j + 1) + 2 + i;

			/*
			*       __
			*       \  |
			*        \ |
			*         \|
			*/
			_indicesArr[nf++] = (_uStepsNum + 1)*j + 1 + i;
			_indicesArr[nf++] = (_uStepsNum + 1)*(j + 1) + 2 + i;
			_indicesArr[nf++] = (_uStepsNum + 1)*j + 2 + i;
		}
	}
}

void EarthGrid3d::generalSphereLines()
{
	double ustep = 1 / (double)_uStepsNum, vstep = 1 / (double)_vStepNum;
	int np = 0;

	_linesNum = (_vStepNum*_uStepsNum)*6;
	_linesArr = new float[_linesNum];

	// 绘制24条经线
	for (int i = 0; i < _uStepsNum; i+= _uStepsNum/24)
	{
		for (int j = 0; j < _vStepNum; j++)
		{
			Point3d pt = getSpherePoint(ustep*i, vstep*j, 1.01f);
			_linesArr[3 * np] = pt.x;
			_linesArr[3 * np + 1] = pt.y;
			_linesArr[3 * np + 2] = pt.z;
			np++;

			Point3d pt2 = getSpherePoint(ustep*i, vstep*(j+1), 1.01f);
			_linesArr[3 * np] = pt2.x;
			_linesArr[3 * np + 1] = pt2.y;
			_linesArr[3 * np + 2] = pt2.z;
			np++;
		}
	}

	// 绘制1条纬线（赤道）
	for (int i = 0; i < _uStepsNum; i++)
	{
		int j = _vStepNum / 2;
		//for (int j = 0; j < _vStepNum; j++)
		{
			Point3d pt = getSpherePoint(ustep*i, vstep*j, 1.01f);
			_linesArr[3 * np] = pt.x;
			_linesArr[3 * np + 1] = pt.y;
			_linesArr[3 * np + 2] = pt.z;
			np++;

			Point3d pt2 = getSpherePoint(ustep*(i+1), vstep*j, 1.01f);
			_linesArr[3 * np] = pt2.x;
			_linesArr[3 * np + 1] = pt2.y;
			_linesArr[3 * np + 2] = pt2.z;
			np++;
		}
	}
}

4、运行结果（3d整体地球）

4.1 opengl / glut / c++ （3d）

在这里插入图片描述

4.2 opengl / glfw / glad / stb_image / c++ （3d）

在这里插入图片描述

4.3 opengl / glfw / glad / stb_image /c++ （3d，天空盒，高度贴图）

在这里插入图片描述

4.4 opengl / win32/ glew / FreeImage / c++ (3d，大气层)

在这里插入图片描述

4.5 opengl / freeglut / glew / FreeImage / c++ (3d，法线贴图)

在这里插入图片描述

4.6 opengl / glfw / glad / stb_image / freetype / c++ (2d/3d，加载geojson)

在这里插入图片描述
在这里插入图片描述

4.7 opengl / glfw / glad / FreeImage / c++ (3d，白天层/黑夜层/云层)

在这里插入图片描述
在这里插入图片描述

4.8 opengl / glut / gl3w / c++ (3d，太阳系)

在这里插入图片描述

在这里插入图片描述

4.9 opengl / glut / glew / glm / openssl / c++ (3d，在线地震数据显示)

GeoJSON 是一种用于对各种地理数据结构进行编码的格式。 GeoJSON 对象可以表示几何、要素或 特征。GeoJSON 使用 JSON 标准。 GeoJSONP 源使用相同的 JSON 响应，但 GeoJSONP 响应 包装在函数调用中，eqfeed_callback。

  • 在线地震数据的请求结果json如下：
{
  type: "FeatureCollection",
  metadata: {
    generated: Long Integer,
    url: String,
    title: String,
    api: String,
    count: Integer,
    status: Integer
  },
  bbox: [
    minimum longitude,
    minimum latitude,
    minimum depth,
    maximum longitude,
    maximum latitude,
    maximum depth
  ],
  features: [
    {
      type: "Feature",
      properties: {
        mag: Decimal,
        place: String,
        time: Long Integer,
        updated: Long Integer,
        tz: Integer,
        url: String,
        detail: String,
        felt:Integer,
        cdi: Decimal,
        mmi: Decimal,
        alert: String,
        status: String,
        tsunami: Integer,
        sig:Integer,
        net: String,
        code: String,
        ids: String,
        sources: String,
        types: String,
        nst: Integer,
        dmin: Decimal,
        rms: Decimal,
        gap: Decimal,
        magType: String,
        type: String
      },
      geometry: {
        type: "Point",
        coordinates: [
          longitude,
          latitude,
          depth
        ]
      },
      id: String
    },
    …
  ]
}

在这里插入图片描述
中国地震台网中心：
http://news.ceic.ac.cn/index.html?time=1674619089

  • 显示最近三天内的地震情况：
    在这里插入图片描述

  • 显示最近一个月内的地震情况：
    在这里插入图片描述

4.10 opengl / fltk / glew / curl / openssl / c++ / geojson (3d，在线地震数据显示2)

  • 显示最近小时内的地震情况：
    在这里插入图片描述
  • 显示最近一个月内的地震情况：
    在这里插入图片描述

5、运行结果（3d瓦片地球）

5.1 opengl / glfw / glew / curl / proj4 / gdal / stb_image / c++ （3d，瓦片贴图）

  • 加载卫星影像图
    在这里插入图片描述

  • 加载行政地图
    在这里插入图片描述

  • 加载shp(shapefile)文件和显示：
    在这里插入图片描述

  • 绘制自定义的图形元素（线、面等）
    在这里插入图片描述

  • 切换地球底图的瓦片图源

5.2 opengl / glfw / glad / boost::asio / proj4 / stb_image / c++ （3d，瓦片贴图）

在这里插入图片描述

6、运行结果（2d瓦片地球）

6.1 opengl / sdl / boost::asio / c++ （2d，瓦片贴图）

在这里插入图片描述

以上章节所有地球代码，均在VS2017开发环境编译通过。

结语

如果您觉得该方法或代码有一点点用处，可以给作者点个赞，或打赏杯咖啡；╮(￣▽￣)╭
如果您感觉方法或代码不咋地//(ㄒoㄒ)//，就在评论处留言，作者继续改进；o_O???
如果您需要相关功能的代码定制化开发，可以留言私信作者；(✿◡‿◡)
感谢各位童鞋们的支持！( ´ ▽´ )ﾉ ( ´ ▽´)っ！！！

爱看书的小沐
关注
  • 41
    点赞
  • 踩
  • 103
    收藏
    觉得还不错? 一键收藏
  • 27
    评论
专栏目录
【UE4】C++：代码中实现实时三维绘画
wwj897570的博客
02-20 804
【UE4】C++：代码中实现实时三维绘画
实现三维图形绘制
06-23
Visual studio和DirectX的编程接口，学会编写窗口，实现三维图形的绘制等功能。
C++实现利用OpenGL实现三维绘图
08-08
C++利用OpenGL实现三维绘图，包括详细的代码注释
c++调用python画3维图形
小鱼儿飞丫飞
03-30 1089
文章目录python plot.pyc++程序编写CMakeLists编译运行 python plot.py #coding=utf-8 import matplotlib as mpl from mpl_toolkits.mplot3d import Axes3D import numpy as np import matplotlib.pyplot as plt def plot(): ...
Qt/C++编写地图应用/离线地图下载/路径规划/轨迹回放/海量点/坐标转换
最新发布
Qt/C++视频监控/推流/物联网/大屏系统/自定义控件/UI定制/输入法
07-03 1939
这个地图组件写了很多年了，最初设计的不叫粗糙，最开始只是为了满足项目需要，并没有考虑太多拓展性，比如最初都是按照百度地图写死在代码中，经过这几年大量的现场实际应用，以及大量的用户提出的改进意见，逐渐萌生了彻底重新编写对应地图相关的代码，比如基类子类的设计，各种功能接口通过js函数交互，而不是一堆代码写在网页中。
Visual C++6.0画三维立体图形
热门推荐
七分熟少女
07-22 1万+
在画三维立体图之前，主要是先要找到二维和三维的对应关系，这是转化的关键。 关键代码如下: S[i].x=P[i][1]+sqrt(2)/3.0*(-P[i][0]); S[i].y=P[i][2]+sqrt(2)/3.0*(-P[i][0]); 一、建立边表结构,在同一屏幕上完成三视图和正等轴测投影图 //消隐之前 //编译环境：Visual C++ 6.0，EasyX_201902...
【小沐学GIS】基于C++绘制三维数字地球Earth（OpenGL、glfw、glut）第二期
爱看书的小沐
01-14 2942
地球仪绘制是通过使用OpenGL技术来实现的。首先，需要利用OpenGL绘制球体的方法来创建一个球体模型。可以使用球体的顶点和索引数组来定义球体的几何形状。然后，利用OpenGL的纹理映射功能将一张地球的图片作为纹理映射到整个球体的表面上。通过将纹理坐标映射到球体的顶点上，可以实现将图片贴在球体上的效果。╮(￣▽￣)╭如果您感觉方法或代码不咋地//(ㄒoㄒ)//，就在评论处留言，作者继续改进；o_O???如果您需要相关功能的代码定制化开发，可以留言私信作者；(✿◡‿◡)！！
【小沐学GIS】基于C++QT绘制三维数字地球Earth（OpenGL）
爱看书的小沐
02-13 4371
Qt（官方发音 [kju:t]，音同 cute）是一个跨平台的 C++ 开发库，主要用来开发图形用户界面（Graphical User Interface，GUI）程序，当然也可以开发不带界面的命令行（Command User Interface，CUI）程序。Qt 经常被当做一个 GUI 库，用来开发图形界面应用程序，但这并不是 Qt 的全部；
通用的三视图和正等轴测图的绘制（C++实现）
hanmo22357的博客
05-06 3107
//实验名称：三视图和正轴测图的绘制 //功 能：对于用户输入的有序空间坐标点数组，绘制这些点所组成的三维立体图形的三视图和正等轴测图（函数可以通用，本题中为了演示选择了一个固定的三维物体） //编译环境：Visual Studio 2022，EasyX_20220116 #include<iostream> #include<graphics.h> #include<conio.h> using namespace std; const int gap =
用C++实现WebGIS
C++与Web
11-06 6748
用C++实现的BS架构的WebGIS，Demo是某公司的办公室地图，显示了会议室，工位，打印机和卫生间的布局，共4层；推而广之，同样适用于大型商场，医院，学校，工厂，小区及其他需要可视化的环境
使用OpenGL实现绘制三维图形.rar_thisks1_visual c++ opengl_三维图形绘制_使用OpenGL实现
09-23
在本项目中，我们将探讨如何在Visual C++ 6.0环境下利用OpenGL来绘制三维图形。 首先，你需要安装支持OpenGL的开发环境，如在Windows上安装Visual C++ 6.0，因为它是许多程序员喜欢的老版本IDE，虽然现在有更新的...
VC++ 三维图形绘制实例
11-22
VC++ 三维图形绘制实例
C++做的电子地图
01-05
C++ 电子地图 有兴趣的可以看看！附可运行的文件 ！
Opengl 绘制的地球
12-21
利用opengl绘制的地球 Data下为纹理图片，必须和执行文件在同一目录
opengl c++ 三维图形
05-20
用c++编写opengl的程序，要求实现三维图形的旋转以及其各种功能
OpenGL绘制三维地形图
01-17
根据高程数据，利用OpenGL绘制三维地形，高程数据为txt文件
VC++使用OpenGL绘制三维图形的实例.rar
03-08
在本实例中，我们将深入探讨如何使用VC++集成开发环境结合OpenGL库来创建并实现三维图形的绘制。OpenGL是一个跨语言、跨平台的编程接口，主要用于渲染2D、3D矢量图形，广泛应用于游戏开发、科学可视化以及工程设计等...
openGL.zip_OpenGL 球体_c++语言画球_opengl界面_三维球体 opengl
09-21
本项目是基于OpenGL库，利用C++编程语言实现的一个三维球体的绘制，旨在展示如何在OpenGL环境下创建一个简单的图形界面。 在C++中，使用OpenGL进行图形渲染通常涉及以下几个关键步骤： 1. **初始化OpenGL上下文**:...
openGL实现三维点云显示-C++代码类资源
10-01
opengl实现三维点云的显示，放大，缩小，旋转，多角度浏览等功能。 opengl实现三维点云的显示，放大，缩小，旋转，多角度浏览等功能。
编写使用glfw glad绘制卧室三维场景的完整C++代码
06-10
好的，下面是一个使用glfw glad绘制卧室三维场景的完整C++代码，代码中包含了一些基本的场景渲染和交互功能，可以作为参考。 ``` #include <glad/glad.h> #include <GLFW/glfw3.h> #include <iostream> void framebuffer_size_callback(GLFWwindow* window, int width, int height); void processInput(GLFWwindow *window); // 窗口大小 const unsigned int SCR_WIDTH = 800; const unsigned int SCR_HEIGHT = 600; int main() { // 初始化glfw glfwInit(); glfwWindowHint(GLFW_CONTEXT_VERSION_MAJOR, 3); glfwWindowHint(GLFW_CONTEXT_VERSION_MINOR, 3); glfwWindowHint(GLFW_OPENGL_PROFILE, GLFW_OPENGL_CORE_PROFILE); // 创建窗口 GLFWwindow* window = glfwCreateWindow(SCR_WIDTH, SCR_HEIGHT, "Bedroom Scene", NULL, NULL); if (window == NULL) { std::cout << "Failed to create GLFW window" << std::endl; glfwTerminate(); return -1; } glfwMakeContextCurrent(window); // 初始化glad if (!gladLoadGLLoader((GLADloadproc)glfwGetProcAddress)) { std::cout << "Failed to initialize GLAD" << std::endl; return -1; } // 设置视口大小 glViewport(0, 0, SCR_WIDTH, SCR_HEIGHT); // 注册窗口大小改变的回调函数 glfwSetFramebufferSizeCallback(window, framebuffer_size_callback); // 渲染循环 while (!glfwWindowShouldClose(window)) { // 处理输入 processInput(window); // 渲染指令 glClearColor(0.2f, 0.3f, 0.3f, 1.0f); glClear(GL_COLOR_BUFFER_BIT); // 交换缓冲区 glfwSwapBuffers(window); glfwPollEvents(); } // 清理 glfwTerminate(); return 0; } void framebuffer_size_callback(GLFWwindow* window, int width, int height) { glViewport(0, 0, width, height); } void processInput(GLFWwindow *window) { if (glfwGetKey(window, GLFW_KEY_ESCAPE) == GLFW_PRESS) glfwSetWindowShouldClose(window, true); } ``` 这个代码只是一个简单的框架，你需要在这个基础上添加你自己的3D场景渲染代码。希望对你有所帮助！
评论 27
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值